3线性方程组迭代法.docVIP

桂林电子科技大学 数学与计算科学学院实验报告 实验室： 实验日期：2010年 4月17日 学号 0800730121 批次 3 姓名 彭兆拓 成绩 课程 名称 数值分析 实验项目 名 称 线性方程组迭代法 指导 教师 李光云 一，实验目的 通过实验,熟悉线性方程组的列选主元高斯消去解法,.LU分解法,Jacobi迭代法,Gauss-Seidel迭代法和SOR迭代法 二，实验原理 求解大型线性方程组的方法一般有直接法和迭代法.迭代法和直接法相比,具有程序简单,占用存储单元少等优点,适用于求解高阶,特别是零系数比较多的稀疏线性方程组. 1. 雅可比迭代法 雅可比迭代法属于迭代法求解线性方程组,它的求解思想是由线性代数方程组改写,建立迭代格式: 当给出一个初始向量,可根据该迭代格式得到一个向量的迭代序列,当迭代序列收敛于时候,这就是方程组的解. 编写雅克比迭代法的matlab程序,应利用该方法的矩阵形式,充分发挥matlab的优势.记 则 , 则雅可比迭代格式可表示为 2. Gauss-Seidel迭代法 在雅可比迭代法的第i个公式中利用,已经算出来的替换公式中的,得到Gauss-Seidel迭代法 其矩阵形式可表示为 ,或者 3. SOR迭代法 为了达到更快的收敛速度,在Gauss-Seidel迭代法中引入松弛因子,其矩阵形式为 SOR迭代法的松弛因子应在区间内取值,当时,称为完全松弛因子, 当时,称为超松弛因子,当时,称为次松弛因子,如何选取松弛因子使得收敛速度加快是讨论SOR迭代法的一个重要内容,目前只对少数模型问题,才有确定使得收敛速度打到最优的松弛因子的理论公式. 三，使用仪器，材料 计算机一台;Mathematica软件. 四，实验内容与步骤 1.用雅克比和高斯-塞德尔迭代法计算一下各题,求出各个方法的迭代矩阵,比较两个方法的收敛速度. (1) (2) 2.请用SOR法计算以上两题,并且找到比雅克比和高斯-塞德尔迭代法收敛速度更快的松弛因子. 五，实验过程原始记录（数据，图表，计算等） 1、 （1）高斯-塞德尔迭代法： （2）雅克比迭代法： （2）高斯-塞德尔迭代法： 从实验数据得出，高斯-塞德尔迭代法要比雅克比迭代法的收敛速度要快 2、 （1）A={{3,-1,1},{3,6,2},{3,3,7}}; b={1,0,4}; x0={0.,0.,0.}; e=10^(-4); w=1; 松弛因子w=1时结果运算如下： （2）题的结果如下图所示： 六，实验结果分析或总结 总体感觉此两次实验难度比较难，由于时间问题，很多知识都无法掌握，加上实验无指导书，对于一些不知道的问题还是无从下手。有不少问题还是不清楚是如何得来的，对于一些关键词还是没能理解。 7