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

前推回代法潮流计算如图所示一个5节点的配电网系统，S1=S2=0.2+j0.2(MV·A)，S3=S4=S5=0.4+j0.4(MV·A)，假定所有负荷均为恒功率负荷，节点1为参考节点相角为0o。计算电网个节点电压及支路功率分布，收敛条件为610。M123451+j1Ω1+j2Ω2+j2Ω2+j1Ω10kVS1S2S4S5S3程序设计如下：clc;z=[01+1i000;1+1i01+2i2+2i0;01+2i000;02+2i002+1i;0002+1i0];%各节点之间的阻抗s=[0.2+0.2i0000;00.2+0.2i000;000.4+0.4i00;0000.4+0.4i0;00000.4+0.4i];%上三角首端功率，下三角末端功率v1=[1010101010];v2=[100000];%节点电压v=v2-v1;n=length(z);m=0;%循环次数whilemax(abs(v))=10^-6fori=n:-1:2forj=4:-1:1if(ij&&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首端功率endendendfori=2:nforj=1:nif(ij&&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末端电压endendendv=v2-v1;m=m+1;v1=v2;endmsv2计算结果