秋 实验6蒙特卡洛模拟monto carlo simulation精品_第1页
秋 实验6蒙特卡洛模拟monto carlo simulation精品_第2页
秋 实验6蒙特卡洛模拟monto carlo simulation精品_第3页
秋 实验6蒙特卡洛模拟monto carlo simulation精品_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、1UsingMonteCarlosimulation,writeaMatlabprogramtocalculateanapproximationto冗byconsideringthenumberofrandompointsselectedinsidethequartercircleQ:x2+y2=1,x0,y0wherethequartercircleistakentobeinsidethesquareS:0 x1,0y1解：counter=0;N=100000;fori=1:Nx=rand(1,1);y=rand(1,1);ifxA2+yA21counter=counter+1;endend

2、p=4*counter/N结果：p=3.14482.UsingMonteCarlosimulation,writeaMatlabprogramtocalculatethevolumeofanellipsoid(椭球体)解：counter=0;N=1000;fori=0:Nx=-4*sqrt(2)+8*sqrt(2)*rand(1,1);y=-8+16*rand(1,1);z=-8*sqrt(2)+16*sqrt(2)*rand(1,1);ifxA2/2+yA2/4+zA2/8=0&d=0.1&d=0.2&d=0.4&d=0.7&d=0.5&p(j+1)=0.5&p(j+2)=0.5t=1;en

3、dendift=1counter=counter+1;endendchance=counter/N结果：chance=0.36765.(Optional)UsingMonteCarlosimulation,writeaMatlabprogramtocalculatethevolumetrappedbetweenthetwoparaboloids(抛物体)z=8一x2一y2andz=x2+3y2Notethatthetwoparaboloidsintersect(相交)ontheellipticcylinderx2+2y2=4解：x=-10:0.1:10;y=x;X,Y=meshgrid(x,y);Z=8-X.A2-Y.A2;mesh(X,Y,Z)holdonz=X.A2+3*Y.A2;mesh(X,Y,z)ezplot(X.A2+2*Y2-4)holdoffcounter=0;N=10000;fori=1:Nx=-2+4*rand(1,1);0-5-5y=-sqrt(2)+2*sqrt(2)*rand(1,1);z=-8+16*rand(1,1);ifxA2+2*yA28-xA2-yA2&zx、2+3T、2counter=counter+1;endendcount

人人文库> 全部分类> 图纸下载 > 毕业设计

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论