gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-35/Ex3523.m

    %例35-23：QR分解法求解下列方程组。
A=[3,1,-4,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; %系数矩阵A
b=[13,-9,6,0]';  %矩阵B
[Q,R]=qr(A);  %QR分解
x=R\(Q\b)     % 求解
[Q,R,E]=qr(A); %  采用QR的第2种格式分解
x1=E*(R\(Q\b))   % 求解