clear;clc;
%================================================%
%==============定義衛(wèi)星軌道參數(shù)===================%
%================================================%
a=6880.140e+3;                                                                  %%軌道半長軸
e=0.003;                                                                        %%軌道扁率
%Orbit_i=97.423*pi/180;                                                          %%軌道傾角
Orbit_i=0*pi/180; 
Upw=-60*pi/180;                                                                 %%升交點(diǎn)赤經(jīng)
Nearw=160*pi/180;                                                               %%近地點(diǎn)幅角
Ta=0;                                                                           %%假定為起始時刻(衛(wèi)星飛經(jīng)近地點(diǎn)時刻)，為絕對時間，這里取相對時刻為0
T0=0;                                                                           %%假設(shè)地球自轉(zhuǎn)過春分點(diǎn)的時刻
Mui=3.986e14;                                                                   %%為常數(shù)，用于計算衛(wèi)星運(yùn)行周期     
Tp=2*pi*sqrt(a^3/Mui);                                                          %%衛(wèi)星運(yùn)行周期                                                                     %%成像的點(diǎn)目標(biāo)數(shù)目
Tt=2000;                                                                        %%定義當(dāng)前時刻               
Ntarget=7;
%==============================================%
%===========衛(wèi)星在慣性坐標(biāo)系下的坐標(biāo)=============%
%==============================================%
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);              %%計算衛(wèi)星在軌道坐標(biāo)系下的坐標(biāo)
Rs_xyz=Aov*Orbit_xyz;
hold on;
plot3(Rs_xyz(1),Rs_xyz(2),Rs_xyz(3),'r*'); 
                                                            %%天線的相位中心相對于衛(wèi)星星體坐標(biāo)系