共计 6690 个字符，预计需要花费 17 分钟才能阅读完成。
在前面的文章中，我们详细介绍了低轨卫星的TLE两行参数每个字段的含义、利用TLE参数计算得到卫星位置的精度、TLE参数的误差对地面低轨卫星终端天线对星的影响，以及TLE参数到底是怎么观测、数据融合得到的。前面的文章列表如下：
- [解析低轨卫星轨道星历的TLE两行参数](https://mp.weixin.qq.com/s/R_EXqja6z8qpsrzAlj24tA)
- [低轨卫星轨道星历TLE参数的精度有多高？](https://mp.weixin.qq.com/s/6pl77j4JHOW-i0wklWQY4A)
- [低轨卫星轨道参数 TLE 到底是怎么产生的？](https://mp.weixin.qq.com/s/CVV_CmhBBOudCuFHYkc14Q)
- [低轨卫星终端天线的对星，仅仅采用 TLE 参数为什么不行？](https://mp.weixin.qq.com/s/pWAVIt_jgDyWkVGlK3HicQ)
我们知道在GNSS卫星导航领域，GNSS卫星（比如北斗、GPS）也会播发轨道星历，地面的GNSS接收机接收后可解算卫星位置。那么，**低轨卫星的TLE参数**与**GNSS导航卫星的广播星历**之间有哪些本质差异？它们的生成方式、数学模型、坐标基准、精度水平及其随时间退化特性如何？
今天，我们就来回答这些问题。
---
### 一、基本概念：什么是TLE参数与GNSS广播星历
#### **1.1 TLE（Two-Line Element Set，两行轨道要素集）**
TLE是一种由北美防空司令部（NORAD）维护的标准化轨道数据格式，用于描述人造地球卫星在某一时刻的平均轨道状态。它包含以下核心参数：
- 轨道倾角、升交点赤经、偏心率、近地点幅角、平近点角
- 平均运动速率（决定半长轴）
- B*阻力系数（经验摄动项）
TLE由全球雷达与光电观测站网观测数据融合，经SGP4/SDP4模型拟合生成。其输出的是**经长期摄动平均化处理后的“平均轨道根数（Mean Orbital Elements）”**，而非某时刻的瞬时几何。
TLE不是描述“真实物理轨道”，而是驱动一个**半解析轨道预测模型（SGP4）**，实现快速估算任意时刻的**近似瞬时位置/速度**。
#### **1.2 GNSS广播星历（Broadcast Ephemeris）**
导航卫星（如北斗、GPS）通过无线电播发轨道参数组，由分布全球的监测站实时测量卫星伪距，主控站解算后上传至卫星再向用户广播。
- 目标：支持米级定位，要求**空间信号测距误差（SISRE） ≤ 2 m**
- 形式：一组基于开普勒轨道模型 + 摄动改正项的15–16个参数包
- 坐标系：直接对应地心地固坐标系ECEF（CGCS2000 / WGS-84）
广播星历反映的是某一参考时刻的**瞬时轨道状态（Osculating Elements）的精密拟合值**，适用于短时间内高精度外推。
---
### 二、参数类型差异：平均根数 vs. 瞬时根数
这是两者最根本的数学区别。
#### **2.1 瞬时根数（Osculating Elements）——GNSS星历的本质**
- **原理**：在某一特定时刻，假设所有的干扰力（摄动）突然消失，卫星将仅受地球中心引力作用。此时卫星的位置和速度所对应的那个开普勒椭圆，就是“瞬时轨道”。
- **特点**：它随着卫星的运动而**不停变化**。每一秒钟，因为受力情况微弱改变，对应的瞬时椭圆都在抖动。
- **例子**：GPS/北斗卫星播发的广播星历，本质上就是对一段时间内瞬时轨道的高度精准拟合。
下表根据《北斗卫星导航系统公开服务信号接口控制文件 B1I (3.0版)》整理的星历参数。
| 参数名称 | 符号 | 比特数 | 最小步长（LSB） | 单位 | 物理意义 |
| ------------------------ | -------------- | ------ | -------------------------- | ---------- | ----------------------- |
| 星历参考时间 | toe | 17 | 8 秒 | s | 参数生效起始时刻（UTC） |
| 长半轴平方根 | √A | 32 | 2-19 米1/2 | √m | 决定轨道大小 |
| 偏心率 | e | 32 | 2-33 | —— | 轨道扁率 |
| 近地点幅角 | ω | 32 | 2-31 × π 弧度 | rad | 轨道指向 |
| 平均运动改正量 | Δn | 16 | 2-43 × π 弧度/秒 | rad/s | 对名义 n 的微调 |
| 参考时间平近点角 | M0 | 32 | 2-31 × π 弧度 | rad | 参考时刻卫星相位 |
| 升交点赤经（参考时间） | Ω0 | 32 | 2-31 × π 弧度 | rad | 轨道平面方向 |
| 升交点赤经变化率 | Ω̇ | 24 | 2-43 × π 弧度/秒 | rad/s | 固体潮+重力场致旋进 |
| 轨道倾角（参考时间） | i0 | 32 | 2-31 × π 弧度 | rad | 轨面倾斜程度 |
| 轨道倾角变化率 | IDOT | 14 | 2-43 × π 弧度/秒 | rad/s | 季节性摄动 |
| 纬度幅角余弦调和改正振幅 | Cuc | 18 | 2-31 × π 弧度 | rad | 轨道偏心率周期性修正 |
| 纬度幅角正弦调和改正振幅 | Cus | 18 | 2-31 × π 弧度 | rad | 同上 |
| 轨道半径余弦调和改正振幅 | Crc | 18 | 2-6 米 | m | 轨道半径调制 |
| 轨道半径正弦调和改正振幅 | Crs | 18 | 2-6 米 | m | 同上 |
| 轨道倾角余弦调和改正振幅 | Cic | 18 | 2-31 × π 弧度 | rad | 倾角微小振荡 |
| 轨道倾角正弦调和改正振幅 | Cis | 18 | 2-31 × π 弧度 | rad | 同上 |
> 📝 说明：所有角度相关参数的实际存储是相对于 π 的分数，例如 LSB = 2-31π 表示一个比特单元代表 ≈ 1.455 × 10-9 弧度。
---
#### **2.2 平均根数（Mean Elements）——TLE的核心**
- **原理**：通过解析方法（如德普里特变换或李变换）对瞬时根数进行时间域平均，去除短周期摄动项（如由 J2 引起的日周期振荡），保留长期演化的趋势部分。
- **特点**：平滑、稳定，适合构建长期传播模型。
- **应用**：TLE中的7个动力学参数均为**平均轨道根数**，需输入SGP4模型计算，才能还原出近似的瞬时位置。
> 📘 举例：就像用一条平滑曲线拟合每日气温波动数据——虽然看不到每小时的变化，但能较好预测季节走向。
TLE不是描述“真实物理轨道”，而是驱动一个**半解析轨道预测模型（SGP4）**，实现快速估算任意时刻的**近似瞬时位置/速度**。
TLE包含7个核心参数：倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均运动速率，以及一个低轨卫星特有的——**B*阻力系数**。这个B*系数非常关键。低轨卫星在300–1500公里的高度运行，大气虽然稀薄，但“拖拽”效应明显。B*就是为了拟合这种效应而存在的。
| 参数名称 | 符号 | 说明 |
| :------------------------------- | :------- | :----------------------------------------------------------- |
| **历元 (Epoch)** | T0 | 该组数据的确切时间点（**UTC时间**）。所有预测都是以此为基准进行外推。 |
| **轨道倾角 (Inclination)** | i | 轨道平面与地球赤道平面的夹角（0°~180°）。 |
| **升交点赤经 (RAAN)** | Ω | 轨道平面在赤道面上的投影与春分点方向的夹角。 |
| **偏心率 (Eccentricity)** | e | 描述轨道的形状（LEO 卫星通常接近 0，即接近圆）。 |
| **近地点幅角 (Arg. of Perigee)** | ω | 从升交点到近地点之间的夹角，决定轨道在平面内的指向。 |
| **平近点角 (Mean Anomaly)** | M | 描述历元时刻卫星在轨道上的位置。它是一个**计算用的角度**，需要通过开普勒方程转换为**真近点角**才代表卫星在真实椭圆轨道上的实际角度位置。 |
| **平均运动 (Mean Motion)** | n | 卫星每天绕地球飞行的圈数（决定了半长轴 a，**圈数越多，轨道越低**）。 |
| **B* 阻力系数 (B-Star)** | B* | **关键模型参数**，它是SGP4/SDP4轨道模型中使用的**经验阻尼项**，综合反映了大气阻力、太阳光压等多种摄动力对卫星轨道变化的总体影响，通过模型拟合得到。 |
---
### 三、数学模型差异：简化分析模型 vs. 精密数值拟合
### 3.1 TLE参数，SGP4模型
SGP4 是一个**简化**的、**半解析**的数学模型，用于表示近地轨道 (LEO) 卫星的轨道。它基于重力扰动和一些主要的非保守力（如大气阻力和太阳光压）的影响，但对这些影响进行了**近似和简化**。
SGP4 在模型设计时就考虑到了简化，以便于快速计算和传播轨道。它通过引入一些经验公式和系数来近似复杂的摄动。
核心是半解析方法，使用初始轨道元素（TLE）通过数值积分模拟卫星位置。TLE包括6个经典轨道元素（半长轴、偏心率、倾角等）加上B星（拖曳系数）和均近点角。
假设主要摄动为地球非球形引力（J2-J4项）、大气拖曳和大太阳辐射压，不考虑更高阶效应。
#### 3.2 导航星历：精密的积分与拟合
GNSS导航星历则运行一套标准化的轨道外推算法。 该算法基于开普勒轨道运动模型，并加入了二阶谐项改正。地面主控站会用高性能计算机进行数值积分，解算复杂的摄动方程。然后，将这些复杂的数值解，通过一套标准化的“开普勒+摄动改正”公式包起来，变成播发给用户的15参数，采用的模型比 SGP4 **更复杂、更精细**。
由于广播星历的有效期通常较短（1–2小时），这套拟合模型能保证在有效期内位置误差呈线性增长且幅度极小。
GNSS接收机需严格遵循ICD中定义的计算次序。 这包括平近点角计算、开普勒方程求解、升交点角距修正等步骤，就能算出米级精度的位置。
### 3.3 对比表格
| 特性 | SGP4 模型 (TLE) | GNSS 星历模型 |
| :----------- | :--------------------------------------- | :----------------------------------------------------------- |
| **模型类型** | 简化半解析模型 | 精确数值积分或高精度解析模型 |
| **复杂性** | 较低 | 很高 |
| **主要应用** | LEO 轨道跟踪, 空间碎片监测, 业余无线电 | 精密导航定位, 大地测量, 科学研究, 精密定轨 |
| **计算速度** | 快 | 慢 |
| **数据量** | 小 (TLE 文件) | 大 (RINEX 等星历文件) |
| **轨道类型** | 主要适用于 LEO | 适用于 LEO, MEO, GEO 等 |
| **摄动考虑** | 简化地球引力, 经验大气阻力, 简化其他摄动 | 高阶地球引力, 精确月地/日地引力, 精确大气阻力, 太阳光压, 相对论效应, 等 |
---
### 四、坐标基准差异：TEME vs. ECEF
这是一个极易引发工程错误的关键点。
| 项目 | TLE/SGP4 输出 | GNSS广播星历 |
| -------------- | ------------------------------------------------------------ | ---------------------------------------------------------- |
| 原生坐标系 | 近似惯性系（传统称 TEME: True Equator, Mean Equinox）
Z轴为协议地极方向，X轴为协议春分点（GCRS框架近似） | 地心地固系（ECEF）
如 CGCS2000（中国）、WGS-84（美国） |
| 是否随地球旋转 | ❌ 否（固定于惯性空间） | ✅ 是（与地球同步旋转） |
| 用户是否需转换 | ✅ 是！必须经过坐标转换→ 转为ECEF，才能获得经纬度 | ❌ 否！直接可用经纬高表示 |
> 🔧 转换风险提示：若忽略地球自转角速度（7.292 × 10-5 rad/s），时间对齐误差1秒，低轨卫星经度偏差可达 **约460米**。务必确保UTC时间同步精度优于1毫秒。
---
### 五、位置精度表现与随时间退化特性
#### 5.1 绝对位置精度的量化对比
根据实测数据，TLE 在历元时刻（即数据刚生成的时刻）的典型绝对位置误差在 **1 公里至 3 公里** 之间。这一误差水平由观测站的几何分布和观测噪声决定。
相比之下，北斗广播星历的空间信号用户测距误差（SISRE）通常优于 **2 米**。在定位精度要求较高的场景下，两者的量级差异达到了三个数量级。
| 数据源 | 典型空间位置误差（RMS） | 数据基础 |
| ---------------------- | ----------------------- | -------------------------------------- |
| TLE（NORAD/CelesTrak） | **1 km ~ 3 km** | 受观测几何、数据融合质量、轨道高度影响 |
| 北斗B1I广播星历 | **≤ 2 m**（SISRE指标） | 主控站精密跟踪+高频更新 |
#### 5.2 精度随时间的漂移
时间是精度的敌人。
**TLE 的误差**随时间呈非线性增长。对于星链（Starlink）这类 550 公里高度的卫星，大气阻力是误差的主要来源。12 小时后的预报误差可能增加到 5 公里。这是由于 B* 系数无法完全预测大气密度受太阳活动影响而产生的随机波动。
如果超过 24 小时未获取新数据，摄动力的累积效应会导致预测轨道彻底偏离真实的物理路径。
**GNSS广播星历**则通过高频更新来规避这一问题。北斗系统的星历更新频率为每 1 小时一次。GPS 的更新频率为每 2 小时一次。GNSS接收机通过这种“接力”式的数据更新，将轨道拟合误差始终控制在米级阈值内。如果12小时没有更新，轨道拟合误差会达到几十米。
下表总结了两类轨道参数的精度退化特性：
| 参数类型 | 初始精度 | 12小时后误差 | 24小时后误差 | 更新机制 |
|----------|----------|---------------|---------------|----------|
| TLE | 1–3 km | ≈ 5 km | > 10 km（不可靠） | 每8–12小时人工/自动更新 |
| GNSS广播星历 | ≤ 2 m | < 5 m | > 30 m（失效） | 每1–2小时自动播发新参数 |
---
### 六、总结对比表
为了直观理解，我们把两类数据的核心参数做成了表格。
| **特性** | **TLE (低轨卫星常用)** | **导航广播星历 (北斗/GPS)** |
| ------------ | ------------------------------------------------------------ | ------------------------------------------------------------ |
| **根数类型** | **平均根数**：7个核心：倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均运动速率、**B* 阻力系数**。 | **瞬时根数**：15–18个参数：6个开普勒根数 + **6个谐项改正数**（Cuc, Cus等）+ 3个变化率（İ, Ω̇, Δn）。 |
| **数学模型** | SGP4（解析模型），基于“平均根数” | 开普勒 + 谐项改正（拟合模型）：基于“瞬时根数” |
| **坐标系** | TEME（准惯性系），需要经过转换才能得到地心地固坐标系中的经纬度 | ECEF（如 CGCS2000 / WGS-84）：地心地固坐标系，直接对应地面经纬度 |
| **初始精度** | 1–3 公里（历元时刻） | 1–2 米（星历参考时刻） |
| **更新频率** | 8–12小时（比如，星链卫星的星历数据（包含 TLE）通常**每 8 小时更新一次**） | 1 小时（北斗），2小时（GPS） |
| **典型用途** | 粗对星、多普勒预偏、可见性预测 | GNSS定位、授时 |
---
### 总结
本文详细探讨了低轨卫星TLE参数与GNSS卫星广播星历之间的差异，包括它们的生成方式、数学模型、坐标基准、精度水平及其随时间退化特性。TLE参数主要用于低轨卫星的轨道预测，其精度较低且随时间退化较快，而GNSS广播星历则提供高精度的卫星位置信息，适用于高精度定位和授时。通过对比，我们可以更好地理解这两种轨道参数的特点和适用场景，从而在实际应用中做出合理的选择。
**码字不易，大家觉得文章不错，就安排一下一键三连吧：点赞、分享、推荐。**
**大家有业务咨询，技术交流的，可以添加个人微信 15701068503。**
**“通信与导航”公众号主页有技术文章的分类合集，感兴趣的朋友可以查询。**
