gaussseidel.m

gauss-seidel算法实现对方程组的求解

function  [x,k]=gaussmethod(A,b,x0,N,emg)
%A是线性方程组的左端矩阵
%B是右端向量
%x0是迭代初始值
%N表示迭代次数上限，
%emg表示控制精度
%用gauss-seidel迭代法求 线性方程
%K表示迭代次数
%x表示用迭代法求得的线性方程组的近似解
n=length(A);
x1=zeros(n,1);
x2=zeros(n,1);
x1=x0;
r=max(abs(b-A*x1));
k=0;
while r>emg
  for i=1:n
     sum=0;
for j=1:n  
    if j>i 
      sum=sum+A(i,j)*x1(j);
elseif j<i
sum=sum+A(i,j)*x2(j);
end
x2(i)=(b(i)-sum)/A(i,i);
end
r=max(abs(x2-x1));
x1=x2;
k=k+1;
if k>N
disp('迭代失败，返回');
return;
end
end
x=x1;
end