版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数学实验报告实验序号:5日期:2004年5月7日班级B班姓名学号实验名称用Monte-Carlo方法计算定积分(随机投点法和平均值法)问题的背景:不少统计问题,如计算概率,各阶矩等,最后都归结为定积分的近似计算问题(特别是高维积分).这一方法是求解数学物理,工程技术及经济管理近似的数值常用方法.实验目的:1.学习和掌握Excel的有关命令.2.进一歩理解大数定律.3.掌握随机模拟的思路.4.比较两种方法的精度.(n=100,1000,10000)实验原理与数学模型:随机投点法设定积分其中在[0,1]上连续,且有.设服从[0,1;0,1]上的均匀分布,则可知服从[0,1]上的均匀分布,且独立.因为,所以,定积分的值就是事件的概率.这样,根据大数定律,我们可以用频率来近似概率.具体做法如下:先产生2上均匀的随机数:.然后的次数,最后的估计值,即注:为了计算一般的定积分,我们可以这样来解决.设,其中有界,不妨设.令则有我们只要求即可.(2)平均值法设其中在[0,1]上连续,再设服从[0,1]上的均匀分布,则的数学期望为,所以估计的值就是估计的数学期望的值.由辛钦大数定律,可以用的观察值的平均值去估计的数学期望.具体做法如下:先产生上均匀的随机数:.然后对每个计算,最后的估计值,即实验所用软件及版本:MicrosoftofficeExcel2010主要内容(要点):(1)用随机投点法和平均值法计算定积分:(2)比较两种方法的精度.(n=100,1000,10000.对每一个n重复做5次),本次实验报告只做了n=1000的情况实验过程:(含解决方法和基本步骤,主要程序清单及异常情况记录等)随机投点法:估计定积分令当.(1)从主菜单中选择【工具|数据分析】,在【数据分析】对话框中滚动列表框,选择“随机数发生器”选项,如图9-1(2)重重两次可得1000对随机数,作变换,然后进行比较,如图所示,最后进行数数,得k=389,s=0.389具体如下所示.估计定积分下面我们估计,从上面不难得到S2=0.389,所以=3.028实验过程:(含解决方法和基本步骤,主要程序清单及异常情况记录等)平均值法:(1)从主菜单中选择【工具|数据分析】,在【数据分析】对话框中滚动列表框,选择“随机数发生器”选项,如图所示:估计定积分因为所以我们首先产生(-1,1)的均匀随机数.其它类似.如图所示估计定积分令.由上表可以算得:=0.401所以,J2=2*c2=0.741755.从上述结果可以看出,误差较大,必须加大容量.实验结果与实验总结(体会):通过本实验的两种方法的使用,将定积分用两种不同的方法计算出来,并得到不同的结果。该两种方法是基于不同的思想原理和步骤进行的。我们在学习和使用这两种方法的过程中必须明确它们各自的原理与不同点,也要求客观的评价和使用这两种方法。当然,作为数学软件的使用,它们都是各有千秋的。
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 急诊PCI手术并发症的护理预防与配合
- 关于销售绩效考核制度
- 合规审计规范财务核算制度
- 审计质量五级控制制度
- 员工手册绩效考核制度
- 四大停用国企审计制度
- 保安公司绩效考核制度
- 审计违反财务制度
- 医院二甲审计制度
- 审计环境宿舍管理制度
- 智慧路灯综合解决方案
- 《大学生心理健康》教案-自我意识课件
- 《春季健康饮食》课件
- 500字作文标准稿纸A4打印模板-直接打印
- 生物化学英文版课件:Chapter 6 Enzyme catalysis
- 23J916-1:住宅排气道(一)
- 慢性病健康管理规范
- 检验检测机构质量手册程序文件质量记录合集(依据2023年版评审准则)
- 冀教版(冀人版)科学六年级下册全册教案
- 国际贸易理论与实务习题答案汇总(王峰第三版)第1-16章+实务案例题
- GB/T 26121-2010可曲挠橡胶接头
评论
0/150
提交评论