蒙特卡洛(MC)方法的原理（续）

一个实例

用概率的方法计算y=x^2在[0,1]区间上的积分值。

 

x=0:0.01:1;
y=x.^2;plot(x,y);


Num =100000; % n的数目


RandData=rand(2,Num); %根据随机点数，产生随机的(x,y)散点,不明白可以试试
Below=find(RandData(1,:).^2>RandData(2,:));%寻找位于曲线下的散点
Result = length(Below)/length(RandData);%最终结果的表示


Result
% How to display the final results

 

Matlab显示的结果为：

 

                        Result = 0.3335

 

如果您跑这段代码的话，最后的结果可能会不同！这也正是该方法的一个特点，她的结果误差是随机，但是只要我们合理的控制条件，这种方法还是

很不错的哦!!!

 

 

总结

感觉自己写的这篇感悟实在是太水了，但毕竟是自己的理解，恳请各位网友批评指正，我以后也会多多努力的！欢迎大家和我交流。

我的邮箱是1184611233@qq.com