数值分析上机实验报告八..docVIP

实验报告八 题目： 线性方程组的迭代法 摘要： 对于工程技术中产生的大型稀疏矩阵方程组，利用迭代法解是最合适的。 前言：（目的和意义） 掌握Jacobi迭代法、Gauss-Serdel迭代法、SOR方法的基本原理和应用。 数学原理： Jacobi迭代法，对于线性方程组Ax=b，如果A为非奇异方阵，即 则可将A分解为A=D-L-U，其中D为对角阵。其元素为A的对角元素，L与U为A的下三角阵和上三角阵： 于是Ax=b转化为： 与之对应的迭代公式为： Gauss-Serdel迭代法，在Jacobi迭代过程中，计算？已经得到，不必再用？，即原来的迭代公式可以改进为，于是得到Gauss-Serdel迭代公式： SOR方法，全称为逐次超松弛迭代法，它是Gauss-Serdel迭代法的一种加速方法，是解大型稀疏矩阵方程组的有效方法之一，逐次超松弛迭代公式为： 其中 程序设计： 本实验采用Matlab编写主程序如下 Jacobi迭代法： A=[5,2,1;-1,4,2;2,-3,10]; MAXTIME=50; eps=1e-4; [n,m]=size(A); x=zeros(n,1); y=zeros(n,1); k=0; disp(迭代过程X的值情况如