摘 要 数学领域中,定积分计算问题应用广泛,经典的定积分数学定义方法可直接用于求解定积分,但是,对于函数解析式未知的情况下,传统的数学定义方法无法进行定积分计算,而蒙特卡罗方法对函数解析式不进行限制,其以概率方法进行近似计算从而逐渐逼近定积分理论值。本文针对函数解析式已知与未知的两种情况,分别以定积分数学定义方法和蒙特卡罗方法进行定积分计算,并从算法收敛速度以及计算结果精确度两方面对算法进行评测。实验结果表明,定积分数学定义法收敛速度快,计算精度高,但是普适性低,对于函数解析式未知情况下无法进行计算;而蒙特卡罗方法尽管收敛速度较慢,但是普适性极高,且函数解析式未知情况下,效果更优。

关键词 微积分;定积分;蒙特卡罗方法;收敛速度

中图分类号 O1 文献标识码 A 文章编号 2095-6363(2017)07-0003-02

1 概述

微积分[ 1 ]是数学领域的一个基础学科,是高等数学中研究函数微分、积分以及有关概念与应用的数学分支,其研究范畴包含3个方面:微分、积分以及微分与积分两者之间的关系。

若( )f x是[a,b]上的连续函数,并且有’( )( )F xf x=,则( )( )( )baf x dxF bF a∫=?。也就是说,一个定积分的值就是原函数积分上限的值与原函数在积分下限的值的差值,即牛顿-莱布尼兹公式计算定积分。其表明对于图形无限细分再累加成为可能,并可将其转化为对积分的计算,揭示了积分与微分本质的关系,因此牛顿-莱布尼兹公式又称微积分基本定理。

然而计算定积分[ 2 ]的数学定义方法以及牛顿-莱布尼兹公式方法都仅限于函数( )f x解析式已知的情况,对于( )f x未知解析式的情况下,无法进行定积分求解。

蒙特卡罗(Monte-Carlo)[3]方法是20世纪40年代中期由于科学技术的发展和电子计算机的发明,被提出的一种以概率理论[4]为指导的一类极其重要的数值计算方法,是以随机抽样为主要手段,使用随机数(或伪随机数)解决数值计算问题的方法,又称统计模拟方法。蒙特卡罗方法是一种重要的利用计算机模拟的近似计算方法,主要用于解决确定性的数学问题(如计算定积分)和随机性问题(如扩散问题),广泛应用于各个领域。对于定积分计算领域因函数( )f x解析式未知而无法运用定义方法进行定积分计算的难题,蒙特卡罗方法以统计模拟方法进行定积分计算。

2 蒙特卡罗方法计算定积分

定积分就是求解函数( )f x在区间[a,b]上图线下方包围的面积,即在Oxy坐标平面上,曲线( )f x与直线xa=、xb=以及x轴围成的曲边梯形的面积值(确定的实数值)。其数学定义为:若函数( )f x在区间[a,b]上连续,以平行于y轴的直线分割图象为无数个矩形,而后累加区间[a,b]上的矩形。具体方法如下:

对于定积分计算问题,若函数( )f x解析式已知,则可以传统的数学定义方法可直接求解定积分,亦可以蒙特卡罗方法以概率的近似计算方法求解定积分;但是对于函数( )f x解析式未知的情况,传统的数学定义方法无法进行定积分计算,而以概率论进行近似计算的蒙特卡罗方法对函数解析式o限定,依旧适用于定积分计算。

3 实验结果

针对函数( )f x解析式已知与未知的两种情况,分别以数学定义方法与蒙特卡罗方法进行定积分计算,实验结果如下。

4 结论

综上所述,两种方法在不同情况下收敛效果不同。定积分数学定义方法适用于给定函数解析式的定积分求解,其收敛速度快,计算结果精度高,然而普适性较低,难以实现函数解析式未知的定积分求解问题;而蒙特卡罗方法尽管收敛速度较慢,精度较低,实验结果误差并非单调递减趋势,但是普适性强,对函数解析式未知的定积分计算依然适用。

参考文献

[1]同济大学数学编委组.微积分:上册[M].3版.北京:高等教育出版社,2009.

[2]马晓涛,马华.定积分计算中的几个常用方法[J].高等数学研究,2005,8(6):36.

[3]尹增谦等.蒙特卡罗方法及其应用[J].物理与工程,2002,12(3):45-49.

[4]马振华.现代应用数学手册:概率统计与随机过程卷[M].北京:清华大学出版社,2000:123-130.

[5]宫野.计算多重积分的蒙特卡罗方法与数论网络法[J].大连理工大学学报,2001,4(11):20-23.

