矩阵的正交-三角分解有四种方法：Householder变换,Givens变换，Gram-Schmidt正交化变换、修正的Gram-Schmidt正交化变换，这四种分解的算法在高惠璇的《统计计算》（或者是统计计算方面的书）里面有，老师要求我们根据算法写出R程序，可是，我算法还不是很明白，不懂写...请高手们赐教，不管是给出Ｒ程序、或者是矩阵分解算法方面的指教，都可以。谢谢了！

自己写呗 每个方法都是几句话就够了 。。。

回复 第1楼 的 ieieya：这是gram-schmidt分解的两个方法，normal是一个矩阵，你参考一下自己再改改。 算法参考维基。

<br />
mean.normal = colMeans(normal)<br />
#### gram-schmidt general method<br />
U = normal.center = sweep(normal,2,mean.normal)<br />
for(j in 2:nvar)<br />
  for(i in 1:(j-1))<br />
    U[,j] = U[,j] - sum(normal.center[,j]*U[,i])/sum(U[,i]*U[,i])*U[,i]</p>
<p>#### gram-schmidt numerical stability method<br />
    U = normal.center = sweep(normal,2,mean.normal)<br />
    for(j in 2:nvar)<br />
      for(i in 1:(j-1))             #gram-schmidt numerical stability method<br />
        U[,j] = U[,j] - sum(U[,i]*U[,j])/sum(U[,i]*U[,i])*U[,i]<br />
</p>