当前位置: 首页 > news >正文

如何理解泊松分布

文章目录

    • 一、引例——鲸鱼研究
    • 二、泊松分布


一、引例——鲸鱼研究

有生态学家对生活在北冰洋水域的鲸鱼进行了跟踪研究,他们利用一台水下无人机来探测鲸鱼数量,这是近十天的数据:

第1天第2天第3天第4天第5天第6天第7天第8天第9天第10天
1010141211158976

最后四天鲸鱼数目的连续减少趋势令人担忧,水质恶化吗?。

注意到每天探测到的鲸鱼不会是同时发生的,而是分散在一天的不同时间段。

这事实上是一个随机事件的时间分布,于是我们将一天时间划分为12段,那么第一天的数据如下:

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

那么对于这组数据,鲸鱼出现的概率 p = 5/6

我们试图把鱼群数量转换成时间上的二项分布,那么 12 个 时间段就转化成了一个二项分布,一天内出现10头鲸鱼的概率如下:
Cnk⋅pk⋅(1−p)n−k=C1210⋅(56)10⋅(16)2=0.2961C_n^k \cdot p^k \cdot (1-p)^{n-k} = C_{12}^{10} \cdot (\frac{5}{6})^{10} \cdot (\frac{1}{6})^2 = 0.2961 Cnkpk(1p)nk=C1210(65)10(61)2=0.2961
这正是 拉普拉斯 得意门生 西莫恩·德尼·泊松(1781~1840) 的思路。

前面我们选取的是第一天的数据,我们现在选取第二天的数据,有14头鲸鱼,12个时间段有些不够用,于是我们划分成24个时间段。

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

则一天中出现14头鲸鱼的概率:
pi=1424=712则Cnk⋅pk⋅(1−p)n−k=C2414⋅(712)14⋅(512)10=0.1634\begin{align} & p_i = \frac{14}{24} = \frac{7}{12} & 则 C_n^k \cdot p^k \cdot (1-p)^{n-k} = C_{24}^{14} \cdot (\frac{7}{12})^{14} \cdot (\frac{5}{12})^{10} = 0.1634 \end{align} pi=2414=127Cnkpk(1p)nk=C2414(127)14(125)10=0.1634
但是我们上面的思考有一个漏洞就:一个时间段是有可能有多头鲸鱼出现的

这说明什么?时间段太大了!我们直接把一天划分成无数个时间段,单独每个时间段都趋于无穷小,这样无论两头鲸鱼的时间挨得有多近,都能区分开来,则有:
limn→∞Cnk⋅pk⋅(1−p)n−k,其中:n代表时间分段数,p代表任意时间段内鲸鱼出现的概率p的估计不再以某一天数据为准,而是选取平均值,即p=μn\begin{align} & lim_{n \rightarrow \infty} C_n^k \cdot p^k \cdot (1-p)^{n-k},其中:\\ & n代表时间分段数,p代表任意时间段内鲸鱼出现的概率 \\ & p的估计不再以某一天数据为准,而是选取平均值,即 p = \frac{\mu}{n} \end{align} limnCnkpk(1p)nk,其中:n代表时间分段数,p代表任意时间段内鲸鱼出现的概率p的估计不再以某一天数据为准,而是选取平均值,即p=nμ

