示例：实数矩阵的 Cholesky 因式分解
使用
Cholesky 函数执行实数矩阵的
Cholesky 因式分解。
|
|
要在进行布尔值比较时避免逻辑不匹配，请在计算选项下拉列表中启用约等。
|
1. 定义正定实方阵 M。
3. 应用
rank 函数以确保
M 为满秩矩阵。
4. 设置自变量 p 和 u 来控制启用/禁用旋转和上界/下界因式分解。
5. 使用 Cholesky 函数对矩阵 M 执行默认因式分解 - 带有旋转和下界的因式分解。
| 默认函数 Cholesky(M) 等同于 Cholesky(M,1,0)。 |
6. 显示 P10T x M x P10 = L10 x L10T。
该关系在逻辑上是正确的。
7. 使用 Cholesky 函数对矩阵 M 执行因式分解 - 不带有旋转和下界的因式分解 (默认)。
| 不指定自变量 u (如在 Cholesky(M, 0) 中) 与将其设置为 0 （如在 Cholesky(M, 0, 0) 中） 是等效的。 |
| 启用旋转时返回的下界矩阵 L10 与禁用旋转时返回的下界矩阵 L00 不相等。 该关系在逻辑上是错误的。 |
8. 显示 M = L00 x L00T。
该关系在逻辑上是正确的。
9. 使用 Cholesky 函数对矩阵 M 执行因式分解 - 带有旋转和上界的因式分解。
10. 显示 P11T x M x P11 = U11T x U11。
该关系在逻辑上是正确的。
11. 使用 Cholesky 函数对矩阵 M 执行因式分解 - 不带有旋转和上界的因式分解。
12. 显示 M = U01T x U01。
该关系在逻辑上是正确的。
