作业帮 > 综合 > 作业

求矩阵LU分解的matlab代码

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/05/09 08:42:20
求矩阵LU分解的matlab代码
.m文件欢迎!
function[L U flag]= LU_Decom(A)
[n,m]=size(A);
if n~=m
error('The rows and the columns of the matrix A must be equal!')
return;
end
flag='OK'; L=eye(n);U=zeros(n);
for k=1:n
for j=k:n
z=0;
for q=1:k-1
z=z+L(k,q)*U(q,j);
end
U(k,j)=A(k,j)-z;
end
if abs(U(k,k))< eps
flag='failure';return;
end
for i=k+1:n
z=0;
for q=1:k-1
z=z+L(i,q)*U(q,k);
end
L(i,k)=(A(i,k)-z)/U(k,k);
end
end