常微分方程的数值解法-多库文档.PDF

第四章 常微分方程的数值解法 Numerical Solutions to Ordinary Differential Equations 概述 一阶常微分方程初值问题: d y = f x , y 求函数y=y(x)满足： P r o b l e m d x ( ) y (x) =f(x, y(x)) y (x 0 ) = y 0 f(x,y)在[a,b]上连续,且满足Lipschitz条件: L,y ,y , s.t. |f(x,y )-f(x,y )| L|y -y | 1 2 1 2 1 2 则初值问题Problem I有唯一解 实际工程技术、生产、科研上会出现大量的微分方程问题 很难得到其解析解，有的甚至无法用解析表达式来表示， 因此只能依赖于数值方法去获得微分方程的数值解。 微分方程的数值解法： p不求y=y(x)的精确表达式,而求离散点x ,x , x 处的函数值 0 1 n p设Problem I 的解y(x)的存在区间是[a,b]，初始点x =a，取 0 [a,b] 内的一系列节点x , x ,...,x 。a= x x x =b 0 1 n 0 1 n 一般 用等距步长。 a b x x x ... x x 0 1 2 n-1 n p用数值方法，求得y(x)在每个节点x 的值y(x ) 的近似 k k 值，用y 表示，即y ≈y(x )，这样y , y ,...,y 称为微分 k k k 0 1 n 方程的数值解。 p求y(x)————求y , y ,...,y 0 1 n @方法： 用步进式和递推法 将[a,b]n等分 a= x x x =b 步长h= b-a x =a+kh 0 1 n n k y = y(x ) 0 0 y = g(h, x y , y , y , ...y ) n+1 n, n n-1 n-2 n-m 计算过程： y0 y1 y2 ynm ynm1 yn yn1 怎样建立递推公式？ Taylor公式 数值积分法 4.1 欧拉公式 @思想:用向前差商近似代替微商.