ex3522.m

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

%例35-22：LU分解法求解下列方程组。
A=[3,1,-4,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; %系数矩阵A
b=[13,-9,6,0]';  %矩阵B
[L,U]=lu(A);  %LU分解
x=U\(L\b)   % 求解
[L,U ,P]=lu(A); %  采用LU的第2种格式分解
x=U\(L\P*b)   % 求解