LU分解的MATLAB实现

LU分解的基本MATLAB实现，有DEMO的演示和可输入参数的代码。

A = [145;41836;31630];
n = size（A1）;
U = zeros（n）;
L =zeros（n）;
P=0;
for i =1:n
    L（ii） =1;
end;
for i =1:n
    U（1i）=A（1i）;
end;
for i=1:n
    L（i1） = A（i1）/U（11）;
end;

i=2;
while i<=n
    for j=i:n
        k=1;
        while k<=j-1
            P=L（ik）*U（kj）+P;
            k=k+1;
        end;
        U（ij）=A（ij）-P;
        P= 0;
    end;
     for j = i:n
        for k = 1:j-1
            P=L（jk）*U（ki）+P;
        end;
        L（ji）=（A（ji）-P）/U（ii）;
        P= 0;
    end;
    i=i+1;
end;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        558  2014-10-08 10:38  LU\Demo.m

     文件        768  2014-10-08 10:42  LU\LU.m

     文件        125  2014-10-08 10:44  LU\readme.txt

     目录          0  2014-10-08 10:44  LU

----------- ---------  ---------- -----  ----

                 1451                    4