算法大全第15章_常微分方程的解法.pdf

第十五章 常微分方程的解法 建立微分方程只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并 加以检验。如果能得到解析形式的解固然是便于分析和应用的，但是我们知道，只有线 性常系数微分方程，并且自由项是某些特殊类型的函数时，才可以肯定得到这样的解， 而绝大多数变系数方程、非线性方程都是所谓“解不出来”的，即使看起来非常简单的 dy 2 2 方程如 y +x ，于是对于用微分方程解决实际问题来说，数值解法就是一个十 dx 分重要的手段。 §1 常微分方程的离散化 下面主要讨论一阶常微分方程的初值问题，其一般形式是 ⎧dy ≤ ≤ ⎪ f (x ,y ) a x b ⎨dx (1) ⎪ y (a) y 0 ⎩ 在下面的讨论中，我们总假定函数 f (x ,y ) 连续，且关于 y 满足李普希兹(Lipschitz)条 件，即存在常数L ，使得 | f (x ,y ) −f (x ,y ) |≤L | y −y | 这样，由常微分方程理论知，初值问题(1)的解必定存在唯一。 所谓数值解法，就是求问题（1）的解y (x ) 在若干点 a x x x Lx b 0 1 2 N 处的近似值 y n (n 1,2,L,N ) 的方法，y n (n 1,2,L,N ) 称为问题（1）的数值解， hn x n +1 −xn 称为由xn 到xn +1 的步长。今后如无特别说明，我们总取步长为常量h 。 建立数值解法，首先要将微分方程离散化，一般采用以下几种方法： （i ）用差商近似导数 y x y x ( n +1 ) − ( n ) 若用向前差商 代替y (x n ) 代入（1）中的微分方程，则得 h y (xn +1 ) −y (xn ) ≈f (xn , y (xn )) (n 0,1,L) h 化简得 y (x ) ≈y (x ) +hf (x ,y (x )) n+1 n n n 如果用y (x n ) 的近似值y n 代入上式右端，所得结果作为y (x n +1 ) 的近似值，记为y n +1 ， 则有 y n +1 y n +hf (xn ,y n ) (n 0,1,L) （2 ） 这样，问题（1）的近似解可通过求解下述问题 y y +hf (x ,y ) (n 0,1,L) ⎧ n +1 n n n ⎨ （3 ） y 0 y (a) ⎩ 得到，按式（3 ）由初值y 可逐次算出y ,y ,L。式（3 ）是个离散化的问题，称为差 0 1 2 分方程初值问题。