首先介绍一个特殊的方阵J∈R2n×2nJ \in R^{2n \times 2n}JR2n×2n,其定义如下:
J=[OnIn−InOn]J=\begin{bmatrix} O_n & I_n\\ -I_n & O_n \end{bmatrix}J=[OnInInOn]
其中On∈Rn×nO_n\in R^{n\times n}OnRn×n是零矩阵;In∈Rn×nI_n \in R^{n \times n}InRn×n是单位矩阵。

那么这个矩阵具有下列的性质:
JT=−JJ^T = -JJT=J
J−1=JTJ^{-1} = J^TJ1=JT
JTJ=I2nJ^TJ = I_{2n}JTJ=I2n
JTJT=−I2nJ^TJ^T = -I_{2n}JTJT=I2n
J2=−I2nJ^2 = -I_{2n}J2=I2n
detJ=±1det J = \pm 1detJ=±1

哈密顿矩阵(Hamiltonian matrix)定义

一个矩阵A∈R2n×2nA\in R^{2n \times 2n}AR2n×2n叫哈密顿矩阵如果JAJAJA 是对称的,也就是说:
JA=(JA)T⇒ATJ+JA=0JA = (JA)^T \Rightarrow A^TJ+JA = 0JA=(JA)TATJ+JA=0,其中JJJ就是上面介绍的特殊矩阵。
H={A∈R2n×2n∣ATJ+JA=0}\mathcal{H} =\{A \in R^{2n\times 2n}|A^TJ+JA = 0\}H={AR2n×2nATJ+JA=0}2n×2n2n\times 2n2n×2n的哈密顿矩阵的集合。

以下三条性质是等价的:

  1. AAA是一个哈密顿矩阵
  2. A=JSA=JSA=JS,其中S=STS=S^TS=ST
  3. (JA)⊤=JA(JA)^\top = JA(JA)=JA

定理:
A,B∈HnA,B \in \mathcal{H}^nA,BHn,那么下面定理是对的:

  1. A+B∈HnA+B \in \mathcal{H} ^nA+BHn
  2. αA∈Hn\alpha A \in \mathcal{H}^nαAHn
  3. [A,B]∈Hn[A,B] \in \mathcal{H}^n[A,B]Hn,其中[A,B]=defAB−BA[A,B]\overset{\underset{\mathrm{def}}{}}{=}AB-BA[A,B]=defABBA

结果:(H,[⋅,⋅])(\mathcal{H}, [\cdot, \cdot])(H,[,])是个李代数(Lie algebra)。

定理

A∈HnA\in \mathcal{H}^nAHnpA(x)p_A(x)pA(x)为矩阵AAA的特征多项式,那么:

  • pA(x)=pA(−x)p_A(x) = p_A(-x)pA(x)=pA(x)
  • 如果pA(c)=0p_A(c) = 0pA(c)=0,那么pA(−c)=pA(cˉ)=pA(−cˉ)=0p_A(-c) = p_A(\bar c)=p_A(-\bar c)=0pA(c)=pA(cˉ)=pA(cˉ)=0,其中c∈Rc\in RcR

代数Riccati Equation

首先定义不变子空间(invariant subspace)
由向量v1,,...vkv_1,,...v_kv1,,...vk张成的线性空间ν\nuν叫做矩阵AAA的不变子空间如果对于任意v∈ν,Av∈νv \in \nu, Av \in \nuvν,Avν
考虑下面的代数Riccati Equation (ARE):
0=R(x)=F+ATX+XA+XA−XGX(1)0=R(x)=F+A^TX+XA+XA-XGX (1)0=R(x)=F+ATX+XA+XAXGX1
其中A,F,G,X∈Rn×nA,F,G,X\in R^{n \times n}A,F,G,XRn×nF,GF,GF,G是对称矩阵,XXX是未知的对称矩阵。
定义H=[AGF−AT]H =\begin{bmatrix}A & G\\ F& -A^T \end{bmatrix}H=[AFGAT],令[UV]\begin{bmatrix} U\\ V \end{bmatrix}[UV]
HHH的不变子空间的一个向量,
也就是说
[AGF−AT][UV]=[UV]Z,      (2)\begin{bmatrix} A & G\\ F & -A^T \end{bmatrix}\begin{bmatrix} U\\ V \end{bmatrix}=\begin{bmatrix} U\\ V \end{bmatrix}Z , \; \; \; (2)[AFGAT][UV]=[UV]Z,(2)
其中Z∈Rn×n,λ(Z)⊂λ(H)Z\in R^{n\times n}, \lambda(Z) \subset \lambda(H)ZRn×n,λ(Z)λ(H)
假定UUU是非奇异的,那么从 (2)(2)(2)的第一行得到:
AU+GV=UZ↔U−1AU+U−1GV=ZAU+GV = UZ \leftrightarrow U^{-1}AU+U^{-1}GV=ZAU+GV=UZU1AU+U1GV=Z
将其插入第二行得到的方程:
FU−ATV=VZ=VU−1AU+VU−1GVFU-A^TV=VZ=VU^{-1}AU+VU^{-1}GVFUATV=VZ=VU1AU+VU1GV
以上方程等于
0=F−ATVU−1−VU−1A−VU−1GVU−10=F-A^TVU^{-1}-VU^{-1}A-VU^{-1}GVU^{-1}0=FATVU1VU1AVU1GVU1
X:=−VU−1X:= -VU^{-1}X:=VU1,我们发现XXX正好是(1)(1)(1)的解。

定理

考虑哈密顿矩阵HHH,假设该矩阵没有特征值在虚轴上,其不变子空间是χ−=Im[X1X2]⊂R2n×n\chi_{-} = Im \begin{bmatrix} X_1\\ X_2 \end{bmatrix} \subset R^{2n \times n}χ=Im[X1X2]R2n×n
这是一个由稳定特征值对应的特征向量张成的空间。
如果X1X_1X1是可逆的,那么X=−X2X1−1X=-X_2X_1^{-1}X=X2X11是代数Ricatti方程的解,并且A−GXA-GXAGX是稳定的。

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