drawsatelliteorbit.m

来自「绘制GPS导航系统中卫星轨道的仿真程序」· M 代码 · 共 44 行

function DrawSatelliteOrbit
%绘制卫星轨道

a=26560;   %卫星椭圆轨道的长半轴
e=0.02;    %椭圆的偏心率
E=[0:0.1:2*pi];
x=a*(cos(E)-e);
y=a*sqrt((1-e^2))*sin(E);
z=0*E;

DtoR=2*pi/360;
A1=[32.8   92.8   152.8   212.6    272.8   332.8];    %卫星星座数据（即星座的平面投影）
for k=1:6
    A=A1(k)*DtoR;     %升交点的经度
    B=55*DtoR;        %轨道倾角
    C=pi/100;
    %总共有6个卫星轨道平面
    R3=[cos(A)   -sin(A)   0;
        sin(A)    cos(A)   0;
        0         0        1];
    
    R1=[1       0       0;
        0    cos(B)   -sin(B);
        0     sin(B)   cos(B)];
    
    R2=[cos(C)   -sin(C)   0;
        sin(C)    cos(C)   0;
        0         0        1];
    
    L1=length(E);
    R312=R3*R1*R2;
    Ans=R312*[x;y;z];     %将卫星轨道坐标系转换到地球赤道坐标

    plot3c(Ans(1,:),Ans(2,:),Ans(3,:),k);
    hold on;
    axis equal;
    axis off;
    
end