计算导纳矩阵.docx

1、计算导纳矩阵 /*计算导纳矩阵*G11=z12r/(z12r*z12r+z12m*z12m)+k*k*z13r/(z13r*z13r+z13m*z13m)+z14r/(z14r*z14r+z14m*z14m);B11=-z12m/(z12r*z12r+z12m*z12m)-k*k*z13m/(z13r*z13r+z13m*z13m)-z14m/(z14r*z14r+z14m*z14m)+y140+y120; G22=z12r/(z12r*z12r+z12m*z12m)+z24r/(z24r*z24r+z24m*z24m);B22=-z12m/(z12r*z12r+z12m*z12m)-z24m

2、/(z24r*z24r+z24m*z24m)+y240+y120;G33=z13r/(z13r*z13r+z13m*z13m);B33=-z13m/(z13r*z13r+z13m*z13m); G44=z14r/(z14r*z14r+z14m*z14m)+z24r/(z24r*z24r+z24m*z24m;B44=-z14m/(z14r*z14r+z14m*z14m)-z24m/(z24r*z24r+z24m*z24m)+y240+y140;G12=G21=-z12r/(z12r*z12r+z12m*z12m);B12=B21=z12m/(z12r*z12r+z12m*z12m);G13=G3

3、1=-k*z13r/(z13r*z13r+z13m*z13m);B13=B31=k*z13m/(z13r*z13r+z13m*z13m);G14=G41=-z14r/(z14r*z14r+z14m*z14m);B14=B41=z14m/(z14r*z14r+z14m*z14m);G23=G32=0.0;B23=B32=0.0; G24=G42=-z24r/(z24r*z24r+z24m*z24m);B24=B42=z24m/(z24r*z24r+z24m*z24m);G34=G43=0.0;B34=B43=0.0;for(i=1;i5;i+)for(j=1;j5;j+)printf(%f+%f

4、j,Gij,Bij);printf();printf(n);/形成节点导纳矩阵/计算各节点不平衡量loop1: printf(迭代次数k1=%dn,k1);for(i=1;i5;i+)floata=0,b=0;for(j=1;j5;j+)a+=Gij*ej-Bij*fj;b+=Gij*fj+Bij*ej; Pi=Psi-(ei*a+fi*b);/计算有功功率的增量Qi=Qsi-(fi*a-ei*b);/计算无功功率的增量 V32=V3s*V3s-e3*e3; printf(有功功率增量P1=%f,P1);printf(,); printf(有功功率增量P2=%f,P2);printf(,);

5、printf(有功功率增量P3=%f,P3); printf(无功功率增量Q1=%f,Q1);printf(,); printf(无功功率增量Q2=%f,Q2);printf(,); printf(电压增量V32=%f,V32); printf(n);/*形成雅克比矩阵*for(j=1;j4;j+)if(1=j)floatc=0,d=0;intm;for(m=1;m5;m+)c+=G1m*em-B1m*fm;d+=G1m*fm+B1m*em;J1*N-1j*N-1=-c-G1j*e1-B1j*f1;J1*N-1j*N=-d+B1j*e1-G1j*f1;J1*Nj*N-1=d+B1j*e1-G1

6、j*f1;J1*Nj*N=-c+G1j*e1+B1j*f1;else J1*N-1j*N-1=-G1j*e1-B1j*f1;J1*Nj*N=G1j*e1-B1j*f1;J1*N-1j*N=B1j*e1-G1j*f1;J1*Nj*N-1=B1j*e1-G1j*f1; /*计算修正方程*for(i=1;iM;i+)Lii=1;Li1=Ji1/U11; for(n=2;nM;n+) for(j=n;jM;j+)sigma1=0;for(s=0;s=n-1;s+)sigma1+=Lns*Usj;Unj=Jnj-sigma1; for(i=n;iM;i+)sigma2=0;for(s=0;s=n-1;s

7、+)sigma2+=Lis*Usn;Lin=(Jin-sigma2)/Unn; b1=P1;b2=Q1;b3=P2;b4=Q2;b5=P3;b6=V32;for(i=1;iM;i+)sigma1=0;for(n=1;n=i-1;n+)for(i=1;i=1;i-)sigma2=0; for(n=i+1;nM;n+)sigma2+=Uin*xn;xi=(yi-sigma2)/Uii;xe1=-x1;xe2=-x3;xe3=-x5;xf1=-x2;xf2=-x4;xf3=-x6;printf(节点电压:n);for(i=1;i4;i+)ei+=xei;fi+=xfi; for(i=1;i4;i+)

8、printf(e%d=,i);printf(%f,ei);printf(,);for(i=1;i4;i+)printf(f%d=,i);printf(%f,fi);printf(,); %给定节点电压初值e = 1 1 1.1 1.05;f = 0 0 0 0 ;%给定PQ节点有用和无功功率，以及PV节点有功和电压P1 = -0.3;Q1 = -0.18;P2 = -0.55;Q2 = -0.13;P3 = 0.5;V3 = 1.1;%手动输入节点导纳矩阵G= 1.042093, -0.588235, 0, -0.453858; -0.588235, 1.069005, 0, -0.48076

9、9; 0 , 0 , 0, 0; -0.453858, -0.480769, 0, 0.934627 ;B= -8.242876, 2.352941, 3.666667, 1.891074; 2.352941, -4.727377, 0, 2.403846; 3.666667, 0 ,-3.333333, 0; 1.891074, 2.403846, 0, -4.261590 ; maxP = 0;maxQ = 0;maxV = 0; I = 0,0; 0,0; 0,0; 0,0 ;for v=1:15 %迭代次数for n=1:4 I(1,1)=I(1,1)+G(1,n)*e(n)-B(1,

10、n)*f(n); I(1,2)=I(1,2)+G(1,n)*f(n)+B(1,n)*e(n);endfor n=1:4 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);endfor n=1:4 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);endfor n=1:4 I(4,1)=I(4,1)+e(n)*e(n)+f(n)*f(n); I(4,2)=0;end H =;N =;M =;

11、L =;R =;S =;J =; %求不平衡量P1 = -0.30-e(1)*I(1,1)-f(1)*I(1,2);Q1 = -0.18-f(1)*I(1,1)+e(1)*I(1,2);P2 = -0.55-e(2)*I(2,1)-f(2)*I(2,2);Q2 = -0.13-f(2)*I(2,1)+e(2)*I(2,2);P3 = 0.50-e(3)*I(3,1)-f(3)*I(3,2);V3 = 1.102-I(4,1); %分块计算雅各比矩阵元素for m=1:3 for n=1:3 if(m=n) H(m,m) = B(m,m)*e(m)-G(m,m)*f(m)-I(m,2); N(m

12、,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); R(m,m) = -2*f(m); S(m,m )= -2*e(m); 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); R(m,m) = 0; S(m,m )= 0; end endend%确定雅克比矩阵各元素J= H(