电分前推回代法潮流计算matlab程序设计

前推回代法潮流计算

如图所示一个5节点的配电网系统，S1=S2=0.2+j0.2(MV·A)，S3=S4=S5=0.4+j0.4(MV·A)，假定所有负荷均为恒功率负荷，节点1为参考节点相角为0o 。计算电网个节点电压及支路功率分布，收敛条件为610ε-<。

3

S4

程序设计如下：

clc;

z=[ 0 1+1i 0 0 0;

1+1i 0 1+2i 2+2i 0;

0 1+2i 0 0 0;

0 2+2i 0 0 2+1i;

0 0 0 2+1i 0 ];%各节点之间的阻抗

s=[ 0.2+0.2i 0 0 0 0 ;

0 0.2+0.2i 0 0 0 ;

0 0 0.4+0.4i 0 0 ;

0 0 0 0.4+0.4i 0 ;

0 0 0 0 0.4+0.4i ];

%上三角首端功率，下三角末端功率

v1=[10 10 10 10 10 ];

v2=[10 0 0 0 0 ];%节点电压

v=v2-v1;

n=length(z);

m=0;%循环次数

while max(abs(v))>=10^-6

for i=n:-1:2

for j=4:-1:1

if (i>j&&z(i,j)~=0)

s(i,j)=sum(s(:,i))-s(i,j);%计算ij末端功率

s(j,i)=s(i,j)+((real(s(i,j))).^2+...

(imag(s(i,j))).^2)./(v1(i).^2).*z(i,j);

%计算ji首端功率

end

end

end

for i=2:n

for j=1:n

if (i>j&&s(i,j)~=0)

v2(i)=sqrt((v2(j)-((real(s(j,i)).*real(z(i,j))+... imag(s(j,i)).*imag(z(i,j)))./v2(j))).^2+...

((real(s(j,i)).*imag(z(i,j))-...

imag(s(j,i)).*real(z(i,j)))./v2(j)).^2);

%计算ij末端电压

end

end

end

v=v2-v1;

m=m+1;

v1=v2;

end

m

s

v2

计算结果