蒙特卡洛方法 matlab 实现

已有 12594 次阅读 2010-10-22 14:12 |个人分类:MATLAB|系统分类:科研笔记|关键词:蒙特卡洛方法 matlab

蒙特卡罗(Monte Carlo)方法,也称为计算机随机模拟方法,是一种基于"随机数"的计算方法。
假设我们有个y=x^2的表达式,如何用MC方法求得函数在[0,1]区间的定积分呢?
定积分可以用面积来求解,也就是通过求箭头下的面积
即:正方形面积为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,:))%x轴为生成的随机数,y轴为随机数的平方


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


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

本文引用地址: http://blog.sciencenet.cn/blog-316653-375888.html   转载请注明来自科学网博客,并请注明作者姓名。

蒙特卡洛方法 matlab 实现相关推荐

  1. 2021-04-09 随机模拟—蒙特卡洛方法 Matlab代码实现

    随机模拟-蒙特卡洛方法 Matlab代码实现 蒙特卡洛方法 蒙特卡洛方法(Monte Carlo method),也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出 ...

  2. 欧式期权matlab编码蒙特卡罗方法,基于MATLAB的金融工程方法与实践第七章 基于蒙特卡洛方法的期权定价.ppt...

    蒙特卡洛的优缺点 MCMC方法的优点 分布假设更一般,描述市场因素可能变化的统计分布既可以是正态.对数正态的,也可以是带跳的扩散分布.t分布等. 随机生成风险因素的各种各样的未来假想情景,可在模型中融 ...

  3. Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动

    Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...

  4. matlab mcl,wsn定位蒙特卡洛方法mcl的matlab

    wsn定位蒙特卡洛方法mcl的matlab clear;clc;%初始化工作Ns = 20;Nn = 200;Vmax = 20;Xrange = 200;Yrange = 200;tr = 50;s ...

  5. 蒙特卡洛方法生成的电动汽车充电需求曲线 matlab平台建模

    蒙特卡洛方法生成的电动汽车充电需求曲线 matlab平台建模 ID:6919672622070241

  6. matlab 蒙特卡洛工具箱_量化投资_轻松实现MATLAB蒙特卡洛方法建模

    1 目录 * MATLAB随机数的产生 - Uniform,Normal & Custom distributions * 蒙特卡洛仿真 * 产生股票价格路径 * 期权定价 - 经典公式 - ...

  7. wsn定位matlab仿真,WSN定位蒙特卡洛方法MCL的MATLAB实现源码

    WSN定位蒙特卡洛方法MCL的MATLAB实现源码 clear; clc; %初始化工作 Ns = 20; Nn = 200; Vmax = 20; Xrange = 200; Yrange = 20 ...

  8. matlab蒙特卡洛方法求积分,matlab-蒙特卡洛法估计积分值

    <matlab-蒙特卡洛法估计积分值>由会员分享,可在线阅读,更多相关<matlab-蒙特卡洛法估计积分值(6页珍藏版)>请在人人文库网上搜索. 1.西安交通大学实验报告课程: ...

  9. 在我方某前沿防守地域 matlab,蒙特卡洛方法模拟小例子

    例 在我方某前沿防守地域,敌人以一个炮排(含两门火炮)为单位对我方进行干扰和破坏.为躲避我方打击,敌方对其阵地进行了伪装并经常变换射击地点. 经过长期观察发现,我方指挥所对敌方目标的指示有50%是准确 ...

最新文章

  1. python hex 补0_Python保留前导零的二进制到十六进制转换
  2. Oracle 11g Dataguard 物理备库配置(四)之broker snapshot standby测试
  3. XGB模型训练报错 terminate called after throwing an instance of ‘std::bad_alloc‘ what()
  4. 培智计算机教学论文,【培智数学论文】_培智数学教学论文
  5. C语言一趟冒泡交换最小值,C语言单链表冒泡排序为啥以下代码实现不了?
  6. python+selenium之自动生成excle,保存到指定的目录下
  7. Jupyter notebook内使用(添加)虚拟环境
  8. VC实现文件拖拽获取文件名
  9. Ubuntu 命令技巧
  10. 手把手教你搭建Java 开发环境
  11. SecureCRT下载安装与破解
  12. php导入json文件_php读取json文件的实例(本地与远程),笔记
  13. 基于JavaWeb的会议室预约管理系统
  14. 【Tableau server日常维护9.0】window server2016 Tableau server保姆教学安装
  15. 启动“附近的人”功能,你有兴趣吗?
  16. 农夫山泉市值超5000亿港元 域名还在优化中?
  17. Unity3D Behavior Designer 行为树1 认识
  18. 2021 年最佳开源软件榜单,Python 高性能框架 FastAPI 上榜了!
  19. 在内嵌RAM中运行程序 No Algorithm found for 错误
  20. 【Linux】【CentOS】CentOS 网络管理

热门文章

  1. servlet容器,web容器,spring容器,springmvc容器的区别
  2. iphone阅读html,iphone下Safari查看网页源代码
  3. 2021年4月19日 深圳头条后台开发实习面试(一面)(含总结)
  4. 大学生云计算期末考试题
  5. 区块链项目ICO评估模型
  6. 【如何开发小程序】自己如何开发小程序?
  7. MinIO 文件服务器简单搭建
  8. Dotween简单安装设置
  9. LVGL-gui_user.c
  10. 同步、异步 阻塞、非阻塞