7常微分方程数值解法.pdfVIP

第 7 章 常微分方程数值解法 基本概念 1. 一阶常微分方程的初值问题 y (x) f (x, y(x )), x (a ,b) y (a) y0 注：若 f 在 D = { a x b , | y|+ } 内连续，且满足 Lip 条件： L 0，使 | f ( x – y1) – f ( x ，y2)| L| y 1 – y2| 则的连续可微解 y( x) 在 [ a，b] 上唯一存在。 2. 初值问题的数值解 称的解 y ( x) 在节点 xi 处的近似值 yi y ( xi ) a x1 x2 ... xn = b. 为其数值解，方法称为数值方法。 i 注：① 考虑等距节点 : x = a + ih ， h = ( b – a)/ n. ② 从初始条件 y ( a) = y0 出发，依次逐个计算 y1 ，y2 ，…， yn 的值，称为步进法。 两种：单步法、多步法。 ③ 二阶常微分方程 y( x) = f ( x ，y( x ) ，y( x )) 可设为一阶常微分方程组的初值问题： 引进新的未知函数 z(x) = y(x) ，则 y ( x ) z(x) z ( x ) f ( x, y (x ), z(x )) 其初始条件为： y (a) y0 z(a) y0 称为一阶微分方程组的初值问题，方法类似。 ④ 边界问题，常用差分方法解。 初值问题数值解法的构造及其精度 构造方法 对于可借助 Taylor 展开（导数法） 、差商法、积分法实现离散化来构造求积公式： 1. 设 y C[ a，b] 将 y( xi +1) = y( xi +h) 在 xi 处展开 2 h y(x ) y (x ) hy (x ) y ( ) i 1 i i 2 2 h y( x ) hf (x , y (x )) y ( ) [ xi ，xi +1] i i i 2 ( i ) + ( i ， i ) 其中 i ( i ). y x +1 yi hf x y y y x 称 i = i + ( i ， i ). = 0 ， 1，2，... ， – 1 y +1 y hf x y i n 为 Euler 求解公式， （Euler 法） y (x ) y(x ) 2. 用差商来表示： i 1 i ( , ( )). y xi y