发布时间 : 星期六 文章题目19：基于MATLAB的电力系统复杂潮流计算更新完毕开始阅读
雅克比矩阵元素J(k)i按式(7-21)和(7-22)计算
增大节点号i=i+10）0）P（=0.0；Q（5=0.0 5(k)(k)解修正方程式,由PiQi和J(k)i(k)根据式（2-3）得各节点功率的不平衡量为： (k)计算电压修正量Vi和i0）（0）?Q?P（=0.5000；2=1.1000 20）0）求出?Q（?P（=-0.3750；Vi=0.0750 ,33i0）0）?P（=-0.4000；?Q（4=-0.0500 4是k=k+1(k)Vi(k+1)=Vi+(k+1)i(k)Vi(k)i=(k)i+(k)(k)=-0.6000；?Q=-0.1000 I??0.3000?jP0.9000?jQ根据公式：I??aI??0.0750?j0.2250Ui (0)I4?03.3牛拉法潮流计算的流程图 (0)I5?0结束3.4 MATLAB算法的计算程序 %开始 clc clear
disp('节点总数为：'); N=5
disp('平衡节点为：'); 1
disp('PQ节点为：'); JD=[2,3,4,5] e=[1.06 1 1 1 1]; f=[0 0 0 0 0]; P1=0; Q1=0; P2=-0.2; Q2=-0.2; P3=0.45; Q3=0.15; P4=0.4; Q4=0.05; P5=0.6;
?P
Q5=0.1;
G=[6.2500,-5.0000,-1.2500,0,0;-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.75,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];%形成电纳矩阵。 disp('节点电导矩阵G为:'); disp(G)
disp('节点电纳矩阵B为:'); disp(B) k=0; for v=1:7
I=[0,0;0,0;0,0;0,0;0,0]; for n=1:5
I(1,1)=I(1,1)+G(1,n)*e(n)-B(1,n)*f(n); I(1,2)=I(1,2)+G(1,n)*f(n)+B(1,n)*e(n); end for n=1:5
I(2,1)=I(2,1)+G(2,n)*e(n)-B(2,n)*f(n); I(2,2)=I(2,2)+G(2,n)*f(n)+B(2,n)*e(n); end for n=1:5
I(3,1)=I(3,1)+G(3,n)*e(n)-B(3,n)*f(n); I(3,2)=I(3,2)+G(3,n)*f(n)+B(3,n)*e(n); end for n=1:5
I(4,1)=I(4,1)+G(4,n)*e(n)-B(4,n)*f(n); I(4,2)=I(4,2)+G(4,n)*f(n)+B(4,n)*e(n); end for n=1:5
I(5,1)=I(5,1)+G(5,n)*e(n)-B(5,n)*f(n); I(5,2)=I(5,2)+G(5,n)*f(n)+B(5,n)*e(n); end H=[]; N=[]; M=[]; L=[]; 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),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);