第五章数值积分精选文档.ppt
《第五章数值积分精选文档.ppt》由会员分享,可在线阅读,更多相关《第五章数值积分精选文档.ppt(37页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第五章数值积分本讲稿第一页,共三十七页 在在a,b上取上取 a x0 x1 xn b,做,做 f 的的 n 次插值次插值多项式多项式 ,即得到,即得到Ak由由 决定,决定,与与 无关。无关。节点节点 f(x)插值型积分公式本讲稿第二页,共三十七页误差误差本讲稿第三页,共三十七页例:例:对于对于a,b上上1次插值,有次插值,有考察其代数精度。考察其代数精度。f(x)abf(a)f(b)梯形公式解:解:逐次检查公式是否精确成立逐次检查公式是否精确成立代入代入 P0=1:=代入代入 P1=x:=代入代入 P2=x2:代数精度代数精度=1本讲稿第四页,共三十七页抛物型求积公式:二次插值求积公式:本讲稿
2、第五页,共三十七页Simpson公式本讲稿第六页,共三十七页若若某某个个求求积积公公式式所所对对应应的的误误差差R f 满满足足:R Pk=0 对对任任意意 k n 阶阶的的多多项项式式成成立立,且且 R Pn+1 0 对对某某个个 n+1 阶阶多多项项式式成成立立,则则称此求积公式的称此求积公式的代数精度代数精度为为 n。代数精度:代数精度:本讲稿第七页,共三十七页怎样验证代数精度:怎样验证代数精度:本讲稿第八页,共三十七页注:形如 的求积公式至少有 n 次代数精度 该公式为插值型(即:)本讲稿第九页,共三十七页思考:思考:代数精度是否是越高代数精度是否是越高越好?越好?本讲稿第十页,共三十
3、七页梯形公式的误差梯形公式的误差本讲稿第十一页,共三十七页定理证明本讲稿第十二页,共三十七页辛普森辛普森(Simpson)(Simpson)求积公式的误差:求积公式的误差:本讲稿第十三页,共三十七页思考:思考:结论是什么?结论是什么?怎么办?怎么办?本讲稿第十四页,共三十七页复合求积:复合求积:高次插值有高次插值有Runge 现象,故采用分段低次插值现象,故采用分段低次插值 分段低次合成的分段低次合成的 复合求积公式。复合求积公式。复合梯形公式:复合梯形公式:在每个在每个 上用梯形公式:上用梯形公式:本讲稿第十五页,共三十七页=Tn/*中值定理中值定理*/怎么办?怎么办?本讲稿第十六页,共三十
4、七页44444=Sn注:为方便编程,可采用另一记法:令 n=2n 为偶数,这时 ,有 复合复合Simpson公式:公式:本讲稿第十七页,共三十七页复化求积例:复化求积例:本讲稿第十八页,共三十七页复化求积例:复化求积例:两种方法谁好?两种方法谁好?本讲稿第十九页,共三十七页给定精度给定精度 ,如何取,如何取 n?通常采取将区间通常采取将区间不断对分不断对分的方法,即取的方法,即取 n=2k上例中上例中2k 68 k=7注意到区间再次对分时注意到区间再次对分时可用来判断迭代是否停止。本讲稿第二十页,共三十七页5.2 5.2 高斯型高斯型积分积分构造具有构造具有2n+1次代数精度的求积公式次代数精
5、度的求积公式将节点将节点 x0 xn 以及系数以及系数 A0 An 都作为待定系数。都作为待定系数。令令 f(x)=1,x,x2,x2n+1 代入可求解,得到的公式具代入可求解,得到的公式具有有2n+1 次代数精度。这样的节点称为次代数精度。这样的节点称为Gauss 点点,公,公式称为式称为Gauss 型求积公式型求积公式。本讲稿第二十一页,共三十七页例:在两点数值积分公式中,如果积分点也作为未知量,在两点数值积分公式中,如果积分点也作为未知量,则有则有4 4个未知量个未知量,可以列出可以列出4 4个方程个方程:(在:(在-1,1-1,1为例)为例)可解出:可解出:数值积分公式数值积分公式具有
6、3阶代数精度,比梯形公式1阶代数精度高本讲稿第二十二页,共三十七页推广:推广:n加权加权GaussGauss积分公式积分公式权函数本讲稿第二十三页,共三十七页例:例:求求 的的 2 点点 Gauss 公式。公式。解:设解:设 ,应有,应有 3 次代数精度。次代数精度。+101100)()()(xfAxfAdxxfx代入代入 f(x)=1,x,x2,x3不是线性方程组,不易求解。本讲稿第二十四页,共三十七页 x0 xn 为为 Gauss 点点 与与任任意意次次数数不不大于大于n 的多项式的多项式 P(x)(带权)正交(带权)正交。定理求 Gauss 点 求w(x)本讲稿第二十五页,共三十七页证明
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 数值 积分 精选 文档
限制150内