lu.m

来自「matlab 多个计算方程租的算法汇集 主要讲的是一般方法 还有正定矩阵的解法」· M 代码 · 共 31 行

function [l,u]=lu(a)
n=length(a);
l=eye(n);
u=zeros(n);
for i=1:n
    u(1,i)=a(1,i);
end
for i=2:n
    l(i,1)=a(i,1)/u(1,1);
end

for r=2:n
    %%%%
    for i=r:n
        uu=0;
        for k=1:r-1
            uu=uu+l(r,k)*u(k,i);
        end
        u(r,i)=a(r,i)-uu;
    end
    %%%%
    for i=r+1:n
        ll=0;
       for k=1:r-1
                ll=ll+l(i,k)*u(k,r);
            end
            l(i,r)=(a(i,r)-ll)/u(r,r);
        end
     %%%%
    end