数值计算方法 用逐次松弛法求方程组的解..docVIP

中北大学理学院 课 程 设 计 题目：用逐次松弛法求方程组的解 课程：数值计算方法 成员： 1107014124 董强 1107014126 李迎 1107014128 冯梦文 题目 用逐次松弛法求方程组的解 具 体 内 容 已知方程组： 取初值，要求。 要求 1.认真读题，了解问题的数学原形； 2.选择合适问题求解的数值计算方法； 3.设计程序并进行计算； 4.对结果进行解释说明 采用方法 及结果 说明 根据题目要求用松弛法求解给出方程组的解。 1、松弛法是在高斯迭代法的基础上为提高收敛速度，采用加权平均而得到的新算法 再由与加权平均得，最终超松弛法的迭代矩阵为 2、容易判断给出的系数矩阵是对称正定矩阵，所以只要选择的加速因子w只要满足0w2, 则解方程组的SOR（超松弛法）方法收敛。 3、对于最佳松弛因子的选择方法，则是根据不同的加速因子在相同的系数矩阵下迭代次数最少的为最优，具体的迭代次数见附录，最后选择最佳加速因子w=1.3，得到结果为，并且所得的精度满足题目的要求。 成员 董强 李迎 冯梦文 已知方程组： 取初值，要求。 一、问题分析 求解这个方程组，从题目中我们可以看出，有三种方法，雅可比迭代法、高赛德尔迭代法、逐次超松弛法。但是迭代发收敛速度太慢，就会增加计算量，而失去使用价值。逐次超松弛(Successive?Over?Relaxation)迭代法，简称SOR迭代法，它是在GS法基础上为提高收敛速度，采用加权平均而得到的新算法。它是大型解稀疏矩阵方程组的有效方法之一，具有计算公式简单，程序设计容易，占用计算机内存少等优点，但需要选择加速因子。 超松弛迭代法的理论基础 逐次超松弛(Successive Over Relaxation)迭代法，简称SOR迭代法，它是在GS法基础上为提高收敛速度，采用加权平均而得到的新算法，设解方程(7.1.3)的GS法记为　　　(1)再由与加权平均得 这里ω＞0称为松弛参数，将()代入则得　　　 () 该法称为SOR迭代法，[WTBX]ω＞0称为松弛因子，当ω=1时()式即为法，将(写成矩阵形式，则得　　　 即　　 从而 k=0、1、2、、、、、于是得SOR迭代的矩阵表示　 加速因子w的选择 对于SOR法，松弛因子的选择对于收敛速度影响较大，关于最优松弛因子W的研究较为复杂，对此我们选用不同的松弛因子比较其收敛速度。 由SOR法收敛，则0W2 因此我们在（0，2）中选取不同的松弛因子并比较得出近似的最优松弛因子。 二、问题求解 取 W=1.1代入程序进行计算 取W=1.2 代入程序计算 取w=1.3带入程序计算 结果说明及分析： 不同松弛因子的不同松弛因子的迭代次数 松弛因子 迭代次数 松弛因子 迭代次数 0.1 339 1.1 19 0.2 179 1.2 14 0.3 120 1.3 12 0.4 88 1.4 15 0.5 69 1.5 19 0.6 55 1.6 25 0.7 45 1.7 35 0.8 36 1.8 51 0.9 39 1.9 96 1 28 注：w=1时，是高斯赛德尔迭代法。 最后对于不同的w取值，我们进行计算并列表比较，可以发现不同w取值迭代收敛速度不同，当w=1.3时，迭代收敛速度最快。 最终方程组的解为与精确解接近，且满足题目所要求的精度 附录：计算程序 #includeiostream #includemath.h #includeiomanip using namespace std; #define kk 50 //定义最大方程元数 int n,i,c,j,ll,hh,gg,mm,f=1; double A[kk][kk],x[kk][kk],b[kk],y[kk],a[kk],z[kk],m,nn,d,e=1,w,fff ; void main() { cout　　　**********************************************************************endl; cout　　　*********数值计算方法之《逐次超松弛法求解线性方程组》*****************endl; cout　　　***************************制