%?現需要在高度為600~700KM范圍內，回歸周期為30天的太陽同步軌道上布置一顆L波段的合成孔徑雷達衛星。
%?1.?請設計合適的軌道參數，并仿真對應軌道下一個運行周期內星下點軌跡圖（緯度和經度兩個圖）；
%?2.?計算多普勒中心頻率，多普勒調頻率與緯度的關系圖（設姿態角為0°，天線視角20°~60°，推薦35°）；
%?3.?按高斯投影方式，將星下點軌跡圖投影到地圖坐標上。

clc
close?all
clear?all

%%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?變量定義?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%

%?常量定義，具體計算見報告（表1、表2）
PI?=?3.1415926;
R_EARTH?=?6371;%?地球半徑（Km）
W_EARTH?=?7.29e-5;%?地球自轉速度（rad）
GGC?=?3.986e5;%?地心引力常數
L_WAVE?=?0.24;%?L波段?波長（m）
LAA?=?35?*PI/180;%?衛星天線視角
D?=?30;%?回歸周期（天）
N?=?444;%?衛星在一個回歸周期內運行的周期數
T?=?97.30?*60;%?衛星運行周期
N_T?=1;%?實驗觀測周期數，可延長觀察時間

%?軌道六根數
a?=?7007.45;%?軌道長半軸（Km）
e?=?0.001;%?軌道偏心率
i?=?97.9?*?PI/180;%?軌道傾角（rad）
w?=?90?*?PI/180;%?軌道近地點幅角（rad）
RAAN?=?0;%?升交點赤經

%?軌道真近心角f（rad），由以下公式計算
Ts?=?5;%?采樣時間間隔（s）
t?=?[0:Ts:fix（N_T*T）];%?采樣時間點
tp?=?T/4;%?衛星經過近地點的時刻，近1/4圓軌
n?=?sqrt（GGC/a^3）;%?衛星平運動速度（rad）
M?=?n?*?（t?-?tp）;%?衛星平近點角（rad）
f?=?M?+?（2*e?-?e^3/4）*sin（M）?+?1.25*e^2?*?sin（2*M）?+13/12?*?e^3?*?sin（3*M）;%?軌道真近心角（rad）

%%?%?%?%?%?