《地球物理计算方法》第2章 数值积分—数值微分.pptVIP

《地球物理计算方法》

第2章数值积分—数值微分*1、差商公式h-0的极限，数值微分：用函数值的线性组合近似函数在某点的导数值。差商比较示意图―――――中点方法*中点格式表达成为误差R误差R正比与h2。从截断误差角度看，h越小越好；从舍入误差角度看，h不能太小（f(a+h)与f(a-h)太接近）；所以实际计算时h选择一个合适值。*2、中点方法的加速中点格式表达成为误差R分析步长减半为h/2时*分析步长为h/2时误差R分析步长为h时误差R*分析步长为h和为h/2时微分值的来估计误差故有其误差事后估计得到：*如果忽略高次项：整理得到（一次校正加速公式）：可以逐步二分h，逐步加速*例：求函数f(x)=ex在x=1处的导数值。（近似值要求三位有效数字）解：中点公式：hG(h)G1(h)G2(h)G3(h)0.83.017652.7159172.7182852.718280.42.791352.7181372.7182760.22.736442.7192670.12.722813、数值求导公式的设计方法*令它对于成立等式，得如1、代数精度来构造**2、插值多项式方法来构造*插值多项式余项依据插值多项式余项定理求导得到：*对于任意的x，其项