[原创]Faruto老师蒙特卡罗方法投点图程序的优化

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

    李洋老师书上那个代码是在循环里面画图，这样没循环一次就得调用一次plot函数，大大减慢了运行速度，改进之前，运行时间19.348197秒。这个时间是很长的，可以稍微更改代码如下，运行时间0.264268 秒。（刚学MATLAB，也只能发些这些技术含量不高的帖子，多赚点经验，各版友见谅我水帖啊）
        代码如下：

1. clear
   
2. clc
   
3. TestNum=1e4;
   
4. %生成[-1,1]*[-1,1]上服从均匀分布的随机数,rand为[0,1]上均匀分布，
   
5. %2*(a-0.5）的方式将起改变为[-1,1]上的均匀分布
   
6. X=2*(rand(TestNum,2)-0.5);
   
7. CircleNum=0;
   
8. figure
   
9. hold on
   
10. circle=[];
    
11. other=[];
    
12. otherNum=0;
    
13. for i=1:TestNum
    
14. if X(i,1)^2+X(i,2)^2<=1
    
15. CircleNum=CircleNum+1;
    
16. circle(CircleNum,:)=X(i,:);
    
17. else
    
18. otherNum=otherNum+1;
    
19. other(otherNum,:)=X(i,:);
    
20. end
    
21. end
    
22. plot(circle(:,1),circle(:,2),'r.');
    
23. plot(other(:,1),other(:,2),'bo');
    
24. %计算（圆周率）Pi值
    
25. SPi=4*CircleNum/TestNum;

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享2 收藏2 回帖

关键词：蒙特卡罗方法 蒙特卡罗 Far 蒙特卡 circle 蒙特卡罗 程序 李洋

厉害支持一下

加油加油

可以把你改进的东西慢慢积累，整理成一个大帖子 我好加精华

fantuanxiaot 发表于 2015-1-19 18:49
可以把你改进的东西慢慢积累，整理成一个大帖子 我好加精华

好的

学习了。

不错