nalupad.m

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

function a=nalupad(a)
% 用途：求可逆方阵的Doolittile分解,紧凑格式。
% 格式：a=nalupad(a) 其中a为可逆方阵，返回下三角矩阵部分为l，上三角矩阵部分为u。
n=length(a);
a(2:n,1)=a(2:n,1)/a(1,1);
for k=2:n
    a(k,k:n)=a(k,k:n)-a(k,1:k-1)*a(1:k-1,k:n);
    a(k+1:n,k)=(a(k+1:n,k)-a(k+1:n,1:k-1)*a(1:k-1,k))/a(k,k);
end