原创力文档- 1、本文档共10页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
蒙特卡罗方法
随机模拟方法，也称为蒙特卡罗（MonteCarlo）方法，是一种基于“随机数”的计算方法．它以概率统计理论为基础，通过对研究的问题或系统进行随机抽样，然后对样本值进行统计分析，进而得到所研究问题或系统的某些具体参数、统计量等．在实际生活中，大量问题包含着随机性因素，这时进行随机模拟是一种有效的方式．随着模拟次数的增多，其精度也逐渐增高．由于需要大量反复的计算，一般用计算机来完成．蒙特卡罗方法在金融工程学，宏观经济学，计算物理学（如粒子输运计算、量子热力学计算、空气动力学计算）等领域应用广泛。
?图8-4正方形的内切圆
???解在一个面积为4的长方形内随机投点，该点落在此区域的概率为该区域的面积与长方形的面积比值，即p=S:4,可得S=4p．?
MATLAB命令如下：p=rand(10000,2);x=2*p(:,1)-1;y=2*p(:,2);N=find(y=2-x.^2y.^3=x.^2);M=length(N);S=4*M/10000%p=M/10000plot(x(N),y(N),b.)运行后在命令行窗口得到：S=2.1204可知区域面积为2.1204，这与精确解接近，图形如图8-5所示.图8-5蒙特卡罗投点区域
例8.14设某团队有50个人组成，试确定在一年中该团队至少2个人相同的概率．?
MATLAB命令如下：functionp=fun814(n)%n表示实验总次数s=0;fori=1:n%循环n次a=zeros(1,365);b=ceil(rand(1,50)*365);%随机生成1~365日期j=1;c=1;
whilej=50c==1k=b(j);a(k)=a(k)+1;%a(k)存放生日为b(j)的人数ifa(k)=2c=0;%找到有2人生日为同一天endj=j+1;endifc==0s=s+1;%统计2人生日相同的实验次数endendp=s/n;%计算概率
在命令行窗口调用p=fun814(10000)，得到p=0.9686．即团队有50个人组成，在一年中该团队至少2个人相同的概率是0.9686，这与理论解接近．
文档评论（0）