蒙特卡罗(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

蒙特卡罗算法的matlab实现相关推荐

  1. 数模算法 — 蒙特卡罗算法及Matlab案例

    本文为博主编写文章,未经博主允许转载,转载请注明出处:https://blog.csdn.net/weixin_42604241/article/details/97382158 蒙特卡罗算法(方法) ...

  2. fcm算法的MATLAB实现,FCM算法的matlab程序(初步)

    FCM算法的matlab程序 1.采用iris数据库 iris_data.txt 5.1 3.5 1.4 0.2 4.9 3 1.4 0.2 4.7 3.2 1.3 0.2 4.6 3.1 1.5 0 ...

  3. 2018-4-8蚁群算法---包子阳《智能优化算法以及Matlab实现》第五章

    资料来源: <智能优化算法以及matlab实现>包子阳  余继周 编著 第五章-----蚁群算法 是一种元启发式优化算法(自己理解:就是作为群体的单位个体也就是元,在里面充当着随机的选择搜 ...

  4. matlab dfp法,DFP算法及Matlab程序.docx

    DFP算法及Matlab程序 作业二 用DFP算法求解,取,.一.求解:求迭代点x1令,得的极小值点,所以得:于是,由DFP修正公式有下一个搜索方向为求迭代点x2令,得的极小值点于是得:,所以:,因H ...

  5. matlab整定串级pid,PID算法在Matlab串级控制中的应用

    PID算法在Matlab串级控制中的应用 自114 1112002039 陈艳 前言:这个专题是由王娟老师给我们授课,我感觉收获挺大的,尤其是matlab仿真软件的使用,为我以后的实验课打下良好的基础 ...

  6. OpenMP 编程实例(蒙特卡罗算法)

    有关clock()函数 1,clock()函数在头文件#include<time.h>中 2,clock()函数的返回值类型为clock_t.clock_t其实是long,即长整形. cl ...

  7. PSO-LSSVM算法及其MATLAB代码

    挺完整的一篇博客,这里转载记录一下. 原文链接:PSO-LSSVM算法及其MATLAB代码 一.PSO 1.概念 粒子群优化算法(PSO:Particle swarm optimization)是一种 ...

  8. matlab hist函数_算法工匠MATLAB专训营:Matlab绘图,小试牛刀

    作者 | 蔡老师 仿真秀专栏作者 首发 | 仿真秀平台 导读:正文之前,我在此详细说明一下,因为本文包含的程序太难得,网上肯定找不到这样的程序.随着讲课的越来越深入,我给出的程序会越来越实用,接近于实 ...

  9. matlab音频基频的提取,(620512681) 自相关基频提取算法的MATLAB实现

    第31卷总第80期 西北民族大学学报(自然科学版) V01.31.No.4 1 0年1 2 0 2月 Journal of Nonhw铭t University for Nationalities(N ...

最新文章

  1. SharePoint服务器如果需要安装杀毒软件, 需要注意什么?
  2. ios 登录 java 后台,IOS苹果登录sign in with apple后端校验
  3. halcon/c++接口基础 之 析构函数和Halcon算子
  4. LeetCode 1014. 最佳观光组合
  5. python格式化输出区别_python格式化输出的区别
  6. 关于STL的一些笔记
  7. C++ string类的说明
  8. 预训练新范式!为什么Prompt会更有效?
  9. OSPF建立邻居、邻接关系 学习笔记
  10. synchronized关键字,Lock接口以及可重入锁ReentrantLock
  11. Java实现简易版的【图书管理系统】
  12. 2004-2021年数据库系统工程师软考中级题目及答案
  13. mdt 计算机名_配置 MDT 部署共享规则
  14. Python菜鸟快乐游戏编程_pygame
  15. 广东百望税控盘初始化设置
  16. 矩阵的 Jordan 标准型
  17. 深度学习 tensorflow 三维矩阵乘法(batch 迭代必须搞懂的矩阵乘法,维度增加)
  18. Matlab 中一些符号的含义
  19. Arduino 单位数码管实验
  20. WMS云端转型 势在必行

热门文章

  1. Java解P2678 [NOIP2015 提高组] 跳石头,有图有注释,通俗易懂
  2. LaTex引用中文论文
  3. 详解APK静态分析引擎`quark-engine`的5大功能
  4. LiveQing视频云平台服务-接入一个或多个LiveNVR实现云端管理及直播
  5. android设置左右声道音量,android控制左右声道实例
  6. 赫兹的单位换算_hz是什么单位(频率和赫兹的换算)
  7. 微信公众号开启服务器配置流程及注意事项
  8. matlab怎么画loess拟合,matlab smooth loess
  9. Unity----VR摄像机(浅谈)
  10. 天池比赛-02-用Apriori算法进行商品关联分析