ImageVerifierCode 换一换
格式:PPT , 页数:39 ,大小:1.91MB ,
资源ID:6944261      下载积分:10 金币
快捷下载
登录下载
邮箱/手机:
温馨提示:
快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。 如填写123,账号就是123,密码也是123。
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

加入VIP,免费下载
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.docduoduo.com/d-6944261.html】到电脑端继续下载(重复下载不扣费)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: 微信登录   QQ登录   微博登录 

下载须知

1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。
2: 试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。
3: 文件的所有权益归上传用户所有。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 本站仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

版权提示 | 免责声明

本文(Newton-Cotes公式.ppt)为本站会员(j35w19)主动上传,道客多多仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知道客多多(发送邮件至docduoduo@163.com或直接QQ联系客服),我们立即给予删除!

Newton-Cotes公式.ppt

1、第四章 数值积分 (Numerical Integration), f(x)的原函数F(x)为初等函数,本章主要讨论如下形式的一元函数积分,在微积分里,按Newton-Leibniz公式求定积分,要求被积函数f(x), 有解析表达式;,一、数值积分的必要性,1 引言 /*Introduction*/,例如函数:,考虑一个实际问题: 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的.,1、f(x)的原函数F(x)不能用初等函数表示,假若要求波纹瓦长4英尺,每个波纹的高度(从中心 线)为1英寸,且每个波纹以近似2英寸为一个周 期. 求制做一块波纹瓦所需铝板的长度L. 这个问题就是要求

2、由函数f(x)=sin x给定的曲 线,从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为:,上述积分称为第二类椭圆积分,它不能用普通方法来计算.,例如函数,并不复杂,但它的原函数却十分复杂:,2. 有些被积函数其原函数虽然可以用初等函数表示成有限形式,但表达式相当复杂,计算极不方便.,3. f(x)没有解析表达式,只有数表形式:,这些都说明,通过原函数来计算积分有它的局限性,因而,研究关于积分的数值方法具有很重要的实际意义.,在几何上可解释为由x=a, x=b, y=0和 y=f(x) 所围 成的曲边梯形的面积.积分计算之所以有困难, 就是因为这个曲边梯形有一条边y=

3、f(x)是曲的.,(1)求积公式的概念,积分值,二、求积公式及其代数精度,依据积分中值定理,对于连续函数f(x) ,在a,b内 存在一点,使得,称f()为区间a,b的平均高度. 问题在于点的具体位置一般是不知道的.这样,只要对平均高度f()提供一种算法,相应地便获得一种数值求积 方法.,左矩形公式: I(f)(b-a)f(a) 右矩形公式: I(f)(b-a)f(b) 中矩形公式: I(f)(b-a)f(a+b)/2,如果简单地选取区间a,b的一个端点或区间中 点的高度作为平均高度,这样建立的求积公式分 别是:,作为平均高度f()的近似值.,此外,众所周知的梯形公式:,I(f)(b-a)f(a

4、)+f(b)/2,和 Simpson公式:,I(f)(b-a)f(a)+4f(a+b)/2)+f(b)/6,则分别可以看作用 a, b, c=(a+b)/2三点高度的加 权平均值,f(a)+f(b)/2,和,f(a)+4f(c)+f(b)/6,更一般地,取区间a,b内n+1个点 xi,(i=0,1,2,n) 处的高度f(xi) (i=0,1,n)通过加权平均的方法 近似地得出平均高度f(),这类求积方法称为 机械求积:,或写成:,数值积分公式,求积系数,求积节点,记,称为求积公式余项(误差).,构造或确定一个求积公式,要讨论解决的问题有,(ii) 求积公式的误差估计和收敛性,为了构造数值求积公

5、式,需要提供一种 判定求积方法精度高低准则,称为数值求积公式,(i) 确定求积系数Ak和求积节点xk ;,定义中的条件(i),(ii)等价于:,具有m次代数精度,如果它满足如下两个条件:,(i)对所有次数 m次的多项式Pm(x),有,(ii)存在m+1次多项式Pm(x) ,使得,(2)求积公式的代数精度,定义1,称求积公式,1. T 公式的代数精度,2. S公式的代数精度,因此S-公式具有三次代数精度。,在积分区间a,b 上取n+1个节点xi,i=0,1,n,作 f(x)的n次代数插值多项式(拉格朗日插值公式):,则有,为插值余项,于是有,三、插值型求积公式,由节点决定,与f(x) 无关。,称

