資源簡介
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:
評論
共有 條評論