主要内容

Qrinsert.

将列或行插入QR分解

句法

[Q1,R1] = qrinsert(q,r,j,x)
[Q1,R1] = qrinsert(q,r,j,x,'col')
[Q1,R1] = qrinsert(q,r,j,x,'行')

描述

[Q1,R1] = qrinsert(q,r,j,x)返回矩阵的QR分解A1,在哪里A1a = q * r用专栏X之前插入A(:,j)。如果一种具有N列和j = n + 1, 然后X在最后一列之后插入一种

[Q1,R1] = qrinsert(q,r,j,x,'col')是相同的qrinsert(q,r,j,x)

[Q1,R1] = qrinsert(q,r,j,x,'行')返回矩阵的QR分解A1,在哪里A1a = q * r额外的行,X,之前插入A(J,:)

例子

a =魔术(5);[q,r] = qr(a);j = 3;x = 1:5;[Q1,R1] = qrinsert(Q,R,J,X, '行')Q1 = 0.5231 0.5039 -0.6750 0.1205 0.0411 0.0225 0.7078 -0.6966 0.0190 -0.0788 0.0833 -0.0150 0.0308 0.0592 0.0656 0.1169 0.1527 -0.9769 0.1231 0.1363 0.3542 0.6222 0.63980.2104 0.3077 0.1902 0.4100 0.4161 -0.7264 -0.0150 0.3385 0.4500 0.4961 -0.6366 0.1761 0.0225 R1 = 32.4962 26.6801 21.4795 23.8182 26.0031 19.9292 0 12.4403 2.1340 4.3271 0 0 24.4514 11.8132 3.9931 0 0 0 20.2382 10.3392 0 0 0 0 0 16.1948 0 0 0 0

返回有效的QR分解,虽然可能不同于

A2 = [A(1:J-1,:);X;a(j:结束,:)];[Q2,R2] = QR(A2)Q2 = -0.5231 0.5039 0.6750-0.1205 0.05 0.011-0.0.6966 -0.0190 0.0788 0.0.0308 0.0.0.0833 -0.0150.0656 -0.1169 0.15231 0.1169 -0.1231 0.1169 -0.1231 0.1363 -0.0.1231 0.1163 -0.3542-0.6222 0.6398 0.2104 0.1231 0.1169 -0.1231 0.1163 -0.011231 0.1163 -0.0639 -6222 0.11231 0.1169 -0.1231 0.1169 -0.1231 0.1163 -0.6398 0.11.0.3077 0.1902 -0.4100 -0.4161 -0.7264 -0.0150 -0.3385 0.4500 -0.4961 0.6366 0.1761 0.0225 R2 = -32.4962 -26.6801 -21.4795 -23.8182 -26.0031 0 19.9292 12.4403 2.1340 4.3271 0 0 -24.4514 -11.8132 -3.9931 0 0 0 -20.2382 -10.33920 0 0 0 16.1948 0 0 0 0 0

算法

Qrinsert.功能插入值的值X进入jth列(行)的R.。然后它使用一系列Givens旋转来释出非零元素R.在对角线上和下方jth列(行)。[1]

参考资料

[1] Golub,Gene H.和Charles F.范贷款。矩阵计算。第四届。巴尔的摩,MD:约翰霍普金斯大学出版社,2013年,第6.5.2-6.5.3节,第335-338页。

也可以看看

||

在R2006A之前介绍