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

《电分前推回代法潮流计算matlab程序设计.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、前推回代法潮流计算如图所示一个5节点的配电网系统，S1=S2=0.2+j0.2(MV·A)，S3=S4=S5=0.4+j0.4(MV·A)，假定所有负荷均为恒功率负荷，节点1为参考节点相角为0o。计算电网个节点电压及支路功率分布，收敛条件为。程序设计如下：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];%上三角首端功率，下三角末端功率

2、v1=[1010101010];v2=[100000];%节点电压v=v2-v1;n=length(z);m=0;%循环次数whilemax(abs(v))>=10^-6fori=n:-1:2forj=4:-1:1if(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首端功率endendendfori=2:nforj=1:nif(i>j&&s(i,j)~

3、=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计算结果