《地球物理计算方法》第三章 微分方程差分法——Euler及其改进方法.pptVIP

《地球物理计算方法》

第3章微分方程差分法——Euler及其改进方法*问题的提出常微分方程初值问题：*(1)只有少数问题如能求出y=f(x)的解析表达式，对大多数微分方程要求出函数的准确表达式，计算量大、甚至不可能；(2)不需要完全得到y=f(x)，有时只需得到在某些点处的函数近似值即可。问题：求解函数yk=f(xk)*初值问题差分方法的特点：

步进式——求解过程顺着节点排列的次序一步一步地向前推进。只要给出从已知信息yn，yn-1，yn-2，…计算yn+1的递推公式。求解的核心：设计差分格式。消除导数，离散化1、Euler方法*Euler格式：微分的离散化——差商代替导数在点xn列出一阶方程

将xk处的导数y’(xk)近似地用差商表示（向前差分格式） *几何解释（又称为欧拉折线法）ynxn=boy0yxa=x0x1x2x3y3y1y2*按泰勒公式的展开：局部截断误差：欧拉公式局部截断误差:*若局部截断误差为O(hp+1)，则称该数值方法的精度为p阶的。精度：欧拉格式的局部截断误差:1阶精度*2、隐式Euler方法隐式Euler方法：向后差商公式。 xn xn+1特点：非显式格式，计算比较困