- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第章 14.1 数值积分的基本概念
14.1.1 求积公式
设是的函数,对于定积分:,如果被积函数是用解析式表示的连续函数,从理论上讲可以使用牛顿-莱布尼兹 (Newton-Leibniz) 公式,即:
(14.1.1)
但如果:
找不到用初等函数表示的原函数,如等等。
虽然找到了原函数,但因表达式过于复杂而不便于计算。
因此,建立有效的数值积分的方法就很有必要了。从上所述,数值求积公式应该避免用原函数表示,最好能由被积函数的直接决定,这种做法是否可行呢?
图.1.1所示上与轴之间的面积。
图14.1.1 积分示意图
由定积分的定义:
(14.1.2)
:
(14.1.3)
式为求积节点,而,它仅与节点的选取有关,而的具体形式上式便称为求积公式。
这类数值积分的方法通常称为机械求积法,它利用一些离散点上的函数的值作线性组合而得出积分的近似值。于是,求积分的问题转化为计算被积函数在节点处函数值的问题了。对于形如(14.1.3)的求积公式,关键在于确定求积节点和相应的求积系数
14.1.2
求积公式与精确值的差为求积公式的余项,也称为:
(14.1.4)
,而对于m+1以上的多项式的却不能精确等于定积分的数值,则称该求积公式的代数精度为m。
任何一个求积公式的代数精度至少为0,即取,求积公式精确地等于:
(14.1.5)
上式是求积系数应满足的基本条件,可以用它来检验一个求积公式的系数的正确性。一般地,要求(14.1.3)的代数精度为,只要令它对于都能精确成立,也就是说要求:
(14.1.6)
14.2 插值型求积公式
14.2.1 插值型求积公式
及函数在这些节点上的值,按前面所讲的插值法,。这样, (14.2.1)
式中,求积通过插值基函数积分得出,即:
(14.2.2)
由插值余项定理,可得到插值型求积公式的为:
(14.2.3)
与变量有关,且。
如果求积公式(14.2.1)是插值型的,由(14.2.3)式知,对于次数小于的多项式,其余项等于0,所以,此时求积公式的代数精度至少为。反之,如果的,则它必定是。
14. 内插求积公式
由插值公式,对任一函数(包括离散形式的离散数据),可用一个次多项式对其进行插值,当插值多项式为拉格朗日多项式时,即:
(14.2.4)
(14.2.5)
其中:
(14.2.6)
(14.2.7)
14.2.3 牛顿-柯特斯公式
牛顿-柯特斯 (Newton-Cotes) 公式:等距节点的内插求积公式。将积分区间等分,得牛顿牛顿-柯特斯公式:
(14.2.8)
其中:
(14.2.9)
这种等距节点的内插求积公式通常称为牛顿-柯特斯公式。下面先介绍几个牛顿-柯特斯公式的特殊形式。
14.24 梯形公式
在牛顿-柯特斯公式中取得梯形公式为:
(14.2.10)
梯形公式是用梯形面积近似代替曲边梯形的面积。
梯形公式的误差估计为:
,() (14.2.11)
图14.2.1 梯形公式示意图
14.25 抛物形公式
在牛顿-柯特斯公式,取得抛物形公式 (又称辛普生Simpson公式):
(14.2.12)
式中:
,
文档评论(0)