微分方程数值解习题课概要.doc

微分方程数值解习题课纲要 微分方程数值解习题课纲要 PAGE / NUMPAGES 微分方程数值解习题课纲要 微分方程 初值问题数值解 习题课 1 一、应用向前欧拉法和改良欧拉法求由以下积分 x e t 2 y dt 0 所确立的函数 y 在点 x =0.5,1.0,1.5 的近似值。 解：该积分问题等价于常微分方程初值问题 y e x2 y(0) 0 此中 h=0.5 。其向前欧拉格式为 y y he ( ih ) 2 i 1 i y0 0 2 改良欧拉格式为 yi 1 yi h (e ( ih)2 e ( i 1)2 h2 ) 2 y0 0 将两种计算格式所得结果列于下表 i xi 向前欧拉法 yi 改良欧拉法 yi 0 0 0 0 1 0.5 0.5 0.44470 2 1.0 0.88940 0.73137 3 1.5 1.07334 0.84969 3 二、应用 4 阶 4 步阿达姆斯显格式求解初值问题 y x y 1 0 x 0.6 y(0) 1 取步长 h=0.1. 解： 4 步显式法一定有 4 个起步值， y0 已知，其余 个 y1, y2 , y3 用 4 阶龙格库塔方法求出。 此题的信息有： 步长 h=0.1 ；结点 xi ih 0.1i (i 0,1, ,6) ； f ( x, y) x y 1, y0 y(0) 1 4 经典的 4 阶龙格库塔公式为 5 y y i h ( k 2k 2 2k 3 k 4 ) i 1 6 1 k1 f (xi , yi ) xi yi 1 k 2 f (x h , y hk1 ) x y 0.05k 1.05 i 2 i 2 i i 1 k 3 f ( x h , y i hk2 ) x i y 0.05k 2 1.05 i 2 2 i k4 f ( xi h, yi hk3 ) xi yi 0.1k3 1.1 算得 y1 1.0048375, y2 1.0187309, y3 1.0408184 6 4 阶 4 步阿达姆斯显格式 y y h (55 f i 59 f i 1 37 f i 2 9 f i 3 ) i 1 i 24 yi 1 5.9 yi 1 3.7 yi 2 0.9 yi 3 0.24i 3.24) (18.5 yi 24 由此算出 y4 1.0703231, y5 1.1065356,y6 1.1488186 7 三、用 Euler 方法求 y ex y x 1, 0 x 1 y 0 1 问步长 h 应当怎样选用，才能保证算法的稳固性？ 解：此题 f x, y ex y x 1 f y x, y ex 0,0 x 1 此题的绝对稳固域为 1 h 1 hex 1 得 0 hex 2 ， 故步长应知足 0 he 2, 0 h 0.736 8 四、 求梯形方法 yk 1 yk h [ f ( xk , yk ) f (xk 1 , yk 1 )] 2 的绝对稳固域。 证明：将 Euler 公式用于试验方程 y y ，获得 yk 1 yk h [ ykyk 1 ] 2 整理 9 1 h yk 1 (1 h) yk 2 2 设计算 yk 时有舍入偏差 k , k 0,1,2, ，则有 1 h k 1 (1 h ) k 2 2 据稳固性定义，要想 k 1 k ，只须 1 h 1 h 2 2 10 所以方法绝对稳固域为复平面 h 的整个左半平面 （？），是 A-稳固的。 11 五、对初值问题 y y y(0) 1 0 x 1 证明：用梯形公式 y n 1 y n h [ f ( x n , y n ) f ( x , y )] 2 n 1 n 1 求得的数值解为 2 n yn h 2 h 并证明当步长 h 0 时， yn 收敛于该初值问题的精准 解 yn e x 12 证明：由梯形公式，有 y y n h [ f ( x , y ) f ( x n 1 , y n 1 )] y n h [ y n y n 1 ] n 1 2 nn 2 整理，得 y n 1 2 h y 2 h n 由此递推公式和初值条件，有 13 yn 2 h n y0 2 h n 2 h 2 h 14 x [0,1] ，则有在区间 0, x 0,1 上有 x xn nh ，步长 h x n ，由前面结果有 lim yn lim 2 h n lim 1 2h 2 h 2 h n n h 0 2 h 2x lim 1 2h 2 h 2h e x h 0 2 h x h 由 x 的随意性，得所证。 15 六、关于微分方程 y f ( x, y)