跳到正文
格致开物MATHWIKI

代数黎卡提方程

代数黎卡提方程(algebraic Riccati equation,ARE)是一类含未知矩阵二次项的矩阵方程,在最优控制和状态估计中出现。控制中的解不只需要满足代数等式,还通常要选择使闭环稳定的那一支;“求到一个实根”不足以完成控制设计。

连续时间控制方程

给定实矩阵 A∈ℝn×n、B∈ℝn×m,以及 Q=Q𝖳⪰0、R=R𝖳≻0,连续时间控制ARE为 A𝖳P+PA−PBR−1B𝖳P+Q=0,P=P𝖳. 二次项 PBR−1B𝖳P 使它不同于线性的李雅普诺夫方程。若 A−BR−1B𝖳P 全部特征值实部为负,称 P 为稳定化解。[1]

把 V=x𝖳Px 代入 x˙=Ax+Bu,展开并配平方: x𝖳Qx+u𝖳Ru+V˙=(u+R−1B𝖳Px)𝖳R(u+R−1B𝖳Px). 等式成立的条件正是ARE。右边非负,且在 u=−R−1B𝖳Px 时为零,由此连接到线性二次调节器的最优性证明。方程中的每个项都由同一代价和动力学产生,而非凭经验拟合的参数。

一个方程的两个根

取 A=B=Q=R=1,方程化为 2p−p2+1=0,两根为 p+=1+2、p−=1−2。闭环系数为 1−p:正根给出 −2,负根给出 +2。负根满足代数方程,却让状态增长,也不是非负代价的值函数。

二次残差2p-p平方+1在1正负sqrt2处为零,正根稳定而负根不稳定
横轴是未知数p,曲线是方程残差。只有右侧根产生稳定闭环;残差为零本身不能选定正确分支。

标准存在唯一性定理需要 (A,B) 可稳定、(Q1/2,A) 可检测。在这些条件下,存在唯一半正定稳定化解。这里“唯一”限定在该类解中,并非说所有对称代数根只有一个。数值求解后应同时检查对称性、方程残差和闭环特征值。

离散时间形式

对 xk+1=Axk+Buk,代价为 ∑k≥0(xk𝖳Qxk+uk𝖳Ruk),对应方程是 P=Q+A𝖳PA−A𝖳PB(R+B𝖳PB)−1B𝖳PA. 它由把下一步代价 (Ax+Bu)𝖳P(Ax+Bu) 加到当前代价,再对 u 配平方得到。控制增益为 K=(R+B𝖳PB)−1B𝖳PA,稳定化条件改为 A−BK 的特征值模小于一。[2]

变式。 离散标量 A=B=Q=R=1 给出 p=1+p−p2/(1+p),即 p2−p−1=0。非负根为 φ=(1+5)/2,反馈 k=φ/(1+φ)=1/φ,闭环系数 1−k=1/φ2 在单位圆内。它与连续时间的根 1+2 不同,不能只换导数符号而沿用同一方程。

参考资料

  1. ↑ Russ Tedrake,Underactuated Robotics, Linear Quadratic Regulators,Basic Derivation,代数Riccati方程及最优反馈。
  2. ↑ Russ Tedrake,Underactuated Robotics, Linear Quadratic Regulators,Discrete-time LQR。