- 1、本文档共30页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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 差分格式的建立
考
您可能关注的文档
- 西南大学网络教育2016建筑结构CAD论述.doc
- 二年级数学上册14.9说课.doc
- 二年级下册第三单元_图形的运动(一)说课.doc
- 二年级下册数学第三单元图形的运动备课说课.doc
- 第九章二重积分部分概述.ppt
- 西南交大晶体位错论述.doc
- 二年级下册数学第五单元教学设计说课.doc
- 第九章二重积分习题课_概述.ppt
- 第九章反应器热稳定性概述.ppt
- 二年级下册语文目标答案2016.3.3说课.doc
- 招生宣传与社交媒体策略.pptx
- 2025-2030橙汁产业规划专项研究报告.docx
- 2025-2030散装盐行业市场现状供需分析及重点企业投资评估规划分析研究报告.docx
- 2025-2030板鞋行业发展分析及投资价值研究咨询报告.docx
- 2025-2030数码文化设备产业市场深度调研及发展趋势与投资研究报告.docx
- 2025-2030气体燃料内燃机行业发展分析及投资价值研究咨询报告.docx
- 2025-2030氯化锌干电池行业市场现状供需分析及重点企业投资评估规划分析研究报告.docx
- 2025-2030核磁共振设备市场前景分析及投资策略与风险管理研究报告.docx
- 2025-2030棉花行业风险投资运行分析及运作模式与投融资研究报告.docx
- 2025-2030殷瓦钢行业市场深度调研及前景趋势与投资研究报告.docx
文档评论(0)