gs.m

有关matlab r2007的实用教程

function GS=GS(A,b,x0,eps)
%该函数文件是通过Gauss-Seidel选代格式求解线性方程组的近似解
%其中,A为方程组的系数矩阵,b为方程组的右端向量,x0为选代初值
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);    %求解出严格上对角阵
G=(D-L)\U;    %求解出迭代矩阵
f=(D-L)\b;
GS=G*x0+f;
while norm(GS-x0)>=eps
    x0=GS;
    GS=G*x0+f;
end