計(jì)算導(dǎo)納矩陣

1、 /*計(jì)算導(dǎo)納矩陣*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

2、)-z24m/(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=B

3、21=z12m/(z12r*z12r+z12m*z12m); G13=G31=-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/(z

4、24r*z24r+z24m*z24m); G34=G43=0.0; B34=B43=0.0;for(i=1;i<5;i+) for(j=1;j<5;j+) printf("%f+%fj",Gij,Bij);printf(" "); printf("n");/形成節(jié)點(diǎn)導(dǎo)納矩陣/計(jì)算各節(jié)點(diǎn)不平衡量loop1:  printf("迭代次數(shù)k1=%dn",k1); for (i=1;i<5;i+) floa

5、t a=0,b=0; for(j=1;j<5;j+) a+=Gij*ej-Bij*fj;  b+=Gij*fj+Bij*ej;    Pi=Psi-(ei*a+fi*b);/計(jì)算有功功率的增量 Qi=Qsi-(fi*a-ei*b);/計(jì)算無(wú)功功率的增量   V32=V3s*V3s-e3*e3;    printf("有功功率增量P1=%f",P1);  printf(" ,");  prin

6、tf("有功功率增量P2=%f",P2);  printf(" ,");  printf("有功功率增量P3=%f",P3);   printf("無(wú)功功率增量Q1=%f",Q1);  printf(" ,");  printf("無(wú)功功率增量Q2=%f",Q2);  printf(" ,");  printf(&qu

7、ot;電壓增量V32=%f",V32);   printf("n");/*形成雅克比矩陣* for(j=1;j<4;j+) if(1=j) float c=0,d=0; int m; for(m=1;m<5;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;&

8、#160;J1*Nj*N-1=d+B1j*e1-G1j*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;      /*計(jì)算修正方程*for(i=1;i<M;i+) Lii=1; Li1=Ji1/U11; 

9、;  for(n=2;n<M;n+)    for(j=n;j<M;j+)  sigma1=0; for(s=0;s<=n-1;s+) sigma1+=Lns*Usj;Unj=Jnj-sigma1;    for(i=n;i<M;i+)  sigma2=0; for(s=0;s<=n-1;s+) sigma2+=Lis*Usn; Lin=(Jin-sigma2)/Unn;   

10、 b1=P1;b2=Q1;b3=P2; b4=Q2;b5=P3;b6=V32;for(i=1;i<M;i+)  sigma1=0; for(n=1;n<=i-1;n+)for(i=1;i<M;i+)  U1i=J1i; sigma1+=Lin*yn; yi=bi-sigma1;   for(i=M-1;i>=1;i-)  sigma2=0;  for(n=i+1;n<M;n+) sigma2+=Uin*xn; xi=

11、(yi-sigma2)/Uii; xe1=-x1;xe2=-x3;xe3=-x5; xf1=-x2;xf2=-x4;xf3=-x6; printf("節(jié)點(diǎn)電壓:n"); for(i=1;i<4;i+) ei+=xei;  fi+=xfi;   for(i=1;i<4;i+) printf("e%d=",i);  printf("%f",ei);  printf(" 

12、,"); for(i=1;i<4;i+) printf("f%d=",i);  printf("%f",fi);  printf(" ,");  %給定節(jié)點(diǎn)電壓初值e = 1 1 1.1 1.05;f = 0 0 0 0 ;%給定PQ節(jié)點(diǎn)有用和無(wú)功功率，以及PV節(jié)點(diǎn)有功和電壓P1 = -0.3;Q1 = -0.18;P2 = -0.55;Q2 = -0.13;P3 = 0.5;V3 = 1.1;%手動(dòng)輸入節(jié)點(diǎn)導(dǎo)納矩陣G= 1.042093,

13、 -0.588235, 0, -0.453858; -0.588235, 1.069005, 0, -0.480769; 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

14、 ;for v=1:15 %迭代次數(shù)for n=1:4 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);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(

15、4,1)=I(4,1)+e(n)*e(n)+f(n)*f(n); I(4,2)=0;end H =;N =;M =;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); %分塊計(jì)算雅各比矩陣元素for m=1:3 for

16、n=1:3 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); 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

17、(m,n); L(m,n) = H(m,n); R(m,m) = 0; S(m,m )= 0; end endend%確定雅克比矩陣各元素J= H(1,1),N(1,1),H(1,2),N(1,2),H(1,3),N(1,3); M(1,1),L(1,1),M(1,2),L(1,2),M(1,3),L(1,3); H(2,1),N(2,1),H(2,2),N(2,2),H(2,3),N(2,3); M(2,1),L(2,1),M(2,2),L(2,2),M(2,3),L(2,3); H(3,1),N(3,1),H(3,2),N(3,2),H(3,3),N(3,3); R(3,1),S(3,1)

18、,R(3,2),S(3,2),R(3,3),S(3,3) ;%求解節(jié)點(diǎn)電壓的修正量A =;C = P1; Q1; P2; Q2; P3; V3 A = -inv(J)*C; %排序找出最大誤差 maxP=C(1); maxQ=C(2); maxV=C(6); %V3直接賦給電壓誤差最大值maxV for d = 1:2:5 %在P1/P3/P5中找出最大值賦給maxP if(C(d)>maxP) maxP=C(d); end end for t = 2:2:4 %在Q2/Q4中找出最大值賦給maxQ if(C(t)>maxQ) maxQ=C(t); end end %修正節(jié)點(diǎn)電壓值

19、 if(maxP>0.00001)&&(maxQ>0.00001)&&(maxV>0.00001) e(1)= e(1)+A(2) e(2)= e(2)+A(4) e(3)= e(3)+A(6) f(1)= f(1)+A(1) f(2)= f(2)+A(3) f(3)= f(3)+A(5) endend%計(jì)算平衡節(jié)點(diǎn)功率和網(wǎng)絡(luò)中的功率分布S4 = (1.05+0*i)*(G(4,1)-i*B(4,1)*(e(1)-i*f(1) +(G(4,2)-i*B(4,2)*(e(2)-i*f(2) +(G(4,3)-i*B(4,3)*(e(3)-i*f(

20、3) +(G(4,4)-i*B(4,4)*(1.05-0*i) %平衡節(jié)點(diǎn)功率 S12 = (e(1)+i*f(1)*e(1)+i*f(1)-e(2)-f(2)*i*(G(1,2)+i*B(1,2)S13 = (e(1)+i*f(1)*e(1)+i*f(1)-e(3)-f(3)*i*(G(1,3)+i*B(1,3)S14 = (e(1)+i*f(1)*e(1)+i*f(1)-1.05-0*i*(G(1,4)+i*B(1,4)S21 = (e(2)+i*f(2)*e(2)+i*f(2)-e(1)-i*f(1)*(G(2,1)+i*B(2,1)S24 = (e(2)+i*f(2)*e(2)+i*f

21、(2)-1.05-0*i*(G(2,4)+i*B(2,4)S31 = (e(3)+i*f(3)*e(3)+i*f(3)-e(1)-f(1)*i*(G(3,1)+i*B(3,1)S41 = (1.05+0*i)*1.05+0*i-e(1)-f(1)*i*(G(4,3)+i*B(4,3)S42 = (1.05+0*i)*1.05+0*i-e(2)-f(2)*i*(G(4,2)+i*B(4,2)課程設(shè)計(jì)任務(wù)書學(xué)生姓名： 專業(yè)班級(jí)： 指導(dǎo)教師： 工作單位： 題 目: 直角坐標(biāo)下牛頓法潮流計(jì)算 初始條件：下圖所示的簡(jiǎn)單電力系統(tǒng)中，網(wǎng)絡(luò)各元件參數(shù)的標(biāo)幺值如下：z12=0.10+j0.40; y120=y210=j0.01528;z13=j0.3; k=1.1; z14=0.12+j0.50; y140=y410=j0.01920; z24=0.08+j0.40; y240=y420=j0.01413。系統(tǒng)中節(jié)點(diǎn)1、2為PQ節(jié)點(diǎn)，節(jié)點(diǎn)3為PV節(jié)點(diǎn)，節(jié)點(diǎn)4為平衡節(jié)點(diǎn)，已給定P1s+jQ1s=-0.30-j0.18; P2s+jQ2s=-0.55-j0.13; P3s=0.5; V3s=1.10; V4s=1.05。容許誤差。節(jié)點(diǎn)電壓用直角坐