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

《地球物理计算方法》

第2章数值积分—数值微分

21、差商公式h-0的极限，数值微分：用函数值的线性组合近似函数在某点的导数值。

差商比较示意图―――――中点方法

4中点格式表达成为误差R误差R正比与h2。从截断误差角度看，h越小越好；从舍入误差角度看，h不能太小（f(a+h)与f(a-h)太接近）；所以实际计算时h选择一个合适值。

52、中点方法的加速中点格式表达成为误差R分析步长减半为h/2时

6分析步长为h/2时误差R分析步长为h时误差R

7分析步长为h和为h/2时微分值的来估计误差故有其误差事后估计得到：

8如果忽略高次项：整理得到（一次校正加速公式）：可以逐步二分h，逐步加速

9例：求函数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.72281

3、数值求导公式的设计方法??

11令它对于成立等式，得如1、代数精度来构造

12

132、插值多项式方法来构造

14插值多项式余项依据插值多项式余项定理求导得到：

15对于任意的x，其项难以确定：但对于节点xk上：该节点xk的余项：

(1)两点公式设已给出两个节点上的函数