一个矩阵A,A按列正交化QR分解,A=QR,Q为阶列正交矩阵,R为阶上三角阵,用修正的Gram-Schmidt正交化方法或者householder变换来实现,哪位高手给指点一下,编个matlab的子程序,实在感激不尽啊。