2020第16节欧拉方法.ppt

改进型欧拉方法 梯形公式虽然提高了精度，但仍是隐式方法，算法复杂， 运算量大。而在实际计算中只迭代一次，这样建立的预测 — 校正系统称作改进的欧拉公式。先用欧拉折线法得到初步近 似值，再用梯形公式校正。 预测 y eu ( n ? 1) ? y ( n ) ? hf ( x ( n ) , y ( n ) ); h eu 校正 y ( n ? 1) ? y ( n ) ? [ f ( x ( n ) , y ( n ) ) ＋ f ( x ( n ? 1) , y ( n ? 1) )], 2 eu ? y ( n ? 1) ? y ( n ) ? hf ( x ( n ) , y ( n ) ) ? eu 也可写为： ? y c ? y ( n ) ? hf ( x ( n ? 1) , y ( n ? 1) ) ? y ? ( y ? y ) / 2 ( n ? 1) p c ? 第六章 微分方程的初值问题 第一节 欧拉方法 微分方程（组） 微分方程的描述： 1 、自变量： x ， t ， .... 2 、因变量： y ， z ， u ， v 3 、常微分方程： F ( y , y ( n ) ( n -1) ,..., y '', y ', y , x ) ? 0 ( n ) 微分方程组： ? F 1 ( y , y ,..., y '', y ', y , z ,..., z , x ) ? 0 ? ( n ) ( n -1) ( n ) ? F 2 ( y , y ,..., y '', y ', y , z ,..., z , x ) ? 0 ? ... ? ( n ) ( n -1) ( n ) ? ? F m ( y , y ,..., y '', y ', y , z ,..., z , x ) ? 0 ( n ) ( n -1) 线性微分方程的通解 线性微分方程： a n ( x ) y ? ,..., ? a 2 ( x ) y '' ? a 1 ( x ) y ' ? a 0 ( x ) y ? 0,( 齐次方程 ) ( n ) a n ( x ) y ? ,..., ? a 2 ( x ) y '' ? a 1 ( x ) y ' ? a 0 ( x ) y ? f ( x ),( 非齐次方程 ) ( n ) 齐次线性微分方程的通解： Y ( x ) ? c 1 y 1 ( x ) ? c 2 y 2 ( x ) ? ... ? c n y n ( x ) 非齐次线性微分方程的通解： Y ( x ) ? c 1 y 1 ( x ) ? c 2 y 2 ( x ) ? ... ? c n y n ( x ) ? y p ( x ) 常系数线性微分方程的通解 常系数线性微分方程解法： a n y ? ,..., ? a 2 y '' ? a 1 y ' ? a 0 y ? f ( x ) ( n ) 拉普拉斯变换 ： L [ a n y ? ,..., ? a 2 y '' ? a 1 y ' ? a 0 y ] ? L [ f ( x )] ? [ a n s ? ,..., ? a 2 s ? a 1 s ? a 0 ] Y ( s ) ? F ( s ) F ( s ) Y ( s ) ? n 2 1 a n s ? ,..., ? a 2 s ? a 1 s ? a 0 则： y ( t ) ? ? Y ( s ) ds n 2 1 ( n ) 微分方程的初值问题与边值问题 微分方程初值问题： ? ? y ' ? f ( y , x ) ? y ( x 0 ) ? b ? ? y '' ? f ( y ' y , x ) ? y ( x 0 ) ? a ? ? y '( x 0 ) ? b 微分方程边值问题： ? ? y '' ? f ( y ' y , x ) ? y ( x 0 ) ? a ? ? y ( x 1 ) ? b 一阶微分方程的解 李普希兹 ( Lipschitz ) 条件 : 如果定义在区间 D ? f ( t , y 1 ) - f ( t , y 2 ) ? L y 1 - y 2 则称该函数满足李普希兹 ( Lipschitz ) 条件， L 为 Lipschitz 常数。 2 的函数