目录

• 1. 导航电文中的轨道参数概述
• 2. 具体的轨道根数与摄动参数
• 3. 卫星运动的几何关系与三个“近点角”
• 4. 卫星位置计算的核心步骤与公式
• AI 总结

---

1. 导航电文中的轨道参数概述 原片 @ 00:00*

Screenshot-[00:00]

GPS 导航电文（Navigation Message）并不直接提供卫星在每一瞬间的空间坐标 (X, Y, Z)，而是提供了一组轨道根数（Orbital Elements）和相关的常数。用户需要利用这组常数，计算出卫星在任意时刻 t 的空间位置。

1.1 参考时刻的概念

• 轨道参数的时效性：由于存在摄动力（Perturbation），轨道根数并非真正的常量，而是随时间缓慢变化的。因此，给出的轨道根数必须对应一个特定的时刻。
• 参考时刻 (t_{oe})：
  • 视频中提到的“701”实际上是指 t_{oe} (Time of Ephemeris)，即星历参考时刻。
  • 如果是卫星钟差的参考时刻，通常称为 t_{oc}。
  • 这两个参考时刻通常是不同的，允许分别定义。
• 计算目标：利用 t_{oe} 时刻的参数，外推计算出观测时刻 t 的卫星状态。

---

2. 具体的轨道根数与摄动参数 原片 @ 02:40*

Screenshot-[02:57]

导航电文通常提供 6 个基本的轨道根数以及 9 个摄动（修正）参数。

2.1 六个基本轨道根数

虽然经典开普勒轨道根数有 6 个，但在 GPS 广播星历中，为了计算方便和精度，参数形式略有调整：

1. M_0 (Mean Anomaly at Reference Time)：参考时刻 t_{oe} 的平近点角。
   • 注意：这是一个相对于 t_{oe} 的值。
2. \sqrt{A} (Square Root of Semi-Major Axis)：轨道长半轴的平方根。
   • Screenshot-[04:09]
   • 为什么要给平方根？ 为了用户后期计算方便，因为在计算平均角速度时需要用到 A^3 的相关项，给 \sqrt{A} 可以减少开方运算，提高计算效率。
3. e (Eccentricity)：轨道偏心率，描述椭圆轨道的形状。
4. i_0 (Inclination)：参考时刻的轨道倾角。
5. \Omega_0 (Longitude of Ascending Node)：参考时刻的升交点赤经。
   • i_0 和 \Omega_0 共同确定了轨道平面在空间中的位置。
6. \omega (Argument of Perigee)：近地点角距。

2.2 九个摄动（修正）参数

为了修正地球引力非球形、日导致引力等摄动影响，还提供了以下参数：

• \Delta n (Mean Motion Difference)：平均角速度的修正值。
  • Screenshot-[05:51]
• \dot{\Omega} (Rate of Right Ascension)：升交点赤经对时间的变化率。
• \dot{i} (Rate of Inclination)：轨道倾角对时间的变化率。
• 调和改正项 (Harmonic Correction Terms)：
  • 用于对三个参数进行周期性修正：升交点角距 (Argument of Latitude u)、地心距离 (Radius r)、轨道倾角 (Inclination i)。
  • 参数符号通常包含 c (cosine) 和 s (sine)：
    • C_{uc}, C_{us}：升交点角距的余弦/正弦调和改正振幅。
    • C_{rc}, C_{rs}：地心向径的余弦/正弦调和改正振幅。
    • C_{ic}, C_{is}：轨道倾角的余弦/正弦调和改正振幅。

---

3. 卫星运动的几何关系与三个“近点角” 原片 @ 07:12*

为了计算卫星在轨道上的位置，核心是确定卫星与近地点的角度关系。这里涉及三个关键的角度概念。

3.1 真近点角 (f 或 \nu)

• Screenshot-[08:08]
• 定义：卫星与地心连线，和近地点方向之间的夹角。
• 作用：这是描述卫星真实位置的最直接角度，计算坐标时最终需要用到它。
• 困难：真近点角的变化极不均匀，直接根据时间计算非常困难。

3.2 开普勒第二定律（面积律）

• Screenshot-[09:21]
• 卫星在单位时间内扫过的面积是常数。
• 在近地点（Perigee）运动角速度快，在远地点（Apogee）运动角速度慢。
• 这导致无法简单地用线性公式计算真近点角。

