drewsatelliteorbit.m

GPS 的地球卫星运转情况仿真 还可以生成GPS信号

function DrewSatelliteOrbit
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];%卫星星座数据，见图B-2
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;]; %将卫星轨道坐标系转换到地球赤道坐标
x1=Ans(1,:);
y1=Ans(2,:);
z1=Ans(3,:);
plot3c(x1,y1,z1,k);  %调用程序一
hold on; 
axis equal;
axis off;
end