潮流計(jì)算matlab-牛頓拉夫遜法法

1、用 matlab 潮流計(jì)算（牛頓拉夫遜法）% 主程序：d=uigetfile('ieee14.m','Select Data File');if pathname=0error('you must select a valid data file') elsel(dfile);% strip off .meval(d); %打開數(shù)據(jù)文件endglobal n;global m;nb,mb=size(bus);nl,ml=size(line);nSW=0;nPV=0;nPQ=0;for i=1:nb, type=bus(i,6);if type=3

2、, nSW=nSW+1; SW(nSW,:)=bus(i,:);elseif type=2, nPV=nPV+1;PV(nPV,:)=bus(i,:);elsenPQ=nPQ+1;PQ(nPQ,:)=bus(i,:);%節(jié)點(diǎn)重新編號(hào)% 平衡節(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ù)目% 統(tǒng)計(jì) PV 節(jié)點(diǎn)數(shù)目% 統(tǒng)計(jì) PQ 節(jié)點(diǎn)數(shù)目endend;bus=PQ;PV;SW;newbus=1:nb'f=bus(:,1);nodenum=newbus bus(:,1);bus(:,1)=newbus;for i=1:nlfor j=1:2for k=1:

3、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:mif i=j Um(i,j)=bus(i,2);endendendb=dPQ(Y,bus);C=jac(bus,

4、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 次迭代時(shí)雅可比矩陣-',K)fprintf(myf, 'n'); for i=1:(n+m-1)for j=1:(n+m-1)fprintf(myf,'%8.6f', C(i,j); fprintf(m

5、yf, '');endfprintf(myf, 'n');endfprintf(myf,'-第 次迭代時(shí) 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 次迭代時(shí) dx(誤差)-',K) fprintf(myf, 'n');for i=1:(

6、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, 'n');fprintf(myf,' 第 %d 次迭代后相角(角

7、度)',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=bus(n,2)*bus(j,2)*(real(Y(n,j)*cos(bus(n,3)-b

8、us(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)、支路的初始編號(hào)r=zeros(1,mb);for t=1:nfor l=t+1:nif bus(t,1)>

9、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','w');fprintf(myf, '- 節(jié)點(diǎn)導(dǎo)納矩陣-n');

10、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, ' 牛頓拉夫遜法潮流計(jì)算結(jié)果n');fprintf(myf, ' 節(jié)點(diǎn)計(jì)算結(jié)果n');fprintf(myf, '-節(jié)點(diǎn)for l=1:nbfor j=1:mbif

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

12、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, '');endendfprintf(myf, 'n');e

13、ndfclose(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=line(k,5)+j*line(k,6);K=line(k,7

14、);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 時(shí)變壓器支路Y(I,I)=Y(I,I)+Yt+Ym;丫(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 時(shí)變壓器支路Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)

15、=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)=P1;endfor i=1:n-1%計(jì)算PV、 PQ 節(jié)點(diǎn)的無功

