QR算法的Matlab 程序,三种实现矩阵QR分解的算法与程序

To learn, to share, to debate, then comes

progress.

--------------------------------------------------------------------------------------------------

QR分解

若非奇异矩阵A能够分解为一个正交矩阵Q与非奇异上三角矩阵R的乘积，即：

A=QR。则称其为A的QR分解。

实现QR分解的方法有很多种，包括Givens变换法，Householder变换法，Schemit正交化法。具体原理可以参考《矩阵论》(程云鹏，西工大出版)教材，这里仅给出三种实现QR分解的程序。

1.基于Givens变换的QR分解

%***********************************************

%矩阵分析与应用-利用Givens变换实现矩阵的QR分解

%输入：A-NxN的矩阵 输出：Q-正交阵

R-上三角阵

17.4 by Howie

function [Q,R]=md_qrg(A)

% if

size(A,1)~=size(A,2)

% disp('please input a matrix by N x N');

% return

%

end

[N,M]=size(A);

R=zeros(N);

T=eye(N);

B=A;

for

j=1:N-1

Tj=eye(N+1-j