哈密顿矩阵和代数Riccati方程
首先介绍一个特殊的方阵J∈R2n×2nJ \in R^{2n \times 2n}J∈R2n×2n,其定义如下:J=[OnIn−InOn]J=\begin{bmatrix}O_n & I_n\\-I_n & O_n\end{bmatrix}J=[On−InInOn]其中On∈Rn×nO_n\in R^{n\times n}On∈Rn×n是零矩阵;In∈Rn×nI_n
首先介绍一个特殊的方阵J∈R2n×2nJ \in R^{2n \times 2n}J∈R2n×2n,其定义如下:
J=[OnIn−InOn]J=\begin{bmatrix} O_n & I_n\\ -I_n & O_n \end{bmatrix}J=[On−InInOn]
其中On∈Rn×nO_n\in R^{n\times n}On∈Rn×n是零矩阵;In∈Rn×nI_n \in R^{n \times n}In∈Rn×n是单位矩阵。
那么这个矩阵具有下列的性质:
JT=−JJ^T = -JJT=−J
J−1=JTJ^{-1} = J^TJ−1=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}A∈R2n×2n叫哈密顿矩阵如果JAJAJA 是对称的,也就是说:
JA=(JA)T⇒ATJ+JA=0JA = (JA)^T \Rightarrow A^TJ+JA = 0JA=(JA)T⇒ATJ+JA=0,其中JJJ就是上面介绍的特殊矩阵。
令H={A∈R2n×2n∣ATJ+JA=0}\mathcal{H} =\{A \in R^{2n\times 2n}|A^TJ+JA = 0\}H={A∈R2n×2n∣ATJ+JA=0}为2n×2n2n\times 2n2n×2n的哈密顿矩阵的集合。
以下三条性质是等价的:
- AAA是一个哈密顿矩阵
- A=JSA=JSA=JS,其中S=STS=S^TS=ST
- (JA)⊤=JA(JA)^\top = JA(JA)⊤=JA
定理:
令A,B∈HnA,B \in \mathcal{H}^nA,B∈Hn,那么下面定理是对的:
- A+B∈HnA+B \in \mathcal{H} ^nA+B∈Hn
- αA∈Hn\alpha A \in \mathcal{H}^nαA∈Hn
- [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]=defAB−BA
结果:(H,[⋅,⋅])(\mathcal{H}, [\cdot, \cdot])(H,[⋅,⋅])是个李代数(Lie algebra)。
定理
令A∈HnA\in \mathcal{H}^nA∈Hn,pA(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 Rc∈R。
代数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+XA−XGX(1),
其中A,F,G,X∈Rn×nA,F,G,X\in R^{n \times n}A,F,G,X∈Rn×n,F,GF,GF,G是对称矩阵,XXX是未知的对称矩阵。
定义H=[AGF−AT]H =\begin{bmatrix}A & G\\ F& -A^T \end{bmatrix}H=[AFG−AT],令[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)[AFG−AT][UV]=[UV]Z,(2)
其中Z∈Rn×n,λ(Z)⊂λ(H)Z\in R^{n\times n}, \lambda(Z) \subset \lambda(H)Z∈Rn×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=UZ↔U−1AU+U−1GV=Z。
将其插入第二行得到的方程:
FU−ATV=VZ=VU−1AU+VU−1GVFU-A^TV=VZ=VU^{-1}AU+VU^{-1}GVFU−ATV=VZ=VU−1AU+VU−1GV。
以上方程等于
0=F−ATVU−1−VU−1A−VU−1GVU−10=F-A^TVU^{-1}-VU^{-1}A-VU^{-1}GVU^{-1}0=F−ATVU−1−VU−1A−VU−1GVU−1。
令X:=−VU−1X:= -VU^{-1}X:=−VU−1,我们发现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=−X2X1−1是代数Ricatti方程的解,并且A−GXA-GXA−GX是稳定的。
更多推荐

所有评论(0)