发布时间 : 星期四 文章(强烈推荐)两机五节点网络潮流计算方法牛拉法和pq法_电力系统稳态分析毕业论文更新完毕开始阅读
else
H(m,n)=-B(m,n)*E(m)+G(m,n)*F(m); N(m,n)=G(m,n)*E(m)+B(m,n)*F(m); J(m,n)=-B(m,n)*F(m)-G(m,n)*E(m); L(m,n)=G(m,n)*F(m)-B(m,n)*E(m);
JJ(2*m-1,2*n-1)=H(m,n); JJ(2*m-1,2*n)=N(m,n); JJ(2*m,2*n-1)=J(m,n); JJ(2*m,2*n)=L(m,n); end end end JJ
%设由P,Q的改变量组成的8×1矩阵为PQ,由E,F的改变量组成的8×1矩阵为dU for m=1:4
PQ(2*m-1)=dP(m); PQ(2*m)=dQ(m); end
dU=inv(JJ)*PQ';%求逆矩阵
C=max(abs(dU));?s绝对值或复数绝对值 for n=1:4
F(n)=F(n)+dU(2*n-1); E(n)=E(n)+dU(2*n);
end for n=1:5
U(n)=E(n)+(F(n))*i; end k=k+1;
k, dU=dU',PQ,U end
%计算S(5),也就是题目中的S1,即平衡节点功率 for m=1:5
I(m)=Y(5,m)*U(m); end
S(5)=U(5)*sum(conj(I)) Y =
3.7500 -11.2500i -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 0
-2.5000 + 7.5000i 10.8333 -32.5000i -1.6667 + 5.0000i -1.6667 + 5.0000i -5.0000 +15.0000i
0 -1.6667 + 5.0000i 12.9167 -38.7500i -10.0000 +30.0000i -1.2500 + 3.7500i
-1.2500 + 3.7500i -1.6667 + 5.0000i -10.0000 +30.0000i 12.9167 -38.7500i 0
0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 6.2500 -18.7500i
S =
-0.6000 - 0.1000i 0.2000 + 0.2000i -0.4500 - 0.1500i -0.4000 - 0.0500i 0 JJ =
11.2500 3.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500
-3.7500 11.2500 2.5000 -7.5000 0 0 1.2500 -3.7500
-7.5000 -2.5000 33.4000 10.5333 -5.0000 -1.6667 -5.0000 -1.6667
2.5000 -7.5000 -11.1333 31.6000 1.6667 -5.0000 1.6667 -5.0000
0 0 -5.0000 -1.6667 38.9750 12.8417 -30.0000 -10.0000
0 0 1.6667 -5.0000 -12.9917 38.5250 10.0000 -30.0000
-3.7500 -1.2500 -5.0000 -1.6667 -30.0000 -10.0000 38.7500 12.9167
1.2500 -3.7500 1.6667 -5.0000 10.0000 -30.0000 -12.9167 38.7500 k = 1
dU =
-0.1076 0.0093 -0.0473 0.0430 -0.0863 0.0154 -0.0922 0.0141 PQ =
-0.6000 -0.1000 0.5000 1.1000 -0.3750 0.0750 -0.4000 -0.0500 U =
1.0093 - 0.1076i 1.0430 - 0.0473i 1.0154 - 0.0863i 1.0141 - 0.0922i 1.0600 JJ =
11.0516 4.3956 -7.3011 -3.3304 0 0 -3.6505 -1.6652
-5.5956 10.8516 3.3304 -7.3011 0 0 1.6652 -3.6505
-7.7040 -2.9621 33.1594 13.0913 -5.1360 -1.9747 -5.1360 -1.9747
2.9621 -7.7040 -12.5804 33.6083 1.9747 -5.1360 1.9747 -5.1360
0 0 -4.9331 -2.1238 38.3848 16.0299 -29.5988 -12.7427
0 0 2.1238 -4.9331 -16.8888 38.0788 12.7427 -29.5988
-3.6876 -1.6135 -4.9168 -2.1513 -29.5009 -12.9079 38.1553 16.2727
1.6135 -3.6876 2.1513 -4.9168 12.9079 -29.5009