应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 【转帖】蒙特卡罗算法计算圆周率
132/2返回列表上一页12
查看: 3160  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 模拟EPI ，点击这里进行查看

月3399

银虫 (初入文坛)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3):谢谢分享 2010-11-09 14:24:43
一个简单的MATLAB程序：
function pai=pai(N)
   s=0;for i=1:N,
       x=2*rand(2,1)-1;
       if ((x(1))^2+(x(2))^2)<=1,s=s+1;
           plot(x(1),x(2)),hold on;
       end
   end
   pai=4*s/N;
   pai
赞一下(2人) 回复此楼
11楼2010-11-09 14:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobulu420

新虫 (初入文坛)
好东东，谢楼主了
回复此楼
12楼2010-11-11 10:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Samuel_wang

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
161408楼: Originally posted by zyj8119 at 2010-09-11 10:25:56:
这是一个关于如何求pi的matlab小程序，用到了Monte Carlo的思想，我个人觉得所谓MC方法就是用概率的思想去解决问题，主要就是靠计算机生成的伪随机数来对一个随机事件多次仿真得到所求解的近似值。
下面是个改进

请问：在上述代码中，怎么体现循环？好像没有地方用到nloops啊，我看得不是太明白~~~求指教~~~
另外，我若想对所求的pi进行误差分析，程序该怎么改呢？
赞一下 回复此楼
13楼2012-05-22 18:55:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyj8119 的主题更新
132/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