LU分解的目的
LU分解主要应用在数值分析中,用来解线性方程和计算行列式,例如求解AX=b,如果把A分为LU的乘积,L是下三角矩阵,U是上三角矩阵,计算的过程就会非常简便
LU分解代码算例
代码部分
A = [1 0 0 0;
1 4 1 0;
0 1 4 1;
0 0 0 1];
b = [0;3; -6; 0];
%求解三角方阵
U=zeros(4);
L=zeros(4);
for i=1:4
L(i,i)=1;
end
U(1,1)=A(1,1);
for r = 1 : n
for j = r : n%计算U矩阵
if r > 1;
for k = 1 : r-1
U(r,j) = A(r,j) - L(r,k)*U(k,j);
end
end
end
for i = r+1 : n%计算L矩阵
if r < n;
for k = 1 : r-1