本文提出了利用 MATLAB 语言来进行电力系统潮流计算的方法。通过实例, 说明了该 方法编程简便、 运算效率高并符合人们的思维习惯.
Zab=0.54+0.65j;Zbc=0.62+0.5j;Zcd=0.6+0.35j;Zbe=0.72+0.75j;
Zef=1.0+0.55j;Zeg=(0.65+0.35j);Zbh=(0.9+0.5j);
Sb=(0.6+0.45j);Sc=(0.4+0.3j);Sd=(0.4+0.28j);
Se=(0.6+0.4j);Sf=(0.4+0.3j);Sg=(0.5+0.35j);
Sh=(0.5+0.4j);V0=(10);Va=10.5;
V=[10 10 10 10 10 10 10];
E=[1 1 1 1 1 1 1];
while (max(E)>0.001)
B=V;
SScd=[real(Sd)^2+imag(Sd)^2]/V(1)^2*Zcd;
Scd=Sd+SScd;
Sc1=Sd+SScd+Sc;
SSbc=[real(Sc1)^2+imag(Sc1)^2]/V(2)^2*Zbc;
Sbc=SSbc+Sc1;
SSef=[real(Sf)^2+imag(Sf)^2]/V(3)^2*Zef;
Sef=SSef+Sf;
SSeg=[real(Sg)^2+imag(Sg)^2]/V(4)^2*Zeg;
Seg=Sg+SSeg;
Se1=Sef+Seg+Se;
SSbe=[real(Se1)^2+imag(Se1)^2]/V(5)^2*Zbe;
Sbe=SSbe+Se1;
SSbh=[real(Sh)^2+imag(Sh)^2]/V(6)^2*Zbh;
Sbh=SSbh+Sh;
Sb1=Sbc+Sbe+Sbh+Sb;
SSab=[real(Sb1)^2+imag(Sb1)^2]/V(7)^2*Zab;
Sab=SSab+Sb1;
Va=10.5;
disp('[Scd Sbc Sef Seg Sbe Sbh Sab]');
S=[Scd Sbc Sef Seg Sbe Sbh Sab]
Vb=Afun(Sab,Zab,Va);Vc=Afun(Sbc,Zbc,Vb);
Vd=Afun(Scd,Zcd,Vc);Ve=Afun(Sbe,Zbe,Vb);
Vf=Afun(Sef,Zef,Ve);Vg=Afun(Seg,Zeg,Ve);
Vh=Afun(Sbh,Zbh,Vb);
disp('[Vd Vc Vf Vg Ve Vh Vb]');
V=[Vd Vc Vf Vg Ve Vh Vb]
E=V-B;
end
function [v] =Afun(m,n,w)%m为节点负荷Sab,n为线路之间阻抗Zab,w为线路首端电压Va Va1=[real(m)*real(n)+imag(m)*imag(n)]/w;%电压降落纵分量 Va2=[real(m)*imag(n)-imag(m)*real(n)]/w;%电压降落横分量 v=[(w-Va1)^2+Va2^2]^0.5;%线路末端电压 end