第四章数值微分与数值积分介绍.pptVIP

第四章 数值微分与数值积分 　　数值微分是指用数值方法求解表达式复 杂的函数（或表格函数）微分（或导数）的 数值方法。 差商型微分计算公式 差商型微分公式的几何含义 向后差商 与计算向前差商的方法类似，由泰勒展开得向后差商的余项： 中心差商 插值型微分计算公式 利用数据表构造函数f x 的插值多项式Pn x ，并取Pn ‘ x 的值作为f ‘ x 的近似值，这样建立的数值公式f ‘ x ≈Pn’ x 称为插值型微分公式。 若取拉氏插值式，则得 从而有 若求节点xi处的导数，因为 ，所以 即以插值型微分公式估计导数时对应的截 断误差 下面是取等距节点 xi x0+ih,i 0,1,…,n 时常用公式 1.两点公式（n 1） 第二节 数值积分 　　在微积分中用Newton-Leibniz公式计算连续函 数的定积分： 插值型求积公式 解 由 得 最后得到数值积分公式： 下面讨论这种求积公式的误差： 　　这种插值积分的误差就是对插值误差的 积分，我们利用拉式插值的余项公式，得 求积公式的代数精度 　　由此可知当求积公式具有m阶代数精度时，对任意的m阶多项式都有 显然，如果一个求积公式的代数精度越 高，它就能对更多的被积函数f x 准确成 立，具有更好的实际计算意义。 例：确定积分公式中的A0，A1，A2 使积分公式的代数精度最高，并指 出实际的代数精度。 继续计算 第三节 牛顿-柯特斯求积公式 把积分区间[a，b]分成n等分，记步长 为h b-a /n，取等分点xi a + ih i 0,1,...,n 作为数值积分节点，A的计算公式在做变量 替换 x a+th后化简为 记 于是，插值型求积公式可写成 下表中列出n从1到6的牛顿-柯特斯系数，从表中可以看出牛顿-柯特斯系数具有对称性。 梯形公式 记 怎样确定梯形公式的代数精度？ 即： 因为 x-a x-b 在[a,b]上不变号，由积分中值定理可知，在上必存在一点η，使得 辛卜生 Simposn 公式 所以 其几何意义即用过三点的抛物线面积近似代替积分的曲边面积 分别将f x 1，x，x2， x3，代入到I f 和S f 中，可以得到辛卜生公式的代数精度为3，也就是它对于次数不超过三次的多项式成立。 柯特斯 Cotes 公式 例 用梯形公式、辛卜生公式及柯特斯公式计算积分 （结果取5位） 解（1）用梯形公式计算 积分的标准值 计算结果表明，柯特斯公式精度最高， 辛卜生公式次之，梯形公式精度较差。 第四节 复化求积方法 把积分区间分割成若干小区间，在每 个小区间[xi,xi+1]上用梯形积分公式，再将 这些小区间上的数值积分累加起来，称为 复化梯形公式。 由梯形积分公式的分析可知，每一个 子区间[xi,xi+1]上，求积的截断误差为 由于 在[a,b]上的连续性，可知在 a,b 内必存在一点 使 将各子区间的公式迭加，我们得到 下表显示了n 8时复化辛卜生公式中节 点和系数的关系：在每个积分区间上提出因 子h/6后，三个节点的系数分别是1，4，1； 将4个区间的系数按节点的位置累加，可以 看到首尾节点的系数是1，奇数点的系数是 4，偶数点是2。 由误差公式 由于函数的二阶导数在区间[a，b]上变化 不大，所以 ，于是 复化梯形的递推算式 当利用复化梯形公式Tn计算需要提高精度时，每次对前一次的子区间分半，可以充分利用老分点上的函数值，每次只需计算新增加的分点的和。 例 用复化梯形公式的递推算式计算 复化积分的算法描述 复化梯形方法： 对[a，b]上n等分， 则有 龙贝格算法计算元素顺序表 例 用龙贝格方法计算积分 解 按上述步骤计算 由此得到启发，对任给的误差控制量ε 0，要 ，只需 即可。 对区间[a，b] n等分，h b-a /n，记[xi，xi+1/2] 上的中点为xi+1/2，可以得到Tn和T2n的表达式如下： 所以 用这种递推公式计算时，只需要计算新增加的节点上的函数值，原有 的计算结果仍然可以使用。 解（1）在区间[0，1]上使用梯形公式，得 （2）将区间[0，1]二等分，x 1/2是新分点 （3）再将各小区间二等分，x 1/4和x 3/4 是新分点 1.令 2.对 计算 3. 复化辛卜生方法： 1.令 2.对 计算 3. 第五节 龙贝格 Romberg 方法 梯形公式计算方法简单，但是收敛速度慢，为了加快其收敛速度，我们引入了龙贝格方法。 在前面关于梯形公式截断误差的事后估计的讨论中，我们得到 即 再将Tn和T2n的表达式代入，得 　　在实际应用中，将积分区间等分以后构造的插值型积分公式，称为牛顿-柯特斯求积公式。 下面可以看到，牛顿－柯特斯积分系数和积分节点以及积分区间无直接关系，系数固定而易于