第 61 章 逆特征值问题 (IEP)¶
前置：特征值 (Ch06) · 矩阵分析 (Ch14) · 正定矩阵 (Ch16)
本章脉络：从正向问题到逆向问题 \(\to\) 逆特征值问题 (Inverse Eigenvalue Problems) 的一般定义 \(\to\) 结构化约束（对称、Jacobi、Toeplitz）的重要性 \(\to\) 雅可比矩阵 (Jacobi Matrices) 的 IEP \(\to\) 基于谱数据的矩阵构造（Lanczos 方法） \(\to\) 存在性与唯一性判据 \(\to\) 数值算法：交替投影、Newton 法 \(\to\) 应用：结构设计（设计具有特定振动频率的桥梁）、地质勘探（根据震波反推地层硬度）、控制系统的极点配置
延伸：逆特征值问题是线性代数的“设计模型”；它将我们从“观察算子性质”提升到“构造满足特定物理特性的算子”，是连接数学分析与工程综合（Synthesis）的核心纽带
在线性代数的标准课本中，我们学习如何计算给定矩阵的特征值。但在工程实践中，问题往往相反：我们预先设定了一组理想的特征值（如桥梁的避震频率），需要构造出一个具有特定结构（如刚度矩阵）且恰好拥有这些特征值的矩阵。这就是逆特征值问题（Inverse Eigenvalue Problems, IEP）。本章将讨论这种“按需定制”矩阵的代数挑战与构造算法。
61.1 问题的定义与分类¶
定义 61.1 (逆特征值问题)
给定一组标量 \(\Lambda = \{\lambda_1, \lambda_2, \ldots, \lambda_n\}\)，寻找一个属于特定矩阵类 \(\mathcal{S}\) 的矩阵 \(A\)，使得 \(A\) 的谱恰好为 \(\Lambda\)。 - 结构约束 \(\mathcal{S}\)：通常要求矩阵是对称的、三对角的（Jacobi 阵）或具有特定的零元素模式。
61.2 雅可比矩阵的 IEP¶
定理 61.1 (Jacobi IEP 唯一性)
给定 \(n\) 个互异实数 \(\{\lambda_i\}\) 作为整个矩阵的特征值，以及 \(n-1\) 个互异实数 \(\{\mu_j\}\) 作为其 \((n-1)\) 阶主子阵的特征值。若满足交错性质： $\(\lambda_1 < \mu_1 < \lambda_2 < \mu_2 < \cdots < \mu_{n-1} < \lambda_n\)$ 则存在唯一的雅可比矩阵（三对角且次对角元为正）满足这些条件。
61.3 数值求解方法¶
技术：交替投影法
寻找满足两个集合交集的点： 1. 集合 1：具有目标特征值的矩阵（等谱轨道）。 2. 集合 2：具有特定结构的矩阵（线性子空间）。 通过在两个集合间反复投影，可以收敛到一个（局部）最优的构造结果。
练习题¶
1. [基础] 若要求构造一个 \(2 \times 2\) 实对称阵，使其特征值为 \(\{1, 3\}\)。给出一个解。
参考答案
构造方法： 最简单的解是直接取对角阵： \(A = \begin{pmatrix} 1 & 0 \\ 0 & 3 \end{pmatrix}\)。 补充：事实上，任何具有形式 \(Q \operatorname{diag}(1, 3) Q^T\) 的矩阵都是解。
2. [结构约束] 若在第 1 题的基础上增加约束：非对角元必须为 1。求矩阵。
参考答案
计算步骤： 1. 设 \(A = \begin{pmatrix} a & 1 \\ 1 & c \end{pmatrix}\)。 2. 迹的性质：\(a + c = \sum \lambda = 1 + 3 = 4\)。 3. 行列式性质：\(ac - 1 = \prod \lambda = 1 \cdot 3 = 3\)。 4. 解方程组：\(c = 4 - a \implies a(4-a) - 1 = 3 \implies a^2 - 4a + 4 = 0\)。 5. 解得 \(a = 2, c = 2\)。 结论：\(A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}\)。
3. [交错定理] 判定：是否可以构造一个 Jacobi 矩阵 \(A\)，使得 \(A\) 的特征值为 \(\{1, 10\}\)，而其 \(1 \times 1\) 子阵特征值为 \(\{12\}\)？
参考答案
结论：不可以。 理由：根据 Cauchy 交错定理，子阵特征值 \(\mu\) 必须被原阵特征值夹持，即 \(1 \le \mu \le 10\)。由于 \(12 > 10\)，该 IEP 问题在 Jacobi 结构约束下无解。
4. [唯一性] 为什么对于一般矩阵（无结构约束），IEP 的解不唯一？
参考答案
代数直观： 特征值只决定了矩阵的“尺度”，而没有决定矩阵的“方向”。 对于任何正交矩阵 \(Q\)，\(Q \Lambda Q^T\) 都具有相同的特征值。这种旋转自由度导致了解的无穷多性。只有施加了严格的结构约束（如三对角），自由度才会被耗尽从而产生唯一性。
5. [应用] 在桥梁避震设计中，逆特征值问题对应什么任务？
参考答案
解释： 1. 桥梁具有固有的共振频率（特征值）。 2. 如果环境噪声（如风、车辆）的频率接近特征值，桥梁会坍塌。 3. 工程师通过改变桥墩的质量和刚度（矩阵元素），使得最终的动力学矩阵特征值避开危险频率区间。这是一个典型的带约束的 IEP。
6. [计算] 给定特征值 \(\{2, -2\}\)，构造一个迹为 0 且行列式为 -4 的 \(2 \times 2\) 矩阵。
参考答案
构造： \(A = \begin{pmatrix} 0 & 2 \\ 2 & 0 \end{pmatrix}\)。 验证：\(\operatorname{tr} = 0, \det = -4\)。特征方程 \(\lambda^2 - 4 = 0 \implies \lambda = \pm 2\)。
7. [极点配置] 控制理论中的“极点配置”是 IEP 吗？
参考答案
是的。 通过状态反馈 \(u = -Kx\)，闭环系统的特征值变为 \(\sigma(A-BK)\)。 我们的目标是选择增益矩阵 \(K\)（结构约束），使得 \(\sigma(A-BK)\) 位于左半平面的指定位置。
8. [性质] 证明：若 \(A\) 是实对称 Jacobi 矩阵且次对角元不为 0，则其特征值必互异。
参考答案
证明思路： 1. 计算 \(A-\lambda I\) 的 \(n-1\) 阶主子阵。由于次对角元非零，该子阵的秩至少为 \(n-1\)。 2. 这意味着对于任何 \(\lambda\)，\(A-\lambda I\) 的几何重数最大为 1。 3. 由于 Jacobi 阵可对角化，代数重数等于几何重数。 结论：所有特征值只能是单根。
9. [数值算法] 简述 Newton 法求解 IEP 的基本思路。
参考答案
将特征值看作矩阵元素 \(a_{ij}\) 的非线性函数 \(f(a) = \Lambda\)。 利用矩阵微积分中的特征值导数公式（见 Ch42），通过计算 Jacobian 矩阵，对矩阵元素进行迭代修正，直至特征值偏差达到容差范围内。
10. [极限] 随着矩阵阶数增加，构造一个具有特定零模式的 IEP 难度如何？
参考答案
结论：难度呈指数级增加。 判定特定零模式下 IEP 是否有解通常是一个高度非线性的代数几何问题。对于大规模系统，通常不追求精确匹配，而是通过最小二乘 IEP 寻找最接近目标的结构化矩阵。
本章小结¶
逆特征值问题是线性代数的“逆向工程”：
- 从观测到创造：它打破了只能分析现有系统的被动局面，确立了按照目标性能设计物理结构的代数准则。
- 约束的约束：IEP 证明了结构的稀疏性与谱的自由度之间存在深刻的拮抗关系，揭示了物理可行性的数学边界。
- 计算的枢纽：作为连接系统辨识、极点配置与结构优化的纽带，IEP 算法是现代高精度工程设计的核心逻辑引擎。