公式的进一步推导
limn→∞Cnk⋅pk⋅(1−p)n−k=limn→∞Cnk⋅(μn)k⋅(1−μn)n−k=lim⁡n→∞n(n−1)(n−2)⋯(n−k+1)k!⋅μknk⋅(1−μn)n−k=μkk!⋅lim⁡n→∞n(n−1)(n−2)⋯(n−k+1)nk⋅(1−μn)n−k=μkk!⋅lim⁡n→∞nn(n−1)n(n−2)n⋯(n−k+1)n⋅(1−μn)n−k=μkk!⋅lim⁡n→∞1⋅1…1⋅(1−μn)n−k=μkk!⋅lim⁡n→∞(1−μn)n−k=μkk!⋅lim⁡n→∞(1−μn)n⋅lim⁡n→∞(1−μn)−k=μkk!⋅lim⁡n→∞(1−μn)n=μkk!e−μ(等价无穷小)\begin{align} & lim_{n \rightarrow \infty} C_n^k \cdot p^k \cdot (1-p)^{n-k} \\ =& lim_{n \rightarrow \infty} C_n^k \cdot (\frac{\mu}{n})^k \cdot (1-\frac{\mu}{n})^{n-k} \\ =& \lim _{n \rightarrow \infty} \frac{n(n-1)(n-2) \cdots(n-k+1)}{k!} \cdot \frac{\mu^{k}}{n^{k}} \cdot\left(1-\frac{\mu}{n}\right)^{n-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} \frac{n(n-1)(n-2) \cdots(n-k+1)}{n^{k}} \cdot\left(1-\frac{\mu}{n}\right)^{n-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} \frac{n}{n} \frac{(n-1)}{n} \frac{(n-2)}{n} \cdots \frac{(n-k+1)}{n} \cdot\left(1-\frac{\mu}{n}\right)^{n-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} 1 \cdot 1 \dots 1 \cdot\left(1-\frac{\mu}{n}\right)^{n-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} \left(1-\frac{\mu}{n}\right)^{n-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} \left(1-\frac{\mu}{n}\right)^{n} \cdot \lim _{n \rightarrow \infty} \left(1-\frac{\mu}{n}\right)^{-k} \\ =& \frac{\mu^{k}}{k!} \cdot \lim _{n \rightarrow \infty} \left(1-\frac{\mu}{n}\right)^{n} \\ =& \frac{\mu^{k}}{k!} e^{-\mu} (等价无穷小) \\ \end{align} =========limnCnkpk(1p)nklimnCnk(nμ)k(1nμ)nknlimk!n(n1)(n2)(nk+1)nkμk(1nμ)nkk!μknlimnkn(n1)(n2)(nk+1)(1nμ)nkk!μknlimnnn(n1)n(n2)n(nk+1)(1nμ)nkk!μknlim111(1nμ)nkk!μknlim(1nμ)nkk!μknlim(1nμ)nnlim(1nμ)kk!μknlim(1nμ)nk!μkeμ(等价无穷小)

二、泊松分布

我们把上面的μ 换成 λ,就得到了我们熟悉的泊松分布公式:
λkk!e−λ,其中λ代表给定时间段或空间区域内,随机事件平均发生的次数。\frac{\lambda^k}{k!}e^{-\lambda},其中 \lambda 代表给定时间段或空间区域内,随机事件平均发生的次数。 k!λkeλ,其中λ代表给定时间段或空间区域内,随机事件平均发生的次数。

我们利用该公式计算引例中各时间段鲸鱼出现数目的可能性,并用柱状图来表示:

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

我们取显著性水平为 5%,那么从图中可以得知,鲸鱼数目在[5, 15]都是正常范围,我们的样本都在正常范围内。

有时候我们会用用泊松分布近似二项分布,一般在 试验次数n 很大,成功概率 p 很小,且 λ = np 为有限值时,二项分布可以用泊松分布来近似,经过前面的推导,不难明白其中缘由。

http://www.dtcms.com/a/297070.html

相关文章:

  • Vue基础(24)_VueCompinent构造函数、Vue实例对象与组件实例对象
  • UG创建的实体橘黄色实体怎么改颜色?
  • STM32项目实战:正弦波
  • 什么是JSON,如何与Java对象转化
  • QT跨平台应用程序开发框架(11)—— Qt网络编程
  • [NLP]一个完整的 UPF 文件示例
  • Vim 编辑器全模式操作指南
  • 纳米编辑器之Nano 编辑器退出**的详细操作指南
  • 【MacOS】发展历程
  • Linux 中 `chown`、`chgrp` 和 `chmod` 命令详解
  • openGauss数据库在CentOS 7 中的单机部署与配置
  • VMware虚拟出来的centos中设置静态ip
  • fish-speech 在50系列显卡使用 --compile加速兼容
  • Rust + Tauri 开发所需环境清单(以 Windows 为例)
  • 【unitrix】 6.16 非负整数类型( TUnsigned )特质(t_unsingned.rs)
  • [Rust 基础课程]猜数字游戏-获取用户输入并打印
  • 智能问答分类系统:基于SVM的用户意图识别
  • 弹性网:基于神经网络的多组分磁共振弹性成像波反演与不确定性量化|文献速递-医学影像算法文献分享
  • 奥比中光的dabai_dcw2相机彩色对齐方案
  • Android Camera setRepeatingRequest
  • 11. isaacsim4.2教程-Transform 树与Odometry
  • Java面试题(中等)
  • Cartographer安装测试与模块开发(三)--Cartographer在Gazebo仿真环境下的建图以及建图与定位阶段问题(实车也可参考)
  • 融合与智能:AI 浪潮驱动下数据库的多维度进化与产业格局重塑新范式
  • 深入解析Linux匿名管道机制与应用
  • 从数据孤岛到融合共生:KES V9 2025 构建 AI 时代数据基础设施
  • Lua 函数
  • JAVA_THIRTEEN_常用API
  • 星慈光编程虫2号小车讲解第三篇--附件概述
  • ai存在意义的对话