椭圆型偏微分方程实验报告课件资料.pdfVIP

实验报告 实验报告 实验项目名称 椭圆型偏微分方程 实验项目名称 椭圆型偏微分方程 实 验 室 数 学 实 验 室 实 验 室 数 学 实 验 室 所属课程名称 微分方程数值方 所属课程名称 微分方程数值方 实 验 类 型 算 设计 实 验 类 型 算 设计 实 验 日期 2014年6月6 日 实 验 日期 2014年6月6 日 班 级 班 级 学 号 学 号 姓 名 姓 名 成 绩 成 绩 实验概述 ： 实验概述 ： 【实验目的及要求】 【实验目的及要求】 实验 目的是通过分析 Possion 问题并用交替迭代法来求解其次边值问 实验 目的是通过分析 Possion 问题并用交替迭代法来求解其次边值问 题，进一步了解交替迭代法的算法特点——即在矩形区域上的差分格式可以大 题，进一步了解交替迭代法的算法特点——即在矩形区域上的差分格式可以大 大降低计算量。实验要求是利用 Peaceman-Rachford 迭代格式编写出相应的代 大降低计算量。实验要求是利用 Peaceman-Rachford 迭代格式编写出相应的代 码解决 Possion 问题。 码解决 Possion 问题。 【实验原理】 【实验原理】 对于简单的椭圆型偏微分方程 Poission 方程: 对于简单的椭圆型偏微分方程 Poission 方程: Ω 采用正方形网格剖分正方形区域 Ω ,对 x 和 y 方向采用中心差分并记 采用正方形网格剖分正方形区域 ,对 x 和 y 方向采用中心差分并记 则对 Poission 方程离散后差分格式可写成 ; 则对 Poission 方程离散后差分格式可写成 ; 改写为 改写为 由此得Peaceman-Rachford 迭代格式为 由此得Peaceman-Rachford 迭代格式为 其分量形式为 其分量形式为 将以上两步写成矩阵形式，第一步迭代为： 将以上两步写成矩阵形式，第一步迭代为： 第二步迭代为： 第二步迭代为： 这里的 gij 和 gij 分别为 这里的 gij 和 gij 分别为 迭代参数可取为： 迭代参数可取为： 实际上每个迭代步相当于解N − 1 个系数矩阵为三对角阵的N − 1 阶线性 实际上每个迭代步相当于解N − 1 个系数矩阵为三对角阵的N − 1 阶线性 代数方程组，可用追赶法求解。