3.3 偏近点角 (E)

• Screenshot-[11:14]
• 几何定义：
  1. 以轨道中心为圆心，分别以长半轴 a 和短半轴 b 为半径作辅助圆（大圆和小圆）。
  2. 过卫星位置作垂直于长轴的垂线，延长与大圆相交。
  3. 该交点与中心的连线，与近地点方向的夹角即为偏近点角 E。
• 作用：作为中间变量，连接真近点角和平近点角。

3.4 平近点角 (M)

• Screenshot-[12:43]
• 定义：假设一个虚拟点在以长半轴为半径的圆上以匀角速度（平均角速度 n）运动。
• 特点：与时间呈严格的线性关系，计算非常简单。
• 公式：
  M_t = M_0 + n(t – t_0) 其中 n 为平均角速度。

---

4. 卫星位置计算的核心步骤与公式 原片 @ 15:49*

计算思路：时间 t \rightarrow 平近点角 M \rightarrow 偏近点角 E \rightarrow 真近点角 f。

步骤一：计算平均角速度 n_0

• Screenshot-[16:00]
• 利用开普勒第三定律导出：
  n_0 = \sqrt{\frac{\mu}{A^3}}
  • \mu：地球引力常数 (GM)，约为 3.986005 \times 10^{14} m^3/s^2（WGS-84标准）。
  • A：轨道长半轴（通过广播星历中的 \sqrt{A} 平方得到）。

步骤二：修正平均角速度 n

• Screenshot-[18:00]
• 计算得到的 n_0 是参考时刻的值，实际平均角速度需加上摄动改正值 \Delta n：
  n = n_0 + \Delta n

步骤三：计算观测时刻的平近点角 M_k

• Screenshot-[18:55]
• 计算观测时刻 t 与参考时刻 t_{oe} 的时间差 t_k = t – t_{oe}。
• 计算平近点角：
  M_k = M_0 + n \cdot t_k

步骤四：解开普勒方程求偏近点角 E_k

• Screenshot-[20:13]
• 开普勒方程：
  M_k = E_k – e \sin E_k
• 求解方法：这是一个超越方程，无法直接解析求解 E_k。通常采用迭代法（如牛顿迭代法）。
  • 初值设为 E_0 = M_k。
  • 迭代公式：E_{j+1} = M_k + e \sin E_j。
  • 由于 GPS 卫星轨道偏心率 e 很小（约 0.002-0.003），收敛速度非常快，通常迭代 3-4 次即可满足精度要求（如误差小于 10^{-9}）。

步骤五：计算真近点角 f (或 \nu)

• Screenshot-[22:24]
• 利用几何关系将 E_k 转换为真近点角 f_k。视频中提到的公式（半角公式形式）：
  \tan \frac{f_k}{2} = \sqrt{\frac{1+e}{1-e}} \tan \frac{E_k}{2} 或者直接使用反正切函数形式：
  f_k = \arctan \left( \frac{\sqrt{1-e^2} \sin E_k}{\cos E_k – e} \right)

通过以上步骤，确定了卫星在轨道平面内的角度位置，结合后续的轨道幅角和升交点赤经计算，即可最终求得卫星在地球固连坐标系下的三维坐标。

---

AI 总结

本节视频详细讲解了 GPS 全球定位系统中卫星位置的计算原理，核心在于利用导航电文提供的轨道根数（Ephemeris Data）推算任意时刻的卫星位置。

1. 数据基础：导航电文不直接广播坐标，而是广播参考时刻 t_{oe} 的 6 个基本轨道根数（如 \sqrt{A}, e, i_0 等）和 9 个摄动参数。
2. 核心难点：卫星运动遵循开普勒第二定律（面积律），角速度不均匀，无法直接由时间计算真近点角。
3. 解决方案：引入平近点角 M（匀速、线性相关）和偏近点角 E（几何过渡）作为中间变量。
4. 计算流程：
   • 由 \sqrt{A} 和 \Delta n 计算修正后的平均角速度 n。
   • 推算观测时刻的平近点角 M_k。
   • 通过迭代求解开普勒方程 (M = E – e \sin E) 得到偏近点角 E_k。
   • 最后将 E_k 转换为真近点角 f，从而确定卫星在轨道上的具体位置。

这套算法是 GPS 接收机进行定位解算的数学基础，属于研究生复试中卫星导航课程的重点考察内容。