# ESA PB05 式 Brouwer–Lyddane 平均根数轨道传播器（J2–J6 + 日月）

以 Brouwer–Lyddane / PB05 传统的半分析平均根数模型传播近地与 GNSS 区域轨道：EGM96 带谐场 J2–J6 下的轨道平均世俗漂移（Gauss 变分方程数值平均——Delaunay 正则化的数值等价物，并以闭式 J2 率内置交叉校验）、一阶 J2 短周期平均↔密切重构、J3 长周期冻结轨道分析（e = −J3·sin i/(2·J2·p)、ω = 90°）、椭圆倾斜摄动轨道上的日月四极双重平均世俗效应（含 39.23° 临界互倾角的 Celletti/Kozai–Lidov 项）、地面轨迹共振检查（1:1、2:1、3:2、5:2 及小有理公度）、太阳同步倾角求解，以及报告内 RK4 数值传播对平均模型的交叉验证。

> 标准页面: https://elysiatools.com/zh/tools/esa-pb05-pce-pck-reddy-porter-keplerian-equation-averaged-element-propagator

- **分类:** Astronomy

- **关键词:** 平均根数传播器, brouwer lyddane 平均根数, j2 世俗率, 升交点漂移, 冻结轨道, kozai lidov, 日月摄动, 太阳同步轨道, 田谐共振, 平均转密切

## 概述

传播器用 RK4 积分平均根数 (a, e, i, Ω, ω, M)，右端项是对所选深度带谐场的 Gauss 变分率轨道平均——即 Brouwer–Lyddane/PB05 Delaunay 正则化的数值等价，并用闭式一阶 J2 率交叉校验（Ω̇ = −3/2·J2·n·(Re/p)²·cos i，ω̇ = 3/4·J2·n·(Re/p)²·(5cos²i−1)）。日月选项叠加四极阶双重平均第三体项（同时对卫星与摄动轨道平均），能再现 GEO 约 0.75–0.95°/年的经典倾角增长，并携带 Kozai–Lidov 临界互倾角 arccos(√(3/5)) = 39.23°。短周期 J2 偏差由瞬时率积分得到，可加回采样密切根数（平均→密切反变换）。报告还运行同力模型的独立 RK4 笛卡尔数值传播并给出两者偏差、地面轨迹共振分类、冻结轨道（J3 长周期）参数与给定几何的太阳同步倾角解。

## 输入项

- **半长轴 a（公里）** (number): 8000
- **偏心率 e** (number): 0.05
- **倾角 i（度）** (number): 52.5
- **升交点赤经 Ω（度）** (number): 30
- **近地点幅角 ω（度）** (number): 70
- **平近点角 M（度）** (number): 0
- **传播时长（天）** (number): 30
- **带谐引力深度** (select)
- **第三体摄动** (select)
- **加平均→密切重构面板** (checkbox)

## 适用场景

- 需要评估近地轨道（LEO）、中圆轨道（MEO）或静止轨道（GEO）在数月至数年尺度下的升交点漂移与近地点进动时
- 设计太阳同步轨道或 J3 长周期冻结轨道并需验证平根数长期稳定性时
- 分析地球静止轨道在日月引力摄动下的倾角增长率及田谐共振特性时

## 工作原理

- 输入开普勒六根数（半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω、平近点角 M）及传播时长
- 配置摄动模型深度，包括 EGM96 带谐阶数（J2、J2–J4 或 J2–J6）以及日月四极双重平均选项
- 程序通过 RK4 积分轨道平均 Gauss 变分率演化平均根数，叠加一阶短周期扰动重构密切根数，并输出包含图表与数值交叉验证的 HTML 分析报告

## 使用案例

- 太阳同步遥感卫星轨道寿命与节点漂移率匹配设计
- 地球同步通信卫星寿命期内南北位置保持与倾角漂移量估算
- 大椭圆轨道与导航星座长期 Kozai–Lidov 偏心率与倾角振荡筛查

## 常见问题

### 平均根数传播与完全数值积分传播有何区别？

平均根数传播通过滤除轨道周期内的短周期振荡，直接对轨道平均变分率进行大步长积分，计算速度极快且便于观察长期世俗演化趋势。

### 工具支持哪些摄动力模型？

支持 EGM96 模型的 J2 至 J6 阶带谐引力摄动，以及基于四极阶双重平均模型的太阳和月球第三体引力摄动。

### 平均→密切重构面板的作用是什么？

该面板将一阶 J2 短周期几何偏差叠加回平均开普勒根数，恢复出瞬时密切轨道根数，以便与雷达观测或星历数据直接对比。

### 冻结轨道（Frozen Orbit）是如何计算的？

工具利用 J2 与 J3 摄动平衡关系，计算当近地点幅角固定在 90° 或 270° 时的平衡偏心率 e = −J3·sin(i) / (2·J2·p)。

### 最大支持传播多长时间？

工具支持 1 天至 1095 天（约 3 年）的轨道长期演化推演。

## 相关工具

