nalu.m

来自「Matlab数值解算法实现代码」· M 代码 · 共 11 行

function [l,u]=nalu(a)
% 用途：求可逆方阵的Doolittile分解。
% 格式：[l,u]=nalu(a) 其中a为可逆方阵，l返回单位下三角矩阵，u返回上三角矩阵。
n=length(a);
u=zeros(n,n);l=eye(3,3);
u(1,:)=a(1,:);l(2:n,1)=a(2:n,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);
    l(k+1:n,k)=(a(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k))/u(k,k);
end