drawsatelliteorbit.m

来自「该程序用于计算卫星轨道和绘制卫星运行轨道轨迹图。」· M 代码 · 共 33 行

function DrawSatelliteOribt
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;
DotR=2*pi/360;
A1=[32.8 92.8 152.8 212.6 272.8 332.8];
for k=1:6
    A=A1(k)*DotR;
    B=55*DotR;
    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