《第5章数值积分精选PPT.ppt》由会员分享，可在线阅读，更多相关《第5章数值积分精选PPT.ppt（40页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第5章数值积分第1页，本讲稿共40页一、梯形公式一、梯形公式 过a,b两点，作直线方程 用 代替 ，得 第2页，本讲稿共40页 二、辛浦生（二、辛浦生（SimpsonSimpson）公式（抛物线求积公式）公式（抛物线求积公式）把 区间二等分，过a,b和 三点，作抛物线：构造 的过程：利用Lagrange插值法得公式：第3页，本讲稿共40页第4页，本讲稿共40页用 代替 则求得：三、牛顿一柯特斯公式（三、牛顿一柯特斯公式（Newton-CotesNewton-Cotes）区间分法：把区间 n等分，其分点为 过这 个节点，构造一个n次多项式 第5页，本讲稿共40页 其中 ，用 代替被积函数 ，则有
2、 其中 计算系数 ，用变量替换 ，于是 第6页，本讲稿共40页第7页，本讲稿共40页故 引进记号 第8页，本讲稿共40页则 计算牛顿一柯特斯系数只要能给出n就能求出 例：时有 相应的求积公式为 第9页，本讲稿共40页时有：N=4时，相应的求积公式为：第10页，本讲稿共40页 *的 公式，特别称为柯特斯公式。*教材表5.1给出了n=1至6的牛顿一柯特斯系数表。例1试用梯形求积公式，抛物线求积公式和牛顿一柯特斯公式（取n=4）计算定积分(1)梯形求积公式(2)抛物线求积公式 第11页，本讲稿共40页(3)牛顿一柯特斯求积公式（取 ）*积分的准确值 5.2 5.2 梯形求积公式和抛物线求积公式的误差
3、估计梯形求积公式和抛物线求积公式的误差估计定义：对一个一般的求积公式 第12页，本讲稿共40页 若求积公式 为任意一个次数不高于m次的代数多项式时，等号成立，而对 是m+1次多项式时，公式不能精确成立，则我们说求积公式具有m次代数精确度。梯形求积公式具有一次代数精确度。当 在上存在 ，用一次多项式 逼近 有关系式：根据假定 是一次多项式，则，因此，则第13页，本讲稿共40页但当 时 梯形公式具有一次代数精确度牛顿一柯特斯公式的代数精确度至少是n。定义定义：对牛顿一柯特斯公式，当在 存在，用n次插值多项式 逼近 有关系式：假若 是n次多项式，则 第14页，本讲稿共40页因此，牛顿一柯特斯公式的代
4、数精确度至少是n。证证明明:当n为偶数时，牛顿一柯特斯公式的代数精确度至少是n，可以达到记 次多项式为则 次导数为对 两边积分得：第15页，本讲稿共40页令 为正整数，引进变换 ，有令第16页，本讲稿共40页则 也就是说 是一个奇函数，故 即当n为偶数时，牛顿一柯特斯公式对 次多项式精确度成立，因此代数精确度达到 ，抛物求积公式是 时的牛顿一柯特斯公式，代数精确度至少是3,容易证明，抛物线求积公式对四次多项式不能精确成立。取 第17页，本讲稿共40页因为 所以抛物线求积公式的代数精确度是3。定理定理5.15.1若 ，则梯形求积公式有误差估计 证明证明：根据定理2.1有第18页，本讲稿共40页其
5、中 是依赖于x的函数，两边积分得 因为 在 上连续，而 在 上恒小于零，利用积分中值定理，在 上存在一点 ，使第19页，本讲稿共40页因此 定理得证。5.3 5.3 复化公式及其误差估计复化公式及其误差估计 一、复化梯形公式 把区间 n等分，节点第20页，本讲稿共40页，对每个小区间用 梯形公式，则得 第21页，本讲稿共40页 称为复化梯形求积公式，下标 表示将区间 等分，若把区间2n等分，在每个小区间上仍用梯形求积公式，则可得到 与之间有关系式：其中（二）复化抛物线公式第22页，本讲稿共40页 因为抛物线公式用到了区间的中点，所以构造复化抛物线公式时，必须把区间 等分为偶数份。为此令 为正整
6、数，在每个小区间 上用抛物线求积公式其中，因此第23页，本讲稿共40页称为复化抛物线求积公式，是区间等分数（三）复化梯形求积公式的误差估计 定理5.3若 ，则其中 证明证明：把区间 等分，并在每个小区间上直接用梯形求积公式的误差估计（5-9）则 第24页，本讲稿共40页由于 在 上连续，利用连续函数的性质知在 中存在一点 ，使这样，就得到了复化梯形求积公式的截断误差第25页，本讲稿共40页（四）复化抛物求积公式的误差估计定理5.4设 ，则对复化抛物型公式有误差估计 其中 例2 计算积分 要求保证有五位有效数字，问若用复化梯形公式，应取多少？若用复化抛物求积公式计算，n又应取多少？第26页，本讲
7、稿共40页解：由 有 ，因b-a=1，故当 时有根据定理5.3有 的真值具有一位整数，所以若要求积分具有五位有效数字，只要取（根据有效数字误差限）第27页，本讲稿共40页两边取对数并整理：只要 即可，也就是说，把区间 等分为68份就可以满足计算要求。用复化抛物线求积公式计算，由（516）有 两边取对数并整理得 第28页，本讲稿共40页因此 即可，或者说，把 六等分就可以满足。5.4 逐次分半法 我们看到，截断误差随 的增长而减小，但对一个给定的积分问题，选定了某种求积的方法，如何确定适当的 ，使近似值和真值之差在允许范围之内呢？所以自动选取积分步长是在求积过程中，根据精度要求，自动确定 并算出
8、满足精度的近似值。当区间 等分时截断误差为 第29页，本讲稿共40页当区间 等分时，截断误差为 将上面两式相减得到 当 在区间 上连续，并假定 充分大时 第30页，本讲稿共40页有：故此可以利用此公式估计误差(允许误差)来判断近似值 是否已满足要求，下面给出计算过程：（1）取 ，计算下 ：（2）把区间 分割为二份，取 计算 第31页，本讲稿共40页（3）把区间 四等分，计算 一般的计算公式为 在计算序列 过程中，每算出一个新的近似值 时，就用条件 （允许误差）来检验，当上式满足，计算就停止。第32页，本讲稿共40页例例3 3.用复化梯形公式，复化抛物线公式和 的牛顿一柯特斯公式计算积分（见教材
9、）5.5 5.5 数值方法中的加速收敛技巧数值方法中的加速收敛技巧 假设有一个量 ，用一个步长为 的函数 去逼近，与 无关，对于 的截断误差有估计式：（）第33页，本讲稿共40页例例 利用数据表利用数据表xk01/81/43/81/25/83/47/81f(xk)43.938463.764703.506853.200002.876402.460002.265492计算积分计算积分这个问题有明显的答案这个问题有明显的答案取取n=8用复化梯形公式用复化梯形公式第34页，本讲稿共40页取取n=4,用辛卜生公式用辛卜生公式第35页，本讲稿共40页都是与 无关的常数。也就是说 逼近 的阶是 。令 用 来
10、代替，则有现在用 乘以 式的两边后和上式相减，整理得：若 ，用 除等两边有第36页，本讲稿共40页其中都是与h无关的常数，令那么 逼近 的误差的阶为 ，一般地，当我们选定 为满足 的适当正数时，计算公式为：第37页，本讲稿共40页用归纳法容易证明 逼近 的误差为这里 是与 无关的常数5.5 Romberg5.5 Romberg求积法求积法 Romberg求积是在复化求积公式的基础上，应用Richaroson外推构造的一种算法。方法：把区间 几等分，然后用复化梯形求积公式求得近似值 记为 ，再把 等分，求得近似值 记为 ，如此等等第38页，本讲稿共40页可以得到一个序列 ，它逼近真值的速度为 。取 并注意到 ，有利用 和 可以得到第39页，本讲稿共40页一般地 其中 对 来说截断误差的阶是第40页，本讲稿共40页