- 1、本文档共26页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
[牛顿科特斯求积公式的代数精度研究买买提克热木·亚森信息085班2013.4.14
编号
学士学位论文
牛顿-柯特斯求积公式的
代数精度研究
学生姓名:买买提克热木·亚森
学 号:20080101025
系 部:数学系
专 业:信息与计算科学
年 级: 2008-5班
指导教师:阿米娜· 沙比尔
完成日期: 2013年4月20日
摘要
本文主要对求数值积分公式的代数精度进行探讨。首先描述了数值积分的矩形法,梯形法,插值求积公式等求积方法的基本思路和代数精度概念,进行了余项估计。然后重点讨论对牛顿—柯特斯求积公式当的情形及其代数精度.最后用数值例题验证牛顿—柯特斯求积公式的代数精度的重要性。
关键词:数值积分;梯形公式;代数精度;牛顿-柯特斯公式;
目 录
摘要 1
引言 2
1基本概念 4
1.1代数精度的概念 4
1.3 插值型的求积公式 6
1.4 求积公式的余项 7
1.5 求积公式的收敛性和稳定性 9
2. 牛顿—柯特斯公式 10
2.1 柯特斯系数与辛普森公式 10
2.2偶阶求积公式的代数度 18
2.3 辛普森公式的余项 19
参考文献 23
致谢 24
引言
1问题的提出
实际问题当中常常需要计算积分,有些数值方法,如微分方程和积分方程的求解,也都和积分计算相联系.依据人们所熟知的微积分基本定理,对于积分,只要找到被积函数的原函数,便有下列牛顿—莱布尼茨(Newton—Leibniz)公式
但实际使用这种求积方法往往有困难,因为大量的被积函数,比如,等,其原函数不能用初等函数表达,故不能用上述公式计算,即使能求得原函数的积分有时计算也十分困难.
例如对于被积函数,其原函数
计算仍然很困难.另外,当是有测量或数值计算给出的一张数据表时,牛顿—莱布尼茨公式也不能直接运用.因此很有必要研究计算方便,计算量少,精度高而且稳定的数值积分方法.
2我的想法
积分中值定理告诉我们,在积分区间内存在一点,成立
就是说,低为而高为的矩形面积恰等于所求曲边梯形的面积,如图4.1 ,问题在于点的具体位置一般是不知道的,因而难以准确算出的值,我们将称为区间上的平均高度。这样只要对平均高度提供一种算法,相应地便获得一种数值求积方法,即如果我们用两端点“高度”与的算法平均值作为平均高度的近似值,这样导出的求积方法
(1.1)
便是我们所熟悉的梯形公式,几何意义看图4.2 ,而如果改用区间中点的“高度”近似地取代平均高度,则有可导出所谓中矩形公式
(1.2)
更一般地,我们可以在区间上适当取某些节点,然后用的加权平均得到平均高度的近似值,这样构造出的求积公式有
(1.3)
式中称为求积节点,称为求积系数,赤称伴随节点的权.权仅仅与节点的选取有关,而不依赖于被积函数的具体形式.这类数值积分方法通常称为机械求积,其特点是将积分求值问题归结为被积函数值的计算,这就避开了牛顿—莱布尼茨公式需要寻求原函数的困难,很适合在计算机上使用.
有时当被积函数的原函数过于复杂时,也不宜套用积分公式,而用数值积分公式.在微积分中,定积分是Riemann(黎曼)和的极限,它是分割小区间趋于零时的极限,即
在数值积分公式中,只能用有限项的和近似上面的极限,通常由函数在离散点函数值的线性组合形式给出,记
,
在本章中,用表示近似积分值,称为求积节点,称为求积系数,确定中积分系数的过程就是构造数值积分公式的过程.
怎样判断数值积分公式的效果?代数精度是衡量数值积分公式优劣的重要标准之一.
由此可知当具有m阶代数精度时,对于不高于m次多项式都有.
1基本概念
1.1代数精度的概念
定义1 如果某个求积公式对于次数不超过的多项式均能准确地成立,但对于次多项式就不准确成立,则称该求积公式具有次代数精度(或代数精度确度).
记上以为积分节点的数值积分公式为
若满足
而
则称具有m阶代数精度 .
不难验证,梯形公式(1.1)和矩形公式(1.2)均具有一次代数精度.
一般地,欲使求积公式(1.3)具有次代数精度,只要令它对于都能准确成立,这就要求
(1.4)
为简洁起见,这里省略了符号中的上、下标.
如果我们事先选定求积节点,比如以区间的等距分点作为节点,这时取求解线性方程组(1.4)即
文档评论(0)