通过matlab计算卫星位置
作者：互联网
卫星星历是描述卫星运动轨道的信息。也可以说卫星星历就是一组对应某一时刻的轨道参数及其变率。有了卫星星历就可以计算出任意时刻的卫星位置及其速度。GPS卫星星历分为预报星历和后处理星历。预报星历又称广播星历。
GPS广播星历参数共有16个，其中包括1个参考时刻，6个对应参考时刻的开普勒轨道参数和9个反映摄动力影响的参数。这些参数通过GPS卫星发射的含有轨道信息的导航电文传递给用户。
1.星历参考时刻 ：
2.长半轴平方根：
3.偏心率：
4.参考历元下平近点角：
5.近地点角距：
6.轨道倾角：
7.本周初始历元的升交点赤经：
8.平运动差（由精密星历计算得到的卫星平均角速度与按给定参数计算所得的平均角速度之差）：n
9.轨道倾角变化率（弧度/秒）：
10.升交点赤经变化率（弧度/秒）：
11.纬度幅角的余弦调和项改正的振幅（弧度）：
12.纬度幅角的正弦调和项改正的振幅（弧度）：
13.轨道半径的余弦调和项改正的振幅（m）：
14.轨道半径的正弦调和项改正的振幅（m）：
15.轨道倾角的余弦调和项改正的振幅（弧度）：
16.轨道倾角的正弦调和项改正的振幅（弧度）：
以下为matlab程序：
%the homework 1 of chapater 1
%student ：Taylen
%time ：2018/9/25
%程序功能：根据所提供的星历参数,计算此卫星在信号发射时刻t( GPS时间)239050.7223s
%时的时空位置
%GPS广播星历参数共有16个，其中包括1个参考时刻，6个对应参考时刻的开普勒轨道参数
%和9个反映摄动力影响的参数。这些参数通过GPS卫星发射的含有轨道信息的导航电文传递给用户。
%时间参数
pra1 = 244800; %te---星历参考时刻，即星历表参考历元(s)
%开普勒六参数
pra2 = -1.064739758; %M0---按参考历元te计算的平近点角（弧度）
pra3 = 0.005912038265; %e---轨道偏心率
pra4 = 5153.65531; %sqrt(a)长半轴平方根
pra5 = -1.717457876; %w0---近地点角距
pra6 = 0.9848407943; %i0---按参考历元计算的轨道倾角(弧度)
pra7 = 1.038062244; %Ω0---本周初始历元的升交点赤经(弧度）
%轨道摄动九参数
pra8 = 4.249105564e-9; % Δ
标签：cos,轨道,位置,matlab,计算,星历,卫星,sin 来源： https://www.cnblogs.com/blogger-Li/p/11630155.html