数值分析第七章.pptVIP

三峡大学理学院 杜廷松 第7章 数值微积分 §7.1 数值积分问题的提出 计算定积分，但在工程技术和科学研究中，常遇到如下情况： 1) 找 的原函数 相当困难 例如： 等，它们的原函数 都不能用初等函数表达成有限形式. 2) 的 原函数存在，但其表达式相当复杂，而且有时难以给出最后的数值结果.如在积分表中可以查到不定积分的解析式 当积分限 时，原函数中第一项无定义，得不到定积分值； 3) 除一些特殊的无穷积分外，通常很难求无穷积分的值; 4) 被积函数没有有限的解析式，而以表格形式给出. 鉴于以上情况，采用数值的方法计算定积分就很必要且很有效.对函数的微分也一样，以表格形式给出的函数，要求出其导数时，还是要依靠数值微分的方法. 本章介绍用插值原理得到的常用的数值求积公式、它们的误差估计、代数精度.之后介绍数值微分. §7.2 插值型求积公式 一 、插值型求积公式 二 、求积公式的代数精度 一 、插值型求积公式 利用插值多项式来构造数值求积公式，具体作法如下： 在积分区间 上取一组节点 (已知 的函数值，作 的 次插值多项式，则 其中 为n次插值基函数，则得： 公式写成一般形式： (2.1) 若记： (2.2) 则得数值求积公式： (2.4) 形如(2.4)式的求积公式称为机械求积公式，其中 称为求积节点， 称为求积系数.如果公式(2.4)中求积系数 是由(2.2)确定的，则称该求积公式为插值型求积公式. 注意：求积系数 仅与求积节点、求积区间有关而与被积函数 的具体形式无关. 二、 求积公式的代数精度 定义2.1 求积公式（2.4）对于任何次数不大于 的代数多项式 均精确成立，而 对 不精确成立，则称该求积公式具有 次代数精度. 例2.1 梯形公式 代数精度 定理7.2.1 含有 个节点 的插值型求积公式的代数精度至少为 . 例2.2求积公式 已知其余项的表达式为 ，试确定系数 ,使该求积公式具有尽可能高的代数精度，并给出求积公式的余项和代数精度的次数. §7.3 Newton－Cotes公式 一 、Newton－Cotes公式 二 、误差分析 三、 数值稳定性 四 、复合Newton－Cotes公式 一 、Newton－Cotes公式 前面介绍了插值型求积公式及其构造方法，在实际应用时，考虑到计算上的方便，常将积分区间等分之，并取分点为求积节点，这样构造出来的插值型求积公式就称为Newton－Cotes公式.一句话，Newton－Cotes公式是节点等距的插值型求积公式. 将 n等分,其节点为 ,节点 处的函数值为 .则Newton－Cotes公式为： 作变换： ，则插值型求积公式中系数 为 记： (3.1) 则 于是得等距节点的插值型的Newton－Cotes求积式： (3.2) 称为Cotes系数. 当 是次数不超过 次的多项式时， 的插值多项式就是 本身， 柯特斯(Cote