%?現(xiàn)需要在高度為600~700KM范圍內(nèi)，回歸周期為30天的太陽(yáng)同步軌道上布置一顆L波段的合成孔徑雷達(dá)衛(wèi)星。
%?1.?請(qǐng)?jiān)O(shè)計(jì)合適的軌道參數(shù)，并仿真對(duì)應(yīng)軌道下一個(gè)運(yùn)行周期內(nèi)星下點(diǎn)軌跡圖（緯度和經(jīng)度兩個(gè)圖）；
%?2.?計(jì)算多普勒中心頻率，多普勒調(diào)頻率與緯度的關(guān)系圖（設(shè)姿態(tài)角為0°，天線視角20°~60°，推薦35°）；
%?3.?按高斯投影方式，將星下點(diǎn)軌跡圖投影到地圖坐標(biāo)上。

clc
close?all
clear?all

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

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

%?軌道六根數(shù)
a?=?7007.45;%?軌道長(zhǎng)半軸（Km）
e?=?0.001;%?軌道偏心率
i?=?97.9?*?PI/180;%?軌道傾角（rad）
w?=?90?*?PI/180;%?軌道近地點(diǎn)幅角（rad）
RAAN?=?0;%?升交點(diǎn)赤經(jīng)

%?軌道真近心角f（rad），由以下公式計(jì)算
Ts?=?5;%?采樣時(shí)間間隔（s）
t?=?[0:Ts:fix（N_T*T）];%?采樣時(shí)間點(diǎn)
tp?=?T/4;%?衛(wèi)星經(jīng)過(guò)近地點(diǎn)的時(shí)刻，近1/4圓軌
n?=?sqrt（GGC/a^3）;%?衛(wèi)星平運(yùn)動(dòng)速度（rad）
M?=?n?*?（t?-?tp）;%?衛(wèi)星平近點(diǎn)角（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）

%%?%?%?%?%?