factorization_qr.m

来自「matlab基础程序 便于初学者学习」· M 代码 · 共 16 行

% factorization_qr.m
% 矩阵的QR分解示例
A =  [ 1     2     3
       4     5     6
       7     8     9
      10    11    12 ]
[Q,R]=qr(A)
disp('Q的特性：')
detQ=det(Q)
QQ=Q'*Q
disp('求解线性方程组：')
b=[1;3;5;7];
x0=A\b
y=Q'*b;
x=R\y