16、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)=-bus(i,2)*bus(j,2)*(real(Y(i,j)*

17、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(丫(i,i)*bus(i,242;K(i,i)=N(i,i)+2*re

18、al(Y(i,i)*bus(i,2)A2;L(i,i尸-H(i,i)+2*imag(Y(i,i)*bus(i,242;endendendN=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 i=1:(n-1)%形成 delP 矩陣(PQ、 PV 節(jié)點(diǎn) )P

19、1=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)-bus(j,3);Q1=Q1+T;enddelQ(1,i)

20、=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);if (K=0)&(J=0)%普通線路潮流S(I,J尸bus(I,2)A2*(conj(Yt)+conj(Y

21、m)-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);endif K>

22、;0%變壓器支路k>0 時(shí)的潮流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,I);endi

23、f K<0%變壓器支路k<0 時(shí)的潮流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=5&Z

24、t=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)編號(hào)，注入有功，注入無功，(Sn=100MVA)電壓幅值，電壓相位，類型% 類型： 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,1.0,0.0,-0.076

25、,-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;14,1.060,0.0,0

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

27、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,0;10,11,0.04699,

28、0.19797,0.0,0.0438,0;10,14,0.01938,0.05917,0.0,0.0528,0;輸出結(jié)果數(shù)據(jù)1 ：-第 1 次迭代時(shí)雅可比矩陣-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.578554-38.240787-0.000

29、000-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.000000-0.000000-0.00000010.

30、786262-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-0.000000-0.000000-0.000000 1.424005-5.3260553.902050-0.000000-0.000000-0.00000010.36539

31、4-14.7683384.402944-0.000000-0.000000-0.000000-0.0000001.880885-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000 -0.0000003.902050-5.782934-0.000000-0.000000-0.000000-0.0000004.402944-11.362870-0.000000-0.000000-0.000000-0.000000-2.467393-0.000000-0.0000006.959926-0.000000-0.0

32、00000-0.000000-0.000000-0.000000 -0.000000-0.0000001.880885-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-7.6511132.251975-0.000000-0.000000-0.000000-0.000000-2.9468155.399139-0.000000-0.000000-0.000000-0.0000002.489025 -0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000

33、-0.0000002.251975-14.9416222.314963-0.000000-0.000000-0.0000002.48902510.374684-0.000000-0.000000-0.000000-0.000000-4.555697 1.136994-0.000000-0.000000-0.000000-0.000000-0.0000003.029050-0.000000-0.000000-0.0000002.314963-5.344014-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000

34、00-0.0000001.136994 -2.5610001.424005-0.0000005.3460515.427654-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-0.000000-32.727644-0.0000005.047017-0.000000-0.000000-0.0000001.7619051.777691-0.000000-0.000000 -0.000000-0.000000-0.0000005.119505-0.000000-0.000000-0.000000-0.000000-0.000000-0.000

35、000-0.000000-0.0000005.047017-0.000000-10.166523-0.000000-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.0000006.9599265.39913910.374684-0.000000-0.0000003.323549-0.000000-29.479246-0.000000-0.000000-0.000000-0.0000002.594145 5.2

36、68177 -0.000000-0.000000-0.000000-0.000000-0.00000010.786262-0.000000-0.000000-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.000000-0.00000010.608721-6.8409810.000000 0.000000 0.000000 0.000000 0.000000 0.00

37、00000.000000-1.761905-2.005835-0.0000000.000000-0.0000000.000000-37.96863121.5785544.7819431.797979-0.000000-0.000000-0.000000-6.8409819.7061230.000000 0.000000 0.000000 0.000000 0.000000 0.0000000.000000-1.7776910.000000-0.0000000.000000-0.0000000.00000021.578554-31.542421-0.000000-0.000000-0.00000

38、0-0.000000-0.0000000.0000000.0000000.000000 0.0000000.0000000.000000 0.0000000.0000000.0000000.0000000.000000-0.0000000.000000-0.0000000.0000004.781943-0.000000-14.4397249.090083-0.000000-0.000000-0.0000000.0000000.0000000.000000 5.326055-3.902050 0.000000 0.000000 0.000000-1.4240050.0000000.000000-

39、0.0000000.000000-0.0000000.0000001.797979-0.0000009.090083-24.2825063.02905010.365394-0.0000000.0000000.0000000.000000 -3.9020505.782934 -1.8808850.000000 0.0000000.0000000.0000000.000000-0.0000000.000000-0.0000000.000000-0.000000-0.000000-0.00000010.365394-0.000000-14.7683384.4029440.0000000.000000

40、0.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.5941452.2519750.000000-0.

41、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.4240050.000000 0.000000 0

42、.000000 -1.1369942.5610000.0000000.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 3.029050 -0.0000002.314963 -5.344014-第1次迭代時(shí)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-0016.056626e-002-9.219354

43、e-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次迭代時(shí)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.721659e+000-4.249173e

44、-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 次迭代后相角（ 角度 ）44.11250048.9578926

45、8.16611180.81979290.848607114.299041125.877362123.89794698.64381324.34596735.398301128.70329568.1661110.000000-88.46859650.77006215.630499 4.956802 0.000000 0.000000 0.000000 0.000000 0.0000008.3512030.0000000.000000-17.67902420.9626216.9766783.695668-0.000000-0.00000053.574257-0.000000-70.163300-0.

46、0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000-0.0000001.8716490.00000012.117328-29.191523-0.000000-0.000000-0.000000-0.00000015.630499-0.000000-0.000000-0.000000-85.47526245.0445950.0000000.0000000.0000000.000000-0.000000-0.0000000.00000024.800168-6.976678-0.0000003.13630210.1

47、12981-0.0000004.956802-0.000000-0.000000-0.000000 -0.00000045.044595-111.55769648.1013580.0000000.0000000.000000-0.000000-0.0000000.000000-0.000000-3.695668-0.000000-10.112981-24.720611-第 2 次迭代時(shí)雅可比矩陣-8.760031-0.0000008.844924-0.0000000.000000-0.00000013.45494128.512426-0.000000-0.000000 -0.000000 12

48、.548251-0.000000-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.000000-0.000000-0.000

49、00014.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.000000-43.168982 21.

50、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.860379-0.000000-0.

51、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.0000000.0000000.0000007.4

52、04987-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.000000-0.000000-68.68

53、4204-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.9568020.0000000.0000000.000000-12.1173280.000000 0.00000017.294580 0.000000 0.000000 0.0