蒙特卡罗算法的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实现相关推荐
- 数模算法 — 蒙特卡罗算法及Matlab案例
本文为博主编写文章,未经博主允许转载,转载请注明出处:https://blog.csdn.net/weixin_42604241/article/details/97382158 蒙特卡罗算法(方法) ...
- 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 ...
- 2018-4-8蚁群算法---包子阳《智能优化算法以及Matlab实现》第五章
资料来源: <智能优化算法以及matlab实现>包子阳 余继周 编著 第五章-----蚁群算法 是一种元启发式优化算法(自己理解:就是作为群体的单位个体也就是元,在里面充当着随机的选择搜 ...
- matlab dfp法,DFP算法及Matlab程序.docx
DFP算法及Matlab程序 作业二 用DFP算法求解,取,.一.求解:求迭代点x1令,得的极小值点,所以得:于是,由DFP修正公式有下一个搜索方向为求迭代点x2令,得的极小值点于是得:,所以:,因H ...
- matlab整定串级pid,PID算法在Matlab串级控制中的应用
PID算法在Matlab串级控制中的应用 自114 1112002039 陈艳 前言:这个专题是由王娟老师给我们授课,我感觉收获挺大的,尤其是matlab仿真软件的使用,为我以后的实验课打下良好的基础 ...
- OpenMP 编程实例(蒙特卡罗算法)
有关clock()函数 1,clock()函数在头文件#include<time.h>中 2,clock()函数的返回值类型为clock_t.clock_t其实是long,即长整形. cl ...
- PSO-LSSVM算法及其MATLAB代码
挺完整的一篇博客,这里转载记录一下. 原文链接:PSO-LSSVM算法及其MATLAB代码 一.PSO 1.概念 粒子群优化算法(PSO:Particle swarm optimization)是一种 ...
- matlab hist函数_算法工匠MATLAB专训营:Matlab绘图,小试牛刀
作者 | 蔡老师 仿真秀专栏作者 首发 | 仿真秀平台 导读:正文之前,我在此详细说明一下,因为本文包含的程序太难得,网上肯定找不到这样的程序.随着讲课的越来越深入,我给出的程序会越来越实用,接近于实 ...
- matlab音频基频的提取,(620512681) 自相关基频提取算法的MATLAB实现
第31卷总第80期 西北民族大学学报(自然科学版) V01.31.No.4 1 0年1 2 0 2月 Journal of Nonhw铭t University for Nationalities(N ...
最新文章
- SharePoint服务器如果需要安装杀毒软件, 需要注意什么?
- ios 登录 java 后台,IOS苹果登录sign in with apple后端校验
- halcon/c++接口基础 之 析构函数和Halcon算子
- LeetCode 1014. 最佳观光组合
- python格式化输出区别_python格式化输出的区别
- 关于STL的一些笔记
- C++ string类的说明
- 预训练新范式!为什么Prompt会更有效?
- OSPF建立邻居、邻接关系 学习笔记
- synchronized关键字,Lock接口以及可重入锁ReentrantLock
- Java实现简易版的【图书管理系统】
- 2004-2021年数据库系统工程师软考中级题目及答案
- mdt 计算机名_配置 MDT 部署共享规则
- Python菜鸟快乐游戏编程_pygame
- 广东百望税控盘初始化设置
- 矩阵的 Jordan 标准型
- 深度学习 tensorflow 三维矩阵乘法(batch 迭代必须搞懂的矩阵乘法,维度增加)
- Matlab 中一些符号的含义
- Arduino 单位数码管实验
- WMS云端转型 势在必行
热门文章
- Java解P2678 [NOIP2015 提高组] 跳石头,有图有注释,通俗易懂
- LaTex引用中文论文
- 详解APK静态分析引擎`quark-engine`的5大功能
- LiveQing视频云平台服务-接入一个或多个LiveNVR实现云端管理及直播
- android设置左右声道音量,android控制左右声道实例
- 赫兹的单位换算_hz是什么单位(频率和赫兹的换算)
- 微信公众号开启服务器配置流程及注意事项
- matlab怎么画loess拟合,matlab smooth loess
- Unity----VR摄像机(浅谈)
- 天池比赛-02-用Apriori算法进行商品关联分析