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

中北大学理学院 课 程 设 计 题目：用逐次松弛法求方程组的解 课程：数值计算方法 成员： 1107014124 董强 1107014126 李迎 1107014128 冯梦文 题目 用逐次松弛法求方程组的解 具 体 内 容 已知方程组： 取初值，要求。 要求 1.认真读题，了解问题的数学原形； 2.选择合适问题求解的数值计算方法； 3.设计程序并进行计算； 4.对结果进行解释说明 采用方法 及结果 说明 根据题目要求用松弛法求解给出方程组的解。 1、松弛法是在高斯迭代法的基础上为提高收敛速度，采用加权平均而得到的新算法再由与加权平均得 取初值，要求。 一、问题分析 求解这个方程组，从题目中我们可以看出，有三种方法，雅可比迭代法、高赛德尔迭代法、逐次超松弛法。但是迭代发收敛速度太慢，就会增加计算量，而失去使用价值。逐次超松弛(Successive?Over?Relaxation)迭代法，简称SOR迭代法，它是在GS法基础上为提高收敛速度，采用加权平均而得到的新算法。它是大型解稀疏矩阵方程组的有效方法之一，具有计算公式简单，程序设计容易，占用计算机内存少等优点，但需要选择加速因子。 超松弛迭代法的理论基础 逐次超松弛(Successive Over Relaxation)迭代法，简称SOR迭代法，它是在GS法基础上为提高收敛速度，采用加权平均而得到的新算法，设解方程(7.1.3)的GS法记为　　　(1)再由与加权平均得 这里ω＞0称为松弛参数，将()代入则得　　　 () 该法称为SOR迭代法，[WTBX]ω＞0称为松弛因子，当ω=1时()式即为法，将(写成矩阵形式，则得　　　 即　　于是得SOR迭代的矩阵表示　 取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　　　***************************制作人：董强 李迎 冯梦文*******************endl; cout　　　***********************说明：方程最多个数为50*************************endl; cout　　　**********************************************************************endl; //*********************************数据的输入******************************** /*bb:*/cout输入的方程元数endl; cinn; cout请输入方程系数矩阵:endl; for(i=0;in;i++) for(j=0;jn;j++) cinA[i][j]; cout请输入右边向量:endl; for(i=0;in;i++) ci