luzhuigan.m

这是线性议程组的直接解法

function [L,U,x]=LUzhuigan(A,f)
% 对于带状矩阵的追赶法
n=length(A);
a=diag(A)';
b=[diag(A,1)',0];
c=[0,diag(A,-1)'];
for i=1:n
    if i==1
        ap(i)=a(i);
        bp(i)=b(i)/ap(i);
        y(i,1)=f(i,1)/ap(i);
    else
        ap(i)=a(i)-c(i)*bp(i-1);
        if i==n
            bp(i)=0;
        else
            bp(i)=b(i)/ap(i);
        end
        y(i,1)=(f(i,1)-c(i)*y(i-1,1))/ap(i);
    end
end
% 这里 L U还不行！
L=diag(ap);
U=eye(n);
for k=n:-1:1
    if k==n
        x(k,1)=y(k,1);
    else
        x(k,1)=y(k,1)-bp(k)*x(k+1,1);
    end
end
end