潮流计算matlab

1、matlab 潮流计算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;%打开数据文件%节点重新编号% 平衡节点数目% PV 节点数目% PQ 节点数目% nb 为总节点数% 统计平衡节点数目SW(nSW,:)=bus(i,:);elseif type=2,nPV=nPV+1; % 统计 PV 节点数目PV(nPV,:)=bus(i,:);elsenPQ=nPQ+1; % 统计PQ节点数目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); % 形成节点导纳矩阵K=0;%迭代次数初值Kmax=10;%最大迭代次数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%计算相角bus(i,3)=bus(i,3)-dX(i,1);endB=dx(nx,n:mx)*Um; %计算电压差 bus(1:m,2)=bus(1:m,2)-B'%计算电压值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,'第次迭代后节点电压(仅PQ节点)',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) %判断是否达到计算精度 break;endend% 计算功率P1=0;T=0;%计算平衡节点的有功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%计算 PV 节点、平衡节点的无功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;%换回各节点、支路的初始编号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); %计算支路潮流及损耗%将节点导纳矩阵、节点潮流计算结果写入文件output2myf=fopen('output2.dat',&#

10、39;w');fprintf(myf, '- 节点导纳矩阵-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,牛顿拉夫逊法潮流计算结果n');fprintf(myf, '-节点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节点计算结果-节点电压', bus(l,j);-n');节点相角', bus(l,j)*180/pi);', bus(l,j);注入有功功率(P)注入无功功率 (Q)类型-n');endfprintf(myf, ' n'); endfprintf(myf, '- 支路计算结果fprintf(myf,节

12、点 (I)-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);% 根据支路参数建立节点导纳矩阵程序：function Y=y(bus,line)% 目的：根据支路参数建立节点导纳矩阵%输入：节点参数矩阵-bus ；支路参数矩阵-line% 输出：节点导纳矩阵-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)%计算 PQ、 PV 节点的有功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%计算PV、 PQ 节点的无功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%计算 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%计算 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;% 计算 dPQ 的程序：function dPQ=dPQ(Y,bus)%delP-有功偏差量%delQ-无功偏差量%Y- 节点导纳矩阵%bus-节点参数( 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 节点)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 节点 )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;% 计算线路损耗、线路潮流程序：function Pf=loss(bus,line)% 计算线路损耗、线路潮流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% 输入的参数数据：% data for test case%各节点参数：节点编号,注入有功,注入无功 ,(Sn=100MV A) 电压幅值 ,电压相位 ,类型% 类型： 1=PQ 节点 ,2=PV 节点 ,3=平衡节点% (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;%各支路参数：起点编号 ,终点编号 , 电阻,电抗,电导,电纳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;输出结果数据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 次迭代后节点电压（仅 PQ 节点 ）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