这是做卫星通信的课程设计时涉及的卫星星下点的作图程序,matlab编写。

我们对卫星星下点的代码做了点修改，把地球自转也考虑进去了。

卫星星下点的运行轨迹如下：

卫星纬度050100150200

卫星经度250300350400

代码如下：（红色的为修改部分）

clc

close all

clear all

%% % % % % % % % % % % %

义 % % % % % % % % % % % % % % % % % % % % % %

PI = 3.1415926;

W_EARTH = 7.29e-5;% 地球自转速度(rad)

GGC = 3.986e5;% 地心引力常数

N_T =10;% 实验观测周期数，可延长观察时间

% 轨道六根数

a = 20000;% 轨道长半轴（Km）

e = 0.3;% 轨道偏心率

i = 30 * PI/180;% 轨道倾角（rad）

w = 70 * PI/180;% 轨道近地点幅角（rad）

RAAN = 10 * PI/180;% 升交点赤经

T = 2 * PI * sqrt((a^3) / GGC);

% % % % % 变量定