微分方程数值解法实验报告.docxVIP

. . WORD WORD 版本. 微分方程数值解法实验报告 ： 班级： 学号： 一：问题描述 求解边值问题： ? ???u ? 2?e?( x? y) (sin? x cos? y ? cos? x sin? ? ? （? x,y)?G ? (0,1)?(0,1) ? ? ??u ? 0,( x, y) ?? ? ? 其精确解为u(x, y) ? e?( x? y ) sin(?x) sin(?y) 问题一：取步长 h=k=1/64,1/128,作五点差分格式，用 Jacobi 迭代法，Gauss_Seidel 迭代法，SOR 迭代法（w=1.45）。求解差分方 程，以前后两次重合到小数点后四位的迭代值作为解的近似值，比较 三种解法的迭代次数以及差分解u (x, y)(h ? 1/ 64,1/128) 与精确解的精度。 h 问题二：取步长h=k=1/64,1/128,作五点差分格式，用单参数和双参数 PR 法解差分方程，近似到小数点后四位。与SOR 法比较精度和迭代步数。 问题三：取步长h=k=1/64,1/128,作五点差分格式，用共轭梯度法和预处理共轭梯度法解差分方程，近似到小数点后四位。与 SOR 法与PR 法比较精度和迭代步数。 二．实验目的： 分别使用五点差分法（Jacobi 迭代，Gauss_Seidel 迭代，SOR 迭代），PR 交替隐式差分法（单参数，双参数），共轭梯度法，预共轭 梯度法分别求椭圆方程的数值解。 三．实验原理： Jacobi 迭代法设线性方程组  Ax ? b (1) 的系数矩阵A 可逆且主对角元素a ,a ,...,a 均不为零,令 ? 11 22 ? nn 并将A 分解成 D ? diag ? a ,a 11 22 ? ,...,a nn 从而(1)可写成 A ? A ? D ? D (2) Dx ? ?D ? A?x ? b 令 x ? B x ? f 1 1 其中B ? I ? D?1 A, f ? D?1b . (3) 1 1 以1B 为迭代矩阵的迭代法(公式) 以 1 x?k ?1? ? B 1  x?k ? ? f 1  (4) 称为雅可比(Jacobi)迭代法(公式),用向量的分量来表示,(4)为 x ( k ?1 ) ? i ? ? 1 ? b a i ii ?n a i j j ?1 j ?i x ( k ) ? j (5) ? i ? 1,2,...n, k ? 0,1,2,... ? ? 其中x?0 ? ? x?0 ? ,x?0 ? ,...x?0 ? T 为初始向量. 1 2 n Guass-Seidel 迭代法 i由雅可比迭代公式可知,在迭代的每一步计算过程中是用x ?k ?的全部分量来计算x?k ?1? 的所有分量,显然在计算第i 个分量x ?k ?1?时,已 i 经计算出的最新分量x ?k ?1? ,...,x  i?1 ?k ?1? 没有被利用，从直观上看,最新计 1算出的分量可能比旧的分量要好些.因此，对这些最新计算出来的第 1 k ? 1 x?k ?1??x ?k ?1? 次近似 的分量 j 加以利用,就得到所谓解方程组的高斯 —塞德（Gauss-Seidel）迭代法. 把矩阵A 分解成 A ? D ? L ? U (6) 其中D ? diag ? a ,a 其中 11 22 ,...,a nn ?, ? L,?U 分别为 A 的主对角元除外的下 ? ? ? D ? L x ? Ux ? b 即 x ? B 2 其中  x ? f 2 B ? ?D ? L??1U , f B2 2 B 以 2 为迭代矩阵构成的迭代法(公式) ? ?D ? L??1 b  (7) x?k ?1? ? B 2 x?k ? ? f 2 (8) 称为高斯—塞德尔迭代法(公式),用 量表示的形式为 x( k ?1 ) ? 1 ? b ?? i a i ? ? ?i?1 a ij x( k ?1 ) j ?n a ij x( k ) ? j j ?1 j ?i?1 ii ? i ? 1,2, n, k ? 0,1,2,... SOR 迭代 T ? (D ??L)?1 (?R ? (1? ?)D) d ? ?(D ??L)?1 b ? opt ? 2 /(1 ? )u ? cos?h 1 ? ? 1 ? ? 2 u i?1, j u i, j ?1 2u 2u  i, j i, j u ? L i?1, j 1 u ? L i, j ?1 2 L ? L ? L 1 2 迭代格式为： (I ?? k L )u 1 k ? 1 ? (I ?? 2k 2 L )uk 2 ?? b k (I ?