该算法通过 Modified Gram-Schmidt 算法计算矩阵“A = q0*r0”的 QR 分解,并在末尾附加 Ap。 它无需重新计算即可更新现有的 QR 分解。 相对于 qrinsert 的优势在于它不需要“Q 必须是正方形”,即它允许给定的 QR 分解为“经济规模”/“瘦形式”。 输入参数 q0 和 r0 可以是空矩阵,在这种情况下,算法计算矩阵 Ap.mM 是定义一般内积的质量矩阵的新细形式 QR,将其设置为“眼睛”或不提供它可以恢复常规点积。
例子:
A = rand(10,5); B = rand(10,2); [q0, r0] = qr(A,0); % 或 qrappend([],[],A); [q1,r1] = qrappend(q0,r0,B);