Doolittle分解（matlab代码）

最新推荐文章于 2023-10-07 17:15:02 发布

原创最新推荐文章于 2023-10-07 17:15:02 发布 · 5.1k 阅读

21 ·

CC 4.0 BY-SA版权

文章标签：

#matlab

matlab学习专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一个简单的LU分解算法实现过程，该算法将矩阵A分解为一个下三角矩阵L和一个上三角矩阵U的乘积。通过逐步迭代的方式更新L和U矩阵的元素，实现了对任意阶方阵的有效分解。

function [L,U] = myLU(A)
    n=length(A);
    L=eye(n,n);
    U=zeros(n,n);
    U(1,1:n)=A(1,1:n);
    L(2:n,1)=A(2:n,1)./U(1,1);

    for r=2:n
        i=r:n;
        U(r,i)=A(r,i)-L(r,1:r-1)*U(1:r-1,i);
        i=r+1:n;
        L(i,r)=(A(i,r)-L(i,1:r-1)*U(1:r-1,r))/U(r,r);
    end
end