迭代法实验报告_米瑞琪.docxVIP

PAGE PAGE 1 数值实验报告Ⅱ 实验名称 迭代法求解Poisson方程 实验时间 2016年 4 月 10 日 姓名 米瑞琪 班级 信息1303 学号 1309010304 成绩 一、实验目的，内容 1、理解并掌握三类迭代方法（Jacobi, Gauss-Siedel, SOR）迭代方法的构造原理； 2、了解矩阵在matlab中不同的存储方式会导致不同的计算效率，学会采用最高效的方式存储矩阵； 3、学会在计算机上实现迭代法，并比较不同方法的效率与误差； 二、算法描述 （一）Jacobi迭代 对于线性方程组： Au=b (1) 为了构造迭代格式，可将上式改写为： u=Tu+d (2) 假定A的对角元均不为0，将A分裂为： A=D-B (3) 其中D为： D=diag 则(1)式可写为： Du=Bu+b (5) 形成Jacobi迭代： u 为了在迭代格式中不出现B，将B=D-A代入(6)式有： u 为了保证迭代格式的收敛，需要使迭代矩阵的谱半径ρT （二）SOR迭代 对于（1）式中的线性方程组，将A分裂为： A=D-L-U (8) 其中L，U分别为A对应的上三角与下三角矩阵的负值。 由Gauss-Siedel迭代格式的中间步骤： D 引入非零参数ω作为松弛因