《电力系统分析潮流计算大作业.pdf》由会员分享，可在线阅读，更多相关《电力系统分析潮流计算大作业.pdf（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、电力系统分析潮流计算大作业电力系统分析潮流计算大作业（源程序及实验报告）源程序及实验报告）源程序如下：源程序如下：采用直角坐标系的牛顿-拉夫逊迭代function chaoliujisuan()m=3;%m=PQ 节点个数v=1;%v=PV 节点个数P=-0.8055-0.18 0;%P=PQ节点的 P 值Q=-0.5320-0.12 0;%Q=PQ节点的 Q 值PP=0.5;%PP=PV 节点的 P 值V=1.0;%V=PV 节点的 U 值E=1 1 1 1.0 1.0;%E=PQ,PV,V 节点 e 的初值F=0 0 0 0 0;%F=PQ,PV,V 节点 f 的初值G=6.3110-3.
2、5587-2.7523 0 0;-3.5587 8.5587-50 0;-2.7523-57.75230 0;0000 0;0000 0;B=-20.4022 11.38799.174300;11.3879-31.00937154.98890;9.174315-28.7757 04.9889;04.988905.24930;004.98890-5.2493;Y=G+j*B;X=;%X=Xn=m+v+1;%总的节点数FX=ones(2*n-2,1);%F(x)矩阵F1=zeros(n-1,n-1);%F(x)导数矩阵a=0;%记录迭代次数EF=zeros(n-1,n-1);%最后的节点电压矩阵w
3、hile max(FX)=10(-5)for i=1:m%PQ 节点 FX(i)=P(i);%P FX(n+i-1)=Q(i);%Qfor w=1:n FX(i)=FX(i)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w);%PFX(n+i-1)=FX(n+i-1)-F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w);%Qendendfor i=m+1:n-1%PV 节点 FX(i)=PP(i-m);%P FX(n+i-1)
4、=V(i-m)2-E(i)2-F(i)2;%Qfor w=1:n FX(i)=FX(i)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w);%Pendendfor i=1:m%PQ 节点for w=1:n-1if i=wF1(i,w)=-(G(i,w)*E(i)+B(i,w)*F(i);F1(i,n+w-1)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,w)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,n+w-1)=G(i,w)*E(i)+B(i,w)*F(i);e
5、lseF1(i,w)=-G(i,i)*E(i)-B(i,i)*F(i);F1(i,n+w-1)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,w)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,n+w-1)=G(i,i)*E(i)+B(i,i)*F(i);for k=1:nF1(i,w)=F1(i,w)-G(i,k)*E(k)+B(i,k)*F(k);F1(i,n+w-1)=F1(i,n+w-1)-G(i,k)*F(k)-B(i,k)*E(k);F1(n+i-1,w)=F1(n+i-1,w)+G(i,k)*F(k)+B(i,k)*E(k);F1(n+i-
6、1,n+w-1)=F1(n+i-1,n+w-1)-G(i,k)*E(k)+B(i,k)*F(k);endendendendfor i=m+1:n-1%PV节点for w=1:n-1if i=wF1(i,w)=-(G(i,w)*E(i)+B(i,w)*F(i);F1(i,n+w-1)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,w)=0;F1(n+i-1,n+w-1)=0;elseF1(i,w)=-G(i,i)*E(i)-B(i,i)*F(i);F1(i,n+w-1)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,w)=-2*E(i);F1(n+i-1,
7、n+w-1)=-2*F(i);for k=1:nF1(i,w)=F1(i,w)-G(i,k)*E(k)+B(i,k)*F(k);F1(i,n+w-1)=F1(i,n+w-1)-G(i,k)*F(k)-B(i,k)*E(k);endendendendX=inv(F1)*(-FX);for i=1:n-1E(i)=E(i)+X(i);F(i)=F(i)+X(n+i-1);enda=a+1;fprintf(第%d 次迭代后的节点电压分别为:n,a);disp(E+j*F);fprintf(第%d 次迭代后功率偏差P Q 电压偏差V 的平方分别为:n,a);disp(FX);enddisp(收敛后的节
8、点电压用极坐标表示为:);EF=E+j*F;for i=1:n-1 fprintf(%d 号节点电压的幅值为:,i)disp(abs(EF(i);fprintf(%d 号节点电压的相角度数为,i)disp(angle(EF(i)*180/pi);endPPH=0;for i=1:nPPH=PPH+EF(n)*conj(Y(n,i)*conj(EF(i);endfprintf(平衡节点的功率);disp(PPH);运行结果：运行结果：运行结果复制如下：运行结果复制如下：第 1 次迭代后的节点电压分别为:1.0034-0.1306i1.0285-0.1019i1.0339-0.0973i1.000
9、0-0.0017i1.0000第 1 次迭代后功率偏差P Q 电压偏差V 的平方分别为:-0.8055-0.180000.5000-0.37200.24740.38750第 2 次迭代后的节点电压分别为:0.9836-0.1290i1.0126-0.1038i1.0183-0.0995i1.0000-0.0035i1.0000第 2 次迭代后功率偏差P Q 电压偏差V 的平方分别为:0.0512-0.0222-0.04030.0002-0.1012-0.0219-0.0099-0.0000第 3 次迭代后的节点电压分别为:0.9831-0.1290i1.0122-0.1038i1.0180-0
10、.0996i1.0000-0.0035i1.0000第 3 次迭代后功率偏差P Q 电压偏差V 的平方分别为:0.0008-0.0003-0.0005-0.0001-0.0021-0.0004-0.0003-0.0000第 4 次迭代后的节点电压分别为:0.9831-0.1290i1.0122-0.1038i1.0180-0.0996i1.0000-0.0035i1.0000第 4 次迭代后功率偏差P Q 电压偏差V 的平方分别为:1.0e-005*0.0280-0.0083-0.0164-0.0121-0.1085-0.0199-0.0135-0.0005收敛后的节点电压用极坐标表示为:1 号节点电压的幅值为:0.99161 号节点电压的相角度数为-7.47482 号节点电压的幅值为:1.01752 号节点电压的相角度数为-5.85483 号节点电压的幅值为:1.02293 号节点电压的相角度数为-5.58644 号节点电压的幅值为:1.00004 号节点电压的相角度数为-0.2021平衡节点的功率0.4968-10.3280i