潮流計算matlab-牛頓拉夫遜法法(共17頁)

1、精選優(yōu)質文檔-傾情為你奉上用matlab潮流計算(牛頓拉夫遜法)%主程序：dfile,pathname=uigetfile('ieee14.m','Select Data File');if pathname=0 error('you must select a valid data file')else lfile=length(dfile); % strip off .m eval(dfile(1:lfile-2); %打開數(shù)據(jù)文件endglobal n;global m;nb,mb=size(bus); %節(jié)點重新編號nl,ml=size(

2、line); nSW=0; % 平衡節(jié)點數(shù)目 nPV=0; % PV節(jié)點數(shù)目 nPQ=0; % PQ節(jié)點數(shù)目 for i=1:nb, % nb為總節(jié)點數(shù) type=bus(i,6); if type=3, nSW=nSW+1; % 統(tǒng)計平衡節(jié)點數(shù)目 SW(nSW,:)=bus(i,:); elseif type=2, nPV=nPV+1; % 統(tǒng)計PV節(jié)點數(shù)目 PV(nPV,:)=bus(i,:); else nPQ=nPQ+1; % 統(tǒng)計PQ節(jié)點數(shù)目 PQ(nPQ,:)=bus(i,:); end end; bus=PQ;PV;SW; newbus=1:nb' f=bus(:,1);

3、 nodenum=newbus bus(:,1); bus(:,1)=newbus; for i=1:nl for j=1:2 for k=1:nb if line(i,j)=nodenum(k,2) line(i,j)=nodenum(k,1); break end end end endY=y(bus,line); %形成節(jié)點導納矩陣K=0; %迭代次數(shù)初值Kmax=10; %最大迭代次數(shù)eps1=1.0e-10;eps2=1.0e-10;m=nPQ;n=nb;Um=eye(m,m);myf=fopen('output1.dat','w');for K=1:

4、Kmax for i=1:m for j=1:m if i=j Um(i,j)=bus(i,2); end end end b=dPQ(Y,bus); C=jac(bus,Y); dX=Cb' dx=dX' nx,mx=size(dx); for i=1:n-1 %計算相角 bus(i,3)=bus(i,3)-dX(i,1); end B=dx(nx,n:mx)*Um; %計算電壓差 bus(1:m,2)=bus(1:m,2)-B' %計算電壓值 dx(nx,n:mx)=B; fprintf(myf,'-第%d次迭代時雅可比矩陣-',K) fprintf

5、(myf, 'n'); for i=1:(n+m-1) for j=1:(n+m-1) fprintf(myf,'%8.6f', C(i,j); fprintf(myf, ' '); end fprintf(myf, 'n'); end fprintf(myf,'-第%d次迭代時dPQ的誤差-',K) fprintf(myf, 'n'); for i=1:(n+m-1) fprintf(myf,'%8.6e', b(1,i); fprintf(myf, 'n'); e

6、nd fprintf(myf, 'n'); fprintf(myf,'-第%d次迭代時dx(誤差)-',K) fprintf(myf, 'n'); for i=1:(n+m-1) fprintf(myf,'%8.6e', dX(i,1); fprintf(myf, 'n'); end fprintf(myf, 'n'); fprintf(myf,'第%d次迭代后節(jié)點電壓(僅PQ節(jié)點)',K) fprintf(myf, 'n'); for i=1:m fprintf(m

7、yf,'%8.6f', bus(i,2); fprintf(myf, 'n'); end fprintf(myf, 'n'); fprintf(myf,'第%d次迭代后相角(角度)',K) fprintf(myf, 'n'); for i=1:n fprintf(myf,'%8.6f', bus(i,3)*180/pi); fprintf(myf, 'n'); end fprintf(myf, 'n'); if (max(abs(dx)<eps1)&(m

8、ax(abs(b)<eps2) %判斷是否達到計算精度 break; endend%計算功率P1=0;T=0; %計算平衡節(jié)點的有功for j=1:n T=bus(n,2)*bus(j,2)*(real(Y(n,j)*cos(bus(n,3)-bus(j,3)+imag(Y(n,j)*sin(bus(n,3)-bus(j,3); P1=P1+T;end bus(n,4)=P1;for k=m+1:n %計算PV節(jié)點、平衡節(jié)點的無功 Q1=0;T=0; for j=1:n T=bus(k,2)*bus(j,2)*(real(Y(k,j)*sin(bus(k,3)-bus(j,3)-imag

9、(Y(k,j)*cos(bus(k,3)-bus(j,3); Q1=Q1+T; end bus(k,5)=Q1;endbus(:,1)=f; %換回各節(jié)點、支路的初始編號r=zeros(1,mb);for t=1:n for l=t+1:n if bus(t,1)>bus(l,1) r=bus(t,:);bus(t,:)=bus(l,:);bus(l,:)=r; end endendfor i=1:nl for j=1:2 for k=1:nb if line(i,j)=nodenum(k,1) line(i,j)=nodenum(k,2); break end end endend f

10、close(myf);Pf=loss(bus,line); %計算支路潮流及損耗 %將節(jié)點導納矩陣、節(jié)點潮流計算結果寫入文件output2 myf=fopen('output2.dat','w'); fprintf(myf, '-節(jié)點導納矩陣-n'); for k=1:n for j=1:n fprintf(myf,'%8.6f', real(Y(k,j); fprintf(myf, '+i*'); fprintf(myf,'%8.6f', imag(Y(k,j); fprintf(myf, 

11、9; '); end fprintf(myf, 'n'); end fprintf(myf, '-牛頓拉夫遜法潮流計算結果-n'); fprintf(myf, '-節(jié)點計算結果-n'); fprintf(myf, '- 節(jié)點 節(jié)點電壓 節(jié)點相角 注入有功功率（P） 注入無功功率(Q) 類型-n'); for l=1:nb for j=1:mb if j=1|j=6 fprintf(myf, ' %8.1f ', bus(l,j); elseif j=3 fprintf(myf, ' %8.6f &#

12、39;, bus(l,j)*180/pi); else fprintf(myf, ' %8.6f ', bus(l,j); end end fprintf(myf, ' n'); end fprintf(myf, '-支路計算結果-n'); fprintf(myf, '- 節(jié)點(I) 節(jié)點(J) 線路功率S(I,J) 線路功率S(J,I) 線路損耗dS(I,J)-n'); for k=1:nl for j=1:5 if j<=2 fprintf(myf,'%8.1f ', Pf(k,j); fprintf(m

13、yf, ' '); else fprintf(myf,'%8.6f', real(Pf(k,j); fprintf(myf, '+i*'); fprintf(myf,'%8.6f', imag(Pf(k,j); fprintf(myf, ' '); end end fprintf(myf, 'n'); end fclose(myf);%根據(jù)支路參數(shù)建立節(jié)點導納矩陣程序：function Y=y(bus,line)%目的：根據(jù)支路參數(shù)建立節(jié)點導納矩陣%輸入：節(jié)點參數(shù)矩陣-bus；支路參數(shù)矩陣-line

14、%輸出：節(jié)點導納矩陣-Ynb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nl I=line(k,1); J=line(k,2); Zt=line(k,3)+j*line(k,4); if Zt=0 disp('對地支路'); Yt=inf; else Yt=1/Zt; end Ym=line(k,5)+j*line(k,6); K=line(k,7); if(K=0)&(J=0) %普通線路 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt+Ym; Y(I,J)=Y(I,J)-Yt

15、; Y(J,I)=Y(I,J); end if(K=0)&(J=0) %對地支路K=J=0，R=X=0 Y(I,I)=Y(I,I)+Ym; end if K>0 %K>0時變壓器支路 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt/K2; Y(I,J)=Y(I,J)-Yt/K; Y(J,I)=Y(I,J); end if K<0 %K<0時變壓器支路 Y(I,I)=Y(I,I)+Yt+Ym; Y(J,J)=Y(J,J)+Yt*K2; Y(I,J)=Y(I,J)+Yt*K; Y(J,I)=Y(I,J); endend %形成雅克矩陣程序

16、：function J=jac(bus,Y)%形成雅可比矩陣global n;global m;for i=1:(n-1) %計算PQ、PV節(jié)點的有功 P1=0;T=0; for j=1:n T=bus(i,2)*bus(j,2)*(real(Y(i,j)*cos(bus(i,3)-bus(j,3)+imag(Y(i,j)*sin(bus(i,3)-bus(j,3); P1=P1+T; end bus(i,4)=P1;endfor i=1:n-1 %計算PV、PQ節(jié)點的無功 Q1=0;T=0; for j=1:n T=bus(i,2)*bus(j,2)*(real(Y(i,j)*sin(bus

17、(i,3)-bus(j,3)-imag(Y(i,j)*cos(bus(i,3)-bus(j,3); Q1=Q1+T; end bus(i,5)=Q1;endfor i=1:n-1 %計算H for j=1:n-1 if i=j H(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j)*sin(bus(i,3)-bus(j,3)-imag(Y(i,j)*cos(bus(i,3)-bus(j,3); N(i,j)=-bus(i,2)*bus(j,2)*(real(Y(i,j)*cos(bus(i,3)-bus(j,3)+imag(Y(i,j)*sin(bus(i,3)-bus(

18、j,3); K(i,j)=bus(i,2)*bus(j,2)*(real(Y(i,j)*cos(bus(i,3)-bus(j,3)+imag(Y(i,j)*sin(bus(i,3)-bus(j,3); L(i,j)=H(i,j); end endendfor i=1:n-1 %計算H for j=1:n-1 if j=i H(i,i)=bus(i,5)+imag(Y(i,i)*bus(i,2)2; N(i,i)=-bus(i,4)-real(Y(i,i)*bus(i,2)2; K(i,i)=N(i,i)+2*real(Y(i,i)*bus(i,2)2; L(i,i)=-H(i,i)+2*ima

19、g(Y(i,i)*bus(i,2)2; end endendN=N(1:n-1,1:m);K=K(1:m,1:n-1);L=L(1:m,1:m);J=H,N;K,L;%計算dPQ的程序：function dPQ=dPQ(Y,bus)%delP-有功偏差量%delQ-無功偏差量%Y-節(jié)點導納矩陣%bus-節(jié)點參數(shù)（P,Q,U及相角）矩陣global n;global m;delP=zeros(1,n-1);delQ=zeros(1,m);for i=1:(n-1) %形成delP矩陣(PQ、PV節(jié)點) P1=0;T=0; for j=1:n T=bus(i,2)*bus(j,2)*(real(Y

20、(i,j)*cos(bus(i,3)-bus(j,3)+imag(Y(i,j)*sin(bus(i,3)-bus(j,3); P1=P1+T; end delP(1,i)=bus(i,4)-P1;endfor i=1:m %形成delQ矩陣(PQ節(jié)點) Q1=0;T=0; for j=1:n T=bus(i,2)*bus(j,2)*(real(Y(i,j)*sin(bus(i,3)-bus(j,3)-imag(Y(i,j)*cos(bus(i,3)-bus(j,3); Q1=Q1+T; end delQ(1,i)=bus(i,5)-Q1;enddPQ=delP,delQ;%計算線路損耗、線路潮

21、流程序：function Pf=loss(bus,line)%計算線路損耗、線路潮流nl,ml=size(line);Pf=zeros(nl,5);for k=1:nl I=line(k,1); J=line(k,2); Zt=line(k,3)+i*line(k,4); if Zt=0 Yt=inf; else Yt=1/Zt; end Ym=line(k,5)+i*line(k,6); K=line(k,7); if (K=0)&(J=0) %普通線路潮流 S(I,J)=bus(I,2)2*(conj(Yt)+conj(Ym)-bus(I,2)*(cos(bus(I,3)+i*si

22、n(bus(I,3)*bus(J,2)*(cos(bus(J,3)-i*sin(bus(J,3)*conj(Yt); S(J,I)=bus(J,2)2*(conj(Yt)+conj(Ym)-bus(J,2)*(cos(bus(J,3)+i*sin(bus(J,3)*bus(I,2)*(cos(bus(I,3)-i*sin(bus(I,3)*conj(Yt); delS(I,J)=S(I,J)+S(J,I); end if(K=0)&(J=0) %對地支路潮流 J=5; S(I,5)=bus(I,2)2*conj(Ym); end if K>0 %變壓器支路k>0時的潮流 S

23、(I,J)=bus(I,2)2*(conj(Ym+Yt*(1-1/K)+conj(Yt/K)-bus(I,2)*(cos(bus(I,3)+i*sin(bus(I,3)*bus(J,2)*(cos(bus(J,3)-i*sin(bus(J,3)*conj(Yt/K); S(J,I)=bus(J,2)2*(conj(Yt)/K2-bus(J,2)*(cos(bus(J,3)+i*sin(bus(J,3)*bus(I,2)*(cos(bus(I,3)-i*sin(bus(I,3)*conj(Yt/K); delS(I,J)=S(I,J)+S(J,I); end if K<0 %變壓器支路k&

24、lt;0時的潮流 S(I,J)=bus(I,2)2*(conj(Ym+Yt)+bus(I,2)*(cos(bus(I,3)+i*sin(bus(I,3)*bus(J,2)*(cos(bus(J,3)-i*sin(bus(J,3)*conj(Yt*K); S(J,I)=bus(J,2)2*(conj(Yt)*K2+bus(J,2)*(cos(bus(J,3)+i*sin(bus(J,3)*bus(I,2)*(cos(bus(I,3)-i*sin(bus(I,3)*conj(Yt*K); delS(I,J)=S(I,J)+S(J,I); end if J=5&Zt=0 Sp=line(k,

25、1) line(k,2) S(I,5) 0 S(I,5); else Sp=line(k,1) line(k,2) S(I,J) S(J,I) delS(I,J); end Pf(k,:)=Sp;end %輸入的參數(shù)數(shù)據(jù)：% data for test case%各節(jié)點參數(shù)：節(jié)點編號,注入有功,注入無功,(Sn=100MVA)電壓幅值,電壓相位,類型%類型：1=PQ節(jié)點,2=PV節(jié)點,3=平衡節(jié)點% (bus#)( volt )( ang )( p )( q )(bus type)bus=1,1.0,0.0,-0.478,0.039,1;2,1.0,0.0,-0.076,-0.016,1;3,

26、1.0,0.0,0.0,0.0,1;4,1.0,0.0,-0.295,-0.166,1;5,1.0,0.0,-0.09,-0.058,1;6,1.0,0.0,-0.035,-0.018,1;7,1.0,0.0,-0.061,-0.016,1;8,1.0,0.0,-0.135,-0.058,1;9,1.0,0.0,-0.149,-0.05,1;10,1.045,0.0,0.183,0.0,2;11,1.010,0.0,-0.942,0.0,2;12,1.70,0.0,-0.112,0.047,2;13,1.90,0.0,0.0,0.174,2;14,1.060,0.0,0.0,0.0,3; %各

27、支路參數(shù)：起點編號,終點編號,電阻,電抗,電導,電納line = 1,2,0.01335,0.04211,0.0,0.0,0;1,3,0.0,0.20912,0.0,0.0,0;1,4,0.0,0.55618,0.0,0.0,0;1,10,0.05811,0.17632,0.0,0.0340,0;1,11,0.06701,0.17103,0.0,0.0128,0;2,10,0.05695,0.17388,0.0,0.0346,0;2,12,0.0,0.25202,0.0,0.0,0;2,14,0.05403,0.22304,0.0,0.0492,0;3,4,0.0,0.11001,0.0,0.

28、0,0;3,13,0.0,0.17615,0.0,0.0,0;4,5,0.03181,0.08450,0.0,0.0,0;4,9,0.12711,0.27038,0.0,0.0,0;5,6,0.08205,0.19207,0.0,0.0,0;6,12,0.09498,0.19890,0.0,0.0,0;7,8,0.22092,0.19988,0.0,0.0,0;7,12,0.12291,0.25581,0.0,0.0,0;8,9,0.17093,0.34802,0.0,0.0,0;8,12,0.06615,0.13027,0.0,0.0,0;10,11,0.04699,0.19797,0.0,

29、0.0438,0;10,14,0.01938,0.05917,0.0,0.0528,0;輸出結果數(shù)據(jù)1：-第1次迭代時雅可比矩陣-38. 21. 4. 1. -0. -0. -0. -0. -0. 5. 5. -0. -0. -10. 6. -0. -0. -0. -0. -0. -0. -0. 21. -38. -0. -0. -0. -0. -0. -0. -0. 5. -0. 6. -0. 6. -9. -0. -0. -0. -0. -0. -0. -0. 4. -0. -24. 9. -0. -0. -0. -0. -0. -0. -0. -0. 10. -0. -0. -0. -0

30、. -0. -0. -0. -0. -0. 1. -0. 9. -24. 10. -0. -0. -0. 3. -0. -0. -0. -0. -0. -0. -0. -5. 3. -0. -0. -0. 1. -0. -0. -0. 10. -14. 4. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. 3. -5. 1. -0. -0. -0. -0. -0. -0. -0. 4. -11. -0. -0. -0. -0. -0. 6. -0. -0. -0. -0. -0. 1. -2. -0. -0. -0. -0. -0. -0. -0. -0. -0

31、. -7. 2. -0. -0. -0. 5. -0. -0. -0. -0. -0. -0. -0. -2. 2. -0. -0. -0. -0. -0. -0. -0. 2. -14. 2. -0. -0. 10. -0. -0. -0. -0. -0. -0. -0. 2. -4. 1. -0. -0. -0. 3. -0. -0. -0. 2. -5. -0. -0. -0. -0. -0. -0. -0. 1. -0. -0. -0. 1. -2. 5. 5. -0. -0. -0. -0. -0. -0. -0. -32. 5. -0. -0. 1. 1. -0. -0. -0.

32、-0. -0. -0. -0. 5. -0. -0. -0. -0. -0. -0. -0. -0. 5. -10. -0. -0. 2. -0. -0. -0. -0. -0. -0. -0. -0. -0. 6. -0. -0. -0. 6. 5. 10. -0. -0. -0. -29. -0. -0. -0. -0. -0. -0. 3. 2. 5. -0. -0. -0. 10. -0. -0. -0. -0. -0. -0. -0. -0. -0. -10. -0. -0. -0. -0. -0. -0. -0. -0. -0. 10. -6. 0. 0. 0. 0. 0. 0.

33、0. -1. -2. 0. 0. -37. 21. 4. 1. -0. -0. -0. -0. -0. -6. 9. 0. 0. 0. 0. 0. 0. 0. -1. 0. 0. 0. 21. -31. -0. -0. -0. -0. -0. -0. -0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 4. -0. -14. 9. -0. -0. -0. -0. -0. 0. 0. 0. 5. -3. 0. 0. 0. -1. 0. 0. 0. 0. 1. -0. 9. -24. 10. -0. -0. -0. 3. 0. 0. 0. -3. 5. -1. 0

34、. 0. 0. 0. 0. 0. 0. -0. -0. -0. 10. -14. 4. -0. -0. -0. 0. 0. 0. 0. -1. 5. 0. 0. 0. 0. 0. -3. 0. -0. -0. -0. -0. 4. -5. -0. -0. -0. 0. 0. 0. 0. 0. 0. 5. -2. 0. 0. 0. -2. 0. -0. -0. -0. -0. -0. -0. -3. 2. -0. 0. 0. 0. 0. 0. 0. -2. 8. -1. 0. 0. -5. 0. -0. -0. -0. -0. -0. -0. 2. -6. 2. 0. 0. 0. -1. 0.

35、0. 0. -1. 2. 0. 0. 0. 0. -0. -0. -0. 3. -0. -0. -0. 2. -5. -第1次迭代時dPQ的誤差-3.e-0016.e-0020.e+000-2.e-001-9.e-0021.e+0001.e+0002.e+000-1.e-0016.e-002-9.e-001-7.e+0000.e+0003.e-0013.e+0005.e+000-1.e-001-5.e-0022.e+0002.e+0004.e+000-5.e-002-第1次迭代時dx(誤差)-7.e-001-8.e-001-1.e+000-1.e+000-1.e+000-1.e+000-2.e

36、+000-2.e+000-1.e+000-4.e-001-6.e-001-2.e+000-1.e+000-5.e-001-5.e-001-1.e+000-1.e+000-1.e+000-1.e+000-2.e+000-1.e+000-1.e+000第1次迭代后節(jié)點電壓(僅PQ節(jié)點)1.1.2.2.2.2.3.2.2.第1次迭代后相角(角度)44.48.68.80.90.114.125.123.98.24.35.128.68.0.-第2次迭代時雅可比矩陣-88. 50. 15. 4. 0. 0. 0. 0. 0. 8. 8. 0. 0. -17. 20. 6. 3. -0. -0. -0. -0

37、. -0. 53. -70. 0. 0. 0. 0. 0. 0. 0. 8. -0. 1. 0. 12. -29. -0. -0. -0. -0. -0. -0. -0. 15. -0. -85. 45. 0. 0. 0. 0. 0. -0. -0. 0. 24. -6. -0. 3. 10. -0. -0. -0. -0. -0. 4. -0. 45. -111. 48. 0. 0. 0. 13. -0. -0. 0. -0. -3. -0. -10. -24. 28. -0. -0. -0. 12. -0. -0. -0. 54. -73. 18. 0. 0. 0. -0. -0. 0.

38、-0. -0. -0. -0. 10. -30. 19. -0. -0. -0. -0. -0. -0. -0. 27. -42. 0. 0. -0. -0. -0. 14. -0. -0. -0. -0. -0. -0. -35. -0. -0. -0. -0. -0. -0. -0. -0. -0. -36. 20. -0. -0. -0. 15. -0. -0. -0. -0. -0. -0. -0. -43. 21. -0. -0. -0. -0. -0. -0. -0. 18. -66. 18. -0. -0. 28. -0. -0. -0. -0. -0. -0. -0. 22.

39、-72. 0. -0. -0. -0. 18. -0. 0. 0. 11. -29. -0. -0. 0. -0. -0. -0. -0. 2. -0. -0. -0. 14. -14. 6. 6. 0. 0. 0. 0. -0. -0. 0. -35. 4. -0. 0. 5. 6. -0. -0. -0. -0. 0. 0. -0. 7. 0. 0. 0. 0. 0. -0. 0. 0. 5. -12. -0. 0. 4. -0. -0. -0. -0. -0. 0. -0. -0. -0. 1. -0. -0. -0. 18. 16. 31. -0. -0. -0. -68. -0. -

40、0. -10. -0. -0. -0. 3. 7. 12. -0. -0. -0. 24. 0. 0. 0. 0. 0. 0. -0. -0. 0. -24. -0. -0. -0. -0. -0. -0. -0. -0. -0. 33. -20. -6. -3. 0. 0. 0. 0. 0. 0. -1. 0. 0. -97. 50. 15. 4. 0. 0. 0. 0. 0. -12. 17. 0. 0. 0. 0. 0. 0. 0. 1. 0. -10. 0. 53. -99. 0. 0. 0. 0. 0. 0. 0. 6. 0. 3. -10. 0. 0. 0. 0. 0. 0. 0.

41、 0. 0. 15. -0. -121. 45. 0. 0. 0. 0. 0. 3. 0. 10. 27. -28. 0. 0. 0. -12. 0. 0. 0. 0. 4. -0. 45. -125. 48. 0. 0. 0. 13. 0. 0. 0. -10. 30. -19. 0. 0. 0. 0. 0. 0. 0. -0. -0. -0. 54. -80. 18. 0. 0. 0. 0. 0. 0. 0. 0. 12. 0. 0. 0. 0. 0. -12. 0. -0. -0. -0. -0. 27. -64. 0. 0. -0. 0. 0. 0. 0. 0. 0. 29. -21.

42、 0. 0. 0. -8. 0. -0. -0. -0. -0. -0. -0. -62. 20. -0. 0. 0. 0. 0. 0. 0. -22. 40. -0. 0. 0. -17. 0. -0. -0. -0. -0. -0. -0. 18. -113. 18. 0. 0. 0. -2. 0. 0. 0. -14. 16. 0. 0. 0. 0. -0. -0. -0. 18. -0. 0. 0. 11. -36. -第2次迭代時dPQ的誤差-7.e+000-6.e+0003.e+0009.e-001-1.e-001-1.e+001-6.e+000-1.e+0019.e-0011.e

43、+0003.e-001-6.e+0000.e+000-4.e+000-1.e+001-1.e+001-7.e+000-3.e+000-1.e+001-1.e+001-2.e+001-3.e+000-第2次迭代時dx(誤差)-7.e-0017.e-0011.e+0001.e+0001.e+0002.e+0002.e+0002.e+0002.e+0002.e-0014.e-0012.e+0001.e+000-1.e-001-1.e-0012.e-002-1.e-001-1.e-0015.e-0022.e-0012.e-001-7.e-002第2次迭代后節(jié)點電壓(僅PQ節(jié)點)1.1.2.2.2.2.

44、2.2.2.第2次迭代后相角(角度)1.3.-5.-8.-12.-18.-19.-19.-17.8.7.-20.-5.0.-第3次迭代時雅可比矩陣-107. 64. 18. 7. -0. -0. -0. -0. -0. 8. 8. -0. -0. -33. 22. -2. -1. -0. -0. -0. -0. -0. 65. -93. -0. -0. -0. -0. -0. -0. -0. 9. 0. 10. -0. 18. -40. -0. -0. -0. -0. -0. -0. -0. 18. 0. -94. 51. -0. -0. -0. -0. -0. 0. 0. -0. 24. 2

45、. -0. -0. -2. -0. -0. -0. -0. -0. 7. 0. 51. -151. 70. -0. -0. -0. 21. 0. 0. -0. 0. 1. -0. 2. -37. 20. -0. -0. -0. 6. 0. 0. 0. 66. -94. 28. -0. -0. -0. 0. 0. -0. 0. -0. -0. -0. 30. -39. 8. -0. -0. -0. 0. 0. 0. 0. 25. -42. -0. -0. 0. 0. 0. 16. 0. -0. -0. -0. -0. 14. -21. -0. -0. -0. 0. 0. 0. 0. 0. 0.

46、-22. 11. 0. 0. 0. 11. 0. -0. -0. -0. -0. -0. -0. -20. 11. -0. 0. 0. 0. 0. 0. 0. 10. -47. 13. 0. 0. 23. 0. -0. -0. -0. -0. -0. -0. 12. -36. 7. 0. 0. 0. 18. 0. -0. -0. 14. -32. 0. 0. -0. 0. -0. -0. -0. 12. -0. -0. -0. 6. -17. 9. 9. -0. -0. -0. -0. -0. -0. -0. -41. 5. -0. -0. 1. 2. -0. -0. -0. -0. -0. -0. -0. 9.