drawsatelliteorbit.m

来自「标准GPS星座的仿真程序。注：采用24颗卫星」· M 代码 · 共 34 行

function [x1,y1,z1]=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;
    %六个轨道平面
    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];  %将卫星轨道坐标数据转换到地球赤道坐标
    x1=Ans(1,:);
    y1=Ans(2,:);
    z1=Ans(3,:);
    plot3c(x1,y1,z1,k);
    hold on;
    axis equal;
    axis off;
end