function [L,U]=doolittle_decomposition(A)
%{
杜利托尔分解
param A:二维矩阵(方阵)
return L:下三角矩阵
return U:上三角矩阵
%}
rows=size(A,1);
cols=size(A,2);
L=eye(rows,cols);
U=A;
for ii=1:rows
for jj=ii+1:rows%从当前行开始消元
temp=-U(jj,ii)/U(ii,ii);
U(jj,:)=U(jj,:)+temp*U(ii,:);
%对单位阵做另一行变换
L(ii,:)=L(ii,:)-temp*L(jj,:);
end
end
L=L';%取转置
end
杜利托尔分解
最新推荐文章于 2024-09-22 22:02:19 发布