精读笔记
Problem Setting
《Optimal Powered Descent Guidance with Pyramid-Shaped Approach-Angle Constraints》（arXiv preprint / 2026）。这篇论文实际处理的是一个很窄但重要的问题：在经典 powered descent LQ guidance 框架里，显式加入 approach-angle / ground-collision avoidance 的状态路径约束，同时尽量保留解析、连续、闭环可实现的 guidance law。
困难点在于路径约束不是终端约束。无约束 LQ powered descent 的漂亮结构来自双积分动力学、常重力和二次控制代价；一旦要求轨迹始终留在某个 glide-slope cone / pyramid 内，最优轨迹会出现约束激活、切换、降维运动和边界弧，原来的全程多项式/线性控制结构不再直接成立。
以前方法大致卡在两端：Apollo/E-guidance/AAPDG 和 LQ guidance 简洁可实现，但不真正 enforce path constraint；fuel-optimal/convex/direct/indirect path-constrained 方法能处理约束，但通常数值化、开环化、bang-bang 化，工程上不如 LQ guidance 干净。关键矛盾是：几何安全约束要求强约束和切换逻辑，而 onboard guidance 又希望闭式、连续、低计算量。
Motivation
作者的动机不是再做一个更精确的 landing optimizer，而是补上 LQ powered-descent guidance 的一个结构性短板：它能把终端位置/速度做得很漂亮，但对中途撞地、浅角进近没有硬保证。
已有 waypoint/intermediate-point 方法本质上是用人为几何点绕开 path constraint，约束没有进入最优控制问题本身，因此只能部分控制接近方向，也不保证整段轨迹安全。已有 fuel-optimal path-constrained 工作虽然更严格，但控制结构复杂，常出现 bang-bang、多个 subarc、数值 shooting/convex solve，和“连续、解析、可闭环”的工程诉求相冲突。
论文的核心观察是：如果接受 quadratic effort 而非 propellant-optimal，把可行域设计成凸倒金字塔，并在每个面上使用 constraint-aligned coordinates，那么 active path constraint 会变成非常简单的二阶坐标约束。这样路径约束不一定会摧毁解析性，反而可以通过降维 LQ 子问题保留闭式结构。
Core Idea
核心思想是把 approach-angle constraint 从“复杂空间几何限制”重写成“沿约束法向坐标非负”的问题。对每个金字塔面建立局部坐标系 (n,t,r)：n 是面法向，t 是面内垂直边界线方向，r 是边方向。于是面约束是 r_n ≥ 0，边约束是 r_t ≥ 0。轨迹一旦触碰面，就要求 r_n=0、V_n=0、u_n+g_n=0；一旦触碰边，再要求 r_t=0、V_t=0、u_t+g_t=0。每个 active constraint 都固定一个法向控制分量，并把最优控制问题降一维。
这和 prior 的本质区别在于：它不是通过 waypoint 改变参考几何，也不是用通用数值优化器求 path-constrained OCP，而是利用二阶积分器和二次代价的结构，把 constrained arc 变成低维 unconstrained LQ guidance。新的 inductive bias 是“安全边界上的最优运动应沿约束流形滑行，法向只做重力抵消”。这个 bias 很强，也正是解析解能成立的原因。
从 scalability 看，它对任意凸多面体不是完全自动 generalize，但对金字塔/多面锥这类由平面和边组成、且约束激活顺序可枚举的几何约束很自然。它的可扩展性来自几何分解和坐标旋转，而不是来自更强优化算法。
Method
1. 金字塔路径约束：解决的是进近角和避碰的硬约束表达问题。把着陆点作为顶点的倒金字塔定义为可行域，使轨迹必须从内部接近目标。这样 approach-angle 不再是终端方向偏好，而是全程几何管道。
2. Constraint-aligned frame：解决的是约束表达复杂的问题。每个面旋转到 (n,t,r) 后，面/边约束分别变成 r_n ≥ 0、r_t ≥ 0。这个步骤的核心变化是让路径约束的微分阶数和切触条件显式化，而不是依赖一般不等式处理。
3. Active constraint 上的降维动力学：面约束为二阶约束，因此 active 后必须满足 r_n=V_n=0，并令 u_n=-g_n；边约束同理令 u_t=-g_t。这样 3D 问题依次降为 2D 面内问题和 1D 边上问题。这里真正关键的是“约束不是惩罚项，而是直接消去自由度”。
4. PMP + switching conditions：解决的是不同弧之间如何拼接仍保持最优的问题。作者用 Hamiltonian/协态跳跃条件推导切换时刻和协态常数，得到连续控制。这个部分的价值在于说明分段结构不是 heuristic，而是满足必要最优性条件。
5. 闭环形式与自由终端时间：最终控制律可写成依赖当前状态的形式。法向未激活时出现 -2V_n^2/(3r_n) 这类非线性项，激活后退化为重力抵消；沿最终未被约束的 r 方向保持经典 LQ landing law。自由终端时间由四次多项式给出，并证明 active constraint 越多，最优时间越短。
Key Insight / Why It Works
最重要的 insight 是：对双积分器 + 常重力 + 二次控制代价，二阶路径约束一旦 active，最优法向控制不是复杂反馈，而是简单的重力抵消。因为在约束流形上 r_normal=0、V_normal=0，为保持不离开流形必须 u_normal+g_normal=0；而无约束 LQ law 在这些状态条件下也恰好给出同样结果。这使得“到达约束后沿约束走”不是额外假设，而是和无约束最优律一致的结构性质。
真正有效的是几何约束和系统结构的匹配：平面/直线约束是二阶状态约束，动力学是二阶积分器，代价是控制平方。三者对齐后，active constraint 可以被解析消元，剩余切向方向继续套用 LQ landing。换句话说，这篇的核心贡献不是更复杂的 PMP 推导，而是找到一种建模方式，让 path-constrained OCP 退化成一串低维经典 LQ 问题。
最可能的核心贡献：金字塔约束下的分段解析结构、约束激活时的法向重力抵消解释、以及 active constraint 使 optimal final time 变短的解析关系。最可能只是辅助的是完整四面金字塔的仿真实现、Euler angle 参数化、以及切换前短暂 open-loop 的平滑工程处理。
这不是 scaling、retrieval、data coverage 一类工作；本质是 better inductive bias / latent structure exploitation。它强依赖问题的可积结构和约束几何，利用的是解析结构而不是数值算力。若换成非线性重力、推力饱和、质量变化或姿态约束，这个结构会迅速变脆。
Relation To Prior Work
最接近的谱系是 D’Souza 的二次控制能量最优 powered-descent guidance、E-guidance 的 LQ 解释、以及 Guo/Wie 一类 waypoint-ZEM/ZEV guidance。本文可以看作把这条解析 LQ guidance 线路推进到显式路径约束，而不是 fuel-optimal powered descent 那条 bang-bang/convex optimization 线路。
相对 APDG/E-guidance/AAPDG，真正新增的是路径约束进入最优控制问题，而不是只满足终端状态或预设 thrust profile。相对 waypoint/intermediate-point 方法，区别是安全几何不是通过中间点间接诱导，而是通过 state inequality 和 constrained arc 明确 enforce。
相对 fuel-optimal path-constrained work，本文牺牲了 propellant optimality，换取连续控制和闭式结构。这里的本质取舍很清楚：它不是要解决最真实的燃料最优 landing，而是提供一个工程上更平滑、更可实时实现的 quadratic-effort guidance law。
看似新的地方中，PMP、Hamiltonian continuity、二阶状态约束切触条件都不是新理论；实质创新在于把这些标准工具和金字塔进近角几何结合，得到一个完整的解析制导律和切换时间表达。
Dataset / Evaluation
评估是纯仿真，没有真实飞行、硬件在环或高保真 6-DoF 验证。样例基于月球着陆参数，展示无约束解穿越地面，而约束解保持在金字塔内；完整 square pyramid 通过四个平面组合，并对多组初始位置和速度扰动测试。
这些实验基本支持两个有限 claim：在理想点质量常重力模型下，推导的控制律能满足金字塔路径约束并保持很小终端误差；完整 pyramid 约束可以通过面/边选择逻辑组合起来工作。
但 evaluation 没有真正验证鲁棒部署能力。控制输入假设 unbounded，质量变化没有进入优化，姿态/推力方向/最小最大推力都没有处理；扰动也只是初始条件变化，不是导航误差、建模误差、执行误差或地形不确定性。文中提到 constrained solution 代价更高，这是预期现象，不构成强验证。
另外，complete pyramid simulation 中 guidance law 是基于 anticipated constraint activation 选择的，文中未充分说明该选择在边界情形、多个约束近同时触发、噪声导致错误预测时是否稳定。
Limitation
最核心限制是问题设定过于理想。解析性来自点质量双积分器、常重力、无推力约束、二次加速度代价和固定平面约束的共同作用。任何一个条件被打破，控制律的简单结构都可能不成立。
第二，最优性范围很窄。它是 minimum control effort，不是 minimum fuel；u=T/m 但质量 m(t) 没有作为状态进入动力学和代价，实际推进剂消耗与这里的 J 只有间接关系。对真实 powered descent，推力饱和、节流限制、发动机最小推力、姿态约束通常是主导约束，本文没有处理。
第三，约束激活结构被强假设为 unconstrained → plane → line，且通过 t1=-3r_n/V_n、t2=-3r_t/V_t 这类表达判断。若出现多面同时激活、初始就在约束附近、速度符号异常、噪声下反复穿越边界，切换逻辑可能不稳。文中未充分说明这种混合系统层面的鲁棒性。
第四，闭环法向控制在 active arc 上理论上不需要，因为切触条件应精确满足；但实际仿真又承认数值误差可能破坏切触，并建议用无约束法向闭环修正。这说明解析 constrained law 和实际鲁棒实现之间仍有 gap。切换前一秒改 open-loop 也是工程补丁，最优性和稳定性影响没有严谨分析。
第五，optimal final time 超出 admissible range 时直接 clip 到最近可行值。这个操作保证避碰，但 clipped 后是否仍是原问题全局最优，文中没有充分论证。
Takeaway
- 1. 最值得记住的是：对二阶积分器，平面/线性路径约束 active 后可以看作动力学降维，而不是一般不等式优化；这会把复杂 path-constrained guidance 变回低维 LQ guidance。
- 2. 进近角约束用倒金字塔/多面锥表达，是一个强但实用的 inductive bias：它牺牲几何灵活性，换取可解析性和可验证的安全边界。
- 3. 这篇推动的是“解析 LQ guidance + 显式几何路径约束”这条路线，而不是燃料最优 powered descent。
- 未来真正有价值的方向是把这种结构和推力/姿态/质量约束结合，看看解析结构能保留多少。
一句话总结
这篇论文是解析 LQ powered-descent guidance 向显式几何路径约束的一次结构化扩展，真正贡献在于利用金字塔约束与双积分器动力学的匹配，把 approach-angle constrained landing 化为分段降维的闭式最优制导律。