蒙特卡洛方法 matlab 实现
2017-02-12 15:05
204 查看
蒙特卡罗(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是基于概率的随机模拟方法
转载本文请联系原作者获取授权,同时请注明本文来自杨乾栩科学网博客。
链接地址:http://blog.sciencenet.cn/blog-316653-375888.html
假设我们有个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是基于概率的随机模拟方法
转载本文请联系原作者获取授权,同时请注明本文来自杨乾栩科学网博客。
链接地址:http://blog.sciencenet.cn/blog-316653-375888.html
相关文章推荐
- 蒙特卡洛方法 matlab 实现
- matlab实现蒙特卡洛方法求解线性规划问题
- 最优化方法的Matlab实现
- 计算方法单点截弦matlab实现
- EM方法解高斯混合模型(GMM)Matlab实现
- MATLAB与C++混合编程:动态链接库方法实现混合编程及常见错误解决办法
- 二分类SVM方法Matlab实现
- C语言实现MATLAB 6.5中M文件的方法
- 【MATLAB】画图使用调色板实现自动线条颜色的方法
- VB与MATLAB接口实现方法
- matlab实现的手写体字符截取字符区域的方法
- 计算方法双点截弦法matlab实现
- 大津(Ostu)阈值分割方法在matlab中的实现
- 用牛顿方法解一元非线性方程的根(Matlab实现)
- 二分类SVM方法Matlab实现
- 计算方法中方程的近似解法中二分法matlab实现
- 不用递归的方法计算米粒的个数和每个米粒的面积,matlab实现
- 分水岭算法 matlab的三种实现方法
- VC 与MATLAB接口编程方法与实现(转)
- matlab实现复杂网络中节点重要度评估的节点收缩方法