Matlab 数值积分 差值 插值 数据积分与数值微分.docVIP

Matlab 数值积分 差值 插值 数据积分与数值微分.doc

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 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)

lizhencai0920 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:6100124015000001

1亿VIP精品文档

相关文档