init2.asv

用matlab将弹道坐标系转换为WGS84中的直角坐标

clear; close all; clc;

 
 
 format ('compact');
 format ('long', 'g');
 b=[];
 l=[];
 h=[];
 [t,v,x,y,z]=txt2array();
  
 
for i=2:1:9601
    
    X= x(i);
    Z= y(i);
    Y=0;
    
%东北天发射坐标----WGS-84下的空间直角坐标系的转化
    [x2,y2,z2]=geo22cart(X,Y,Z);
    x6(i)=x2;
    y6(i)=y2;
    z6(i)=z2;
%WGS-84下的空间直角坐标系变为WGS-84下的大地坐标系'
   [b1, l1, h1] = cart2geo(x2, y2, z2);
   b(i) =b1;
   l(i) =l1;
   h(i) =h1;
   
end
  b(1)=39.9583*pi/180;
  l(1)=116.3090*pi/180;
  h(1)=0;
  
  
 b=b';
 l=l';
 h=h';
 [m1,n1]=size(b)
 m9={'motb'};
v9={'v1_m1'};
%因为不能识别m所以得自己手敲
t=0:1:119;
xlswrite('F:\v1-m.xls',t,'A1:A120');
m9={'motb'};
v9={'v1_m1'};
xlswrite('F:\v1-m.xls',t,'A1:A120');
xlswrite('F:\v1-m.xls', m9 ,'B1:B120');
xlswrite('F:\v1-m.xls', v9,'C1:C120');
xlswrite('F:\v1-m.xls',39.9583*pi/180,'D1:D120');
xlswrite('F:\v1-m.xls',116.3090*pi/180,'E1:E120');
xlswrite('F:\v1-m.xls',h,'F1:F120');

t=120:0.01:0.01*(m1-1);
xlswrite('F:\v1-m.xls',t,'A121:A9721');
xlswrite('F:\v1-m.xls', m9 ,'B121:B9721');
xlswrite('F:\v1-m.xls', v9,'C121:C9721');
xlswrite('F:\v1-m.xls',b,'D121:D9721');
xlswrite('F:\v1-m.xls',l,'E121:E9721');
xlswrite('F:\v1-m.xls',h,'F121:F9721');