Matble中雅可比迭代法和高斯塞德尔迭代法代码

 题目

求线性方程组的近似解及相应迭代次数

 

 要求||x^(k+1)-x^(k)||<=0.0001,初值为向量b

前提说明：A为非奇异实矩阵，b为向量，n为矩阵的阶，x0为初值，eps为精度，N为迭代次数

雅可比迭代法

代码

function x=jacobi(A,b,n,eps,N)
x=zeros(n,1);
x0=b;
k=0;
while k

运行结果

高斯——塞德尔迭代法

它所采用的数据和上述一样

function x=Gauss_Seidel(A,b,n,eps,N)
x=zeros(1,n);
x0=b';
k=0;
while k

 运行结果