L NEWPeiWangChaoLiu配电网潮流计算matlab程序_matlab

%形成邻接矩阵 data=xlsread('E:33netdata.xlsx'); a=zeros(33,33); for j=1:32; a(data(j,2)+1,data(j,3)+1)=1; end %%电压赋初值 v=ones(33,1); v(1)=1.05; delt=zeros(33,1); V=v.*cos(delt)+i*v.*sin(delt); %%节点负荷功率 SL=zeros(33,1); for k=2:33 SL(k)=(data(k-1,6)+i*data(k-1,7))/10; end SL(1)=0; %%线路阻抗 Z=zeros(33,1); for m=2:33 Z(m)=(data(m-1,4)+i*data(m-1,5))/10; end Z(1)=0; S=zeros(33,1); %迭代 for dd=1:100 V=v.*cos(delt)+i*v.*sin(delt); U=V; u=abs(U) for j=1:33 if a(j,:)==0 S(j)=SL(j)+conj(SL(j))*SL(j)/(conj(V(j))*V(j))*Z(j); end end for k=33:-1:2 if any(a(k,:)) for j=k:33 if a(k,j)==1 S(k)=S(k)+S(j); end end S(k)=SL(k)+S(k); S(k)=S(k)+Z(k)*conj(S(k))*S(k)/(conj(V(k))*V(k)); end end for n=1:33; for j=n:33 if a(n,j)==1 V(j)=V(n)-(conj(S(j))/V(n))*Z(j); end end end v=abs(V) if max(abs(v-u))<0.0000001 break end %v=abs(V); delt=angle(V); S=zeros(33,1); end