[微分方程OK]资料.pptVIP

第七章 常微分方程的数值解法 科学和工程技术中常常要求解常微分方程。根据实际背景不同,所遇到的问题可分为两类:一类是常微分方程初值问题；一类是常微分方程边值问题。一般地,要找出这两类问题的解析解往往非常困难,甚至是不能的。本章将介绍它们的数值解法。所谓数值解法,就是在没有办法知道未知函数的解析表达式的情况下,我们近似计算未知函数在其定义域中的某些离散点上的函数值。 当然,如果这些离散点在函数的定义域内的分布很密,且相应点的函数值的计算又非常准确,那么就意味着基本上找到了微分方程的解。 函数是事物的内部联系在数量方面的反映,如何寻找变量之间的函数关系,在实际应用中具有重要意义。在许多实际问题中，往往不能直接找出变量之间的函数关系,但是根据问题所提供的情况，有时可以列出含有要找的函数及其导数的关系式。这就是所谓的微分方程，从而得出微分方程模型。 例1 物体冷却过程的数学模型 将物体放置于空气中，在时刻 时，测量的它的温度为 , 10分钟后测量得温度为 。 我们要求此物体的温度 和时间 的关系,并计算20分钟后物体的温度。这里我们假定空气温度保持为 。 注意到热量总是从温度高的物体向温度低的物体传导，因而 。所以温度差 恒正；又因物体将随时间而逐渐冷却,故温度变化速度恒负。故有： (7.1.1) 这里 是比例常数。方程（7.1.1）就是物体冷却过程的数学模型，它含有未知函数 及它的一阶导数,这样的方程称为一阶微分方程。为了解出物体的温度 和时间 的关系,我们要从方程（7.1.1）中解出u=u(t)。 令 ,则有 （7.1.4） 再根据初始条件 可得 ,于是 （7.1.5） 如果k的数值确定了,(7.1.5)就完全决定了温度 和时间 的关系。 根据条件 时， ，得到 由此得到 从而 （7.1.6） 20分钟后物体的温度就是 在找到了变量之间所要满足的微分方程后,还需要找出代表所考虑的问题的初始状态的条件,这就是所谓的初始条件。求一个微分方程满足一定的初始条件的解的问题，称为微分方程的初值问题。 对于一个微分方程的初值问题，通常要讨论如下三个问题： （1）解的存在性，即初值问题是否有解？ （2）解的唯一性，即初值问题的解是否只有一个？ （3）解的稳定性，即初值问题的解关于初值、参数等的连续性、可导性等等。以上三个问题也叫做微分方程的适定性。 不是所有的微分方程模型都可用例题中的方法求解出来。1638年莱布尼茨提出模型 该问题1886年才被数学家刘维尔证明没有解析解。只能借助于本章要介绍的数值解法。 7.2 初值问题数值解法的推导方式及常用解法 我们先考虑常微分方程初值问题。它的数学模型是：求 ，使之满足 （7.2.1） 其中 是已知常数, 是已知函数,且满足条件: （7.2.2）