微分方程数值解测试试题及答案.ppt

整理得(3.36)的差分格式为(3.37) 其中 解矩阵方程AU=K，K由边值条件及非齐次项h2gl,m构成。 例1 求解椭圆型方程第一边值问题 取步长为 h = 1/3。 解 正方形网格为 (1,2) (2,2) (1,1) (2,1) u(1/3,1)=1 u(2/3,1)=1 u(1/3,0)=0 u(2/3,0)=0 u(0,2/3)=2/3 u(0,1/3)=1/3 u(1,2/3)=4/9 u(1,1/3)=1/9 差分方程为 差分方程组为 差分方程组为化解为 令 得矩阵方程AU=K U = A-1K = [ 0.2450 0.1875 0.5775 0.5222 ]T 上机作业（P.161）习题三 第8题，要求： （1）用Jacobi迭代法或Guass-Seidel迭 代法编制程序； （2）画解曲面图。 椭圆型差分方程的迭代法 微分方程数值解 例1 用五点差分格式求解 Laplace方程 在区域 内的近似解，边界值为： 取 。 解 网格点如图所示 U7 U8 U9 U4 U5 U6 U1 U2 U3 u(1,0)=20 u(2,0)=20 u(3,0)=20 u(1,4)=180 u(2,4)=180 u(3,4)=180 u(0,3)=80 u(0,2)=80 u(0,1)=80 u(4,3)=0 u(4,2)=0 u(4,1)=0 ，解得 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 例1 用五点差分格式求解 Poisson方程 在区域 内的近似解，齐次边界条件为： 取 。 解 五点差分格式为 U7 U8 U9 U4 U5 U6 U1 U2 U3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 方程组 ，解得 U = [U1 U2 U3 U4 U5 U6 U7 U8 U9]T=A-1K = [0.1719 0.2187 0.1719 0.2187 0.2812 0.2187 0.1719 0.2187 0.1719]T 例1 在单位正方形区域Ω上解Laplace方程的 Nenmann问题 解 网格节点如图所示 U7 U8 U9 顶点 U4 U5 内点 U6 边界点 U1 U2 U3 矩阵方程为 令 则矩阵方程为 一般二阶线性椭圆型方程的五点差分格式 1、 Jacobi 迭代格式求解 2、 Guass-Seidel 迭代格式求解 例1 利用五点差分格式近似求解Laplcce方程的Dirichlet问题 取步长h=1，试用Jacobi 迭代法求解，迭代误差限ε=0.1。 解 Jacobi 迭代法格式 dirichjac.m lapjacobi 运行结果 迭代次数 k = 89