用蒙特卡罗积分法近似求解定积分的值及举例

在这里插入图片描述

一、背景知识

1、连续随机变量的概率密度函数

对于连续型随机变量的概率密度函数(PDF),其在整个定义域上的积分必须等于1。这是概率密度函数的一个基本属性,它确保了随机变量取任何值的概率之和等于1,符合概率论的公理。

也即:如果x是一个连续型随机变量,其概率密度函数表示为f(x),则必须满足以下条件:

∫ l u f ( x )   d x = 1 ,其中 l 、 u 表示积分的上下界限 \int_{l}^{u} f(x) \, dx = 1,其中l、u表示积分的上下界限 luf(x)dx=1,其中lu表示积分的上下界限

这个积分表示了随机变量x取值在l到u之间所有可能值的概率之和必须等于1。

2、连续随机变量的均匀分布

均匀分布,又称矩形分布,是一种连续概率分布,其概率密度函数(Probability Density Function, PDF)在某个区间内是恒定的,而在该区间外是零。均匀分布的特点是随机变量在定义域内任何子区间上取值的概率是相同的。

对于连续型均匀分布,其数学表达式如下:

f ( x ) = { 1 b − a for  a ≤ x ≤ b , 0 otherwise . f(x) = \begin{cases} \frac{1}{b-a} & \text{for } a \le x \le b, \\ 0 & \text{otherwise}. \end{cases} f(x)={ba10for axb,otherwise.

其中,a 是分布的下限,b 是分布的上限,b > a。

对于连续型均匀分布,期望值(均值)E[X]是区间中点: E [ X ] = a + b 2 E[X] = \frac{a + b}{2} E[X]=2a+b

在这里插入图片描述

二、蒙特卡罗积分法

蒙特卡罗积分法是一种基于随机数的数值积分方法,核心思想是利用随机变量的期望值来近似积分的值,它通过随机抽样生成大量的样本点的方式来近似计算定积分的值,这是因为积分值服从大数定律。

具体来说,对于被积函数h(x),如果h(x)=f(x)p(x),其中p(x)为一个在积分区间的概率密度函数,即满足在积分区间的积分值为1,即:

∫ a b h ( x )   d x = ∫ a b f ( x ) p ( x )   d x ,其中 ∫ a b p ( x )   d x = 1 \int_{a}^{b} h(x) \, dx=\int_{a}^{b} f(x) p(x)\, dx,其中\int_{a}^{b} p(x)\, dx=1 abh(x)dx=abf(x)p(x)dx,其中abp(x)dx=1

则:
∫ a b h ( x )   d x ≈ 1 n ∑ i = 1 n f ( x i ) ,其中 n 为抽样的样本点数量, f ( x ) = h ( x ) / p ( x ) \int_{a}^{b} h(x) \, dx≈\frac{1}{n} \sum_{i=1}^{n}f(x_i),其中n为抽样的样本点数量 ,f(x)=h(x)/p(x) abh(x)dxn1i=1nf(xi),其中n为抽样的样本点数量,f(x)=h(x)/p(x)

蒙特卡罗方法的误差主要取决于样本点的数量。一般来说,随着样本点数量的增加,误差会减小。
在实际应用中,蒙特卡洛积分法的一个关键步骤是生成随机数。这些随机数通常是均匀分布的,但也可以是其他类型的分布,如正态分布、泊松分布等,主要看被积函数h(x)怎么分解成f(x)p(x),使得p(x)在被积区间的积分值为1,如果p(x)是均匀分布则p(x)=1/(b-a),f(x)=(b-a)h(x)。

三、用蒙特卡罗求定积分求 ∫ 1 10 x d x \int_{1}^{10}xdx 110xdx案例

采用p(x)为均匀分布,则p(x)=1/10,f(x)=10x,取n=11, x i x_i xi=1,2,…,10,针对11个采样点,计算10x的值并求和:
在这里插入图片描述
则:
∫ 1 10 x   d x ≈ 1 11 ∑ i = 1 11 10 x i = 550 / 11 = 50 \int_{1}^{10} x\, dx≈\frac{1}{11} \sum_{i=1}^{11}10x_i=550/11=50 110xdx111i=11110xi=550/11=50

如果用正规的积分计算方式,x的积分为x²/2,则 ∫ 1 10 x d x = ( x 2 / 2 ) ∣ 0 10 = 50 \int_{1}^{10}xdx=(x²/2)|^{10}_0=50 110xdx=(x2/2)010=50,可以看到这2种方式计算的积分完全相同。

四、小结

本文介绍了蒙特卡罗积分法近似求解定积分的方法,并举例说明了具体使用方法。

在这里插入图片描述

更多统计学习基础知识请参考专栏《统计学习基础知识》。

更多人工智能基础知识请参考专栏《人工智能基础知识》。

写博不易,敬请支持:

如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!

关于老猿的付费专栏

  1. 付费专栏《https://blog.csdn.net/laoyuanpython/category_9607725.html 使用PyQt开发图形界面Python应用》专门介绍基于Python的PyQt图形界面开发基础教程,对应文章目录为《 https://blog.csdn.net/LaoYuanPython/article/details/107580932 使用PyQt开发图形界面Python应用专栏目录》;
  2. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10232926.html moviepy音视频开发专栏 )详细介绍moviepy音视频剪辑合成处理的类相关方法及使用相关方法进行相关剪辑合成场景的处理,对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/107574583 moviepy音视频开发专栏文章目录》;
  3. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10581071.html OpenCV-Python初学者疑难问题集》为《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的伴生专栏,是笔者对OpenCV-Python图形图像处理学习中遇到的一些问题个人感悟的整合,相关资料基本上都是老猿反复研究的成果,有助于OpenCV-Python初学者比较深入地理解OpenCV,对应文章目录为《https://blog.csdn.net/LaoYuanPython/article/details/109713407 OpenCV-Python初学者疑难问题集专栏目录 》
  4. 付费专栏《https://blog.csdn.net/laoyuanpython/category_10762553.html Python爬虫入门 》站在一个互联网前端开发小白的角度介绍爬虫开发应知应会内容,包括爬虫入门的基础知识,以及爬取CSDN文章信息、博主信息、给文章点赞、评论等实战内容。

前两个专栏都适合有一定Python基础但无相关知识的小白读者学习,第三个专栏请大家结合《https://blog.csdn.net/laoyuanpython/category_9979286.html OpenCV-Python图形图像处理 》的学习使用。

对于缺乏Python基础的同仁,可以通过老猿的免费专栏《https://blog.csdn.net/laoyuanpython/category_9831699.html 专栏:Python基础教程目录)从零开始学习Python。

如果有兴趣也愿意支持老猿的读者,欢迎购买付费专栏。

老猿Python,跟老猿学Python!

☞ ░ 前往老猿Python博文目录 https://blog.csdn.net/LaoYuanPython ░

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/685139.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

家用洗地机哪个牌子好?专家推荐榜单助你挑选最合适的洗地机

随着科技不断发展,智能家居产品逐渐融入我们日常生活中,洗地机作为家庭清洁必备工具,越来越受到消费者青睐,但是面对市面上种类繁多的洗地机,我们如何挑选到适合自己的产品呢?专家推荐榜单助你挑选最合适的…

在vue项目中使用markdown-it回显markdown文本

前言 其实有很多插件都是可以用来回显markdown文本的,这个插件也是其中之一。 文档地址:markdown-it | markdown-it 中文文档 这个文档在vue2和vue3里面都可以使用,所以还是比较推荐的 使用 安装 npm install markdown-it --save 应用 <template><div><…

正邦科技(第10天)

这里写目录标题 任务一任务二任务三 任务一 下位机报上来的十进制数据进行解析&#xff1a; 360170 固定报文&#xff0c;一个F对应一个字节&#xff0c;温度值&#xff0c;湿度值&#xff0c;烟雾浓度值是十进制转16进制&#xff0c;告警状态需要高低位移位&#xff0c;然后再…

【Pycharm】功能介绍

1.Code Reformat Code 格式化代码&#xff0c;可以帮助我们去自动调整空格等&#xff0c;根据python语法规范自动调整 2.Settings 1.创建py文件默认填充模版 3.读写py文件编码格式一致性 顶部代码指定的编码方式作用&#xff1a; 可以保证python2/3解释器在读取文件的时候按…

个人项目———密码锁的实现

布局组件 布局效果 组件绑定 密码锁的实现代码 using TMPro; using UnityEngine; using UnityEngine.UI;public class PasswordPanel : MonoBehaviour {// public Button button;// 所有按键的父物体public Transform buttonPanel;// 输入字符串的文本框public TMP_Text input…

英国树莓派五大天王和你相约上海国际嵌入式展!

6月12日-14日 上海世博展览馆3号馆 H3馆 237展位 树莓派(Raspberry Pi),这个曾经让全球掀起"创客热潮"的小型单板电脑,如今已经成为嵌入式行业不可或缺的一员。作为行业先驱,树莓派基金会正携手团队,亮相2024年6月12日至6月14日在上海举办的 Embedded World上海国…

【Elasticsearch】es基础入门-03.RestClient操作文档

RestClient操作文档 示例&#xff1a; 一.初始化JavaRestClient &#xff08;一&#xff09;引入es的RestHighLevelClient依赖 <!--elasticsearch--> <dependency><groupId>org.elasticsearch.client</groupId><artifactId>elasticsearch-rest…

基于springboot+vue的家乡特色推荐系统

开发语言&#xff1a;Java框架&#xff1a;springbootJDK版本&#xff1a;JDK1.8服务器&#xff1a;tomcat7数据库&#xff1a;mysql 5.7&#xff08;一定要5.7版本&#xff09;数据库工具&#xff1a;Navicat11开发软件&#xff1a;eclipse/myeclipse/ideaMaven包&#xff1a;…

SpringBoot高手之路03-事务传播行为

那么就是 a事务调用了b事务 日志技术 当解散部门的时候,那么就直接进行 操作日志 就是什么时间点吧部门解散 成功失败都需要记录日志 首先一个日志表 那么日志技术,在电商平台,不论是否支付订单,那么都需要保存订单信息 这个时候我们就使用传播事务 传播事务首先是出现在两…

云南区块链商户平台发票助手成品

目录 1 概述2 功能对比3 项目演示图4 核心逻辑4.1智能赋码4.2 解密方法4.3 登录与检测4.4 发票金额大写转换4.5 检查登录是否失效4.6 验证码识别5 演示效果6 项目部署6.1 Web站点部署6.1.1 环境6.1.2 前端6.1.3 后端6.2 Docker部署6.2.1 构建镜像6.2.2 创建容器6.3.3 访问项目域…

混剪素材哪里找?分享几个热门混剪素材下载网站

在短视频和新媒体的世界里&#xff0c;高质量的混剪素材是吸引观众的关键。今天&#xff0c;我将为大家详细介绍几个优秀的素材网站&#xff0c;它们不仅资源丰富&#xff0c;而且完全满足新媒体创作者的需求。这篇文章将帮助你理解如何有效利用这些平台提升你的视频创作。 蛙…

Python中的“点阵字体”

“点阵字体”是个啥&#xff1f;&#xff0c;在python中怎么使&#xff1f;在现在全面高清的 5 G 5G 5G时代&#xff0c;它还有用“武”之地&#xff1f; (笔记模板由python脚本于2024年06月01日 18:44:31创建&#xff0c;本篇笔记适合会基本编程的coder翻阅) 【学习的细节是欢…

BabylonJS 6.0文档 Deep Dive 动画(四):通过动画排序制作卡通片

一种最为直接的方法是为每个动画剪辑&#xff08;Animatin Clip&#xff09;指定开始时间&#xff0c;最终形成一个卡通动画&#xff08;Cartoon&#xff09;。 1. 设计 1.1 概述 动画的脚本如下&#xff1a; 摄像机显示了一栋带门的建筑物。摄像机靠近门并停止。门打开&am…

⾃动化批量管理-Ansible

目录 一、ansible 简介 自动化工具选择 &#xff08;了解&#xff09;​编辑 1、ansible 是什么&#xff1f; 2、ansible 特点 3、ansible 架构图 二、ansible 任务执行 1、ansible 任务执行模式 2、ansible 执行流程 3、ansible 命令执行过程 三、ansible 配置详解 …

Win32和c++11多线程

Win32和c11多线程 一、概念1.线程的特点线程内核对象线程控制块线程是独立调度和分派的基本单位共享进程的资源 2.线程的上下文切换引起上下文切换的原因 3.线程的状态 二、Windows多线程API1.CreateThread创建线程2.获取线程ID3.关闭线程句柄4.挂起线程5.恢复线程6.休眠线程的…

必应bing国内广告账户如何注册推广呢?

作为全球第二大搜索引擎&#xff0c;必应Bing以其庞大的用户基础和精准的定向能力&#xff0c;为企业提供了拓展市场的绝佳平台。对于许多企业来说&#xff0c;必应Bing广告账户的注册与推广流程可能显得复杂而繁琐。此时&#xff0c;您不妨考虑携手云衔科技&#xff0c;共同开…

echaerts图例自动滚动并隐藏翻页按钮

效果图 代码 legend: {itemHeight: 14,itemWidth: 14,height: "300", //决定显示多少个// 通过 CSS 完全隐藏翻页按钮pageButtonItemGap: 0,pageButtonPosition: end,pageIconColor: transparent, // 隐藏翻页按钮pageIconInactiveColor: transparent, // 隐藏翻页按…

神经网络 torch.nn---优化器的使用

torch.optim - PyTorch中文文档 (pytorch-cn.readthedocs.io) torch.optim — PyTorch 2.3 documentation 反向传播可以求出神经网路中每个需要调节参数的梯度(grad)&#xff0c;优化器可以根据梯度进行调整&#xff0c;达到降低整体误差的作用。下面我们对优化器进行介绍。 …

2024年6月6日第十五周六级系列1作文

对于2024年上半年六级作文预测的大纲 一、作文总体格式与要求 六级作文将继续沿用以往的大致格式&#xff0c;主要包括两个主要任务&#xff1a; Task 1: 针对给定话题的简短摘要或观点文章&#xff08;150-200字&#xff09;。Task 2: 就更复杂的话题进行深入分析的作文&…

23 - 每位教师所教授的科目种类的数量(高频 SQL 50 题基础版)

23 - 每位教师所教授的科目种类的数量 考点&#xff1a; 排序和分组 selectteacher_id,count(distinct subject_id) cnt fromTeacher group byteacher_id;