satellitictrace.m

matlab实现的卫星轨道模拟

clear;clc;
%================================================%
%==============定义卫星轨道参数===================%
%================================================%
a=6880.140e+3;                                                                  %%轨道半长轴
e=0.003;                                                                        %%轨道扁率
%Orbit_i=97.423*pi/180;                                                          %%轨道倾角
Orbit_i=0*pi/180; 
Upw=-60*pi/180;                                                                 %%升交点赤经
Nearw=160*pi/180;                                                               %%近地点幅角
Ta=0;                                                                           %%假定为起始时刻(卫星飞经近地点时刻)，为绝对时间，这里取相对时刻为0
T0=0;                                                                           %%假设地球自转过春分点的时刻
Mui=3.986e14;                                                                   %%为常数，用于计算卫星运行周期     
Tp=2*pi*sqrt(a^3/Mui);                                                          %%卫星运行周期                                                                     %%成像的点目标数目
Tt=2000;                                                                        %%定义当前时刻               
Ntarget=7;
%==============================================%
%===========卫星在惯性坐标系下的坐标=============%
%==============================================%
Aov=TransMaov(Upw,Orbit_i,Nearw);
s=0:10:Tp;
L=length(s);
R0=6378.16e3;
for q=1:1:L
    [f,Orbit_xyz,r,s(q)]=Earthorbit(a,e,Tp,s(q),Mui,Ta);
    Rs_xyz=Aov*Orbit_xyz;
    x(q)=Rs_xyz(1);
    y(q)=Rs_xyz(2);
    z(q)=Rs_xyz(3);
end
figure;
grid on;
plot3(x,y,z,'b');
[X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z;
colormap(gray);
grid on,hold on,surf(X,Y,Z),shading interp

[f,Orbit_xyz,r,Tt]=Earthorbit(a,e,Tp,Tt,Mui,Ta);              %%计算卫星在轨道坐标系下的坐标
Rs_xyz=Aov*Orbit_xyz;
hold on;
plot3(Rs_xyz(1),Rs_xyz(2),Rs_xyz(3),'r*'); 
                                                            %%天线的相位中心相对于卫星星体坐标系