数值方法实验三 - 副本.doc

数学与计算科学学院 实 验 报 告 实验项目名称 迭代法解线性方程组 所属课程名称 数值方法A 实 验 类 型 验证型 实 验 日 期 2014.12.03 班 级 信计12- 学 号 201253100 姓 名 成 绩 一、实验概述： 【实验目的】 1）了解Gauss-Seidel迭代法和松弛迭代法的迭代思想； 2）熟悉掌握Gauss-Seidel迭代法和超松弛迭代法的基本原理和基本方法； 3）学会用逐次Gauss-Seidel迭代法和超松弛迭代法解简单的方程组； 4）在上机操作过程中体会数值计算方法的收敛性及收敛速度； 5）掌握判别迭代法收敛的基本方法及其运用 6）培养利用计算机和所学理论知识分析解决实际问题的能力。 【实验原理】 1.Gauss–Seidel迭代法 迭代格式: （若可逆） ， ，，为初始向量（矩阵形式） ，（分量形式） 2.弛迭代法(简称SOR) 迭代格式: 在Seidel迭代公式中， =， 。 若在修正项前加上一个参数w ，便得到松弛法的迭代公式： ，，k=0, 1, 2, … 将上式写成矩阵形式得 ， k=0, 1, 2, … 其中w 叫做松弛因子, 当w 1时叫做超松弛, w 1时叫做低松弛。 注:（1）w =1时就是Seidel迭代法； （2）适当选取松弛因子w , 可加快收敛速度； （3）关于迭代格式（3），引进收敛速度概念： 【实验环境】 1.硬件环境 2.软件环境 （1） （2）VC++ 6.0， 二、实验内容： 【实验方案】 1．Gauss-Seidel迭代法方案： 用Gauss-Seidel迭代法求解下列方程组，当满足时结束迭代。 (1)输入系数矩阵A，右端项b，矩阵阶数n。 最大误差限，最大迭代步数N及初始迭代向量?? (2)置??????? 1)令最大误差? 2)对于????????????? a) ?? b) ??? c)判断是否成立,成立则；（比较本次与上次迭代中的对应的值，寻找最大误差）? 3)判断是否成立? a)成立，找到近似解，则输出迭代步数和近似解，结束程序。?? ? b)不成立，则判断kN是否成立? ⅰ成立，则k=k+1，回到步骤1)? ⅱ否则，输出失败信息，结束程序. 2.松弛迭代法方案： 用松弛迭代法（取）求解下列方程组，当满足时结束迭代，说明迭代的收敛性，并讨论初值的变化对收敛的影响。 (1)输入系数矩阵A，右端项b，矩阵阶数n。 最大误差限，最大迭代步数N及初始迭代向量?? (2)置??????? 1)令最大误差? 2)对于????????????? a) ?? b) ??? c)判断是否成立,成立则；（比较本次与上次迭代中的对应的值，寻找最大误差）? 3)判断是否成立? a)成立，找到近似解，则输出迭代步数和近似解，结束程序。?? ? b)不成立，则判断kN是否成立? ⅰ成立，则k=k+1，回到步骤1)? ⅱ否则，输出失败信息，结束程序. 4）根据运行结果说明迭代的收敛性 5) 讨论初值的变化对收敛的影响。 【实验过程】（实验步骤） 一、实验步骤 1、上机前想好解决问题的思路与方法； 2、上机前写好解决问题的算法； 3、根据算法编写解决问题的程序； 4、上机时，将事先写好的程序输入到VC++的编辑框； 5、调试所输入的程序并逐步改正错误； 6、调试无误后，运行程序并得出最终结果； 7.说明收敛性及初值变化对收敛的影响。 二.收敛性及收敛速度比较 1）SOR迭代法是G-S迭代法的一种加速方法，故前者收敛速度要快些； 2）由G-S迭代法和SOR迭代法求解的程序运行结果可知初值的选取影响迭代格式的收敛快慢； 3）对于SOR迭代法，收敛快慢依初值和松弛因子的变化而不同。 【实验结论】（结果） 1.G-S迭代结果： 2.SOR迭代结果： 1）初值X取（0,0,0）时： 2）初值X取（1,0,0）时: