网站大量收购独家精品文档,联系QQ:2885784924

第十章偏微分方程数值解法说课.doc

  1. 1、本文档共30页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第十章 偏微分方程数值解法 偏微分方程问题,其求解十分困难。除少数特殊情况外,绝 大多数情况均难以求出精确解。因此,近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。 §1 差分方法的基本概念 1.1 几类偏微分方程的定解问题 椭圆型方程:其最典型、最简单的形式是泊松(Poisson)方程 特别地,当时,即为拉普拉斯(Laplace)方程,又称 为调和方程 Poisson方程的第一边值问题为 其中为以为边界的有界区域,为分段光滑曲线, 称为定解区域,,分别为,上的已知连续函数。 第二类和第三类边界条件可统一表示为 其中为边界的外法线方向。当时为第二类边界条件, 时为第三类边界条件。 抛物型方程:其最简单的形式为一维热传导方程 方程可以有两种不同类型的定解问题: 初值问题 初边值问题 其中,,为已知函数,且满足连接条件 边界条件称为第一类边界条件。 第二类和第三类边界条件为 其中,。当时,为第二类边界条件, 否则称为第三类边界条件。 双曲型方程: 最简单形式为一阶双曲型方程 物理中常见的一维振动与波动问题可用二阶波动方程 描述,它是双曲型方程的典型形式。方程的初值问题为 边界条件一般也有三类,最简单的初边值问题为 1.2 差分方法的基本概念 差分方法又称为有限差分方法或网格法,是求偏微分方程定 解问题的数值解中应用最广泛的方法之一。 它的基本思想是:先对求解区域作网格剖分,将自变量的连 续变化区域用有限离散点(网格点)集代替;将问题中出现的连 续变量的函数用定义在网格点上离散变量的函数代替;通过用网 格点上函数的差商代替导数,将含连续变量的偏微分方程定解问 题化成只含有限个未知数的代数方程组(称为差分格式)。如果 差分格式有解,且当网格无限变小时其解收敛于原微分方程定解 问题的解,则差分格式的解就作为原问题的近似解(数值解)。 因此,用差分方法求偏微分方程定解问题一般需要解决以下问题: (1)选取网格; (2)对微分方程及定解条件选择差分近似,列出差分格式; (3)求解差分格式; (4)讨论差分格式解对于微分方程解的收敛性及误差估计。 下面,用一个简单的例子来说明用差分方法求解偏微分方程 问题的一般过程及差分方法的基本概念。 设有一阶双曲型方程初值问题。 选取网格: -2h -h 0 h 2h 3h 首先对定解区域作网格剖分,最简单 常用一种网格是用两族分别平行于轴与轴的等距直线 ,将分成许多小矩形 区域。这些直线称为网格线,其交点称为网格点,也称为节点, 和分别称作方向和方向的步长。这种网格称为矩形网格。 对微分方程及定解条件选择差分近似,列出差分格式: 如果用向前差商表示一阶偏导数,即 其中。 方程 在节点处可表示为 其中。由于当足够小时,在式 中略去,就得到一个与方程相近似的差分方程 此处,可看作是问题的解在节点处的近似值。同初值条件 结合,就得到求问题的数值解的差分格式。 式 称为差分方程的截断误差。 如果一个差分方程的截断误差为,则称差分方 程对是阶精度,对是阶精度的。显然,截断误差的阶数 越大,差分方程对微分方程的逼近越好。 若网格步长趋于0时,差分方程的截断误差也趋于0,则称 差分方程与相应的微分方程是相容的。这是用差分方法求解偏微 分方程问题的必要条件。 如果当网格步长趋于0时,差分格式的解收敛到相应微分方 程定解问题的解,则称这种差分格式是收敛的。 §2 椭圆型方程第一边值问题的差分解法 本节以Poisson方程为基本模型讨论第一边值问题的差分方法。 2.1 差分格式的建立 考

文档评论(0)

金不换 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档