数值积分本章主要内容机械求积的方法梯形公式simpson.DOCVIP

第二章 数值积分 本章主要内容： 机械求积的方法，梯形公式，simpson公式，求积公式代数精度，牛顿-柯特斯求积公式，精度分析，复化求积思想，复化梯形公式，复化辛普森公式，变步长梯形法，梯形的递推化算法，龙贝格加速算法，高斯型求积公式，数值微分。 教学目的及要求： 使学生了解数值求积公式构造的思想及方法，了解求积公式的代数精度等概念，掌握梯形的递推化算法，龙贝格算法，了解基本数值积分算法的数值实现，会求低阶的高斯型求积公式。 教学重点： 求积公式的代数精度、收敛性和稳定性的概念，复积求积思想以及梯形的递推化算法。 教学难点： 龙贝格算法的实现，复化求积公式的余项估计。 教学方法及手段： 通过定积分中值定理，引入定积分数值计算的关键是提供“平均高度”的一个有效算法，而给出一般数值求积公式的定义。由定积分的几何意义，引入梯形公式，通过分析梯形公式是被积函数的线性插值函数的积分，给出辛普森公式以及构造求积公式的一个方法，即利用被积函数的插值函数的积分，构造插值型求积公式的方法。在此基础上，自然地给出等距节点的柯特斯求积公式。通过“龙格现象”引入复化求积思想，给出复化梯形公式，复化辛普森公式，讨论它们的余项以及关系。通过复化梯形公式与复化辛普森公式的关系，给出龙贝格求积的思想以及相应的计算方法。通过构造最高精度的求积公式，给出高斯型求积公式定义及其构造方法。 在教学中，还要注重对各种求积算法的介绍与上机实践。 教学时间： 本章的教学的讲授时间为10学时，实验学时4学时。 教学内容： 2.1 机械求积 2.2 Newton-Cotes公式 2.3 Guass公式 2.4 复化求积法 2.5 Romberg加速算法 2.7 数值微分 §2.1 机械求积 一 数值求积基本思想 计算定积分，可应用牛顿-莱布尼兹公式 ，其中是在[a,b]上的原函数 但在实际计算中，会遇到两个问题： 1 在[a,b]区间上不具有有限形式的原函数。如，，被积函数，在[0,1]上没有有限形式的原函数； 2 被积函数由表格形式给出。 因此，我们有必要讨论定积分的数值计算问题。下面，我们来探讨一下如何实现定积分的近似计算。 当（）时，定积分在几何上代表着以[a,b]区间为底，以为曲边的曲边梯形面积。定积分中值告诉我们，存在着一点，使得。这表明，该曲边梯形面积等于底宽为b-a，高为的矩形面积。如果选取得不恰当，那么，这就给出定积分的一个近似计算公式。可见，定积分近似计算关键是寻找。我们可将理解为被积函数在[a,b]上的平均高度。而定积分近似计算的关键是为平均高度提供一种有效的算法。 如果我们取，那么 (4.1) 在几何上，它表示用梯形面积近似曲边梯形面积。若记，则，通常称这一公式为梯形公式。 对上述公式，我们换一个思路来分析。作的线性插取函数 那么，，，而 这给了我们启示：可以用被积函数的拉氏插值多项式来近似地代表被插值函数，而获得定积分近似公式。 如果我们利用数据点 ，， 其中，作拉格朗日插值多项式 对在[a,b]上积分，有 类似地， 故 (4.2) 若记，则，称该公式为辛普森公式。这里的平均高度为。 分析定积分近似计算的梯形公式，辛普森公式，不难发现这样一个事实：它们都是被积函数值的线性组合，且组合系数均与被积函数无关。 更一般地，我们可以在区间[a,b]上适当选取节点，然后用加权平均得到平均高度的近似值，这样构造出的求积公式具有下列形式： (4.3) 其中称为求积节点，称为求积系数，且仅与节点的选取有关，与被积函数无关。 这类数值积分公式通常称为机械求积公式，其特点是将积分求值问题归结为函数值的计算，这就避开了牛顿-莱布尼兹公式需要寻求原函数的困难。不难发现，构造形如(4.3)式的数值积分公式比较随意，那么，我们需要建立起评价数值积分公式“好坏”的标准。 二 代数精度 数值积分方法是近似方法，我们自然希望求积公式能对“尽可能多”的函数准确成立，这就提出了所谓代数精度的概念。 【定义1】如果某个求积公式对于次数不超过的多项式均能准确成立，但对于次多项式不准确成立，则称该求积公式具有次代数精度。 【引例1】讨论梯形公式，辛普森公式的代数精度。 令,,准确成立。 再令,,准确成立。 再令,,不准确成立。 故梯形公式具有1次代数精度。 令,,准确成立。 再令,,准确成立。 再令, 准确成立。 再令,,而 也是准确成立。 再令,,而 显然不准确成立。 故辛普森公式具有3次代数精度。 一般地，欲使求积公式(4.3)具有m次代数精度，只要令它对于都能准确成立，这就要求 三 插值型求积公式 设给定一组节点