资源简介
A = QR; 分解后Q的维数和A相同。和MATLAB 自带的QR不同。在有的文献会用到这种QR分解。
代码片段和文件信息
% 对任意的矩阵A,进行另一种QR分解,A为m*n维的; Q为m*n维列正交的矩阵,R为n*n维上三角矩阵的。A = Q*R
function [Q R] = QrDecomposition( A )
[m n] = size(A);
Q = zeros( m n );
R = zeros( n n );
%% 第一列
R(11)= norm( A( : 1) );
Q(: 1) = A( : 1)/R( 1 1 );
%% 第二列到第n列
for k = 2: n
for j = 1:
评论
共有 条评论