《精》第四章 数值微分与数值积分.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。 例 求 ，用复化梯形和复化辛卜生求积公式的分点应取多少？（计算要求5位有效数字） 由复化梯形误差公式得到： 由复化辛卜生截断误差公式，有 由误差公式 由于函数的二阶导数在区间[a，b]上变化 不大，所以 ，于是 复化梯形的递推算式 当利用复化梯形公式Tn计算需要提高精度时，每次对前一次的子区间分半，可以充分利用老分点上的函数值，每次只需计算新增加的分点的和。 例 用复化梯形公式的递推算式计算 复化积分的算法描述 复化梯形方法： 对[a，b]上n等分， 则有 龙贝格算法计算元素顺序表 例 用龙贝格方法计算积分 解 按上述步骤计算 上式表明Tn( f )的误差大约是T2n( f )误差的4倍 由此得到启发，对任给的误差控制量ε0，要 ?????????????????? ，只需 即可。 对区间[a，b] n等分，h=(b-a)/n，记[xi，xi+1/2] 上的中点为xi+1/2，可以