本文提出了利用 ＭＡＴＬＡＢ 语言来进行电力系统潮流计算的方法。通过实例， 说明了该 方法编程简便、 运算效率高并符合人们的思维习惯.

 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