微分方程边值问题地数值方法.docVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h SEQ MTSec \r 1 \h SEQ MTChap \r 1 \h 微分方程边值问题的数值方法 本部分内容只介绍二阶常微分方程两点边值问题的的打靶法和差分法。 二阶常微分方程为 当关于为线性时，即，此时 GOTOBUTTON ZEqnNum765521 REF ZEqnNum765521 \! (1.1)变成线性微分方程 对于方程 GOTOBUTTON ZEqnNum765521 REF ZEqnNum765521 \! (1.1)或 GOTOBUTTON ZEqnNum705174 REF ZEqnNum705174 \! (1.2)，其边界条件有以下3类： 第一类边界条件为 当或者时称为齐次的，否则称为非齐次的。 第二类边界条件为 当或者时称为齐次的，否则称为非齐次的。 第三类边界条件为 其中，当或者称为齐次的，否则称为非齐次的。微分方程 GOTOBUTTON ZEqnNum765521 REF ZEqnNum765521 \! (1.1)或者 GOTOBUTTON ZEqnNum705174 REF ZEqnNum705174 \! (1.2)附加上第一类，第二类，第三类边界条件，分别称为第一，第二，第三边值问题。 1 打靶法介绍 下面以非线性方程的第一类边值问题 GOTOBUTTON ZEqnNum765521 REF ZEqnNum765521 \! (1.1)、 GOTOBUTTON ZEqnNum578060 REF ZEqnNum578060 \! (1.3)为例讨论打靶法，其基本原理是将边值问题转化为相应的初值问题求解。 【原理】假定，这里为解在处的斜率，于是初值问题为 令，上述二阶方程转化为一阶方程组 原问题转化为求合适的，使上述初值问题的解在的值满足右端边界条件 这样初值问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)的解就是边值问题 GOTOBUTTON ZEqnNum765521 REF ZEqnNum765521 \! (1.1)、 GOTOBUTTON ZEqnNum578060 REF ZEqnNum578060 \! (1.3)的解。而对给定的，求 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)的初值问题可以用欧拉方法、龙格-库塔方法等初值问题的数值解法求解。 理论上是隐含的连续函数，如果已知，要使得 GOTOBUTTON ZEqnNum468364 REF ZEqnNum468364 \! (1.8)成立，可以通过求非线性方程 GOTOBUTTON ZEqnNum468364 REF ZEqnNum468364 \! (1.8)的零点来得到合适的，这可用任何方程求根的方法，例如牛顿法、或者其它迭代法。 实际上，是很难找到的，因此必须寻找满意的离散解数值解。 下面叙述打靶法的计算过程：（这里为允许误差，的修改使用线性插值方法） Step 1：先设，求解初值问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)，得到； 若，则为问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)的满意的离散解，结束； Step2： 若时，令，求解初值问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)，得到； 若，则为问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)的满意的离散解，结束； 否则转Step3； Step3：由线性插值得到一般计算公式 Step4： 令，求解初值问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)，得到； 若，则为问题 GOTOBUTTON ZEqnNum415675 REF ZEqnNum415675 \! (1.7)的满意的离散解，结束； 否则转Step3。 这个过程好比打靶，为子弹发射率，为靶心，当时则得到解，故称打靶法。 【例1】用打靶法求解非线性两点边值问题 要求误差。精确解为 。 【解】：首先将原问题化成初值问题