ex3523.m

来自「数字信号处理的MATLAB实现」· M 代码 · 共 8 行

%例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))   % 求解