(简缩版)10计算方法第十节常微分方程初值问题的数值解法(2011-2012学年).pdf

第十章 常微分方程初值问题的数值解法 10.0 引言 在工程和科学计算中，所建立的各种常 微分方程的初值或边值问题，除很少几类 的特殊方程能给出解析解，绝大多数的方 程是很难甚至不可能给出解析解的，其主 要原因在于积分工具的局限性。因此，人 们转向用数值方法去解常微分方程，并获 得了相当大的成功，讨论和研究常微分方 程的数值解法是有重要意义的。 计算方法 一阶常微分方程初值问题的一般形式是: ⎧y ′ f (x , y ),a =x b ⎨ (1) ⎩y (x 0 ) y 0 D {(x , y ) a =≤x ≤b ,c ≤ y ≤d } 称 f (x ,y ) 在区域 D 上对 y 满足 Lipschitz 条件是指: ∃L 0s .t . ∀x ∈[a ,b ], y , y ∈[c ,d ] 1 2 f (x ,y 1 ) −f (x ,y 2 ) ≤L y 1 −y 2 计算方法 已经知道：若f (x ,y ) 在区域 D 上连续,且对 y 满足 Lipschitz 条件,则初值问题(1)在[a,b]上存在唯一的连续 可微解. 所谓微分方程的数值解法，就是寻求方程（1）的解y (x ) 在一系列离散节点 x x x x 1 2 n n+1 上的近似值 y , y , , y , y , 1 2 n n+1 相邻两个节点之间的距离h xn+1 =−xn 称为步长。今后总 设步长为定数。这样就有xn x0 =+nh 计算方法 10.1 初值问题的Euler方法 考虑常微分方程初值问题 ⎧y ′ f (x , y ) ⎨ ⎩y (x 0 ) y 0 计算方法 这样，就得到著名的 Euler 在 x 点列出方程 n 格式： ′ y (x n ) f (x n , y n ) 并令 ⎧y n+1 y n =+hf (xn , y n ) ⎨ ( ) ′ y (x n+1 ) −y (x n ) y (x n +h ) −y (x n ) ⎩y x 0 y 0 y (x n ) h h 代入方程就得到 y (x n+1 ) −y (x n ) ≈f (x n ,y (x n )) h 即