6、为插值型求积公式,其中,形如 的求积公式至少有 n 次代数精度 ,求积系数满足:,误差,定理1,推论1,该公式为插值型即:,2 Newton-Cotes 公式/*Newton-Cotes Rule*/,一、Newton-Cotes公式及数值稳定性,二、复化求积公式及误差估计, 取节点为等距分布:,令,Cotes系数,注:Cotes 系数仅取决于 n 和 k, 可查表得到. 与 f (x)及区间a,b均无关。,由此构造的插值型求积公式称为Newton-Cotes公式, 此时求积系数,(1)Newton-Cotes 公式,一、Newton-Cotes公式及数值稳定性,则,求积公式变为,称为n阶闭型

7、Newton-Cotes求积公式.,记,变为,Newton-Cotes公式的误差为:,与x有关,注意: Cotes系数只与j和n有关,与f(x)和积分区间a,b无关,且满足:,只需验证当n为偶数时,Newton-Cotes公式对f(x)=xn+1的余项为零.,由于f(x)=xn+1,所以f(n+1)(x)=(n+1)! .,引进变换t=u+n/2,因为n为偶数,故n/2为整数,于是有,据此可断定R(f)=0,因为上述被积函数是个奇函数.,定理3,当阶数n为偶数时, Newton-Cotes公式至少 具有n+1次代数精度.,证明,(2)Newton-Cotes公式的数值稳定性,设计算Cj(n)没

8、有误差,中间计算过程中的舍入误差 也不考虑,则在式In(f)的计算中,由j引起的误差为,现在讨论舍入误差对计算结果产生的影响.,设用公式,近似计算积分,时,其中计算函数值f(xj)有误差j (j=0,1,2,n).,故en是有界的,即由j引起的误差受到控制,不超过的(b-a)倍,保证了数值计算的稳定性.而当n7时,Cj(n)将出现负数,保证数值稳定性. 因此高阶公式不宜采用,有实用价值的仅仅是几种低阶的求积公式.,将随n增大,因而不能,如果Cj(n)都是正数,并设,则有,n = 1:,Trapezoidal Rule,/* 令 x = a+th, h = ba, 用中值定理 */,代数精度 =

9、 1,n = 2:,Simpsons Rule,代数精度 = 3,n = 4: Cotes Rule, 代数精度 = 5,高次插值有Runge 现象,高阶Newton-Cotes公式会出现数值不稳定,低阶Newton-Cotes公式有时又不能满足精度要求.解决这个矛盾的办法是将积分区间a,b分成若干小区间,在每个小区间上用低阶求积公式计算,然后将它们加起来,这就是复化求积方法.,二、复化求积公式及误差估计,(1)复化梯形公式:,在每个 上用梯形公式:,= Tn,/*中值定理*/,复化梯形公式积分法,(2)收敛性,由上述的误差估计式可知,当f(x) C2a,b 时,只 要h0时 数列Tn(f)

10、I(f),且收敛速度为二阶 O(h2).,但是f(x)C2a,b条件相对苛刻, 现假定f(x)在a,b上Riemann可积,讨论复化求积公式的收敛性,(3)复化 Simpson 公式,复化Simpson公式积分法,每个子区间上的误差估计式为,将n个子区间的误差相加得,由闭区间上连续函数的介值性质可知在a,b上至少存在一点,使,可见,当f(x)有四阶导数时,复化Simpson公式具有4阶收敛.,(4)误差估计,将区间0,1划分为8等分,应用复化梯形法求得,=0.9456909,例1,解,应用复化Simpson法计算,得,比较上面两个结果T8和S4,它们都需要提供9个点上的函数值工作量基本相同,然而精度却差别很 大.,同积分的准确值I(f)=0.9460831比较,复化梯形法的结果T8=0.9456909只有两位有效数字, 而复化,Simpson法的结果S4=0.9460832却有六位有效数字.,=0.9456909,

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报