- 1、本文档共57页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值积分 计算方法课件及实验 培训讲解.ppt
2. Simpson公式及其余项(n=2) Cotes系数为 求积公式为 上式称为Simpson求积公式,也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有 次代数精度。 3 3. Cotes公式及其余项(n=4) Cotes系数为 求积公式为 上式称为Cotes求积公式,也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度。 注:n?8时,Cotes系数出现负数,会引起误差增大,计算不稳定。 因此,在实际应用中一般不使用高阶Newton-Cotes公式,而是采用低阶复合求积法(下节)。 Cotes系数表: n Ck(n) 1 2 3 4 5 8 1/2 1/2 1/6 4/6 1/6 1/8 3/8 3/8 1/8 7/90 16/45 2/15 16/45 7/90 19/288 25/96 25/144 25/144 25/96 19/288 ………… 989/28350 5888/28350 -928/28350 10496/28350 -4540/28350 … 第三节 复合求积公式 高次插值有Runge 现象,故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式: 在每个 上用梯形公式: = Tn 由介值定理知: 使 即有: 余项: ? 复合 Simpson 公式: 4 4 4 4 4 = Sn 注:为方便编程,可采用另一记法:令 n’ = 2n 为偶数, 这时 ,有 例1: 分别利用复合梯形公式和复合Simpson公式计算积分: 积分的相对精确值为 解:设 =0 步长 h=1/8。 =0运算量基本相同 ?复合求积法的余项和收敛阶: 复合梯形( Trapezoid )公式的余项: 复合辛甫生(Simpson)公式的余项: 复合柯特斯(Cotes)公式的余项: * 计 算 方 法 第四章 数值积分 孙成立 §第四章 数值积分 §4.1 机械求积公式 §4.2 Newton_Cotes公式 §4.3 变步长求积公式及其加速 收敛技巧 §4.1 机械求积公式 第1节 引言 第2节 数值积分的基本方法 第3节 代数精度法 第4节 插值求积法 §4.1.1 引言 定积分的计算可用著名的牛顿-莱布尼兹公式来计算: 其中 F(x) 是 f (x) 的原函数之一,可用不定积分求得. 被积函数 f (x) 是用函数表格提供; f(x) 极为复杂,求不出原函数; 大量函数的原函数不容易或根本无法求出. 只能运用数值积分, 求积分近似值 . 问题 其中, 称为积分节点, 称为求积系数。 §4.1.2 数值积分的基本方法 就是在区间 [a , b] 内取 n+1 个点 利用被积函数 f (x) 在这 n+1 个点的函数值的某一种线性组合来近似作为待求的定积分. 求积分 此积分的几何意义相当于如图所示的曲边梯形的面积。 用 f (x) 的零次多项式 来近似代替 于是有 简单算例 解 左矩公式 推广: 右矩公式 中矩公式 用 f(x) 的一次多项式 来近似代替 , 于是, 梯形公式 推广: 来近似代替 ,于是, 特别地:当 有 用 f (x) 的二次插值多项式, 推广: Simpson公式 §4.1.3 代数精度法 为了使一个求积公式能对更多的积分具有较好的实际 计算意义,就要求它对尽可能多的被积函数都准确地成立. 因此定义代数精度的概念: 若积分 的数值积分公
文档评论(0)