原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第2章
椭圆型方程的差分法
偏微分方程的差分法
本章我们介绍偏微分方程定解问题的数值方法。
由于偏微分方程定解问题是表述自然与工程技术领域中各种现
象最重要的数学工具之一，应用十分广泛，因此，偏微分方程定解
问题的数值求解的讨论就构成了微分方程数值方法的主要内容。
我们主要讨论三类典型方程(椭圆、抛物和双曲型)定解问题的
差分法。
内容一般围绕以下几个方面的问题：
1.差分格式的构造（差分逼近问题）；
2.相容性、稳定性及收敛性问题；
3.差分格式的求解问题。
椭圆型方程的有限差分法
（1）一维情形—常微分方程两点边值问题
d ⎛ d u(x) ⎞ d u(x)
− p x +r x +q x u x f x
( ) ( ) ( ) ( ) ( )
⎜ ⎟ a x b
d x ⎝ d x ⎠ d x
u (a) αu (b) β
（2）二维情形—二阶椭圆型方程边值问题
⎛ ⎞ ⎛ ⎞
∂ ∂u ∂ ∂u
− ⎜p ⎟+ ⎜p =+qu ⎟ f ( )
⎜ ⎟ ⎜ ⎟ x y , ∈Ω
∂xx ⎝ ∂y ⎠ ∂ ⎝ ∂y ⎠
( )
u x (y , ) x y , ϕ( )x y , ∈Γ
差分格式构造的基本问题
（1）对求解区域作网格剖分（连续区域离散化）
一维情形：分割成等距或不等距小区间（单元）；
二维情形：分割成均匀或不均匀矩形、三角形等；
（2）构造逼近微分方程定解问题的差分格式（连续
算子的离散化）
采用的方法：
① 直接差分化；② 积分插值法(有限体积法）；
§1 差分逼近的基本概念
考虑二阶常微分（椭圆）方程两点边值问题：
2
d u x
Lu =− ( )+q x u x fa x x b （2.1.1 ）
2 ( ) ( ) ( )
dx
u a α u b β （2.1.2 ）
( ) ( )
其中f (x),q (x)≥0为区间[a,b]上的连续函数,α, β 为给定常数 。
现将区间[a,b]离散化为N等份，分点为x =a+ih,i=0,1,…,N
i
b −a
h N 。从而得到关于I=[a,b]的一个网格剖分。 xi称为网格节点,
h称为步长。 再将方程（2.1.1 ）在节点xi处离散化。
为此，对于充分光滑的解u (x)在x x 处进行Taylor展开，得
i
文档评论（0）