% 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