• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

蒙特卡洛方法 matlab 实现

原作者: [db:作者] 来自: [db:来源] 收藏 邀请
蒙特卡罗(Monte Carlo)方法,也称为计算机随机模拟方法,是一种基于"随机数"的计算方法。
假设我们有个y=x^2的表达式,如何用MC方法求得函数在[0,1]区间的定积分呢?
定积分可以用面积来求解,也就是通过求箭头下的面积
为了衔接方便,照顾新手,给出作图程序吧
x=0:0.01:1;y=x.^2;plot(x,y);


MC方法实现非常简单,通过下面的代码就可以
staus=10;
for i=1:4  %4次模拟
point=staus.^i; %模拟的随机点数
RandData=rand(2,point); %根据随机点数,产生随机的(x,y)散点,不明白可以试试   %scatter(RandData(1,:),RandData(2,:))
Below=find(RandData(1,:).^2>RandData(2,:));%寻找位于曲线下的散点
Outcome(i)=length(Below)/length(RandData);%最终结果的表示
end
Outcome =

    0.3000    0.3600    0.3180    0.3311
从Outcome看,通过不断增加随即点数,结果越与真实值相符
当散点数为10^4时,所得图见下
BelowData=RandData(:,Below);
hold on
scatter(BelowData(1,:),BelowData(2,:))


如果我们选取的散点数为10^5,则定积分值为0.3335,所得图形见下(代码略,同上)


以上是对MC方法最简单的理解,不过思想上是融会贯通的,适合新手学习。所以很明确,MC是基于概率的随机模拟方法 。

鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
Delphi中的published发布时间:2022-07-18
下一篇:
DELPHI开发LINUX包发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap