昆明理工数值大作业线性方程组的迭代法介绍.doc

工科研究生数值分析上机实验 实验报告 题目： 线性方程组的迭代法 学院： 专业： 年级： 201 级 学生姓名： 学号： 任课教师： k大于所给的迭代次数，对迭代法敛散性的意义； 4、体会初始解x0 ,松弛因子的选取，对计算结果的影响。 二、问题简述： 分别选用Jacobi 迭代法，Gauss-Seidol迭代法和SOR方法计算下列线性方程组的解。 方程组一: 方程组二(设对称正定阵系数阵线方程组)： 方程组三(三对角形线性方程组 )： 本题目用到的三种基本解线性方程组的迭代法简述： 设有线性方程组：Ax=b x为n维向量 A为行列标分别为i,j的元素aij素构成的n阶系数矩阵 Jacobi 迭代法： 将A分解为三部分： A=-- =D-L-U. Jacobi 迭代法之迭代公式如下： x(k+1)=D-1(L+U)x(k)+D-1b Gauss-Seidol迭代法迭代公式: x(k+1)=(D-L)-1Ux(k)+(D-L)-1b SOR方法迭代公式:() x(k+1)=(D-ωL)-1((1-ω)D+ωU)x(k)+(D-ωL)-1ωb 其中ω为SOR方法权因子。 根据实验要求设计计算流程图： Jacobi 迭代法： 否 是 否 是 Gauss-Seidol迭代法： 否 是 否 是 SOR方法： 否 是 否 是 五、实验内容： 本次实验共使用三种迭代方法采用不同精确度和相关因子对三种不同类型的线性方程组进行了上机编程计算。除了在不同精度下采用三种迭代方法进行计算不同方法同一方程的迭代次数不同，不同精度同一方程的迭代次数也不同。其中SOR方法因权因子选取不同而迭代次数不同，本实验采用同一精确值分别对三个方程对应选取不同的三组权因子进行了计算。并且采用二分法逼近的方法寻找了相对最佳权因子。 六、计算结果： 通过编程和迭代得到以下计算结果 Jacobi 迭代法 方程一： 方程一用Jacobi 迭代法 方程二： 方程三： 分析： Gauss-Seidol迭代法： 方程一： 方程一用Gauss-Seidol迭代法：算时，迭代系数矩阵谱半径大于1，故迭代法发散。 方程二： 方程三： 分析： SOR 方法 方程一： 方程一用SOR 迭代法计算时,以0.001，从(0,2)选取权因子ω ,迭代系数矩阵谱半径均大于1，选权因子为0，=1，且随着ω增大程增大趋势，故判断迭代法发散。从而方程一用此三种方法均无法计算！！！！ 方程二： 方程三： 选取迭代向量初值： 　 w=0.8 w=0.9 w=1.0 　 　 e=e-1 e=e-2 e=e-3 e=e-1 e=e-2 e=e-3 e=e-1 e=e-2 x1 2.008 2.0184 2.0018 2.0015 2.012 2.0012 2.0009 1.9995 x2 0.8832 1.0287 1.0018 1.0013 1.0125 1.0007 1.0002 1.0063 x3 -3.087 -2.9782 -2.9996 -3 -2.9936 -3.0008 -3.001 -3.0052 x4 0.0499 0.0047 -0.0016 -0.0016 -0.006 -0.0021 -0.0019 -0.0076 x5 1.0634 0.9802 0.9971 0.9975 0.9832 0.9976 0.9982 0.9909 x6 -1.9766 -2.0371 -2.003 -2.0024 -2.0175 -2.0017 -2.0011 -2