给定一个具有独立列的矩阵 A(不一定是正方形),我能够应用 Gram-Schmidt 迭代并使用 Matlab 的函数 qr
为其列空间(以正交矩阵 Q 的形式)生成正交基
A=[1,1;1,0;1,2]
[Q,R] = qr(A)
然后
>> Q(:,1:size(A,2))
ans =
-0.577350269189626 -0.000000000000000
-0.577350269189626 -0.707106781186547
-0.577350269189626 0.707106781186547
您可以验证列是否正交
Q(:,1)'*Q(:,2) equals zero and
norm(Q(:,1)) equals norm(Q(:,2)) equals 1
给定一个具有独立列的矩阵(如 A),R 中是否存在产生(Gram-Schmidt)正交矩阵 Q 的函数? R 的
qr
函数不会产生正交 Q。
最佳答案
qr
有效,但它使用独特的约定并生成一个 qr
对象,您可以使用 qr.Q
和 qr.R
进一步操作该对象:
> A
[,1] [,2]
[1,] 1 1
[2,] 1 0
[3,] 1 2
> A.qr <- qr(A)
> qr.Q(A.qr)
[,1] [,2]
[1,] -0.5773503 -5.551115e-17
[2,] -0.5773503 -7.071068e-01
[3,] -0.5773503 7.071068e-01
> qr.R(A.qr)
[,1] [,2]
[1,] -1.732051 -1.732051
[2,] 0.000000 1.414214
这是你想要的输出吗?
关于r - Gram-Schmidt 正交化,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/3238242/