计算方法-4.1数值微分.pptVIP

华长生制作 * * 4.1 数值微分 * 微积分中，函数的导数是通过差商的极限定义，表示函数在某点的瞬时变化率，即平均变化率的极限。其几何意义为曲线的斜率。 * 数值计算中，函数以离散点值给出或函数的表达式非常复杂时，函数的导数用差商近似表示。 向前差商 近似取 余项为 * 4.1.1 数值微分 * 向后差商 中心差商 * 类似用泰勒展开 * 步长h越小，精度越高 ？ * 数值微分的精度与步长 7.3891 3.0042 2.7456 2.7210 2.7183 2.7156 2.6912 2.4596 1 2 1.100 1.010 1.001 1 0.999 0.990 0.900 0 例： * * 0.4763 0.0047 0.0017 0.0183 3.1946 2.7230 2.7200 2.7000 0.9525 0.1407 0.0117 0.0183 4.6708 2.8590 2.7300 2.7000 1 0.1 0.01 0.001 步长h由1减少到0.01时，两种方法计算的误差都显著减小，当h由0.01变为0.001时，误差反而增大。 用函数的差商近似导数，会产生截断误差，且截断误差与|h|（步长）有关，当|h|越小，则截断误差也会小；但是截断误差只是误差的一部分，此外还有舍入误差，而数值微分恰好对舍入误差非常敏感，它有时会随着|h|的缩小而增大。 * 设定最佳步长？ 方法：通常用事后估计方法选取步长h，若记D(h), D(h/2)分别为步长h, h/2 的差商计算公式，给定误差限ε，当| D(h) -D(h/2) | ε时，h/2就是合适的步长。 或以| D(hi) -D(hj) |最小确定步长。 * * 4.1.2 插值型求导公式 --------(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不要取得太小 * * 4.1.3 低阶插值型求导公式的分段构造 由于高次插值易产生Runge现象,数值微分一般采用分段低次插值公式,常见的就是分段两点、三点和五点公式 1.分段两点求导公式 对于任取的相邻两点 由两点公式有 --------(10) 称(10)式为分段两点公式 * * * * 2.分段三点求导公式 对于任取的相邻三点 由三点公式,有 --------(11) 称(11)式为分段三点公式 --------(12) 3.分段五点求导公式 由五点公式 * * 看一个实例: 已知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世纪的(相对)年增长率 * * 解: 先用精度较高的分段五点公式求出节点处的导数值 则 求增长率必须先求导数 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 * * * 将节点处的增长率作 三次