【整理】卫星真近点角、平近点角与偏近点角之间的计算

limbo_250

已于 2023-07-18 08:27:40 修改

阅读量1.5w

点赞数 29

文章标签： gps

于 2023-07-17 11:31:50 首次发布

本文链接：https://blog.csdn.net/qq_36107898/article/details/131720155

版权

文章详细阐述了GPS卫星轨道计算的关键步骤，包括卫星向径与偏近点角、真近点角的关系，平近点角的计算，以及如何通过平近点角推导出偏近点角和真近点角，运用了开普勒定律和几何关系。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

阅读了范伟，王贵文，刘发发的文章《GPS卫星轨道计算方法的研究》，将其中关于卫星真近点角、平近点角与偏近点角之间的计算步骤整理如下。

本文图片来源为参考文献截图。

0.概念

轨道根数：6个用于确定轨道位置及卫星位置的参数

卫星过近地点时刻 $t_{0}$ ：字面含义

真近点角 $\theta$ ：观测时刻t卫星和地球连线与近地点形成的夹角

平近点角M：假设卫星以平均角速度在圆轨道上运行（该圆轨道与椭圆轨道共中心，轨道半径等于椭圆轨道长半轴长度），观测时刻t该虚拟卫星和地球连线与近地点形成的夹角

偏近点角E：如图所示，过观测时刻t卫星所在位置P作长半轴垂线，交点为R，延长RP与圆轨道相交于点Q，Q点和圆心连线与近地点形成的夹角

图1 GPS卫星轨道平面图

1.卫星向径r与偏近点角E的关系

如上图，以地球所在位置为原点O建立直角坐标系。

设观测时刻卫星的坐标为P(x，y)，容易知道:

$x= a\cos E-c=a\left ( \cos E-e\right )$ (1)

$y=b\sin E=a\sqrt{1-e^{2}}\sin E$ (2)

因此向径： $r=\sqrt{x^{2}+y^{2}}=a\left (1-e\cos E \right )$ （3）

2.卫星向径r与真近点角θ的关系

由图1容易看到 $\cos \theta =\frac{OR}{OP}=\frac{O_{1}R-O_{1}O}{OP}=\frac{a\left ( \cos E-e \right )}{r}$ (4)

由（3）式知 $\cos E=\frac{1}{e}\left ( 1-\frac{r}{a} \right )$ ，代入（4）得：

$r=\frac{a\left ( 1-e^{2} \right )}{1+e\cos \theta }$ (5)

3.计算平近点角M

根据平近点角的定义容易知道

$M=M_{0}+n\left ( t-t_{oe} \right )$ (6)

其中： $t_{oe}$ 为星历参考时刻，由广播星历给出；

$M_{0}$ 为参考时刻的近点角；

$n$ 为观测时刻的平均角速度，由卫星的平均角速度经摄动改正计算得出。

4.通过平近点角M计算偏近点角E

根据开普勒定律，卫星矢径单位时间内扫过的面积恒定，即：

$\frac{\pi ab}{T}=\frac{\Delta A}{\Delta t}$ (7)

其中 $\Delta A$ 为卫星矢径在 $\Delta t$ 时间内扫过的面积，下面计算 $\Delta A$ ：

当卫星转过一个微小角度 $d\theta$ 时，矢径扫过的面积为

$dA=\frac{r^{2}}{2}d\theta$ (8)

将 $\theta =0$ 对应的观测时刻记为 $t_{0}$ ，那么 $t_{0}$ 至 $t$ 时刻卫星扫过的面积为：

$\Delta A=\int_{0}^{\theta }\frac{r^{2}}{2}d\theta$ (9)

对应的 $\Delta t=t-t_{0}$ (10)

将（5）式代入（9）中，并记 $p=a\left ( 1-e^{2} \right )$ ，可得：

$\Delta A=\frac{p^{2}}{2}\int_{0}^{\theta }\frac{1}{\left ( 1+e\cos \theta \right )^{2}}d\theta$ (11)

结合（7）式，有：

$\frac{\pi ab}{T}\left (t-t_{0} \right )=\frac{p^{2}}{2}\int_{0}^{\theta }\frac{1}{\left ( 1+e\cos \theta \right )^{2}}d\theta$ (12)

下面计算（12）式中的积分：

由式（3）（5）知：

$a\left (1-e\cos E \right )=\frac{a\left ( 1-e^{2} \right )}{1+e\cos \theta }$ （13）

式（13）两端求微分，整理后得：

$\frac{1}{1-e^{2}}\frac{\sin E}{\sin \theta }dE=\frac{d\theta }{\left ( 1+e\cos \theta \right )^{2}}$ (14)

由（2）式， $\sin E=\frac{y}{b}$ ，从图1可以看到， $\sin \theta =\frac{y}{r}$

结合（3）式 $r=a\left (1-e\cos E \right )$ ，(14)式可化为：

$\frac{1}{1-e^{2}}\frac{a\left ( 1-e\cos E \right )}{b}dE=\frac{d\theta }{\left ( 1+e\cos\theta \right )^{2}}$ (15)

将（15）代入（12）式，则：

$\frac{\pi ab}{T}\left (t-t_{0} \right )=\frac{p^{2}}{2}\int_{0}^{E}\frac{1}{1-e^{2}}\frac{a\left ( 1-e\cos E \right )}{b}dE$ （16）

而：

$\frac{p^{2}}{2}\int_{0}^{E}\frac{1}{1-e^{2}}\frac{a\left ( 1-e\cos E \right )}{b}dE=\frac{p^{2}}{2}\int_{0}^{E}\frac{a^{2}}{pb}\left ( 1-e\cos E \right )dE=\frac{ab}{2}\left ( E-e\sin E \right )$ （17）

将（17）代入（16），有：

$\frac{\pi}{T}\left ( t-t_{0} \right )=\frac{E-e\sin E}{2}$ (18)

根据平近点角的定义，易知：

$M=\frac{2 \pi }{T}\left ( t-t_{0} \right )$ ，代入（18）可得：

$M=E-e\sin E$ (19)

5.根据偏近点角E计算真近点角θ

由（1）（2）可知：

$\tan \theta =\frac{y}{x}=\frac{a\sqrt{1-e^{2}}\sin E}{a\left ( \cos E-e \right )}$ (20)

于是： $\theta =\arctan \frac{\sqrt{1-e^{2}}\sin E}{\cos E-e}$ (21)

参考文献：

范伟, 王贵文, 刘发发. 2015. GPS卫星轨道位置计算方法的研究. 山西师范大学学报(自然科学版).