庞特里亚金极值原理（Pontryagin's Maximum Principle）

本文作者：天疆说

本站地址：https://cislunarspace.cn

定义

庞特里亚金极值原理（PMP）由 Pontryagin 等 1962 提出，给出最优控制必须满足的一组必要条件。它取代经典变分法，能够处理带不等式约束（如有界推力）的控制集。给定系统 x˙=f(x,u,t)\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x},\mathbf{u},t) 和增广性能指标，PMP 引入协态 λ\boldsymbol{\lambda} 与 Hamilton 函数 H=L+λ ⁣⊤fH=L+\boldsymbol{\lambda}^{\!\top}\mathbf{f}，断言最优控制在每一时刻使 HH 在容许控制集上取极大（俄语学派原始约定）或极小（工程界，常用于燃料/时间问题），两者仅差 HH 与 λ0\lambda_0 的整体符号，最优轨迹相同，见 Pontryagin 等 1962、Bryson 与 Ho 1975 和 Betts 2010[548][217][198]。

一阶必要条件全集

对 Bolza 问题 J=Φ(x(tf),tf)+∫t0tfL dtJ=\Phi(\mathbf{x}(t_f),t_f)+\int_{t_0}^{t_f}L\,dt，终端约束 Ψ=0\boldsymbol{\Psi}=0，PMP 给出五条条件，见 Ch. 1[260][678]：

  1. 状态方程 x˙=∂H/∂λ=f(x,u∗,t)\dot{\mathbf{x}}=\partial H/\partial\boldsymbol{\lambda}=\mathbf{f}(\mathbf{x},\mathbf{u}^*,t)。
  2. 协态方程 λ˙=−∂H/∂x\dot{\boldsymbol{\lambda}}=-\partial H/\partial\mathbf{x}（详见 协态变量与共轭方程）。
  3. 极值条件（最优性）：u∗\mathbf{u}^* 在容许集 UU 上逐点极小化 HH。无界控制对应 ∂H/∂u=0\partial H/\partial\mathbf{u}=0。有界控制极值在 UU 的边界取得。
  4. 横截条件 λ(tf)=∂Φ/∂x∣tf+(∂Ψ/∂x) ⁣⊤ν\boldsymbol{\lambda}(t_f)=\partial \Phi/\partial\mathbf{x}|_{t_f}+(\partial\boldsymbol{\Psi}/\partial\mathbf{x})^{\!\top}\boldsymbol{\nu}，其中 ν\boldsymbol{\nu} 为终端约束乘子。
  5. 参数条件：终端时刻 tft_f 自由时 H(tf)+∂Φ/∂t+ν ⁣⊤∂Ψ/∂t=0H(t_f)+\partial\Phi/\partial t+\boldsymbol{\nu}^{\!\top}\partial\boldsymbol{\Psi}/\partial t=0。对时间无关问题退化为 H(tf)=0H(t_f)=0，若 HH 自治则 H≡0H\equiv 0。

这五条共同构成一阶必要条件。联立求解即间接法，得到以 (x,λ)(\mathbf{x},\boldsymbol{\lambda}) 为变量的 TPBVP。

切换函数与 bang-bang 结构

当控制以线性形式进入 HH（例如推力幅值 T∈[0,Tmax]T\in[0,T_{max}]），HH 的极小化要求 TT 取边界值。定义标量切换函数 ρ(t)\rho(t) 作为 HH 中 TT 的系数：ρ<0\rho<0 时 T=TmaxT=T_{max}，ρ>0\rho>0 时 T=0T=0，由此导出 bang-bang / bang-off-bang 结构，见 Lawden 1963 与 Bang-bang 控制[452]。若 ρ≡0\rho\equiv 0 持续一段有限时间则为奇异弧，TT 由高阶条件决定。

推力方向的最优值为 −λv/∥λv∥-\boldsymbol{\lambda}_v/\|\boldsymbol{\lambda}_v\|，引出引燃矢量 p=−λv\mathbf{p}=-\boldsymbol{\lambda}_v，用于评估与改进脉冲转移，见 Lawden 1963 与 引燃矢量。

与变分法的衔接

PMP 是变分法（calculus of variations, CoV）的现代形式。对光滑无约束问题，经典 Euler-Lagrange 方程

ddt∂L∂q˙−∂L∂q=0\frac{d}{dt}\frac{\partial L}{\partial\dot{\mathbf{q}}}-\frac{\partial L}{\partial\mathbf{q}}=0

是作用量 ∫L dt\int L\,dt 取极值的必要条件。Legendre 变换给出 Hamilton 形式 q˙=∂H/∂p,  p˙=−∂H/∂q\dot{\mathbf{q}}=\partial H/\partial\mathbf{p},\;\dot{\mathbf{p}}=-\partial H/\partial\mathbf{q}。PMP 把变分法推广到带界控制、不等式约束与非光滑动力学的情形，Euler-Lagrange 方程是 u=q˙\mathbf{u}=\dot{\mathbf{q}} 不受约束时的特例。

离散 Euler-Lagrange 方程是 DMOC（离散力学与最优控制）的离散对应物：把作用量替换为离散和、施加离散 Lagrange-d'Alembert 原理导出离散必要条件，保持连续问题的辛结构，见 Marsden 与 West 2001[481]。

应用要点

  • PMP 给出的是必要而非充分条件。充分性需要凸性或附加二阶判据（Legendre-Clebsch、共轭点，见 共轭点与二阶最优性）。
  • 对 LL 关于 u\mathbf{u} 线性的最小燃料问题，光滑极值条件退化，PMP 的边界规则才是决定控制的关键。
  • 自治 Hamilton 函数沿最优轨迹为常数，这是检验 TPBVP 数值精度的常用判据。

相关概念

完善页面
最近更新:
Contributors: cislunarspace, ouyangjiahong