chase.m

来自「线性方程组的数值解法」· M 代码 · 共 23 行

function perchase(A,B)
%追赶法，参数说明(系数矩阵A,常数项B)
d=max(size(A));
L=zeros(d,d);         %LU分解
U=zeros(d,d);
Y=zeros(d,1);           %设置中间过渡方程组
U(1,1)=A(1,1);
L(1,1)=1;
for i=2:d
    L(i,i)=1;
    L(i-1,i)=A(i-1,i)/U(i-1,i-1);
    Y(1,1)=B(1);
    Y(i,1)=B(i)-Y(i-1,1)*L(i-1,i);
    U(i,i)=A(i,i)-L(i-1,i)*A(i,i-1);     
end

X=zeros(d,1);                   %求解
X(d,1)=Y(d,1)/U(d,d);
for n=d-1:-1:1
    X(n,1)=(Y(n,1)-A(n+1,n)*X(n+1,1))/U(n,n);
end

Result=X