数学方程.pptVIP

数学物理方程;第八章 有限差分法; 前面各章我们介绍了求解数学物理问题显式解(严格解)的几种方法.注意到前面所研究问题的方程是简单的，且所讨论的区域边界是规则的.然而对于生活中的许多问题，其方程是较复杂的，区域边界是不规则的，对于这种问题，要想求出其显式解是很困难的.因此试着求其近似解(数值解).近似解法有摄动法、变分法、模拟法等，本章将介绍一种最常用的近似解法—有限差分法.;§8.1 有限差分法的构造;类似地， 还可以近似为;(8.4);对于偏微分形式，写成差商形式与常微分形式是类似地.因而，对于一个偏微分方程我们可以通过这种近似把它写成差分方程的形式.如二维拉普拉斯方程;当然，这种近似所带来的误差总是存在的，但只要它所产生的误差不超出我们限制的范围，那么这种方法对我们所讨论的问题就是可行.因此，这种解法一个重要的步骤是作误差估计.如对(8.5)的误差估计为：设函数u=u(x,y)充分光滑，利用泰勒公式展开到四阶导数项，可得;缘烂孩境花眨氨钮祖呕虱疼蛊蜡海尿重峙雄弘稳浇矩识升荫凳译碾打鹊展数学方程数学方程;其中ξ=x+θ1Δx，η=x-θ2Δx (0θ1,θ21);因此，用方程(8.5)近似二维拉普拉斯方程，其截断误差(由用泰勒级数的前有限项代替无穷项而造成)为(Δx)2+ (Δy)2的数量级.若令Δx=Δy=h，则截断误差为h2数量级.;近似，其截断误差为 O(Δt)+O((Δx)2).;§8.2 调和方程的差分格式;用有限差分法求解调和方程狄利克雷问题;(1)取步长 (n0整数)，用x=ih(i=0,1,…,n)及y=jh(j=0,1,…,n)这些平行于坐标轴的直线在区域上画出网格，直线的交点成为节点，其坐标为(x,y)=(ih, jh)，简记为(xi,yi)(i,j=0,1,…,n).近似处理的第一步是：不直接去求解u=u(x,y)，而是求在节点上解得近似值uij(i,j=0,1,…,n).注意到当步长h越小，网格分布就越窄，从而就会获得关于解u(x,y)越多的信息.因此，就将求区域Ω上的一个未知函数的问题化为求在节点上解的近似值问题，此问题只须求有限个;未知数.从而有限差分法将一个无限维问题转化为一个有限维问题.;况选取适当的格式，然后用迭代法求解.;整理后可得;将 作为解的一次近似，右端四个函数若有落在边界节点上，则其值以相应的边界条件中已知函数值 代入.如此迭代k次之后，可得; 时，此序列必收敛于差分方程(8.11)的解.因此当k充分大时， 就给出了所要求的近似值.一般地，对充分大的k，当相邻两次迭代解 间的误差小于某个预先给定的适当小的控制数 ，就可结束迭代过程.;§8.3 热传导方程的差分格式;用有限差分法求解一维热传导方程的初边值问题;首先，在求解区域Ωt(0≤x ≤l,t≥0)上作矩形网格.在x轴上以步长Δx=l/J(J为正整数)把[0,l]区间J等分，关于各分点做平行于t轴的网格线.在t轴以步长Δt作平行于x轴的网格线.直线的交点(xi, tj)=(iΔx，jΔt)称为网格节点.;在网格的内节点(xi, tj)处分别用;其中 表示解u(x,t)在节点的近似值. (8.14)中的初值条件与边值条件各自代为边界节点上的初值条件和边值条件.;求解定解问题(8.13)就转化为求方程组(8.14)在条件(8.15)下的解.记;由(8.16)可以清楚看出：它的解可以按t增加的方向逐排求出，即利用初边值条件(8.15)可以求出第0排上的值 ；利用(8.16)可以求出第一排上的值 ；然后由 和边值条件 ，再利用(8.16)在j=2的情形计算出 .如此逐步进行下去，可以求出所有内节点处的值 .;§8.4 波动方程的差分格式;用有限差分法求解一维波动方程的初边值问题;作两族平行线;Δx;近似替代 ，并记 ，则 (8.17)中的微分方程可化为 ;(8.17)式中的初边值条件各自化为在边界节点上的初边值条件;(8.19)下求解代数方程组(8.18).计算方法与热传导方程一样，将(8.19)中的初始条件及;以此类推，可得定解问题(8.17)的近似解.;作业 习题六 第101页第1、3、4题