列正交化-施密特正交化方法

a1=[1 1 0 0];
a2=[1 0 1 0];
a3=[-1 0 0 1];
A=[a1',a2',a3];
u=A;
[m,n] = size(u);
y = u;
y(:,1) = y(:,1)/norm(y(:,1));
for k = 2:n   
    y(:,k)=u(:,k);
end
for k = 2:n   
    for j=k:n
        y(:,j)=y(:,j)-y(:,k-1)'*y(:,j)*y(:,k-1);
    end
    p1 = y(:,k)/norm(y(:,k));
    y(:,k) = p1;
  
end

对矩阵A的列向量，按照顺序逐步正交化，采用施密特方法。

疑问1.施密特正交化首选第1列，再选2列，这样正交，如果选择的顺序不一样，得到的结果是不是也不一样？

         回答：是的，可以选择0向量较多的向量作为第1个向量，再选择0向量第2多的向量作为第2个向量，这样计算量小。

        2.用一个矩阵，采用不同的顺序正交化，最后的结果之间有什么关系，只是列向量产生的空间相同吗？

         回答：应该只是列向量产生的列空间相同，还可以使用MATLAB函数Orth，得到列空间相同的正交矩阵。

      

A =

     1     1    -1
     1     0     0
     0     1     0
     0     0     1

结果y

y =

    0.7071    0.4082   -0.2887
    0.7071   -0.4082    0.2887
         0         0.8165    0.2887
         0         0               0.8660

y'*y

ans =

    1.0000    0.0000   -0.0000
    0.0000    1.0000    0.0000
   -0.0000    0.0000    1.0000

 

 参考文献 蔡改香：关于施密特正交化的一点注释与应用