《数学建模公选课课件》第八讲 数值积分.pptVIP

误差估计 每个子区间上的误差估计式为 复化Simpson公式的几何解释 将n个子区间的误差相加得 由闭区间上连续函数的介值性质可知在[a,b]上至少存在一点?，使 可见,当f(x)有四阶导数时,复化Simpson公式具有4阶收敛. * * 第八讲 数值积分 引言 2、积分是微分方程、概率论等的基础； 1、积分是高等数学中的一个最基本的运算，在实际问题中有着重要应用； 3、解决积分问题的最主要方法之一是Newton—Leibniz公式： 其中函数F是被积函数f在区间[a,b]上的一个原函数。 4、有很多函数“积不出来”，即使形式非常简单，例如： 5、即使能够求出原函数，但原函数的计算也很复杂，例如： 它的一个原函数是 6、f(x)没有解析表达式，只有数表形式: 8.5 5 8 4 6 3 4.5 2 4 1 f(x) x 这些都说明,通过原函数来计算积分有它的局限性； 事实上，在实际问题当中，存在大量Newton-Leibniz公式无法解决的积分。多数情形下，只能利用特定方法给出积分的近似值。 ——数值积分 8.1 实例及其数学模型 卫星轨道长度 问题： 人造地球卫星轨道可视为平面上的椭圆。我国第一颗人造地球卫星近地点距地球表面439km，远地点距地球表面2384km。已知地球半径为6371km，求该卫星的轨道长度。 x y O a 卫星 b c r s1 s2 卫星轨道示意图 常识： 卫星轨道是椭圆，地球处在椭圆的一个焦点上。 数学模型 （1）椭圆轨道的长半轴为a，短半轴为b，焦距为c； （2）近地点、远地点离地球表面分别为s1km、s2km； （3）地球半径为r km； 由图可知： 椭圆周长计算公式 设椭圆的参数方程为： 则椭圆周长为： 它称为椭圆积分，无法用解析方法计算。 8.2 数值积分公式、截断误差及代数精度 一、一般数值积分公式 假设已知函数f在节点xk处的函数值为f(xk)，其中k=0,1,2, …，则 称为数值积分公式。 xk称为求积节点，Ak（与函数f无关）称为求积系数。 构造或确定一个求积公式，要讨论解决的问题有 （3）求积公式的误差估计和收敛性。 （2）为了构造求积公式，需要提供一种判定求积方法精度高低准则； （1）确定求积系数Ak和求积节点xk； 三、数值积分公式的截断误差(求积余项) 二、数值积分公式的思想起源—— 积分的定义 四、数值积分公式的代数精度 定义1：称求积公式 具有m次代数精度,如果它满足如下两个条件: （1）对所有次数≤m的多项式 ,有 （2）存在m+1次多项式 ,使得 定义1中的条件（1）（2）等价于： 8.3 三个简单的数值积分公式 梯形公式 Simpson公式 其中， Cotes公式 一、三个数值积分公式 二、截断误差（求积余项）与代数精度 可以证明：梯形公式、Simpson公式、Cotes公式分别具有1次、3次和5次代数精度 8.4 插值型求积公式 一、公式的导出 在积分区间[a,b]上取n+1个节点xi，i=0,1,2,…,n, 构造f(x)的n次代数插值多项式（拉格朗日插值公式）： 其中 则有 其中， 为插值余项。 于是有 取 由 节点 决定， 与 f(x) 无关。 记 称(*)式为插值型求积公式,其中求积系数Ak由下 式确定. （*） 二、插值型求积公式的几何解释 三、插值型求积公式的截断误差 由 可知，插值型求积公式的截断误差为 四、插值型求积公式的代数精度 定理：具有n+1个求积节点的数值积分公式是插值型求积公式的充分必要条件为该公式至少有n次代数精度。 8.5 一类特殊的插值型求积公式—— Newton—Cotes公式 取节点为等距分布： 由此构造的插值型求积公式称为Newton-Cotes公式,此时求积系数 令 Cotes系数 注：Cotes 系数仅取决于 n 和 k，可查表得到。与 f (x) 及区间[a, b]均无关。 并且有 一、公式的导出 即， 其中， 并且有 上述数值积分公式被称为n阶Newton—Cotes公式 Newton-Cotes公式的误差为: Newton—Cotes公式的代数精度： 二、Newton—Cotes公式的截断误差与代数精度 当阶数n为偶数时, Newton-Cotes公式至少具有n+1次代数精度。 三、Newton-Cotes公式的数值稳定性 假设函数值 有误差 (j=0,1,2,…,n)。 的计算中,由 不考虑中间计算过程中的舍入误差，则在 引起的误差为 如果Cj(n)都是正数,并设 则有 故en是有界的,即由εj引起的误差受到控