求matlab LU 分解程序要写出详尽的程序及注释

来源:学生作业帮助网 编辑:作业帮 时间:2024/11/06 07:52:59

求matlab LU 分解程序要写出详尽的程序及注释
求matlab LU 分解程序
要写出详尽的程序及注释

求matlab LU 分解程序要写出详尽的程序及注释
matlab有多种LU分解程序下面算一种:
function [L,U]=myLU(A)
%实现对矩阵A的LU分解,L为下三角矩阵
A
[n,n]=size(A);
L=zeros(n,n);
U=zeros(n,n);
for i=1:n
L(i,i)=1;
end
for k=1:n
for j=k:n
U(k,j)=A(k,j)-sum(L(k,1:k-1).*U(1:k-1,j)');
end
for i=k+1:n
L(i,k)=(A(i,k)-sum(L(i,1:k-1).*U(1:k-1,k)'))/U(k,k);
end
end
用法,在控制台输入
A=[1 2 3 -4;-3 -4 -12 13;2 10 0 -3;4 14 9 -13];