数值分析实验报告验证线性方程组迭代法实验.doc

1．实验项目名称：验证线性方程组迭代法实验 2．实验目的和要求： 1．掌握常用的几种迭代格式；2．掌握迭代收敛性及误差估计 3．实验使用的主要仪器设备和软件： 硬件：计算机 软件：Windows 和MATLAB 7.0 4．实验的基本理论和方法：（将实验所涉及的基础理论、算法原理列出） Jacobi迭代法majacobi.m，Gauss-Seidel迭代法maseidel.m 5.实验内容与步骤：（程序清单以附录的形式给出） 迭代法的收敛速度 实验题目：用迭代法分别对，解方程组，其中 (1) 选取不同的初值和不同的右端向量b，给定迭代误差，用两种迭代法计算，观测得到的迭代向量并分析计算结果给出结论； (解题提示： format rat %输出格式是分数形式，在产生矩阵之前用了命令； %X = diag(v,k) 以向量v的元素作为矩阵X的第k条对角线元素，当k=0时，v为X的主对角线；当k>0时，v为上方第k条对角线；当k<0时，v为下方第k条对角线。) (2) 取定初值和右端向量b，给定迭代误差，将A的主对角元成倍放大，其余元素不变，用Jacobi迭代法计算多次，比较收敛速度，分析计算结果并给出结论。 （解题提示：分别求出两种迭代式的迭代矩阵，计算它们的谱半径，依据谱半径的大小判断收敛性，根据谱半径的大小，选择谱半径最小的。） （1） Jacobi迭代法 实验数据记录： 在Command window 中输入： >> format rat a1=(-1/5)*ones(18,1); a2=(-1/3)*ones(19,1); a3=4*ones(20,1); a11=diag(a1,-2); a22=diag(a2,-1); a33=diag(a3); A=a11+a22+a33+a11'+a22'; b=ones(20,1) ; format long x=majacobi(A,b) 实验结果： k= 10 x = 0.29345365537489 0.31947759682957 0.33661094883147 0.33944509786893 0.34056201203451 0.34080423363128 0.34088213096119 0.34090126219855 0.34090686062291 0.34090821432777 0.34090821432777 0.34090686062291 0.34090126219855 0.34088213096119 0.34080423363128 0.34056201203451 0.33944509786893 0.33661094883147 0.31947759682957 0.29345365537489 Gauss-Seidel迭代法 在Command window 中输入： >> x=maseidel(A,b) 实验结果： k= 7 x = 0.29345360167372 0.31947755437792 0.33661092356511 0.33944509104599 0.34056202187334 0.34080425766134 0.34088216758613 0.34090131198082 0.34090692672885 0.34090830040647 0.34090832079406 0.34090698302684 0.34090139312093 0.34088226271525 0.34080435952441 0.34056212659272 0.33944519651199 0.33661102813733 0.31947765291123 0.29345368914947 实验结果分析：n=20时用Jacobi迭代法求解得到k= 10，用Gauss-Seidel迭代法进行求解得到k= 7。因为7<10。所以用Gauss-Seidel迭代法进行求解比用Jacobi迭代法来的快，迭代的次数比较少。因此Gauss-Seidel迭代法用来解本题比教好。 Jacobi迭代法 实验数据记录： 在Command window 中输入： >> format rat a1=(-1/5)*ones(198,1); a2=(-1/3)*ones(199,1); a3=4*ones(200,1); a11=diag(a1,-2); a22=diag(a2,-1