数值微分与数值积分.pptVIP

第七章　数值微分与数值积分 复合梯形公式： 即分段用Simpson公式再求和可得复化Simpson公式. 2. 复化Simpson公式 x0 x2 x f(x) x4 h h xn-2 h xn …... h x3 x1 xn-1 4 4 4 4 4 复化Simpson公式： 例：若用复化Simpson公式计算积分 问积分区间要等分为多少份才能保证计算结果 有4位有效数字？ 如果用复化梯形公式计算，则由误差公式 注意到区间再次对分时 注意到区间再次对分时 3 .3 逐次分半算法 递推公式 （一）梯形公式的逐次分半算法 （二）Simpson公式的逐次分半算法 为Simpson公式的逐次分半算法的事后误差估计公式， 实际计算时常以 为停步准则。 梯形公式算法7.1 Simpson公式算法7.2 §4　Romberg求积公式 4.1　Richardson外推法 特别地，若 4.2　Romberg求积公式 由Euler－Maclaurin公式(《数值逼近》黄友谦，李岳生) Romberg方法计算过程 结论： §5　Gauss型求积公式 （一） Gauss型求积公式 解：记 * * 1.1 差商型求导公式 §1 数值微分 向前差商公式 向后差商公式 中心差商公式 1.2 插值型求导公式 适用于求节点处的导数值， 以插值多项式的导数 作为函数导数的近似，即 常用公式 　　 1.两点公式　 2.三点公式　　 误差分别为 i-2 i-1 i i+1 i+2 抛物线 1.3 利用样条插值函数求数值微分 且有误差估计 数值微分公式 §2 Newton－Cotes求积公式 2.1 数值积分的基本思想： 用插值多项式的积分近似函数的积分： 插值型求积公式： 求积公式的 截断误差为 2.2 Newton－Cotes公式 此求积公式称为n阶Newton-Cotes公式 称为Cotes系数，它与积分区间、被积函数均无关， 只与积分区间的等分数n有关. n阶N-C公式的截断误差为 x f(x) L(x) y x0 x1 x f(x) x2 h h L(x) x0 x1 x f(x) x2 h h L(x) x3 h 则称该求积公式具有m次代数精确度。 2.3 误差估计 (一) 求积公式的代数精确度 结论： ∴　抛物线公式具有3次代数精确度 （二）梯形公式与Simpson公式的误差估计 则　 §3　复化求积公式 基本思想：高次插值有Runge 现象， 用分段低次插值函数近似被积函数求积分近似值。 1.复化梯形公式 x0 x1 x f(x) x2 h h x3 h h x4 * *