安徽工程大学课程设计说明书
J=[];
P2=P2-e(2)*I(2,1)-f(2)*I(2,2); %有功功率的不平衡量 Q2=Q2-f(2)*I(2,1)+e(2)*I(2,2); %无功功率的不平衡量 P3=P3-e(3)*I(3,1)-f(3)*I(3,2); Q3=Q3-f(3)*I(3,1)+e(3)*I(3,2); P4=P4-e(4)*I(4,1)-f(4)*I(4,2); Q4=Q4-f(4)*I(4,1)+e(4)*I(4,2); P5=P5-e(5)*I(5,1)-f(5)*I(5,2); Q5=Q5-f(5)*I(5,1)+e(5)*I(5,2); for m=2:5 for n=2:5 if(m==n)
H(m,m)=-B(m,m)*e(m)+G(m,m)*f(m)+I(m,2); N(m,m)=G(m,m)*e(m)+B(m,m)*f(m)+I(m,1); M(m,m)=-G(m,m)*e(m)-B(m,m)*f(m)+I(m,1); L(m,m)=-B(m,m)*e(m)+G(m,m)*f(m)-I(m,2); 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); M(m,n)=-N(m,n); L(m,n)=H(m,n); end end end
J=[H(2,2),N(2,2),H(2,3),N(2,3),H(2,4),N(2,4),H(2,5),N(2,5);M(2,2),L(2,2),M(2,3),L(2,3),M(2,4),L(2,4),M(2,5),L(2,5);H(3,2),N(3,2),H(3,3),N(3,3),H(3,4),N(3,4),H(3,5),N(3,5);M(3,2),L(3,2),M(3,3),L(3,3),M(3,4),L(3,4),M(3,5),L(3,5);H(4,2),N(4,2),H(4,3),N(4,3),H(4,4),N(4,4),H(4,5),N(4,5);M(4,2),
21
安徽工程大学课程设计说明书
L(4,2),M(4,3),L(4,3),M(4,4),L(4,4),M(4,5),L(4,5);H(5,2),N(5,2),H(5,3),N(5,3),H(5,4),N(5,4),H(5,5),N(5,5);M(5,2),L(5,2),M(5,3),L(5,3),M(5,4),L(5,4),M(5,5),L(5,5)]; disp('雅克比矩阵J:'); disp(J); A=[];
C=[P2;Q2;P3;Q3;P4;Q4;P5;Q5] A=J\\C;%解修正方程式 disp('第M次修正方程的解A:'); disp(A);
f(2)=f(2) +A(1,1);
e(2)=e(2) +A(2,1); %计算新值 f(3)=f(3) +A(3,1); e(3)=e(3) +A(4,1); f(4)=f(4) +A(5,1); e(4)=e(4) +A(6,1); f(5)=f(5) +A(7,1); e(5)=e(5) +A(8,1);
disp('各点的电压实部e(单位:V)为(节点号从小到大排列):'); disp(e)
disp('各点的电压虚部f单位:V)为(节点号从小到大排列):'); disp(f); u=e+f*i;
disp('节点电压的第C(k)次近似值:'); disp(u); k=k+1;
22
安徽工程大学课程设计说明书
disp('迭代次数:'); disp(k); end for m=1:5
I(m)=(G(1,m)+B(1,m)*i) *u(m); end
disp('平衡节点的功率');
S1=u(1)*sum(conj(I))%计算平衡节点的功率 for m=1:5 for n=1:5
S(m,n)=u(m)*(conj(u(m))-conj(u(n)))*conj(-(G(m,n)+B(m,n)*i));%计算各支路功率 end end
disp('各支路功率');disp(S) %结束
3.5 MATLAB的计算结果
节点总数为： N = 5 平衡节点为： ans = 1 PQ节点为： JD =
2 3 4 5 节点电导矩阵G为:
6.2500 -5.0000 -1.2500 0 0
23
安徽工程大学课程设计说明书
-5.0000 10.8340 -1.6670 -1.6670 -2.5000 -1.2500 -1.6670 12.9170 -10.0000 0 0 -1.6670 -10.0000 12.9170 -1.2500 0 -2.5000 0 -1.2500 3.7500 节点电纳矩阵B为:
-18.7500 15.0000 3.7500 0 0 15.0000 -32.5000 5.0000 5.0000 7.5000 3.7500 5.0000 -38.7500 30.0000 0 0 5.0000 30.0000 -38.7500 3.7500 0 7.5000 0 3.7500 -11.2500 雅克比矩阵J:
Columns 1 through 6
33.4000 10.5340 -5.0000 -1.6670 -5.0000 -1.6670 -11.1340 31.6000 1.6670 -5.0000 1.6670 -5.0000 -5.0000 -1.6670 38.9750 12.8420 -30.0000 -10.0000 1.6670 -5.0000 -12.9920 38.5250 10.0000 -30.0000 -5.0000 -1.6670 -30.0000 -10.0000 38.7500 12.9170 1.6670 -5.0000 10.0000 -30.0000 -12.9170 38.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500 2.5000 -7.5000 0 0 1.2500 -3.7500 Columns 7 through 8 -7.5000 -2.5000 2.5000 -7.5000 0 0 0 0 -3.7500 -1.2500 1.2500 -3.7500 11.2500 3.7500
24
安徽工程大学课程设计说明书
-3.7500 11.2500 C = 0.1000 0.7000 0.5250 0.3750 0.4000 0.0500 0.6000 0.1000
第M次修正方程的解A: 0.0473 0.0847 0.0863 0.1123 0.0922 0.1136 0.1076 0.1183
各点的电压实部e(单位:V)为(节点号从小到大排列): 1.0600 1.0847 1.1123 1.1136 1.1183 各点的电压虚部f单位:V)为(节点号从小到大排列): 0 0.0473 0.0863 0.0922 0.1076 节点电压的第C(k)次近似值: Columns 1 through 4
1.0600 1.0847 + 0.0473i 1.1123 + 0.0863i Column 5 1.1183 + 0.1076i
25
1.1136 + 0.0922i