?? satellitictrace.m
字號:
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)系
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -