ts_d.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 24 行

function T=TS_D(A,b,x0,eps)
%该函数文件是通过双步选代格式求解线性方程组的近似解
%其中，A为方程组的系数矩阵，b为方程组的右端向量，xO为选代初值
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    error
    return
end
D=diag(diag(A));    %求对角阵
L=-tril(A,-1);    %求解出严格下对角阵
U=-triu(A,1);    %求解出严格上对角阵
B1=(D-L)\U;
f1=(D-L)\b;
B2=(D-U)\L;
f2=(D-U)\b;
T=B1*x0+f1;
T=B2*x0+f2;
while norm(T-x0)>=eps
    x0=T;
    T=B1*x0+f1;
    T=B2*x0+f2;
end