蒙特卡洛随机模拟的MATLAB实例解析纪录
蒙特卡罗(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 随机模拟型综合评价概述 2 基于随机模拟的综合评价步骤 2.1 问题描述 2.2 自主优势量矩阵计算 2.3 具体步骤 2.4 综合评价计算 3 Python实现 1 随机模拟型综合评价 ...
- 2021-04-09 随机模拟—蒙特卡洛方法 Matlab代码实现
随机模拟-蒙特卡洛方法 Matlab代码实现 蒙特卡洛方法 蒙特卡洛方法(Monte Carlo method),也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出 ...
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...
- matlab优化设计例程,优化设计Matlab_实例解析
<优化设计Matlab_实例解析>由会员分享,可在线阅读,更多相关<优化设计Matlab_实例解析(9页珍藏版)>请在人人文库网上搜索. 1.1,单级直齿圆柱齿轮减速器的优化设 ...
- matlab数值模拟的例子,第12章 MATLAB 数值模拟实例解析课件
第12章 MATLAB 数值模拟实例解析课件 采用比母材金属熔点低的金属材料作钎料,将焊件和钎料加热到高于钎料熔点.低于母材溶化温度,利用液态钎料润湿母材 采用比母材金属熔点低的金属材料作钎料,将焊件 ...
- matlab 模拟心电图,使用原生js+canvas实现模拟心电图的实例
从2015年2月转行进入IT行业,到现在也有将近两年的时间了,从最开始的java到现在的前端,前进的路上一直靠自己摸索,一路走到现在,前端大神是绝对谈不上的,最多算一只刚入门的菜鸟. 从最开始的懵懵懂 ...
- 随机网络图 matlab,MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity)...
MATLAB实例:构造网络连接图(Network Connection)及计算图的代数连通度(Algebraic Connectivity) 1. 图的代数连通度(Algebraic Connecti ...
- 基于拟蒙特卡洛模拟法的随机潮流计算matlab程序
电力系统随机潮流计算中常采用模拟法,该方法原理简单.使用方便,能够精确地模拟实际物理过程,但是简单的蒙特卡洛模拟法收敛速度很慢,要得到精确的结果需要以大量的计算时间为代价.本章在此基础上提出了基于拟蒙 ...
- 基于matlab的储能选址定容程序 采用蒙特卡洛随机算法,非粒子群算法
基于matlab的储能选址定容程序 采用蒙特卡洛随机算法,非粒子群算法 图中每个点代表一种配置方式,红点为最优方式 程序稳定运行,每条语句均有详细注释,0基础可看懂 有对应文献 ID:27120675 ...
最新文章
- win8下notepad++无法设置文件关联
- 搞怪的html代码,这个恶搞网页代码是肿么写出来的?
- halcon 圆环类缺陷检测的一种方法(极坐标变换法)
- 防止DISCUZ根域名跳转到forum.php的方法
- MySQL不走联合索引_mysql group by 多列优化思路?为什么不走联合索引?
- oracle Instance status: READY–lsnrctl status|start|stop
- 漫步数理统计二十一——变换:随机向量
- Bootstrap3 Font Awesome 字体图标
- [转载]SIFT(尺度不变特征变换)算法小结
- vconsole在手机真机模拟器上显示console控制台
- quartus仿真28:JK触发器实现的脉冲分配器(分析)
- Linux 中 ext、ext2、ext3、ext4 文件系统介绍与区别
- python抽奖小程序_python实现简单的抽奖小程序,抽奖的内容从文件里面读取
- 计算机内存不足图片,电脑无法显示图片说内存不足
- mybatis 插入insert对象
- bom event周期_BOM-event事件
- linux系统文件制作编译器,制作mipsel-linux交叉编译工具
- 计算机桌面无法解锁,电脑win10系统锁屏后解锁却无法进入桌面的解决方法
- 旧时王谢堂前燕,飞入寻常百姓家
- 如何在uni-app中使用vant组件?最细教学,还看不懂你来捶我
热门文章
- 金山现任CEO张宏江将退休 西山居CEO继任
- win10查看计算机系统版本,Win10版本怎么看?Win10系统版本查看方法
- doris数据库环境搭建报错(invalid cluster id. ignore)
- 大学生英语竞赛大揭秘准备攻略
- WebGIS第九课:智慧校园项目(1)
- DS-5 csat简单使用-dump memory
- Vmware VDI\桌面虚拟化\虚拟化技术\IT 管理
- 多少秒算长镜头_电影中什么是长镜头画面(9个经典长镜头)
- 新手小白如何用linux云服务器搭建wordpress个人网站
- Revit建模:楼板三维视图旋转模型如何使图案跟着旋转