正面教材分享 http://blog.sciencenet.cn/u/wdlang 70%的以色列人是无神论者,不过他们都相信上帝给了他们那块土地。这个世界经不起思考

博文

计算方法:梯形求积公式及其标度律

已有 7605 次阅读 2017-12-16 10:43 |个人分类:计算方法|系统分类:教学心得

梯形求积公式是最简单的数值积分工具。下面的程序利用积分

来计算ln2, 并研究其误差随采样点增加而减少的标度律。

很明显是平方律:

这个平方律是比较差的。我们往往需要配合Richardson插值法(也就是romberg方案 )来获得高精度结果。

但是,在某些时候,梯形公式会有极高的精度,这时标度律会是指数律,即误差随采样点增加指数减小。比如,我们可以考虑利用下面的表达式计算整数阶的bessel函数

相关程序为

从下图可以看到,对于二阶(一型)bessel函数在0.5处的值,只需要14个采样点,就可以达到机器精度!

这背后是有原因的。这里积分是对一个周期函数作一个周期内的积分。在这种情况下,fourier理论可以用上了。积分的准确值是零分量,数值误差来自某些高阶分量。偏偏这里的被积函数又有一个很好的性质,即如果往复平面解析延拓开去,其在全平面解析。按照paley-wiener定理的精神,这种函数的傅里叶分量会指数减小。

作业: 利用梯形公式计算第一和第二型椭圆积分



https://wap.sciencenet.cn/blog-100379-1089963.html

上一篇:德国地铁广告
下一篇:彩虹的数学与物理
收藏 IP: 218.85.66.*| 热度|

2 杨正瓴 徐令予

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 19:05

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部