MATLAB_数值计算_高斯-赛德尔(Gauss-Seidel)迭代的矩阵形式

在雅可比迭代时，每次迭代只用到前一次的迭代值，而在高斯-赛德尔迭代时，每次迭代充分利用当前最新的迭代值。
Gauss-Seidel迭代法百度百科
wiki百科

function[x,k]=gaussseidel(A,b,x0,ep,Nmax)
%Gauss-Seidel迭代
%x0为迭代初值
%ep为精度，默认1e-5
%Nmax为迭代次数上限以防止发散，默认为500
n=length(b);
if nargin<5 %naigin函数输入参数数目
    Nmax=500;
end
if nargin<4 
    ep=1e-5;
end
if nargin<3 
    x0=zeros(n,1);%创建零数组
end
x0=sparse(x0);b=sparse(b);A=sparse(A);%sparse创建稀疏矩阵
x=x0;x0=x0+2*ep;x0=sparse(x0);
k=0;
A1=tril(A);%返回矩阵A的下三角部分
iA1=inv(A1);%计算逆矩阵
while norm(x0-x,inf)>ep &&k<Nmax%norm求解矩阵范数
k=k+1;
x0=x;
x=-iA1*(A-A1)*x0+iA1*b;
end
x=full(x);
if k==Nmax
warning('迭代上限');
end

迭代法算法简单，程序设计容易，可以节省计算机存储空间，但是它存在着收敛性和收敛速度等问题。
rank计算矩阵的秩
diag提取矩阵对角线元素
tril提取下三角阵
triu提取上三角阵
luLU分解