第 53A 章 辛矩阵与 Hamilton 矩阵¶
前置：内积空间 (Ch08) · 矩阵方程 (Ch20) · 经典力学基础
本章脉络：从欧几里得几何到辛几何 \(\to\) 辛形式 (Symplectic Form) 与标准反对称阵 \(J\) \(\to\) 辛矩阵 (Symplectic Matrices) 的定义与群结构 \(\to\) 辛矩阵的谱性质（特征值的配对性） \(\to\) Hamilton 矩阵的定义与李代数结构 \(\to\) Hamilton 矩阵与 Riccati 方程的联系 \(\to\) 辛分解（辛 QR、辛 Schur 分解） \(\to\) 应用：哈密顿力学中的相空间保体积性、最优控制中的状态-协态方程、量子光学
延伸：辛矩阵是描述“面积保持”变换的算子；它揭示了物理系统在能量守恒下的内在代数对称性，是连接经典力学、偏微分方程和控制理论的高级代数语言
在线性代数中，我们熟悉保持内积不变的正交矩阵。但在物理学（特别是哈密顿力学）中，更重要的是保持一种“斜对称内积”不变的变换。辛矩阵（Symplectic Matrices）正是这类算子的代数刻画。它们描述了相空间中面积的守恒，并引出了一类具有完美对称特征值结构的哈密顿矩阵（Hamiltonian Matrices）。本章将探索这一深刻影响物理演化与最优控制的特殊矩阵类。
53A.1 辛形式与辛矩阵¶
定义 53A.1 (标准辛阵 \(J\))
定义 \(2n \times 2n\) 的分块矩阵 \(J\) 为： $\(J = \begin{pmatrix} 0 & I_n \\ -I_n & 0 \end{pmatrix}\)$ 满足 \(J^2 = -I\) 且 \(J^T = -J\)。它定义了辛内积 \([x, y] = x^T J y\)。
定义 53A.2 (辛矩阵)
方阵 \(M \in M_{2n}\) 称为辛矩阵，如果它保持辛形式不变： $\(M^T J M = J\)$ 所有 \(2n\) 阶辛矩阵构成的集合称为 辛群 \(Sp(2n)\)。
53A.2 Hamilton 矩阵¶
定义 53A.3 (Hamilton 矩阵)
方阵 \(H \in M_{2n}\) 称为 Hamilton 矩阵，如果 \(JH\) 是对称的，即： $\((JH)^T = JH \iff H^T J + JH = 0\)$ 性质：若 \(M(t)\) 是辛矩阵且 \(M(0)=I\)，则其导数 \(\dot{M}(0)\) 必为 Hamilton 矩阵。
53A.3 谱性质与配对性¶
定理 53A.1 (特征值配对)
- 辛矩阵：若 \(\lambda\) 是其特征值，则 \(1/\lambda\) 也必是特征值。
- Hamilton 矩阵：若 \(\lambda\) 是其特征值，则 \(-\lambda\) 也必是特征值。 物理意义：这种配对性反映了物理系统中的稳定与不稳定模态的某种对称平衡。
练习题¶
1. [基础] 计算 \(2 \times 2\) 辛矩阵的行列式值。
参考答案
结论：1。 证明：由 \(M^T J M = J\) 两边取行列式： \(\det(M^T) \det(J) \det(M) = \det(J)\)。 由于 \(\det(J) = 1\)，得 \((\det M)^2 = 1 \implies \det M = \pm 1\)。 更为深刻的理论证明（Pfaffian）表明，辛矩阵的行列式必然为 +1。这对应了相空间体积守恒。
2. [Hamilton] 判定 \(H = \begin{pmatrix} A & B \\ C & -A^T \end{pmatrix}\) 是否为 Hamilton 矩阵（其中 \(B, C\) 对称）。
参考答案
验证步骤： 1. \(JH = \begin{pmatrix} 0 & I \\ -I & 0 \end{pmatrix} \begin{pmatrix} A & B \\ C & -A^T \end{pmatrix} = \begin{pmatrix} C & -A^T \\ -A & -B \end{pmatrix}\)。 2. 取转置：\((JH)^T = \begin{pmatrix} C^T & -A^T \\ -A & -B^T \end{pmatrix}\)。 3. 由于 \(B, C\) 对称， \(C^T=C, B^T=B\)。 结论：\((JH)^T = JH\)，故该分块形式是 Hamilton 矩阵的标准构造。
3. [计算] 验证 \(J = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix}\) 本身是否是辛矩阵。
参考答案
计算： \(J^T J J = (-J) J J = (-J)(-I) = J\)。 结论：是的。\(J\) 既是辛矩阵，也是 Hamilton 矩阵（在 2 维情形下）。
4. [谱性质] 若 Hamilton 矩阵的一个特征值为 \(2+3i\)，求其余三个必然存在的特征值。
参考答案
利用配对性： 1. 既然是实矩阵，必有共轭根：\(2-3i\)。 2. 根据 Hamilton 谱性质，必有相反数根：\(-(2+3i) = -2-3i\)。 3. 同理，相反数的共轭：\(-2+3i\)。 结论：特征值以 \(\{ \pm \lambda, \pm \bar{\lambda} \}\) 的四重奏形式出现。
5. [保面积性] 证明辛算子在 \(\mathbb{R}^2\) 平面上保持平行四边形面积不变。
参考答案
证明： 1. 面积可以由辛形式 \(x^T J y\) 给出。 2. 变换后 \((Mx)^T J (My) = x^T (M^T J M) y\)。 3. 由辛矩阵定义，该式等于 \(x^T J y\)。 结论：变换前后的辛内积相等，故有向面积保持不变。
6. [Riccati] Hamilton 矩阵如何用于求解代数 Riccati 方程 (ARE)？
参考答案
代数纽带： ARE 的解 \(X\) 对应于其关联 Hamilton 矩阵 \(H\) 的稳定不变子空间。 通过求出 \(H\) 的特征向量并进行分块组合，可以直接构造出 ARE 的正定解。这避免了非线性迭代，是现代控制理论的标准求法。
7. [逆矩阵] 证明：若 \(M\) 是辛矩阵，则其逆矩阵 \(M^{-1}\) 也是辛矩阵。
参考答案
证明： 1. 已知 \(M^T J M = J\)。 2. 两边左乘 \((M^T)^{-1}\)，右乘 \(M^{-1}\)： 3. \(J = (M^T)^{-1} J M^{-1} = (M^{-1})^T J M^{-1}\)。 结论：满足辛矩阵定义。
8. [李代数] 证明：两个 Hamilton 矩阵的交换子 \([H_1, H_2]\) 仍是 Hamilton 矩阵。
参考答案
证明思路： 这是因为 Hamilton 矩阵构成了辛群 \(Sp(2n)\) 的李代数 \(\mathfrak{sp}(2n)\)。李代数对交换子运算封闭。
9. [计算] 写出 \(4 \times 4\) 标准辛阵 \(J\) 的逆矩阵。
参考答案
结论：\(-J\)。 理由：由于 \(J^2 = -I\)，故 \(J(-J) = I\)。在物理上，这对应于辛旋转的逆向操作。
10. [应用] 简述“辛积分器”（Symplectic Integrator）在轨道计算中的优势。
参考答案
理由： 1. 传统的数值积分（如龙格-库塔法）会随着时间积累能量误差，导致行星轨道漂移。 2. 辛积分器在每一步都强制满足 \(M^T J M = J\)。 3. 这保证了系统的哈密顿量（总能量）在长期演化中几乎守恒，能够长时间保持轨道的闭合性与稳定性。
本章小结¶
辛矩阵与 Hamilton 矩阵揭示了动力系统的深层几何秩序：
- 面积的守护者：辛矩阵通过保持斜对称形式，确立了相空间流动的不可压缩性，是统计力学刘维尔定理的代数根源。
- 谱的对称美：Hamilton 矩阵特征值的正负配对性，完美刻画了保守系统在平衡点附近的局部稳定性，为寻找稳定流形提供了天然基准。
- 控制的桥梁：通过与 Riccati 方程的深刻联系，辛代数实现了从抽象物理对称性到工程最优控制解的跨越，是现代系统科学最强有力的代数支柱。