数学建模算法第15章 常微分方程的解法.docxVIP

- - - - 第十五章 常微分方程的解法 建立微分方程只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并加以检验。如果能得到解析形式的解固然是便于分析和应用的，但是我们知道，只有线性常系数微分方程，并且自由项是某些特殊类型的函数时，才可以肯定得到这样的解， 而绝大多数变系数方程、非线性方程都是所谓“解不出来”的，即使看起来非常简单的 方程如 dy ? y 2 ? x2 ，于是对于用微分方程解决实际问题来说，数值解法就是一个十 dx 分重要的手段。 §1 常微分方程的离散化 下面主要讨论一阶常微分方程的初值问题，其一般形式是 ? dx? dy ? f ( x, y ? dx ? ?? y(a) ? y0 a ? x ? b  (1) 在下面的讨论中，我们总假定函数 f ( x, y) 连续，且关于 y 满足李普希兹(Lipschitz)条件，即存在常数 L ，使得 | f ( x, y) ? f ( x, y) |? L | y ? y | 这样，由常微分方程理论知，初值问题(1)的解必定存在唯一。所谓数值解法，就是求问题（1）的解 y( x) 在若干点 a ? x0 ? x1 ? x2 ? L? xN ? b 处的近似值 yn (n ? 1,2,L, N ) 的方法， yn (n ? 1,2,L, N ) 称为问题（1）的数值解， hn ? xn ?1 ? x