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;
  
t=1:1:9601  
figure(1)
plot(t,h)