目录

• 1. 升交点赤经与地球自转的关系
• 2. 升交点经度 \Omega_k 的计算原理
• 3. 卫星位置计算步骤详解（基于广播星历）
  • 3.1 计算平均运动角速度
  • 3.2 计算平近点角与偏近点角
  • 3.3 计算真近点角与升交距角
  • 3.4 摄动改正与二体问题修正
• 4. 坐标转换：从轨道平面到地固坐标系
• 5. 极移与地球自转不均匀性的讨论
• 6. 精密星历（SP3）与内插计算
• AI 总结

---

1. 升交点赤经与地球自转的关系 原片 @ 00:00*

在计算 GPS 卫星位置时，理解参考系和时间参数至关重要，特别是关于 升交点赤经（Right Ascension of Ascending Node, RAAN） 的处理。

• \Omega_0 的特殊性：在广播星历提供的 6 个轨道根数中，唯独 \Omega_0 不是参考时刻 t_{oe} 的轨道根数，而是参考时刻 t_{oe} 的轨道根数减去本周开始时刻的格林尼治恒星时（GAST）。
• 地球自转的影响：
  • 春分点是惯性系的基准，而格林尼治子午线（本初子午线）随地球自转而转动。
  • 地球自转角速度 \dot{\Omega}_e 是一个常数（平均值）。
  • 如果已知本周开始时刻（Week start, t=0）春分点与格林尼治子午线的夹角，且已知地球自转速度，则可以推算出任何时刻两者的相对位置。

Screenshot-[01:06]

截图说明：视频展示了春分点、格林尼治子午线随时间变化的角度关系，解释了地球自转角速度 \dot{\Omega}_e 如何影响角度计算。

2. 升交点经度 \Omega_k 的计算原理 原片 @ 02:27*

为了方便计算卫星在地固坐标系（ECEF）中的位置，需要求出观测时刻 t 升交点相对于格林尼治子午线的经度 \Omega_k。

• 计算公式推导：
  任何时刻的升交点赤经 \Omega(t) 可以表示为：
  \Omega(t) = \Omega_{t_{oe}} + \dot{\Omega}(t – t_{oe}) 其中 \dot{\Omega} 是升交点赤经的变化率（由地球扁率等摄动因素引起，导航电文已给出）。
• 转化为地固系经度：
  观测时刻的升交点经度（在大地坐标系中）需减去当前的格林尼治恒星时（GAST）。
  \Omega_k = \Omega_{t_{oe}} + \dot{\Omega} \cdot \Delta t – (GAST_{week\_start} + \dot{\Omega}_e \cdot t_{total}) 为了简化计算，导航电文中的 \Omega_0 实际上已经包含了部分项。
• 最终计算公式：
  \Omega_k = \Omega_0 + (\dot{\Omega} – \dot{\Omega}_e)t_k – \dot{\Omega}_e t_{oe}
  • \Omega_0：导航电文给出的参数（参考时刻升交点赤经与本周起始时刻格林尼治恒星时之差）。
  • \dot{\Omega}：升交点赤经变化率。
  • \dot{\Omega}_e：地球自转角速度。
  • t_k：归化时间（观测时刻 – 参考时刻）。

Screenshot-[05:50]

截图说明：图示详细解析了角度关系：\Omega（大Omega）是从春分点到升交点的角度（赤经），而我们需要的是从格林尼治子午线到升交点的角度（经度 \lambda 或 \Omega_k）。图解展示了如何通过减去格林尼治恒星时来实现这一转换。

3. 卫星位置计算步骤详解（基于广播星历） 原片 @ 10:10*

利用导航电文提供的 16 个参数（6个轨道根数 + 9个摄动参数 + 参考时间 t_{oe}），计算卫星在 WGS-84 坐标系下的位置。

3.1 计算平均运动角速度 原片 @ 11:00*

• 首先利用开普勒第三定律计算平均角速度 n_0：
  n_0 = \sqrt{\frac{GM}{A^3}}
  • 注意：导航电文中给出的是 \sqrt{A}（长半轴的平方根），这是为了便于后续计算 A^{3/2}。
• 计算修正后的平均角速度 n：
  n = n_0 + \Delta n
  • \Delta n 是导航电文给出的摄动改正参数。

3.2 计算平近点角与偏近点角 原片 @ 12:49*

• 计算观测时刻的平近点角 M_k：
  M_k = M_0 + n \cdot t_k
  • M_0：参考时刻的平近点角。
  • t_k：归化时间（观测时刻 t – 参考时刻 t_{oe}）。
• 计算偏近点角 E_k（迭代求解开普勒方程）：
  M_k = E_k – e \sin E_k
  • 通常迭代 2-3 次即可收敛。

3.3 计算真近点角与升交距角 原片 @ 14:06*

• 真近点角 f_k (或 \nu_k)：
  利用几何关系或公式由 E_k 求得。
  \tan \frac{f_k}{2} = \sqrt{\frac{1+e}{1-e}} \tan \frac{E_k}{2}
• 升交距角 (Argument of Latitude) \Phi_k：
  \Phi_k = \omega + f_k
  • \omega：近地点角距（Argument of Perigee）。
  • 注意：\omega 和 f_k 都是在轨道平面内定义的角度。

