matlab 提示秩 1,Cholesky 分解的秩 1 更新

cholupdate

Cholesky 分解的秩 1 更新

语法

R1 = cholupdate(R,x)

R1 = cholupdate(R,x,'+')

R1 = cholupdate(R,x,'-')

[R1,p] = cholupdate(R,x,'-')

说明

R1 = cholupdate(R,x)(其中 R

= chol(A) 是 A 的原始 Cholesky 分解)返回 A + x*x' 的上三角 Cholesky 因子，其中 x 是具有合适长度的一个列向量。cholupdate 仅使用 R 的对角线和上三角。R 的下三角将被忽略。

R1 = cholupdate(R,x,'+') 与 R1 = cholupdate(R,x) 相同。

R1 = cholupdate(R,x,'-') 返回 A - x*x' 的 Cholesky 因子。当 R 不是有效的 Cholesky 因子或旧矩阵不是正定矩阵时，将会报告一条错误消息，这样将没有 Cholesky 分解。

[R1,p] = cholupdate(R,x,'-') 将不返回错误消息。如果 p 为 0，则 R1 是 A - x*x' 的 Cholesky 因子。如果 p 大于 0，则 R1 是原始 A 的 Cholesky 因子。如果 p 为 1，则 cholupdate 将失败，因为旧矩阵不是正定矩阵。如果 p 为 2，则 cholupdate 将失败，因为 R 的上三角不是有效的 Cholesky 因子。

示例

A = pasca