drawsatelliteorbit.m

来自「上一个的继续,连续有九个压缩文件,大家可以打开解压后自己慢慢享用,这个都是本人调」· M 代码 · 共 37 行

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.6]; %卫星星座数据
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         0];
    R1=[1         0         0;
        0       cos(B)  -sin(B);
        0       sin(B)  cos(B)];
    R2=[cos(C)  -sin(C)     0;
        sin(C)     sin(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