高建军—第三章-常微分方程差分法1-2016.10.pdfVIP

《地球物理计算方法》 第3章 常微分方程的数值求法 1 第3章 常微分方程的数值求法 3.1 欧拉（Euler ）方法 3.2 改进的欧拉方法 3.3 龙格-库塔（Runge-Kutta ）方法 3.4 亚当姆斯（Adams ）方法 3.5 收敛性与稳定性 3.6 例题选讲 2 问题的提出 常微分方程初值问题： 在许多实际问题中，往往不能直接给出函数的解析表达式 y =f (x) ，但根据已知条件，有时可列出含有待求函数及其 导数的关系式——微分方程。 如：常微分方程及边界条件： y (x ) f (x , y )   , (a ≤ x ≤ b) y x y 0  0 求解函数y=f (x) 3 问题的提出 例如： 一曲线通过点(1 , 2)，且在该曲线上任一点M(x , y )处 的切线斜率为2x ，求该曲线的方程。 设所求曲线为y =y (x) ，由已知条件可得 dy  2x  dx y |x 1 2  （注：1阶方程只需要1个初始条件） 4 问题的提出 (1) 求出一般解： c为积分常数 (2) 已知条件定解： 2 由已知条件：y | =2 → 2=1 +c → c=1 x=1 故函数解析表达式为：y =x2+1 5 问题的提出 问题： (1) 只有少数问题如能求出y =f (x) 的解析表达式，对大多数 微分方程要求出函数的准确表达式，计算量大、甚至不可 能； (2) 不需要完全得到y=f (x) ，有时只需得到在某些点处的函 数近似值即可。 6 问题的提出 办法：数值解法－差分法 利用给定常微分方程及边界条件数值方法解出函数y =f (x)在若 干离散点处近似值的方法， 即在区间[a,b]上有若干离散点： a=x x ...x =b， 0 1 n 用离散化方法求出y ，作为精确值y (x )的近似值： k k y 1，y 2 ，…，y n ， 一般：h =x -x 称为步长。 i+1 i 7 问题的提出 初值问题差分方法的特点： 步进式——求解过程顺着节点排列的次序一步一步地向前 推进