潮流計(jì)算matlab

1、matlab 潮流計(jì)算08292007 電氣 0809 班% 主程序：dfile,pathname=uigetfile('ieee14.m','Select Data File');if pathname=0error('you must select a valid data file')elselfile=length(dfile);% strip off .m eval(dfile(1:lfile-2);endglobal n;global m;nb,mb=size(bus);nl,ml=size(line);nSW=0;nPV=0;nPQ

2、=0;for i=1:nb, type=bus(i,6);if type=3, nSW=nSW+1;%打開數(shù)據(jù)文件%節(jié)點(diǎn)重新編號% 平衡節(jié)點(diǎn)數(shù)目% PV 節(jié)點(diǎn)數(shù)目% PQ 節(jié)點(diǎn)數(shù)目% nb 為總節(jié)點(diǎn)數(shù)% 統(tǒng)計(jì)平衡節(jié)點(diǎn)數(shù)目SW(nSW,:)=bus(i,:);elseif type=2,nPV=nPV+1; % 統(tǒng)計(jì) PV 節(jié)點(diǎn)數(shù)目PV(nPV,:)=bus(i,:);elsenPQ=nPQ+1; % 統(tǒng)計(jì)PQ節(jié)點(diǎn)數(shù)目PQ(nPQ,:)=bus(i,:);endend;bus=PQ;PV;SW;newbus=1:nb'f=bus(:,1);nodenum=newbus bus(:,1);

3、bus(:,1)=newbus;for i=1:nlfor j=1:2for k=1:nbif line(i,j)=nodenum(k,2) line(i,j)=nodenum(k,1); breakendendendendY=y(bus,line); % 形成節(jié)點(diǎn)導(dǎo)納矩陣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:Kmaxfor i=1:mfor j=1:m if i=jUm(i,j

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

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

6、39;,K) fprintf(myf, 'n');for i=1:(n+m-1)fprintf(myf,'%8.6e', dX(i,1); fprintf(myf, 'n');endfprintf(myf, 'n');fprintf(myf,'第次迭代后節(jié)點(diǎn)電壓(僅PQ節(jié)點(diǎn))',K) fprintf(myf, 'n');for i=1:mfprintf(myf,'%8.6f', bus(i,2);fprintf(myf, 'n');endfprintf(myf, &#

7、39;n');fprintf(myf,' 第 %d 次迭代后相角(角度)',K)fprintf(myf, 'n');for i=1:nfprintf(myf,'%8.6f', bus(i,3)*180/pi);fprintf(myf, 'n'); endfprintf(myf, 'n');if (max(abs(dx)<eps1)&(max(abs(b)<eps2) %判斷是否達(dá)到計(jì)算精度 break;endend% 計(jì)算功率P1=0;T=0;%計(jì)算平衡節(jié)點(diǎn)的有功for j=1:nT=b

8、us(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;endbus(n,4)=P1;for k=m+1:n%計(jì)算 PV 節(jié)點(diǎn)、平衡節(jié)點(diǎn)的無功Q1=0;T=0;for j=1:nT=bus(k,2)*bus(j,2)*(real(Y(k,j)*sin(bus(k,3)-bus(j,3)-imag(Y(k,j)*cos(bus(k,3)-bus(j,3); Q1=Q1+T;endbus(k,5)=Q1;endbus(:,1)=f;%換回各節(jié)點(diǎn)、支路的初始編號r=ze

9、ros(1,mb);for t=1:nfor l=t+1:nif bus(t,1)>bus(l,1) r=bus(t,:);bus(t,:)=bus(l,:);bus(l,:)=r;endendendfor i=1:nlfor j=1:2for k=1:nbif line(i,j)=nodenum(k,1)line(i,j)=nodenum(k,2);breakendendendendfclose(myf);Pf=loss(bus,line); %計(jì)算支路潮流及損耗%將節(jié)點(diǎn)導(dǎo)納矩陣、節(jié)點(diǎn)潮流計(jì)算結(jié)果寫入文件output2myf=fopen('output2.dat',&#

10、39;w');fprintf(myf, '- 節(jié)點(diǎn)導(dǎo)納矩陣-n');for k=1:nfor j=1:nfprintf(myf,'%8.6f', real(Y(k,j);fprintf(myf, '+i*');fprintf(myf,'%8.6f', imag(Y(k,j);fprintf(myf, '');endfprintf(myf, 'n');endfprintf(myf,fprintf(myf,牛頓拉夫遜法潮流計(jì)算結(jié)果n');fprintf(myf, '-節(jié)點(diǎn)for

11、l=1:nbfor j=1:mbif j=1|j=6fprintf(myf, ' %8.1felseif j=3fprintf(myf, ' %8.6f elsefprintf(myf, ' %8.6f end節(jié)點(diǎn)計(jì)算結(jié)果-節(jié)點(diǎn)電壓', bus(l,j);-n');節(jié)點(diǎn)相角', bus(l,j)*180/pi);', bus(l,j);注入有功功率(P)注入無功功率 (Q)類型-n');endfprintf(myf, ' n'); endfprintf(myf, '- 支路計(jì)算結(jié)果fprintf(myf,節(jié)

12、點(diǎn) (I)-n');節(jié)點(diǎn) (J)線路功率 S(I,J)線路功率 S(J,I)線路損耗dS(I,J)-n');for k=1:nlfor j=1:5if j<=2fprintf(myf,'%8.1f ', Pf(k,j);fprintf(myf, '');elsefprintf(myf,'%8.6f', real(Pf(k,j);fprintf(myf, '+i*');fprintf(myf,'%8.6f', imag(Pf(k,j);fprintf(myf, '');enden

13、dfprintf(myf, 'n');endfclose(myf);% 根據(jù)支路參數(shù)建立節(jié)點(diǎn)導(dǎo)納矩陣程序：function Y=y(bus,line)% 目的：根據(jù)支路參數(shù)建立節(jié)點(diǎn)導(dǎo)納矩陣%輸入：節(jié)點(diǎn)參數(shù)矩陣-bus ；支路參數(shù)矩陣-line% 輸出：節(jié)點(diǎn)導(dǎo)納矩陣-Ynb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+j*line(k,4);if Zt=0disp('對地支路)；Yt=inf;elseYt=1/Zt;endYm=lin

14、e(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;Y(J,I)=Y(I,J);endif(K=0)&(J=0)%對地支路 K=J=0 ， R=X=0Y(I,I)=Y(I,I)+Ym;endif K>0%K>0 時變壓器支路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J尸丫(J,J)+Yt/KA2;Y(I,J)=Y(I,J)-Yt/K;Y(J,I)=Y(I,J);endif K<0%K<0 時

15、變壓器支路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt*KA2;Y(I,J)=Y(I,J)+Yt*K;Y(J,I)=Y(I,J);endend% 形成雅克矩陣程序：function J=jac(bus,Y)% 形成雅可比矩陣global n;global m;for i=1:(n-1)%計(jì)算 PQ、 PV 節(jié)點(diǎn)的有功P1=0;T=0;for j=1:nT=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;endbus(i,4)=

16、P1;endfor i=1:n-1%計(jì)算PV、 PQ 節(jié)點(diǎn)的無功Q1=0;T=0;for j=1:nT=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;endbus(i,5)=Q1;endfor i=1:n-1%計(jì)算 Hfor j=1:n-1if i=jH(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)=

17、-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);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);endendendfor i=1:n-1%計(jì)算 Hfor j=1:n-1if j=iH(i,i尸bus(i,5)+imag(Y(i,i)*bus(i,242;N(i,i尸-bus(i,4)-real(Y(i,

18、i)*bus(i,242;K(i,i)=N(i,i)+2*real(Y(i,i)*bus(i,2)A2;L(i,i尸-H(i,i)+2*imag(Y(i,i)*bus(i,242;endend endN=N(1:n-1,1:m);K=K(1:m,1:n-1);L=L(1:m,1:m);J=H,N;K,L;% 計(jì)算 dPQ 的程序：function dPQ=dPQ(Y,bus)%delP-有功偏差量%delQ-無功偏差量%Y- 節(jié)點(diǎn)導(dǎo)納矩陣%bus-節(jié)點(diǎn)參數(shù)( P,Q,U 及相角)矩陣global n;global m;delP=zeros(1,n-1);delQ=zeros(1,m);for

19、i=1:(n-1)%形成 delP 矩陣(PQ、PV 節(jié)點(diǎn))P1=0;T=0;for j=1:nT=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;enddelP(1,i)=bus(i,4)-P1;endfor i=1:m % 形成 delQ 矩陣 (PQ 節(jié)點(diǎn) )Q1=0;T=0;for j=1:nT=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

20、)-bus(j,3);Q1=Q1+T;enddelQ(1,i)=bus(i,5)-Q1;enddPQ=delP,delQ;% 計(jì)算線路損耗、線路潮流程序：function Pf=loss(bus,line)% 計(jì)算線路損耗、線路潮流nl,ml=size(line);Pf=zeros(nl,5);for k=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+i*line(k,4);if Zt=0Yt=inf;elseYt=1/Zt;endYm=line(k,5)+i*line(k,6);K=line(k,7);S(I,J尸bus(I,2)A2*(conj(Yt)+

21、conj(Ym)-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);S(J,I)=bus(J,2F2*(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);endif(K=0)&(J=0)%對地支路潮流J=5;S(I,5)=bus(I,242*conj(Ym);endi

22、f K>0%變壓器支路 k>0 時的潮流S(I,J尸bus(I,2F2*(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,2F2*(conj(Yt)/KA2-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,

23、I);endif K<0%變壓器支路k<0 時的潮流S(I,J尸bus(I,2F2*(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,2F2*(conj(Yt)*KA2+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);endif J=

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

25、1.0,0.0,-0.076,-0.016,1;3,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;

26、14,1.060,0.0,0.0,0.0,3;%各支路參數(shù)：起點(diǎn)編號 ,終點(diǎn)編號 , 電阻,電抗,電導(dǎo),電納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.0

27、492,0;3,4,0.0,0.11001,0.0,0.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,

28、0;10,11,0.04699,0.19797,0.0,0.0438,0;10,14,0.01938,0.05917,0.0,0.0528,0;輸出結(jié)果數(shù)據(jù)1 ：- 第 1 次迭代時雅可比矩陣-38.62403321.5785544.7819431.797979-0.000000-0.000000-0.000000-0.000000-0.0000005.346051-0.0000005.119505-0.000000-0.000000-0.000000-0.000000 -0.000000-10.4172586.840981-0.000000-0.000000-0.00000021.57855

29、4-38.240787-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000005.427654-0.000000-0.000000-0.0000006.745496-0.000000-0.000000 -0.0000006.840981-9.429913-0.000000-0.000000-0.0000004.781943-0.000000-24.6582889.090083-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.00000010.7

30、86262-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000 -0.000000 -0.000000 -0.0000001.797979-0.0000009.090083-24.28250610.365394-0.000000-0.000000-0.0000003.029050-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-5.3260553.902050-0.000000-0.000000-0.000000-0.000000-0.000000 1.4240

31、05-0.00000010.365394-14.7683384.402944-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000003.902050-5.7829341.880885-0.000000-0.000000-0.000000-0.000000 -0.000000-0.000000-0.0000004.402944-11.362870-0.000000-0.000000-0.000000-0.000000-0.0000006.959926-0.000000-0.

32、000000-0.000000-0.000000-0.0000001.880885-2.467393-0.000000-0.000000-0.000000-0.000000 -0.000000-0.000000-0.000000-0.000000-0.000000-7.6511132.251975-0.000000-0.000000-0.0000005.399139-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-2.946815-0.0000002.489025 -0.000000-0.00000

33、0-0.000000-0.000000-0.0000002.251975-14.9416222.314963-0.000000-0.00000010.374684-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000002.489025-0.000000-4.555697 1.136994-0.0000003.029050-0.000000-0.000000-0.0000002.314963-5.344014-0.000000-0.000000-0.000000-0.000000-0.000000-0.000

34、000-0.0000001.424005-0.000000-0.0000005.346051-0.0000005.4276541.136994 -2.561000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-32.7276445.047017-0.000000-0.0000001.7619051.777691-0.000000-0.000000-0.000000-0.0000005.119505-0.000000-0.000000-0.000000 -0.000000-0.000000-0.000000-0.00

35、0000-0.000000-0.000000-0.000000-0.0000005.047017-10.166523-0.000000-0.0000002.005835-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000006.745496-0.000000 -0.000000-0.000000-0.000000-0.0000006.9599265.39913910.374684-0.000000-0.000000-0.000000-29.479246-0.000000-0.000000-0.000000-0.000000-

36、0.000000-0.0000003.323549-0.0000002.594145 5.268177 -0.000000-0.00000010.786262-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-10.786262-0.000000-0.000000-0.000000-0.000000-0.000000-0.00000010.608721-0.000000-6.840981-0.000000 -0.0000000.000000 0.000000 0.000000 0.0

37、00000 0.000000 0.0000000.000000-1.761905-2.0058350.0000000.000000-37.96863121.5785544.7819431.797979-0.000000-0.000000-0.000000-6.840981-0.0000009.706123-0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.0000000.000000-1.7776910.0000000.0000000.00000021.578554-31.542421-0.000000-0.000000-0.0000

38、00-0.000000-0.0000000.000000-0.0000000.000000-0.0000000.000000 0.0000000.0000000.000000 0.0000000.0000000.0000000.0000000.0000000.0000000.0000004.781943-0.000000-14.4397249.090083-0.000000-0.000000-0.0000000.000000-0.0000000.000000-0.0000000.000000 5.326055-3.902050 0.000000 0.000000 0.000000-1.4240

39、050.0000000.0000000.0000000.0000001.797979-0.0000009.090083-24.28250610.365394-0.000000-0.0000000.000000-0.0000000.0000003.0290500.000000 -3.9020505.782934 -1.8808850.000000 0.0000000.0000000.0000000.0000000.0000000.000000-0.000000-0.000000-0.00000010.365394-14.7683384.402944-0.0000000.000000-0.0000

40、000.000000-0.0000000.0000000.000000-1.8808855.2044330.000000 0.0000000.0000000.0000000.000000-3.3235490.000000-0.000000-0.000000-0.000000-0.0000004.402944-5.631166-0.000000 -0.000000 -0.0000000.0000000.0000000.0000000.0000000.000000 0.0000005.083169-2.4890250.0000000.0000000.000000-3.204764-2.594145

41、2.2519750.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000000.0000000.0000000.0000000.0000000.000000 0.000000-2.489025 8.894195-1.1369940.0000000.0000002.251975-5.268177-6.3977650.0000002.314963-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000000.0000000.0000000.000000-1.424005

42、0.000000 0.000000 0.000000 -1.1369942.5610000.0000000.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 3.029050 -0.0000002.314963 -5.344014-第 1 次迭代時dPQ 的誤差-0.000000-0.000000-3.822688e-0016.210513e-0020.000000e+000-2.950000e-001-9.000000e-0021.333520e+0001.007177e+0002.034249e+000-1.490000e-001

43、6.056626e-002-9.219354e-001-7.942109e+0000.000000e+0003.667009e-0013.333183e+0005.109282e+000-1.660000e-001-5.800000e-0022.847852e+0002.207175e+0004.213929e+000-5.000000e-002第1次迭代時dx（誤差）-7.699084e-001-8.544764e-001-1.189723e+000-1.410571e+000-1.585607e+000-1.994895e+000-2.196974e+000-2.162427e+000-1

44、.721659e+000-4.249173e-001-6.178169e-001-2.246296e+000-1.189723e+000-5.568104e-001-5.586033e-001-1.299237e+000-1.208867e+000-1.285646e+000-1.499291e+000-2.011550e+000-1.901143e+000-1.488513e+000第 1 次迭代后節(jié)點(diǎn)電壓（僅 PQ 節(jié)點(diǎn) ）1.5568101.5586032.2992372.2088672.2856462.4992913.0115502.9011432.488513第 1 次迭代后相角 （

45、 角度 ）44.11250048.95789268.16611180.81979290.848607114.299041125.877362123.89794698.64381324.34596735.398301128.70329568.1661110.000000-88.46859650.77006215.6304994.9568020.000000 0.000000 0.000000 0.000000 0.0000008.3512030.0000000.000000-17.67902420.9626216.9766783.695668-0.000000-0.000000-0.000000

46、-0.00000053.574257-70.1633000.0000000.0000000.000000 0.000000 0.000000 0.0000000.000000-0.0000001.8716490.00000012.117328-29.191523-0.000000-0.000000-0.000000-0.000000-0.000000-0.00000015.630499-0.000000-85.47526245.0445950.0000000.0000000.0000000.000000-0.000000-0.0000000.00000024.800168-6.976678-0

47、.0000003.13630210.112981-0.000000-0.000000-0.000000-0.0000004.956802-0.00000045.044595-111.55769648.1013580.0000000.0000000.000000- 第 2 次迭代時雅可比矩陣-8.760031-0.0000008.844924-0.0000000.000000-0.00000013.454941-0.000000-0.0000000.000000-0.000000-3.695668-0.000000-10.112981-24.72061128.512426-0.000000-0.

48、000000-0.000000 -0.000000 12.548251-0.000000-0.00000054.962689-73.76120518.7985160.0000000.0000000.000000-0.000000-0.0000000.000000-0.000000-0.000000-0.000000-0.00000010.286004-30.26975019.866386-0.000000-0.000000 -0.000000 -0.000000-0.000000-0.000000-0.00000027.350216-42.1319390.0000000.000000-0.00

49、0000-0.000000-0.00000014.781722-0.000000-0.000000-0.000000-0.000000-0.000000-0.152201-35.701334-0.000000-0.000000 -0.000000 -0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-36.26958720.414751-0.000000-0.000000-0.00000015.854836-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0

50、00000-43.168982 21.053877 -0.000000-0.000000-0.000000-0.000000-0.000000-0.00000018.912486-66.24245918.617657-0.000000-0.00000028.712316-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.00000022.413069 -72.744607 0.293713-0.000000-0.00000018.246829-0.0000000.0000000.00000011.613550-29

51、.860379-0.000000-0.0000000.000000-0.000000-0.000000-0.000000-0.0000002.355263-0.000000-0.0000006.904762-0.000000 14.554342 -14.8093936.537084 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000-35.8518614.723753-0.000000 0.000000 5.396002 6.042147 -0.000000 -0.000000 -0.000000 -0.000000

52、0.0000000.0000007.404987-0.0000000.000000 0.0000000.0000000.000000 0.000000-0.0000000.0000000.0000005.183063-12.588050-0.0000000.0000004.294174-0.000000-0.000000-0.000000-0.000000-0.0000000.000000-0.000000-0.000000 -0.0000001.871649-0.000000-0.000000-0.00000018.91440916.62516731.272979-0.000000-0.00

53、0000-0.000000-68.684204-0.000000-0.000000-10.345614-0.000000-0.000000-0.0000003.718215-0.0000007.001258 12.708639 -0.000000-0.000000 24.800168 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000-0.000000-0.0000000.000000-24.800168-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.00000033.280772-0.000000 -0.000000-20.962621-6.976678 -3.695668 0.000000 0.000000 0.000000 0.000000 0.0000000.233337-1.8791420.0000000.000000-97.16587550.77006215.6304994.95680