数值线性代数第二版徐树方高立张平文上机习题第四章实验报告.docxVIP

第四章上机习题 考虑两点边值问题 d 2 y dy a, 0 a 1 dx2 dx y(0) 0, y(1) 1. 容易知道它的精确解为 x y 1 a1 (1 e ) ax e 为了把微分方程离散化，把 [0,1] 区间 n 等分，令 h=1/n ， xi ih, i 1, , n 1 得到差分方程 yi 1 2yi yi 1 yi 1 yi a, h2 h 简化为 ( h) yi 1 ( 2 h) yi yi 1 ah2 , 从而离散化后得到的线性方程组的系数矩阵为 (2 h) h ( 2 h) h A (2 h) h ( 2 h) 对 1, a 1/ 2, n 100, 分别用 Jacobi迭代法， G-S迭代法和 SOR迭代法求线性方 程组的解，要求有 4 位有效数字，然后比较与精确解得误差。 对 0.1, 0.01, 0.0001, 考虑同样的问题。 解 （1）给出算法： 为解 Ax b ，令 A D L U ，其中 A [ aij ] ， D diag (a11 ,a22 , , ann ) 0 a21 0 L a31 a32 0 ， an1 an2 an,n 10 - 1 - 0 a12 a13 a1n 0 a23 a2n U an 1,n ， 0 0 利用 Jacobi 迭代法， G-S迭代法， SOR迭代法