计算定积分.doc

数学实验报告 实验序号： 日期：2011 年 6 月 28 日 班级 09B 姓名 郭世伟 学号 094080106 实验 名称 用Monte-Carlo方法计算定积分(随机投点法和平均值法) 问题的背景： 不少统计问题,如计算概率，各阶矩等,最后都归结为定积分的近似计算问题(特别是高维积分).这一方法是求解数学物理,工程技术及经济管理近似的数值常用方法. 实验目的： 1. 学习和掌握Excel的有关命令. 2. 进一歩理解大数定律. 3. 掌握随机模拟的思路. 4. 比较两种方法的精度.(n=100,1000,10000) 实验原理与数学模型： （1） 随机投点法 设 其中在[0，1]上连续，且有. 设服从[0,1;0,1]上的均匀分布，则可知服从[0,1]上的均匀分布，且独立.因为 所以，定积分的值就是事件的概率.这样，根据大数定律，我们可以用频率来近似概率. 具体做法如下：先产生2上均匀的随机数：.然后 的次数，最后的估计值，即 注：为了计算一般的定积分，我们可以这样来解决.设 其中有界，不妨设. 令 则有 我们只要求即可. （2） 平均值法 设 其中在[0，1]上连续. 设服从[0,1]上的均匀分布，则的数学期望为 所以估计的值就是估计的数学期望的值.由辛钦大数定律，可以用的观察值的平均值去估计的数学期望.具体做法如下：先产生上均匀的随机数：.然后对每个计算，最后的估计值，即 实验所用软件及版本： 主要内容(要点): (1)用随机投点法和平均值法计算定积分: (2)比较两种方法的精度.(n=100,1000,10000.对每一个n重复做5次) 实验过程：（含解决方法和基本步骤，主要程序清单及异常情况记录等） 问题1（随机投点法）：估计定积分 令 当. （1）从主菜单中选择【工具|数据分析】，在【数据分析】对话框中滚动列表框，选择“随机数发生器”选项，如图9－1 图9－1 （2）重重两次可得100对随机数，作变换，然后进行比较，如图9－2，最后进行数数，得. 图9－2 估计定积分 下面我们估计，从上面不难得到，所以 ＝3.039. 问题2（平均值法）：估计定积分 令 . 由上表可以算得：＝0.401 估计定积分 因为所以我们首先产生 （－1，1）的均匀随机数.其它类似. 如图9－3. 所以，. 从上述结果可以看出，误差较大，必须加大容量. 实验结果与实验总结(体会)： 第一次n=100 k= s= J= 41 0.41 2.663089 第二次n=1000 k= s= J= 443 0.443 2.818215 第三次n=10000 k= s= J= 4220 0.422 2.719498 进一步讨论或展望: 做了试验后我发现投点法和平均值法求的积分值很相近，而做一次是不够的，我们需要多次做后才能发现这个规律，而这个也是检验我们实验是否做着正确的依据。 当做同一n值的多次试验时，只要我们改动其中的随机变量的值，也就是用随机变量发生器重新产生一组数据我们就可以得到一组新的数据，这就会很方便，这是因为Excel2003的优越性，一旦公式写好后只要改变相应值就会引其结果的变化，这在以后我们的学习中会有很大帮助。 教师评语与成绩：