蒙特卡罗方法C语言求定积分,蒙特卡罗方法计算定积分相关推荐

  1. 创建一个矩形类的子类——正方形类,类中定义求面积方法、重写求周长的方法.在主类中,输入一个正方形边长,创建正方形对象,求正方形的面积和周长

    题目: 编写一个应用程序,创建一个矩形类,类中具有长.宽两个成员变量和求周长的方法.再创建一个矩形类的子类--正方形类,类中定义求面积方法.重写求周长的方法.在主类中,输入一个正方形边长,创建正方形对 ...

  2. 蒙特卡罗方法C语言求定积分,邓一硕: 蒙特卡洛方法与定积分计算 | 统计之都 (中国统计学门户网站,免费统计学服务平台)...

    大致看了看,非常好的一份材料,不仅研究了Monte Carlo,还研究了Bayesian,非常值得大家看哦,不过缺点就在于这份PDF没有目录,所以我为这个材料做了一份目录,希望可以帮助大家学习: Mo ...

  3. c语言水仙花两种方法,c语言:3种方法;求出0~999之间的所有“水仙花数”并输出。...

    方法一: #include int main() { int i,j,k,n; printf("水仙花数:",n); for(n=100;n<1000;n++) { i=n/ ...

  4. c语言常用库函数使用方法,c语言常用库函数使用方法及用途

    absread()读磁盘绝对扇区函数 原形:int absread(int drive,int num,int sectnum,void *buf) 功能:从drive指定的驱动器磁盘上,sectnu ...

  5. R语言编程的高效方法

    R语言编程的高效方法 R语言编程的高效方法 高效编程之细节知识点 并行计算 1.查看核数 2.可以使用并行计算的场景 3.Parellel 包 R语言编程的高效方法 学习资料来源: datacamp ...

  6. 计算机显示语言改成英语,如何更改电脑语言 电脑语言改为英文方法 (全文)

    今年evad3rs团队发布的iOS 7完美越狱工具有些特殊,中国版iOS7完美越狱工具中,捆绑有太极助手插件,此插件类似流氓插件,容易导致iOS越狱设备出现白苹果现象(国外版无此插件).目前比较好的解 ...

  7. c语言文件有什么作用,C语言文件的使用方法

    C语言文件的使用方法 C语言文件的使用方法:所谓"文件"是指一组相关数据的有序集合. 这个数据集有一个名称,叫做文件名. 实际上在前面的各章中我们已经多次使用了文件,例如源程序文件 ...

  8. 高等数学:第五章 定积分(1) 定积分的概念

    §5.1  定积分的概念 一.从阿基米德的穷竭法谈起 [引例]从曲线与直线,, 所围图形的面积. 如图:在区间  上插入  个等分点 ,得曲线上点 ,过这些点分别向轴,轴引垂线,得到阶梯形.它们的面积 ...

  9. Matlab在线IDE:MATLAB Online介绍与计算定积分案例

    目录 1.MATLAB Online介绍 功能与特点 命令行窗口和编辑器窗口 适用场景 计费方式 使用方法 2.注册登录 3.计算定积分 1.MATLAB Online介绍 MATLAB Online ...

  10. 用C语言编码定积分,C语言__用六种方法求定积分C语言__用六种方法求定积分.doc...

    C语言__用六种方法求定积分C语言__用六种方法求定积分 描述问题 利用①左矩形公式,②中矩形公式,③右矩形公式 ,④梯形公式,⑤simpson公式,⑥Gauss积分公式求解定积分. 分析问题 2.1 ...

最新文章

  1. Opengl-光照-基本光照-光照贴图(现在告别单调的方块弄个箱子)
  2. MySQL This function has none of DETERMINISTIC, NO SQL...错误1418 的原因分析及解决方法
  3. 搜索推荐系统实战:进化篇
  4. python从sqlserver提取数据_通过Python读取sqlserver数据写成json文件的总结
  5. 《Neo4j全栈开发》_陈韶健
  6. 了解如何使用Vue.js CLI
  7. MFC工作笔记0005---::在vc++中是什么意思
  8. Oracle管理文件OMF (oracle managed files)
  9. Java正则表达式入门概念与范例代码
  10. (原创)javascript进阶(三)
  11. STM32单片机全自动锂电池容量电量检测放电电流电池电压ACS712
  12. 产品经理职责技能和所需证书
  13. Transformer课程 业务对话机器人Rasa 3.x 生成自然语言理解NLU数据
  14. fatal: the remote end hung up unexpectedly问题解决办法
  15. 爬虫(python)—下载技巧
  16. Redis缓存击穿、雪崩、穿透!(超详细)
  17. 元素选择器(type selectors)
  18. 推荐系统论文重要的三个指标——Recall、NDCG、RMSE
  19. Fragment+RadioButton实现点击切换页面效果
  20. 保险成电商角逐互联网金融利器

热门文章

  1. 计算机怎么求地址线,如何计算地址线和数据线
  2. Gd-DTPA钆元素钆元素对比磁共振造影剂 氟-19磁共振/荧光双模态显影剂19F MRI
  3. linux压缩到最小命令,Linux 压缩打包命令详细教程
  4. Java人脸识别考勤系统_GitHub - magician1128/WorkAttendanceSystem: 一个基于opencv、dilb的员工人脸识别考勤系统...
  5. Android-APP内存优化
  6. 求解一元二次方程的根
  7. 微信能转发邮件吗?邮件怎么用微信发送?微信邮件怎么发?
  8. macbook删除windows后合并磁盘分区
  9. 微信引流黑科技:手机浏览器直接唤起微信方案调研
  10. Unity实现音乐播放器