原创力文档- 1、本文档共48页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
Chapter5 常微分方程边值问题的数值解法
§5.1 数学模型
二阶常微分方程的边值问题如下：
控制方程：
其边界条件（Bounding Condition, B.C）可分为如下三类：
（1）第一边值条件: （1.1.1）
（2）第二边值条件： （1.1.2）
（3）第三边值条件： （1.1.3）
其中:第一边值条件给定了在边界处的函数值；第二边值条件给定了在边界处的导数值；第三边值条件为混合边界条件，给出了边界处的函数值和导数值之间的关系。
当方程 GOTOBUTTON ZEqnNum273319 REF ZEqnNum273319 \* Charformat \! (1.1)为线性方程时，其一般形式为：
设其系数、及自由项均为区间内的连续已知函数，故方程 GOTOBUTTON ZEqnNum380199 REF ZEqnNum380199 \* Charformat \! (1.2)为变系数的线性方程，为消除其中的一阶项，处理方法如下：
将 式两端乘以，即得：
由求导公式，则上式可改写为：
(1.3a)
令：，则方程（1.3a）可以新变量表为如下形式：
此形式仍为二阶变系数线性常微分方程，但其中已不显含一阶导数项，而边值条件的形式不变。这表明，对于一般形式的二阶常微分方程，总可以通过上述变量变换的处理将其转化成不含一阶导数项的形式。为此，在下面考虑线性方程的情况时，可以假设方程中不再显含项。
对于方程 GOTOBUTTON ZEqnNum201194 REF ZEqnNum201194 \* Charformat \! (1.4)在多数情况下，无法求得其解析解，故只能采用数值方法求解。
关于边值问题的数值解法有多种，本章只介绍其中比较常用的两种方法：差分法和试射法。
MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h SEQ MTSec \h §5.2 解线性边值问题的差分方法
§5.2.1差分方程的建立
现对线性方程的情况来讨论之，设方程的形式为：
1. 将区间离散化：
取等间距将区间划分为等分，诸结点为：
如此，我们将原在上求解方程 GOTOBUTTON ZEqnNum186661 REF ZEqnNum186661 \* Charformat \! (2.1)解的问题转化为求各结点上近似解值问题。为此，我们首先应将方程 GOTOBUTTON ZEqnNum186661 REF ZEqnNum186661 \* Charformat \! (2.1)中的变量进行离散化处理，其具体作法是：
2. 构造差分方程
对内部结点，利用二阶中心差商来近似原二阶微商，由二阶中心差商的关系式：
其中，将它代入方程 GOTOBUTTON ZEqnNum186661 REF ZEqnNum186661 \* Charformat \! (2.1)中，得在结点处所满足的关系式：
其中，
在上式中，若略去截断误差项（它是间距的二阶小量），则可得原微分方程 GOTOBUTTON ZEqnNum186661 REF ZEqnNum186661 \* Charformat \! (2.1)的近似差分方程为：
这是含有共个未知数的线性方程组，而方程的个数为个，欲使此方程组有唯一解，还需由两个边值条件补充两个方程。
对于第一边值条件，可直接由给定的边值条件给出两个补充方程，即：
将方程 GOTOBUTTON ZEqnNum748986 REF ZEqnNum748986 \* Charformat \! (2.3)与 GOTOBUTTON ZEqnNum135300 REF ZEqnNum135300 \* Charformat \! (2.4)联立即构成了第一边值问题完整的差分方程组。
对于第二和第三边值条件，由于两者均给出了边界处的一阶导数信息，相应于差分方程的形式，对于边界的导数值或表达式，我们亦须用差商来近似表示之。因为我们无法利用区间之外结点的信息，所以在引进两个边界导数的差分近似表达式时，就不能再利用中心差商公式。
若要求的截断误差为一阶的，即，则我们可利用简单的前、后差商表达式来近似原导数边值条件，即有：
若要求的截断误差为二阶的，即，则原导数边值条件的差商近似表示需利用Newton等距插值公式（补图）：
前插公式：
文档评论（0）