数值计算 SOR方法.docVIP

一.SOR迭代法流程图： 开始 设置一维数组首地址 设置二维数组首地址 设置精度e(precision) 设置系数矩阵的阶数N 设置增广矩阵的各值 设置初始向量 设置松弛因子w 设置两解向量的差的范数x_0[i] 设置最大迭代次数并取到MAX中 设置k记录迭代次数 初值是0 Temp=temp+a[i][j]*x_k[j] Temp=temp+a[i][j]*x_0[j] X_k[i]=(x_k[i]-temp)/a[i][j] X_k[i]=(1-w)*x_0[i]+w*x_k[i] X_0[i]=x_k[i]-x_0[i] Y Matrix_category(x_0,n)precision N Y X_0[i]=x_k[i] I=i+1 K=k+1 N Kmax 迭代次数k为，解向量为： 实验名称： 松弛法 实验解题小组成员（班级：09医软（1）班）： 姚飞 参与程序的编写 闫化晴 参与搜集资料与编写程序 余雷 参与搜集资料与后期运行调试 张珊 参与程序的编写 实验内容： SOR迭代法理论： 松弛法是 Gauss －Seidel 迭代 迭代法的一种加速方法.若记 △X = X(K+1) - X(K) = LX(K+1) + UX(K) + f - X(K) 则X(K+1) = X(K) + △X，这样X(K+1) 可以看作是 X(K)加上修正项 △X而得到.若在修正项△X前面添加一个因子= 1,就是Gauss －Seidel 迭代.通过选择可使迭代法收敛的更快.松弛法简称SOR方法，它的计算格式为： 这里称为松弛因子.当 1时称为低松弛迭代，当1 2时称为超松弛迭代. 实验素材及结果： 三、SOR迭代法 例1、用SOR迭代法求解线性方程组： 取初始点，松弛因子，精度要求. 解 ：根据上面程序运行结果如下： 例2、用SOR迭代法求解方程组： 8X1 + X2 - 2X3 = 9, 3X1 - 10X2 +X3 = 19, 5X1 - 2X2 + 20X3 = 72, 取初始点X(0) = (0,0,0)T，松弛因子= 1，精度要求= 10-5. 解 ：根据上面程序运行结果如下： 结束