五点差分格式.pptVIP

微分方程数值解 计算科学系 杨韧 第三章 椭圆型方程的差分格式 §3.1 正方形区域中的Laplace 方程  Dirichlet边值问题的差分模拟 设Ω是 xy 平面中的具有正方形边界 的 一个有界区域，考虑Laplace方程的第一边值 Dirichlet ）问题 网格节点（l ， m） 处的二阶中心差商代替 二阶微商 Laplace方程的五点差分格式(3.6)为 截断误差为O(h2)。 令 则Laplac方程的五点差分格式为（3.8） 即 例1 用五点差分格式求解 Laplace方程 在区域 内的近似解，边界值为： 取 。 解 网格点如图所示 矩阵方程AU=K，K由边界条件所确定，解得 U = [U1 U2 U3 U4 U5 U6 U7 U8 U9]’=A-1K = [55.7143 43.2143 27.1429 79.6429 70.0000 45.3571 112.8571 111.7857 84.2857]T 加密网格，取 h = 0.5 定义向量 为从左