[论文解读] A Method for Fast Diagonalization of a 2x2 or 3x3 Real Symmetric Matrix

M. J. Kronenburg|arXiv (Cornell University)|Jun 26, 2013

Advanced optical system design参考文献 4被引用 3

一句话总结

本文提出了一种快速、符号化的算法，用于通过使用旋转矩阵和欧拉角推导出的闭式公式，对 2×2 和 3×3 实对称矩阵进行对角化，计算其特征值和特征向量。主要贡献在于一种数值稳定且高效的算法，避免了迭代求解器，并通过计算机代数系统进行了符号证明，同时提供了经过验证的 C++ 代码实现。

ABSTRACT

A method is presented for fast diagonalization of a 2x2 or 3x3 real symmetric matrix, that is determination of its eigenvalues and eigenvectors. The Euler angles of the eigenvectors are computed. A small computer algebra program is used to compute some of the identities, and a C++ program for testing the formulas has been uploaded to arXiv.

研究动机与目标

• 开发一种快速、符号化的算法，用于在不使用迭代数值求解器的情况下，计算 2×2 和 3×3 实对称矩阵的特征值和特征向量。
• 利用旋转矩阵和欧拉角，推导出特征值和特征向量的闭式解析表达式。
• 为小型对称矩阵提供一种数值稳定且计算高效的对角化方法，作为标准数值对角化方法的替代方案。
• 通过符号计算和经过测试的 C++ 实现，验证所推导公式的正确性。

提出的方法

• 对于 2×2 矩阵，特征值通过带符号调整的平方根公式计算；特征向量通过计算非对角元素比值的反正切函数确定。
• 对于 3×3 矩阵，特征值通过三次特征方程的三角函数解求解，参数化为 p、q 和 Δ，使用反余弦函数表示。
• 特征向量通过从矩阵元素和特征值推导出的三个旋转角 φ₁、φ₂、φ₃ 重构，使用绕固定轴的旋转矩阵。
• 该方法利用计算机代数程序（Mathematica）验证的符号恒等式，确保所推导表达式的正确性。
• 提供了 C++ 实现以测试和验证解析公式，确保数值稳定性和性能。
• 从旋转角中计算欧拉角，以描述特征向量基相对于标准基的取向。

实验结果

研究问题

• RQ1 能否为 2×2 和 3×3 实对称矩阵的对角化推导出一种避免使用迭代求解器的闭式符号方法？
• RQ2 使用旋转参数时，小型对称矩阵的特征值和特征向量的精确解析表达式是什么？
• RQ3 如何利用从矩阵不变量导出的欧拉角参数化特征向量的方向？
• RQ4 所推导公式在何种符号和数值条件下保持稳定和精确？
• RQ5 所推导的公式能否系统性地验证并实现为可投入生产的 C++ 代码？

主要发现

• 2×2 特征值公式使用带符号调整的平方根表达式推导，确保了正确排序并避免了数值不稳定性。
• 对于 3×3 矩阵，特征值通过三次特征方程的三角函数解计算，其中 Δ = arccos(q / (2√p³)) 确保在 q² ≤ 4p³ 时具有实根。
• 特征向量旋转角 φ₁、φ₂、φ₃ 在符号和 π-周期性意义下唯一确定，该方法保持了特征向量的四重对称性。
• 对退化特征值情况（如 p=0 或 q²=4p³）进行了显式处理，当两个特征值相等时，φ₂ 和 φ₃ 的表达式被简化。
• 符号计算机代数程序确认了所有推导恒等式的正确性，包括向量 g1 和 g2 表达式与矩阵方程的一致性。
• C++ 实现成功验证了分析公式，展示了对 2×2 和 3×3 实对称矩阵进行快速且精确对角化的性能。