边值问题的数值解法在具体求解常微分方程时.PPTVIP

第八章常微分方程数值解法 8.6 边值问题的数值解法 在具体求解常微分方程时，必须附加某种定解条件。定解条件通常有两 种，一种是初始条件，另一种是边界条件。与边界条件相应的定解问题称为 边值问题。本节介绍求解两点边值问题 （8.6.1） （8.6.2） 的数值解法。当 关于 和 是线性时，式（8.6.1）为线性两点边值问题 8.6.1 打靶法 打靶法 的基本原理是将两点边值问题（8.6.1）转化为下列形式的初值问题 （8.6.3） 这里的 为 在 处的斜率。令 ，上述二阶方程可降为一阶方程组 （8.6.4） 因此，边值问题变成求合适的 ，使上述方程组初值问题的解满足原边值问题 的右端边界条件 ，从而得到边值问题的解。这样，把一个两点边值问 题的数值解问题转化为一阶方程组初值问题的数值解问题。方程组初值问题的 所有数值方法在这里都可以使用。问题的关键是如何去找合适的初始斜率的试 探值 。 对给定的 ，设初值问题（8.6.3）的解为 ，它是 的隐函数。 假设 随 是连续变化的，记为 ，于是我们要找的 就是方程 的根。可以用第6章的迭代法求上述方程的根。比如用割线法有 （8.6.5） 这样，可以按下面简单的计算过程进行求解。先给定两个初始斜率 ， 分别作为初值问题（8.6.4）的初始条件。用一阶方程组的数值方法求解它 们，分别得到区间右端点的函数的计算值 和 。如果 或 ，则以 或 作为两点边值问题的解。否则 用割线法（8.6.5）求 ，同理得到 ，再判断它是否满足精度要求 。如此重复，直到某个 满足 ，此时得到 的 和 就是边值问题的解函数值和它的一阶导数值。上述方程 好比打靶， 作为斜率为子弹的发射， 为靶心，故称为打靶法。 值得指出的是，对于线性边值问题（8.6.2），一个简单又实用的方法是用解 析的思想，将它转化为两个初值问题： 求得这两个初值问题的解 和 ，若 ，容易验证 (8.6.6) 为线性两点边值问题（8.6.2）的解。 例 8.7 用打靶法求解线性边值问题 其解的解析表达式为 。 解 先将该线性边值问题转化为两个初值问题 令 ，将上述两个边值问题分别降为一阶方程组初值问题 的打靶法计算值 ，部分点上的计算值、精确值和误差列于表 8-12。 取 h=0.02 ，用经典R-K法分别求这两个方程组解 和 的计算值 和 ，然后按（ 8.6.6 ）得精确解 i y 1 表 8-12 0 0 0 0 0 0 0.2 -0.002407991 0.204007989 0.2016000053 0.2016000000 0.4 -0.006655031 0.432255024 0.4256000080 0.4256000000 0.6 0.019672413 0.709927571 0.7296000083 0.7296000000 0.8 0.145529585 1.064070385 1.2096