- [相关性分析器](https://elysiatools.com/zh/tools/correlation-analyzer): 高级相关性分析工具，计算变量之间的相关系数以衡量它们线性关系的强度和方向。完美适用于统计分析、金融建模、科学研究和数据探索。 功能特点： - 多种相关性方法（皮尔逊、斯皮尔曼、肯德尔） - 相关性矩阵生成 - 统计显著性检验（p值） - 置信区间计算 - 热图可视化 - 散点图矩阵生成 - 缺失值处理策略 - 异常值检测和处理 - 分组分析能力 - 详细统计报告 常见用途： - 金融市场分析和风险评估 - 科学研究和假设检验 - 客户行为和营销分析 - 医疗和医疗数据分析 - 质量控制和流程优化 - 教育绩效评估
- [图片调色板转设计令牌](https://elysiatools.com/zh/tools/image-to-design-tokens): 从图片中用 k-means 聚类提取主色，导出为 CSS 变量 / SCSS 变量 / Tailwind 配置 / JSON 设计令牌，并自动为每个颜色生成命名色阶
- [租 vs 买总拥有成本导师（按揭/税/保险/维护 + NPV 盈亏平衡 + IRR 敏感性龙卷风）](https://elysiatools.com/zh/tools/rent-vs-buy-total-cost-of-ownership-mortgage-tax-insurance-maintenance-irr-tutor): 并排月度模拟：按揭摊销 + 房产税/保险/HOA/维护 vs 租金+租客保险；租客把首付与每月差额按机会成本利率投资；输出逐年 NPV、盈亏平衡年、差额现金流的 IRR、SALT 上限下的逐项抵扣税盾与 ±20% 参数龙卷风敏感性图。
- [局部（管件）水头损失计算器（h = K·v²/2g）](https://elysiatools.com/zh/tools/minor-loss-calculator): 计算管道管件/阀门/弯头的局部水头损失：h_m=K·v²/(2g)，ΔP_m=K·ρ·v²/2。内置 14 种常见管件的典型 K 值（闸阀、截止阀、止回阀、90°/45° 弯头、三通、突然扩大/收缩、出入口等），可直接选择或自定义 K。支持 n 个相同管件串联，输出水头（m/ft）、压降（Pa/kPa/bar）。密度支持 kg/m³、g/cm³、lb/ft³。
- [管道达西摩擦压降计算器（ΔP = f·(L/D)·ρv²/2）](https://elysiatools.com/zh/tools/pipe-pressure-drop-darcy): 计算直管段的达西-韦斯巴赫主要（摩擦）压降：ΔP=f·(L/D)·(ρ·v²/2) Pa；水头损失 h_f=ΔP/(ρ·g)=f·(L/D)·v²/(2g) m。用户直接输入达西摩擦系数 f（非范宁系数）。长度支持 m/km/ft，管径支持 m/cm/mm/inch，密度支持 kg/m³/g/cm³/lb/ft³，内部统一换算为 SI。同时输出 kPa、bar 以及 m/ft 水头。
- [边界值处理器](https://elysiatools.com/zh/tools/data-boundary-processor): 高级边界值处理工具，用于识别和处理数值数据中的最小值和最大值。完美适用于数据验证、范围检查、统计分析和数据预处理。 功能特点： - 多种边界检测方法（绝对值、百分位数、标准差） - 灵活处理策略（裁剪、删除、替换、变换） - 自定义范围验证 - 非对称边界处理 - 批量处理能力 - 综合边界统计 - 数据质量评估 - 可视化边界报告 常见用途： - 数据验证和质量控制 - 传感器数据范围检查 - 金融数据限制执行 - 统计数据预处理 - 机器学习特征工程 - 数据库约束验证
- [水电气账单行项解码器与隐藏费用探测器（CSR/CFL/MHC/HFA rider · 峰谷需量 · kWh/Btu/therm 换算）](https://elysiatools.com/zh/tools/csr-cfl-cfl-mhc-hfa-btu-ton-hours-utility-bill-line-item-validator-and-rgp-peak-explainer): 粘贴账单行项即可解码：客户费/需量费/峰谷电价/输配 rider/可再生能源基金/解耦 rider/容量标签等归类解释，kWh↔Btu↔therm↔Ccf↔m³↔冷吨时 EIA 标准换算，隐藏固定费与年度化测算，负荷率与削峰省钱建议。
- [密度计算器 (ρ = m/V)](https://elysiatools.com/zh/tools/density-calculator): 由三个量中的两个求第三个，附比重与浮沉判断

## 示例

- [Web Python 图像处理示例](https://elysiatools.com/zh/samples/web-image-processing-python): Web Python 图像处理示例，使用 PIL/Pillow 包括读取、保存、缩放和格式转换
- [Android Java 图像处理示例](https://elysiatools.com/zh/samples/android-image-processing-java): Android Java 图像处理示例，包括图像读取保存、缩放和格式转换
- [Android Kotlin 图像处理示例](https://elysiatools.com/zh/samples/android-image-processing-kotlin): Android Kotlin 图像处理示例，包括图像读取保存、缩放和格式转换
- [macOS Objective-C 图像处理示例](https://elysiatools.com/zh/samples/macos-image-processing-objectivec): macOS Objective-C 图像处理示例，包括图像读取保存、图像缩放和格式转换