如何使用矩阵极分解计算器
矩阵极分解计算器将方阵 A 分解为 A = UP,其中 U 为正交矩阵,P 为半正定矩阵。选择矩阵维度(2×2 或 3×3),填入各元素,矩阵极分解计算器即返回 U、P、验证乘积 U·P 以及完整计算过程。
- 选择矩阵维度 — 从下拉菜单中选择 2×2 或 3×3。
- 输入矩阵 A — 填入每个元素,支持小数和负数。
- 查看分解结果 — 矩阵极分解计算器以矩阵形式展示 U、P 和 U·P。
- 检查可逆性 — 工具显示 A 是否可逆,这决定分解的唯一性。
调整元素值,可探索 A 变化时旋转和拉伸分量如何随之改变。
公式与理论 - 矩阵极分解计算器
矩阵极分解计算器实现右极分解:
A = UP
P = √(AᵀA) (对称半正定矩阵)
U = A · P⁻¹ (A 可逆时为正交矩阵)
| 符号 | 含义 |
|---|---|
| A | 输入的 m×m 实矩阵 |
| U | 正交因子:UᵀU = I |
| P | 对称半正定因子 |
| AᵀA | Gram 矩阵;其平方根即为 P |
计算 P = √(AᵀA):
矩阵极分解计算器使用 Denman-Beavers 迭代,从 X₀ = AᵀA、Y₀ = I 出发,迭代执行:
Xₙ₊₁ = ½(Xₙ + Yₙ⁻¹)
Yₙ₊₁ = ½(Yₙ + Xₙ⁻¹)
该迭代以二次速率收敛到 P = √(AᵀA)。得到 P 后,U = A·P⁻¹。
唯一性: 若 A 可逆(行列式非零),极分解唯一;若 A 奇异,P 依然存在但 U 可能不唯一。
假设与限制
当前矩阵极分解计算器支持 2×2 和 3×3 方阵。迭代运行 30 次,对大多数条件良好的矩阵已足够。矩阵求逆过程中绝对值小于 10⁻¹⁰ 的主元被视为数值奇异。
矩阵极分解计算器的应用场景
矩阵极分解计算器适用于线性代数与工程领域的多种应用:
- 计算机图形学 — 将仿射变换矩阵分解为旋转和缩放,用于动画插值(球面线性插值 slerp)。
- 连续介质力学 — 变形梯度 F = RU(右极分解)将刚体旋转与纯拉伸分离。
- 矩阵分析 — 分别研究 U 和 P 的性质,深入理解 A 的谱特性和几何特性。
- 机器人学 — 从含噪声的姿态估计中提取 U,恢复最近的旋转矩阵。
- 数值线性代数 — 以极分解为预处理步骤,辅助计算 SVD。
使用矩阵极分解计算器后,可验证 UᵀU ≈ I(U 为正交矩阵),并通过检查 P 的特征值均非负来确认 P 的半正定性。