3.4 摄动改正与二体问题修正 原片 @ 16:50*

GPS 卫星受摄动影响，轨道不是完美的椭圆。需要计算三个方向的周期性改正项（利用二倍角 2\Phi_k）：

1. 升交距角改正 \delta u_k：
\delta u_k = C_{uc} \cos(2\Phi_k) + C_{us} \sin(2\Phi_k) 2. 向径（距离）改正 \delta r_k：
\delta r_k = C_{rc} \cos(2\Phi_k) + C_{rs} \sin(2\Phi_k) 3. 轨道倾角改正 \delta i_k：
\delta i_k = C_{ic} \cos(2\Phi_k) + C_{is} \sin(2\Phi_k)

• S 结尾的参数对应 sine 项，C 结尾的参数对应 cosine 项。
• 计算改正后的参数：
  • 校正后的升交距角：u_k = \Phi_k + \delta u_k
  • 校正后的向径：r_k = A(1 – e \cos E_k) + \delta r_k
  • 校正后的轨道倾角：i_k = i_0 + \delta i_k + \dot{i} \cdot t_k

Screenshot-[17:42]

截图说明：展示了利用 C_{uc}, C_{us}, C_{rc}, C_{rs}, C_{ic}, C_{is} 进行周期性摄动改正的公式列表。

4. 坐标转换：从轨道平面到地固坐标系 原片 @ 19:28*

• 轨道平面坐标系 (x’_k, y’_k)：
  原点在地球质心（椭圆焦点），x’ 指向升交点，y’ 在轨道平面内垂直于 x’。
  x’_k = r_k \cos u_k y’_k = r_k \sin u_k

Screenshot-[19:56]

截图说明：轨道平面示意图，清晰标注了近地点、升交点、真近点角 f、近地点角距 \omega 以及卫星位置。

• 转换到地固坐标系 (ECEF)：
  需要进行两次旋转：
  1. 绕 X 轴旋转 -i_k（将轨道平面转到赤道面）。
  2. 绕 Z 轴旋转 -\Omega_k（将升交点转到相对于格林尼治子午线的位置）。
  最终坐标 (X_k, Y_k, Z_k) 公式：
  X_k = x’_k \cos \Omega_k – y’_k \cos i_k \sin \Omega_k Y_k = x’_k \sin \Omega_k + y’_k \cos i_k \cos \Omega_k Z_k = y’_k \sin i_k

Screenshot-[25:30]

截图说明：最终的坐标转换公式矩阵，展示了如何将轨道平面坐标 (x’, y’) 转换为大地坐标 (X, Y, Z)。

5. 极移与地球自转不均匀性的讨论 原片 @ 26:01*

• 广播星历的近似性：
  • 上述计算假设地球自转是均匀的（使用平均角速度）。
  • 未考虑极移（Polar Motion）的影响，即瞬时地球坐标系与协议地球坐标系（CTP）之间的微小差异。
  • 计算得到的坐标实际上是瞬时地球坐标系下的坐标。
• 高精度需求：
  • 对于高精度定位，需要考虑地球自转的不均匀性（UT1-UTC）和极移参数（x_p, y_p）。
  • 现代导航电文或精密处理软件会包含这些修正参数，将坐标转换到协议地球坐标系。

6. 精密星历（SP3）与内插计算 原片 @ 29:30*

除了广播星历，IGS 等机构还提供精密星历（Precise Ephemeris）。

• 数据形式：
  • 不是提供轨道根数，而是直接提供离散时间点（例如每 15 分钟一个）的卫星三维坐标 (X, Y, Z) 和速度。
• 内插方法（Interpolation）：
  • 由于观测时刻通常不在 15 分钟的整点上，需要进行内插。
  • 拉格朗日内插法（Lagrange Interpolation） 或其他高阶多项式拟合常被用于计算任意时刻的精密坐标。
  • 注意：对于 15 分钟间隔的数据，简单的线性内插精度不够，通常需要使用高阶多项式（如 9 阶或更高）以保证精度。

Screenshot-[31:31]

截图说明：讲解精密星历（如 SP3 格式）的数据结构，以及如何通过内插算法求得观测时刻的卫星位置。

---

AI 总结

本节课深入讲解了利用 GPS 广播星历计算卫星空间位置的详细数学模型。核心难点在于理解升交点经度 \Omega_k 的计算，它结合了轨道参数、地球自转角速度以及时间参考系的转换。

课程梳理了标准的 7 步计算法：

1. 计算平均角速度 n。
2. 计算平近点角 M_k。
3. 解开普勒方程求偏近点角 E_k。
4. 计算真近点角 f_k 和升交距角 \Phi_k。
5. 加入二阶摄动改正（\delta u, \delta r, \delta i）。
6. 计算轨道平面坐标 (x’, y’)。
7. 通过旋转矩阵转换至地固坐标系 (X, Y, Z)。

最后，课程简要对比了广播星历与精密星历（SP3），指出精密星历采用离散坐标点结合高阶内插（如拉格朗日内插）的方法，适用于更高精度的定位需求，并提及了极移和地球自转不均匀性对高精度解算的影响。