作业帮求matlab LU 分解程序
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/05/05 23:33:39
求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];
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];