解线性方程组的直接方法（1）：杜利特尔LU分解MATLAB实例

最新推荐文章于 2024-07-29 00:30:00 发布

qblslion

最新推荐文章于 2024-07-29 00:30:00 发布

阅读量9.2k

点赞数 5

分类专栏： MATLAB 文章标签： LU分解MATLAB实现

本文链接：https://blog.csdn.net/qblslion/article/details/47043515

版权

MATLAB 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

如果方阵A可以写成A=LR的形式，则称A可以三角分解（LU分解）。

A=LR 被称为A的一个三角分解或者杜利特尔分解。

MATLAB程序演示：

function [L,U]=Doolittle(a)
[n,m]=size(a);e=1e-5;
if n~=m error('矩阵不是方阵');return;end
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))<e error('失败');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

eg1:

eg2: