潮流计算程序及计算结果.docVIP

附表1： 计算机计算潮流程序： %本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号 clear; n=13;%input(请输入节点数：n=); nl=13;%input(请输入支路数：nl=); isb=1;%input(请输入平衡母线节点号：isb=); pr=0.00001;%input(请输入误差精度：pr=); B1=[];%input(请输入由支路参数形成的矩阵： B1=); B2=[];%input(请输入各节点参数形成的矩阵： B2=); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %－－－－－－－修改部分－－－－－－－－－－－－ ym=0; SB=100;UB=220; %ym=input(您输入的参数是标么值？（若不是则输入一个不为零的数值）); if ym~=0 %SB=input(请输入功率基准值:SB=); %UB=input(请输入电压基准值:UB=); YB=SB./UB./UB; BB1=B1; BB2=B2; for i=1:nl B1(i,3)=B1(i,3)*YB; B1(i,4)=B1(i,4)./YB; end disp(B1矩阵B1=); disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB; B2(i,2)=B2(i,2)./SB; B2(i,3)=B2(i,3)./UB; B2(i,4)=B2(i,4)./UB; B2(i,5)=B2(i,5)./SB; end disp(B2矩阵B2=); disp(B2) end % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于低压侧 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 end %求导纳矩阵 disp(导纳矩阵 Y=); disp(Y) %---------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); %PV节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量 end %=================================================================== P=real(S);Q=imag(S); ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; while IT2~=0 IT2=0;a=a+1; for i=1:n if i~=isb %非平衡节点 C(i)=0;D(i)=0; for j1=1:n C