MATLAB程序设计：杜立特尔三角分解法求线性方程组

揽阳_Shadows

于 2023-10-07 17:15:02 发布

阅读量1.5k

点赞数 2

分类专栏：数值分析计算方法MATLAB程序文章标签： matlab 开发语言

本文链接：https://blog.youkuaiyun.com/qq_72375109/article/details/133647574

版权

数值分析计算方法MATLAB程序专栏收录该内容

12 篇文章

订阅专栏

  function [L,U]=LU(a)
    %Doolittle三角分解法
[m,n]=size(a);
if m~=n
    error('输入m的值应与n相同')
end
L=eye(n,n);
%单位下三角阵
U=zeros(n,n);
U(1,:)=a(1,:);
L(:,1)=a(:,1)/U(1,1);
for k=2:n
    U(k,k:n)=a(k,k:n)-L(k,1:k-1)*U(1:k-1,k:n);
    %计算U的第k行元素
    L(k:n,k)=(a(k:n,k)-L(k:n,1:k-1)*U(1:k-1,k))/U(k,k);