func2p2.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 23 行

%func2p2.m用矩阵分解求解方程组
a=[20,2,3;1,8,1;2,-3,15];
[L,U]=lu(a)
b=eye(3,3)
n=3
X=[0,0,0]'
Y=[0,0,0]'
for i=1:n
   Y=inv(L)*b(:,i)      %计算结果，运用LU分解结果
   X=inv(U)*Y
   b(:,i)=X
end

%b =
%    0.0517   -0.0164   -0.0093
%   -0.0055    0.1237   -0.0072
%   -0.0080    0.0269    0.0665
%与直接算逆矩阵的结果一致