5 数值微分.pptVIP

* §5 数值微分 目标 求 向前差商公式 向后差商公式 中心差商公式 例：用中心差商公式求ex在x=1的导数值，h=0.8。 类似的： 两点公式和三点公式的比较图 a-h a a+h 从图象中可以看出，当h足够小时，上述公式可以很好的近似一阶导数值 但h足够小时，分子部分相近的两数相减，舍入误差增大 希望：收敛速度快，精度较高的算法 加速算法是在对上述公式误差进行分析的基础上提出的 ? 泰勒公式 ∈ [ ， ] 其中 f’(a)近似值 f’(a)近似值 类似 截断误差 截断误差 两式相减 f’(a)近似值 ① ② 类似 截断误差 ? 加速公式 ① 4②- ①得： ② 例：用中点加速算法求ex在x=1的导数值，h=0.8起算。 H 0.8 0.4 0.2 0.1 G(h) 3.01765 2.79135 2.73644 2.72281 G1(h) 2.715917 2.718137 2.719267 G2(h) 2.718285 2.718276 G3(h) 2.71828 ≈ ? 插值型求导公式 !!! 类似的可以求出节点处高阶导数 两点公式 三点公式 向前向后差商 中心差商 输入：积分区间，误差限 输出：序列Tn，Sn,Cn,Rn及积分结果 实验二 数值积分 1 2 设计方案计算国土面积（选做） 用龙贝格算法计算 参考书本P71的表2-5 为了计算瑞士国土的面积,首先对地图作了如下测量:以西 向东方向为x轴,由南向北方向为y轴,选择方便的原点,并将 从最西边界到最东边界在x轴上的区间适当地划分为若干 段,在每个分点的y方向测出南边界点和北边界点的y坐标, 数据如表(单位mm): x 7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 y1 44 45 47 50 50 38 30 30 34 y2 44 59 70 72 93 100 110 110 110 x 61.0 68.5 76.5 80.5 91.0 96.0 101.0 104.0 106.5 y1 36 34 41 45 46 43 37 33 28 y2 117 118 116 118 118 121 124 121 121 x 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0 y1 32 65 55 54 52 50 66 66 68 y2 121 122 116 83 81 82 86 85 68 瑞士地图的外形如图(比例尺18mm:40km) 试由测量数据计算瑞士国土的近似面积,并与其精确值 41288平方公里比较 理解插值型求积公式和代数精确度的概念 熟悉两点梯型公式及三点Simpson公式； 熟悉复化梯形公式和复化Simpson公式及截断误差 熟悉梯形递推公式会编程实现Romberg算法 理解高斯公式基本概念，了解两点，三点高斯公式的计算 了解简单的数值微分的方法，如中点法 本章要点 P94. 2. 3(3). 8. 10. 12. 13. 17. 26. 本章作业 * *