%%monte_carlo_ff.m

%被积函数(二重)

function ff=monte_carlo_ff(x,y)

ff=x*y^2;%函数定义处

end

%%monte_carlo.m

%蒙特卡洛计算二重积分

function result=monte_carlo(a,b,c,d,n,m)

%先y后x积分,a是x积分下限,b是x积分上限,c是y积分下限,d是y积分上限,n,m是蒙特卡洛参数

sumxff=0;

for i=1:n

sumyff=0;

xff=a+(b-a)*rand();

for j=1:m

yff=c+(d-c)*rand();

sumyff=sumyff+monte_carlo_ff(xff,yff);

end

aversumyff=sumyff/m;

sumxff=sumxff+(b-a)*aversumyff;

end

result=sumxff/n;

end

%%show.m

clear all;

clc;

format long;

%二重积分的int函数法,用于校对蒙特卡洛法

syms x y;

z=x^2*y^2;

x1=1;

x2=2;

y1=3;

y2=4;

disp(‘int法结果‘)

ans_int=int(int(z,y,y1,y2),x,x1,x2)

%蒙特卡洛法计算结果

disp(‘monte_carlo法结果‘)

ans_monte_carlo=monte_carlo(1,2,3,4,10000,20000)

%二者误差

disp(‘误差‘)

error=abs(ans_monte_carlo-ans_int)

原文:http://www.cnblogs.com/blog-3123958139/p/6128149.html

蒙特卡罗计算积分matlab,matlab下二重积分的蒙特卡洛算法相关推荐

  1. matlab腔内光子寿命,mcFORnp matlab环境下,利用蒙特卡洛模拟光子包在生物组织内的光路传输 271万源代码下载- www.pudn.com...

    文件名称: mcFORnp下载  收藏√  [ 5  4  3  2  1 ] 开发工具: matlab 文件大小: 215 KB 上传时间: 2014-12-29 下载次数: 8 提 供 者: 徐某 ...

  2. matlab碎纸拼接相似函数,基于蒙特卡洛算法构建能量函数的碎纸图片拼接方法

    基于蒙特卡洛算法构建能量函数的碎纸图片拼接方法 [专利摘要]本发明提供了一种基于蒙特卡洛算法构建能量函数的碎纸图片拼接方法,主要涉及双面打印文件的拼接及复原问题,通常由于图片较多,信息量较大,故通常为 ...

  3. python算积分蒙特卡罗_蒙特卡罗计算积分

    作者|Cory Maklin 编译|VK 来源|Towards Datas Science 通常情况下,我们不能解析地求解积分,必须借助其他方法,其中就包括蒙特卡罗积分.你可能还记得,函数的积分可以解 ...

  4. C++ 偏微分数值计算库_蒙特卡罗计算积分 - 人工智能遇见磐创

    作者|Cory Maklin 编译|VK 来源|Towards Datas Science 通常情况下,我们不能解析地求解积分,必须借助其他方法,其中就包括蒙特卡罗积分.你可能还记得,函数的积分可以解 ...

  5. matlab实验求不定积分函数,实验四用matlab计算积分

    实验四用matlab计算积分 实验四 用matlab计算积分 4.1积分的有关理论 定积分:积分是微分的无限和,函数在区间上的积分定义为 其中 从几何意义上说,对于上非负函数,记分值是曲线与直线及轴所 ...

  6. 在matlab环境下实现支持向量机算法

    作者:faaronzheng 转载请注明出处! 关于支持向量机的理论部分我就不进行介绍了,网上有很多关于支持向量机的讲解.如: 支持向量机通俗导论(理解SVM的三层境界)有兴趣的可以自行查找一下.这里 ...

  7. matlab 裕度指令,稳定裕度的计算开环传递函数Matlab命令num.PPT

    稳定裕度的计算开环传递函数Matlab命令num 4.1 s平面和z平面之间的映射 4.2 稳定性分析 4.3 稳态误差分析 4.4 时域特性分析 4.5 频域特性分析 4.6 应用实例 4.1.1 ...

  8. matlab trapz求二重积分,matlab求积分(超详细,含int integral integral2/3 quad trapz

    matlab求积分 matlab求积分函数工具: int 用法1: 格式: int(fun,x,a,b) 功能: 计算定积分 用法2: 格式: int(f,x) 功能: 计算不定积分 注: 使用int ...

  9. 静电场的有限差分法与matlab 仿真课程设计,计算物理和MATLAB课程设计--自激振动系统的MATLAB仿真.doc...

    东北石油大学课程设计任务书 课程 计算物理和MATLAB课程设计 题目 自激振动系统的MATLAB仿真 专业 姓名 学号 主要内容.基本要求.主要参考资料等 主要内容: 研究范?德?波耳(Van de ...

  10. matlab 图像矢量量化,MATLAB环境下基于矢量量化的说话人识别系统(1)

    第21卷第6期湖 北 工 业 大 学 学 报2006年12月 Vol.21No.6 Journal of Hubei Univer sity of Technology Dec.2006 [收稿日期] ...

最新文章

  1. 【UML建模】机房中的UML图
  2. 视觉稿与H5页面之间的终端适配
  3. 对于PS和flash协作做过光字的一点心得
  4. 亮道智能发布自动驾驶环境感知系统测试验证服务|2019 上海车展...
  5. 专题 20 常用调试技术
  6. LeetCode Surrounded Regions(floodfill算法)
  7. Elasticsearch深入11 运维
  8. unity3d游戏数据加密
  9. IDEA实现到JSTL
  10. python要配置几个环境变量_Python环境变量配置
  11. 高速缓冲存储器的功能、结构与工作原理
  12. .NET应用迁移到.NET Core--调查案例
  13. java上三角数组_二维数组的三角填充 两种java实现的方法
  14. (27)css3 2D转换综合案例—卡包效果
  15. Oracle中报ora-01033
  16. Springboot URL Rewrite
  17. 遗传算法详解与MATLAB实现
  18. MAVEN 仓库加速 阿里云镜像
  19. R统计绘图-VPA(方差分解分析)
  20. OpenStack Days China Day 2:厨子与食客的热恋

热门文章

  1. 阿里巴巴国际站全屏分类产品切换效果怎么做?代码生成器产品分类效果切换轮播全屏代码
  2. 计算机硬件系统基本的工作原理是,硬件系统基本组成及工作原理
  3. java取万位的值,Excel中表格数值进行取万位整数的操作方法
  4. 智能商场管理系统(可当做毕设)
  5. 【MySQL学习笔记(十六)】之redo日志超详细讲解
  6. 蒟蒻的第一次豆瓣读书爬虫实践
  7. Oracle 繁简体转换(未解决)
  8. 如何处理计算机显示器故障,电脑显示器常见的故障及解决办法
  9. MD5的认识,建议所有菜菜都看下
  10. matlab设计计算器,MATLAB-计算器的设计.doc