潮流计算的运行结果.doc

系统接线图：程序：%开始clccleardisp(节点总数为：);N=5disp(平衡节点为：);1disp(PQ节点为：);JD=2,3,4,5Y=6.25-18.75i -5+15i -1.25+3.75i 0 0; -5+15i 10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i; -1.25+3.75i -1.667+5i 12.917-38.75i -10+30i 0; 0 -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i;0 -2.5+7.5i 0 -1.25+3.75i 3.75-11.25i disp(导纳矩阵为:Y=);%求出导纳矩阵ei=1.06;1;1;1;1;fi=0;0;0;0;0;P=0;,-0.2;0.45;0.40;0.60;Q=0;-0.20;0.15;0.05;0.10;disp(Y);wucha=0.0001;k=0;G=real(Y);B=imag(Y);N1=4;while wucha0.00001 for m=2:5 pp=0;qq=0; for n=1:5 pp=pp+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); % Pi（0） qq=qq+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)-ei(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); I(m)= (pp-qq*i)/conj(ei(m)+fi(m)*i);%节点注入电流Ii（0） end DP(m)=P(m)-pp; %有功功率的不平衡量 DQ(m)=Q(m)-qq; %无功功率的不平衡量 H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+imag(I(m); N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+real(I(m); J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+real(I(m); L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)-imag(I(m);%公式P131 JJ(2*m-3,2*m-3)=H(m,m);JJ(2*m-3,2*m-2)=N(m,m);JJ(2*m-2,2*m-3)=J(m,m);JJ(2*m-2,2*m-2)=L(m,m); end for m=2:N for n=2:N if n=m else H(m,n)=-B(m,n)*ei(m)+G(m,n)*fi(m); N(m,n)=G(m,n)*ei(m)+B(m,n)*fi(m); J(m,n)=-B(m,n)*fi(m)-G(m,n)*ei(m); L(m,n)=G(m,n)*fi(m)-B(m,n)*ei(m);JJ(2*m-3,2*n-3)=H(m,n);JJ(2*m-3,2*n-2)=N(m,n);JJ(2*m-2,2*n-3)=J(m,n);JJ(2*m-2,2*n-2)=L(m,n);%公式P130 end end end disp(雅克比矩阵JJ:); disp(JJ); for m=2:N DW(2*m-3)=DP(m);DW(2*m-2)=DQ(m);%形成DP和DQ不平衡量的矩阵 end DY=JJDW;%解修正方程式 wucha=max(abs(DY); disp(第M次修正方程的解DY:); disp(DY); for n=2:N fi(n)=fi(n) +DY(2*n-3); ei(n)=ei(n) +DY(2*n-2); %计算各节点电压新值 end u=ei+fi*i; disp(节点电压的第C(k)次近似值:); disp(u); disp(各点的电压实部ei(单位:V)为(节点号从小到大排列):); disp(ei); disp(各点的电压虚部fi(单位:V)为(节点号从小到大排列):); disp(fi); k=k+1; disp(迭代次数:); disp(k);endfor n=1:N1+1 I(n)=Y(1,n)*u(n);%(公式P139)enddisp(平衡节点的功率);S1=u(1)*sum(conj(I)%计算平衡节点的功率for m=1:N1+1 for n=1:N1+1 S(m,n)=u(m)*(conj(u(m)-conj(u(n)*conj(-Y(m,n);%计算m-n各支路功率 S2(n,m)=u(n)*(conj(u(n)-conj(u(m)*conj(-Y(n,m);%计算n-m各支路功率 DS(m,n)=S(m,n)+S2(n,m);%计算各支路功