计算方法(六).pptVIP

* * 插值函数的应用 第 5 章 插值函数的应用 插值方法是一种重要的函数逼近方法，它在 数值微积分和常微分方程数值解中有重要应用． 由 Newton-Leibniz公式，连续函数 在 上的定积分 其中 是 的原函数。 5.1.1 数值求积公式及其代数精度 无能为力。 不能用初等函数表示，即 找不到的原函数； ， ， 没有解析表达式，用表格方式给出时； 大多数的无穷积分，除特殊的无穷积分外。 N-L公式已经 但是大多数实际问题， 常常遇到的困难是： 虽然找到 的原函数，但是它比被积函数复杂的多， 上述的积分就只能利用数值积分公式进行近似计算。 （5-1） 设 是定义在 上的可积函数，考虑带权积分 在 上非负可积，且至多有有限个零点。 其中权函数 所谓数值求积就是用 本节只讨论 的情形。 近似计算 的值。 ． （5-2） 数值求积公式 公式（5-2）称为数值求积公式， 是与 无关的常数，称为求积系数， 其中 上的点 称为求积节点。 求 积 系 数 求积节点 大家熟知第一积分中值定理： 但是的具体位置不可确定。 其几何意义为： 数值积分公式产生的背景 矩形 的面积= 曲边梯形 的面积。 我们可以采用不同的的近似值的方法得到下述数值求积公式： 称为左矩形数值求积公式； 称为右矩形数值求积公式； 称为中矩形数值求积公式； 称为梯形数值求积公式。 （称为步长），将分点 取为插值节点（也是求积节点）， 得到的数值求积公式称为插值型求积公式。 本节采用的逼近函数是 在等距节点上的插值多项式， 进行 等分，令 将 则 可表示为它的Lagrange插值多项式及其余项之和，即 （5-3） 所以 称为 点的Newton-Cotes公式，其中求积系 这样得到的插值型求积公式 (5-6) （5-4） (5-5) 求积余项 (5-7) 标志着求积公式的误差大小。 时的三个公式， 在Newton-Cotes公式中，最常用的是 (5-8) 此时 这就是梯形求积公式： 即 梯形求积公式 此时 这称为Simpson求积公式： (5-9) 进一步可得 Cotes公式 (5-10) Simpson求积公式 Cotes求积公式 练习题 用梯形求积公式和Simpson求积公式计算积分 解： 由梯形求积公式： 由Simpson求积公式： 练习题 用梯形求积公式和Simpson求积公式计算积分 解： 由梯形求积公式： 由Simpson求积公式： 如果某个数值求积公式对比较多的函数能够准确成立，即 那么这个公式的使用价值就较大，可以说这个公式 的精度较高．为衡量数值求积公式的精度，引进代数精度的概念。 如果某个数值求积公式，对于任何次数不超过 次的代数多项式都是精确成立的 但对于 次代数多项式不一定能准确成立，即 则称该求积公式具有 次代数精度． 定义5.1 次代数精度的充要条件是它对 显然，一个数值求积公式具有 这是确定代数精度的最常用方法。 都能准确成立，但对 不能准确成立。 下面求梯形数值求积公式和Simpson数值求积公式的代数精度。 ，我们可得 对于 故梯形数值求积公式具有1次代数精度。 , 我们可得 对于 故Simposon数值求积公式具有3次代数精度。 而 一般的n+1点Newton-Cotes公式的求积余项，有如下定理： 当然也可以通过求积余项估计, 得到代数精度．以下先推导 几个求积余项，进而指出n+1点Newton-Cotes公式的代数精度。 定理5.1 若 其中 是奇数，且 ;若 ，则 其中 ． 定理5.1 是偶数，且 ，则 当 为偶数时， 由于对 次多项式 所以由上述定理可知， 点的Newton-Cotes公式的代数精度为 梯形公式、Simpson公式及Cotes公式的代数精度分别为1，3，5. 当 为奇数时， 点的Newton-Cotes公式