数值分析4.7数值微分资料.pptVIP

华长生制作 下面再看一个用泰勒展开得二阶数值导数的例子。 * 4.7 数值微分 先看一个实例: 已知20世纪美国人口的统计数据为(单位:百万) 年份 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 人口 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 试计算美国20世纪的(相对)年增长率 一、插值型求导公式 --------(1) 对(1)式两边求导,有 --------(2) --------(2) --------(3) (2)式称为插值型求导公式, (3)式为相应产生的误差 由于公式(2)采取的是n次Lagrange插值多项式,而高次 插值会产生Runge现象,因此实际应用中多采用低次插 值型求导公式 二、低阶插值型求导公式 1.两点公式 --------(4) --------(5) (4)(5)式称为带余项的两点求导公式 即 精度1阶 2.三点公式 --------(6) --------(7) --------(8) (6)(7)(8)式称为带余项的三点求导公式 其中(7)式又称为中点公式,其精度稍高 在分段求导公式中有着重要的地位 精度2阶 3.五点公式 --------(9) 组(9)称为带余项的五点求导公式 精度4阶 综合考虑上述三种公式,可知五点公式的精度最高 可发现 且 并且当步长h越小时,误差会越小 但是不是h越小公式越好呢? 因此实际应用中步长h不要取得太小 两点公式和三点公式的比较图 例 设 ，对h=0.01，计算 的近似值。 解 由一阶三点公式有 由五点求导公式有 精确值 。计算结果显然与它们的余项相一致，由五点求导公式 计算所得的结果最精确。 三、低阶插值型求导公式的分段构造 由于高次插值的Runge现象,数值微分一般采用分段低 次插值公式,常见的就是分段两点、三点和五点公式 1.分段两点求导公式 对于任取的相邻两点 由两点公式有 --------(10) 称(10)式为分段两点公式 2.分段三点求导公式 对于任取的相邻三点 由三点公式,有 --------(11) 称(11)式为分段三点公式 --------(12) 3.分段五点求导公式 由五点公式 例: 回到实例(美国人口) 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 解: 则 先用精度较高的分段五点公式求出节点处的导数值 求增长率必须先求导数 x=2.1508 1.3458 1.6158 1.1908 1.2267 2.5000 2.7633 2.3075 2.2692 2.8358 r=dx/dt/x =0.0283 0.0146 0.0152 0.0097 0.0093 0.0166 0.0154 0.0113 0.0100 0.0113 数值微分就是用离散方法即使的近似地求出函数在某点的导数 值.按照Taylor展开原理可得 其中h为一增量，对于一阶导数这几个公式是很实用的。 4.7.2泰勒展开法 * * * * *