雅可比与高斯塞德尔迭代法收敛性MATLAB实验报告.docxVIP

实验雅可比迭代法与高斯-塞德尔迭代法的收敛性和收敛速度

【实验目的】

初步认识非线性问题的迭代法与线性问题迭代法的差别，探讨迭代法及初始值与迭代收敛性的关系。

掌握迭代法求解线性方程组的基本思想和计算步骤

掌握Jacobi迭代法和Gauss-Seidel迭代的MATLAB程序实现方法。

【实验内容】

研究用雅可比迭代法与高斯-塞德尔迭代法解下列方程组Ax=

方程组1、2：

A

方程组3、4：

A

方程组5：

A

【解】:

对1有Jacobi和Gauss-Seidel迭代收敛，2的Jacobi发散，Gauss-Seidel收敛，3的Jacobi和Gauss-Seidel均发散

【程序如下】:

收敛性分析：

雅可比迭代法：

functionM=Jacobi_judge(A,n)

disp(判断Jacobi迭代法是否收敛！！！);

I=eye(n);

D=I.*A;

disp(输出D);

disp(D);

M=I-inv(D)*A;

disp(输出迭代矩阵M);

disp(M);%迭代矩阵

e=eig(M);%求特征值

disp(输出特征值);

disp(e);

fori=1:n%判断谱半径是否小于1，大于1则发散

if(abs(e(i))=1)

disp(迭代发散！！！);

break;

end

if(i==nabs(e(i)