蒙特卡洛方法 matlab 实现
蒙特卡洛方法 matlab 实现
已有 12594 次阅读 2010-10-22 14:12 |个人分类:MATLAB|系统分类:科研笔记|关键词:蒙特卡洛方法 matlab
假设我们有个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,:))%x轴为生成的随机数,y轴为随机数的平方
如果我们选取的散点数为10^5,则定积分值为0.3335,所得图形见下(代码略,同上)
以上是对MC方法最简单的理解,不过思想上是融会贯通的,适合新手学习。所以很明确,MC是基于概率的 随机模拟方法
本文引用地址: http://blog.sciencenet.cn/blog-316653-375888.html 转载请注明来自科学网博客,并请注明作者姓名。
蒙特卡洛方法 matlab 实现相关推荐
- 2021-04-09 随机模拟—蒙特卡洛方法 Matlab代码实现
随机模拟-蒙特卡洛方法 Matlab代码实现 蒙特卡洛方法 蒙特卡洛方法(Monte Carlo method),也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出 ...
- 欧式期权matlab编码蒙特卡罗方法,基于MATLAB的金融工程方法与实践第七章 基于蒙特卡洛方法的期权定价.ppt...
蒙特卡洛的优缺点 MCMC方法的优点 分布假设更一般,描述市场因素可能变化的统计分布既可以是正态.对数正态的,也可以是带跳的扩散分布.t分布等. 随机生成风险因素的各种各样的未来假想情景,可在模型中融 ...
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...
- matlab mcl,wsn定位蒙特卡洛方法mcl的matlab
wsn定位蒙特卡洛方法mcl的matlab clear;clc;%初始化工作Ns = 20;Nn = 200;Vmax = 20;Xrange = 200;Yrange = 200;tr = 50;s ...
- 蒙特卡洛方法生成的电动汽车充电需求曲线 matlab平台建模
蒙特卡洛方法生成的电动汽车充电需求曲线 matlab平台建模 ID:6919672622070241
- matlab 蒙特卡洛工具箱_量化投资_轻松实现MATLAB蒙特卡洛方法建模
1 目录 * MATLAB随机数的产生 - Uniform,Normal & Custom distributions * 蒙特卡洛仿真 * 产生股票价格路径 * 期权定价 - 经典公式 - ...
- wsn定位matlab仿真,WSN定位蒙特卡洛方法MCL的MATLAB实现源码
WSN定位蒙特卡洛方法MCL的MATLAB实现源码 clear; clc; %初始化工作 Ns = 20; Nn = 200; Vmax = 20; Xrange = 200; Yrange = 20 ...
- matlab蒙特卡洛方法求积分,matlab-蒙特卡洛法估计积分值
<matlab-蒙特卡洛法估计积分值>由会员分享,可在线阅读,更多相关<matlab-蒙特卡洛法估计积分值(6页珍藏版)>请在人人文库网上搜索. 1.西安交通大学实验报告课程: ...
- 在我方某前沿防守地域 matlab,蒙特卡洛方法模拟小例子
例 在我方某前沿防守地域,敌人以一个炮排(含两门火炮)为单位对我方进行干扰和破坏.为躲避我方打击,敌方对其阵地进行了伪装并经常变换射击地点. 经过长期观察发现,我方指挥所对敌方目标的指示有50%是准确 ...
最新文章
- python hex 补0_Python保留前导零的二进制到十六进制转换
- Oracle 11g Dataguard 物理备库配置(四)之broker snapshot standby测试
- XGB模型训练报错 terminate called after throwing an instance of ‘std::bad_alloc‘ what()
- 培智计算机教学论文,【培智数学论文】_培智数学教学论文
- C语言一趟冒泡交换最小值,C语言单链表冒泡排序为啥以下代码实现不了?
- python+selenium之自动生成excle,保存到指定的目录下
- Jupyter notebook内使用(添加)虚拟环境
- VC实现文件拖拽获取文件名
- Ubuntu 命令技巧
- 手把手教你搭建Java 开发环境
- SecureCRT下载安装与破解
- php导入json文件_php读取json文件的实例(本地与远程),笔记
- 基于JavaWeb的会议室预约管理系统
- 【Tableau server日常维护9.0】window server2016 Tableau server保姆教学安装
- 启动“附近的人”功能,你有兴趣吗?
- 农夫山泉市值超5000亿港元 域名还在优化中?
- Unity3D Behavior Designer 行为树1 认识
- 2021 年最佳开源软件榜单,Python 高性能框架 FastAPI 上榜了!
- 在内嵌RAM中运行程序 No Algorithm found for 错误
- 【Linux】【CentOS】CentOS 网络管理