伽马(gamma)函数_伽马分布的分布函数怎么求

伽马(gamma)函数_伽马分布的分布函数怎么求一。ΓΓ\Gamma分布指数分布是两次事件发生的时间间隔ΓΓ\Gamma分布是n倍的指数分布即,ΓΓ\Gamma分布表示发生n次(αα\alpha次)事件的时间间隔的概率分布可以直观地认为ΓΓ\Gamma分布是Possion分布在正实数集上的连续化版本Possion(X=k|λ)=λke−λk!Possion(X=k|λ)=λke−λk!Possion(X=k|\lambda…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

一。 Γ \Gamma Γ分布

指数分布是两次事件发生的时间间隔
Γ \Gamma Γ分布是n倍的指数分布
即, Γ \Gamma Γ分布表示发生n次( α \alpha α次)事件的时间间隔的概率分布。

其实 Γ \Gamma Γ分布 就是Possion分布在正实数集上的连续化版本

P o s s i o n ( X = k ∣ λ ) = λ k e − λ k ! Possion(X=k|\lambda )=\frac{\lambda ^ke^{-\lambda} }{k!} Possion(X=kλ)=k!λkeλ

可以看做横坐标是k,纵坐标便是X=k的概率分布: P o s s i o n ( X = k ∣ λ Possion(X=k|\lambda Possion(X=kλ), λ \lambda λ为一个常数,代表单位时间内事件发生的次数。

回忆一下Poisson分布的直观含义:
随机变量X代表出生婴儿的个数,P {X=k} 代表出生k个婴儿的概率,
λ \lambda λ为已知数,代表平均单位时间出生婴儿的个数。
求t时间内出生k个婴儿的概率:P {X=k} ,令t=1,就是
这个公式 P o s s i o n ( X = k ∣ λ ) = λ k e − λ k ! Possion(X=k|\lambda )=\frac{\lambda ^ke^{-\lambda} }{k!} Possion(X=kλ)=k!λkeλ
因此Possion公式的直观意义就是:
已知单位时间内平均出生 λ \lambda λ个婴儿, 得到单位时间内出生k个婴儿的概率。
如果将k看成是一个变量, Possion公式就是单位时间内出生婴儿个数的概率分布。
直观理解,当然是单位时间出生 λ \lambda λ个婴儿的概率最大。

在 P o i s s o n 分 布 中 , λ 是 一 个 已 知 数 , 是 一 个 常 数 , {\color{Red}{ 在Poisson分布中,\lambda 是一个已知数,是一个常数,}} Poissonλ
如 果 我 们 把 λ 看 成 一 个 变 数 , 假 设 是 x {\color{Red}{如果我们把\lambda看成一个变数,假设是x}} λx
那 么 得 到 的 分 布 就 叫 G a m m a 分 布 {\color{Red}{那么得到的分布就叫Gamma分布}} Gamma,显然Gamma比Poisson更高一维的分布。

= > 将 λ 转 为 一 个 连 续 实 数 x {\color{Red}{=> 将 \lambda 转为一个连续实数 x}} =>λx

在Gamma分布的密度中取 α = k + 1 \alpha =k+1 α=k+1,得

G a m m a ( x ∣ α = k + 1 ) = x k e − x Γ ( k + 1 ) = x k e − x k ! Gamma(x|\alpha =k+1)=\frac{x^{k} e^{-x} }{\Gamma (k+1)} = \frac{x^k e^{-x } }{k!} Gamma(xα=k+1)=Γ(k+1)xkex=k!xkex

G a m m a ( x , k ) = x k e − x k ! Gamma(x,k) = \frac{x^k e^{-x } }{k!} Gamma(x,k)=k!xkex

由此可见,Gamma函数是一个关于x和k的二维概率分布。x是单位时间内事件发生的平均次数,k是单位时间内事件发生的某一特定次数,得到类似于下图,可见,它是一个指数分布,k与越接近,概率越大,在k与x相等的地方,概率达最大值。(如果将x固定一个常数,就是Poisson分布。)

所以,Gamma分布与Possion分布在数学形式上是一致的,只是Poisson分布是离散的,Gamma分布是连续的,可以直观的认为Gamma分布是Poission分布在正实数集上的连续化版本。

在这里插入图片描述

二。 Γ \Gamma Γ函数
定义
Γ ( s ) = ∫ 0 + ∞ e − x x s − 1 d x ( s > 0 ) \Gamma (s)=\int_{0}^{+\infty }e^{-x}x^{s-1}dx(s>0) Γ(s)=0+exxs1dx(s>0)

性质

  1. s>0时,此反常积分收敛
  2. Γ ( s + 1 ) = s Γ ( s ) ( s > 0 ) \Gamma(s+1)=s\Gamma(s) (s>0) Γ(s+1)=sΓ(s)(s>0),特别 Γ ( n + 1 ) = n ! \Gamma(n+1)=n! Γ(n+1)=n!
    3)当 s → 0 + s\to0+ s0+时, Γ ( s ) → + ∞ \Gamma(s) \to +\infty Γ(s)+
  3. Γ ( s ) Γ ( 1 − s ) \Gamma(s)\Gamma(1-s) Γ(s)Γ(1s)= π s i n π s ( 0 &lt; s &lt; 1 ) , 则 Γ ( 1 2 ) = π \frac{\pi }{sin\pi s} (0&lt;s&lt;1), 则\Gamma (\frac{1}{2})=\sqrt{\pi} sinπsπ(0<s<1),Γ(21)=π

$\Gamma(n) = (n-1)! $ , Gamma(5+1) = 5! =120
$\Gamma(s) = (s-1)! $ , 5Gamma(5) = 54! =120

三。 Γ \Gamma Γ函数应用
k ! = ∫ 0 ∞ x k e − x d x k!=\int_{0}^{\infty }x^ke^{-x}dx k!=0xkexdx

Γ ( s ) = ∫ 0 ∞ x s − 1 e − x d x \Gamma(s)=\int_{0}^{\infty }x^{s-1}e^{-x}dx Γ(s)=0xs1exdx 中,
作x=u^2的代换可得
Γ ( s ) = 2 ∫ 0 ∞ e − u 2 u 2 s − 1 d u \Gamma (s)=2\int_{0}^{\infty }e^{-u^2}u^{2s-1}du Γ(s)=20eu2u2s1du
再令 t=2s-1,即有
∫ 0 ∞ e − u 2 u t d u \int_{0}^{\infty }e^{-u^2}u^{t}du 0eu2utdu = 1 2 Γ ( 1 + t 2 ) \frac{1}{2}\Gamma(\frac{1+t}{2}) 21Γ(21+t), t>-1
特别,令 s = 1 2 s=\frac{1}{2} s=21, 可得概率论中常用积分

泊松积分
泊松积分: ∫ − ∞ ∞ e − x 2 d x = π \int_{-\infty}^{\infty }e^{-x^2}dx=\sqrt{\pi} ex2dx=π
是一个很重要的结论,在概率论中有重要应用。
∫ 0 ∞ e − x 2 d x = π 2 \int_{0}^{\infty }e^{-x^2}dx=\frac{\sqrt{\pi}}{2} 0ex2dx=2π

∫ − ∞ ∞ e − x 2 d x = π \int_{-\infty}^{\infty }e^{-x^2}dx=\sqrt{\pi} ex2dx=π

(可以用正态分布的公式证明,正态分布公式 Φ ( x ) = 1 2 π ∫ − ∞ ∞ e − t 2 2 d t = 1 \Phi (x)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}e^{-\frac{t^2}{2}}dt=1 Φ(x)=2π
1
e2t2dt=
1

参考链接:
Gamma分布Wiki百科
poisson-gamma-exponential 泊松-Gamma以及指数分布的关系
Gamma distribution in R语言
MATLAB Gamma
神奇的Gamma函数(scipy)

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/200646.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • Linux正确删除软连接[通俗易懂]

    在Linux上删除软连接不要使用rm-rf!!!在Linux上删除软连接不要使用rm-rf!!!在Linux上删除软连接不要使用rm-rf!!!因为如果使用rm-rflinkName的方式,如果不小心在目录后面加了“/”,或者按了Tab键补全,执行之后会删除源目录文件。如果要使用rm-rflinkName的时候一定要注意源、目标文件或目录都不要在后面加…

    2022年4月14日
    195
  • jquery 字符串转数组

    如何把字符串”123,abc,xy,hi”转换成数组或者jQuery数组vararr=’123,abc,xy,hi’.split(‘,’);for(variinarr){alert(arr[i])}出处:http://www.cnblogs.com/xiaomuv587/archive/2012/10/18/2729766.html

    2022年4月14日
    142
  • Vmware Links(转自VMware-land)「建议收藏」

    Vmware Links(转自VMware-land)「建议收藏」这一阵子在专研虚拟机的VSS备份,无意中发现了VMware-land很好的网站,不知道为什么无法访问,难道也被和谐掉了???以下内容是从Google的页面缓存弄出来的,在Google搜索http://vmware-land.com/Vmware_Links.html第一个就是包含了你所知道的和不知道的,关于VMwareESX的方方面面链接地址Backups:Vir…

    2022年10月3日
    6
  • docker连接宿主机mysql_docker搭建mysql主从

    docker连接宿主机mysql_docker搭建mysql主从相信我,看完这一篇,mysql主从复制能遇到的错误在我这里都遇到了,docker能遇到的错误在我这里也遇到了,包括centos的错误如果你还没有用docker搭建过mysql先不要急着搭建集群,可以显示着用docker搭建一个mysql,这里有教程,可以先试着去搭建一个mysql:之前我们使用docker搭建了一个Mysql,那既然一个Mysql我们可以搭建成功,集群还不是分分钟,今天我们就简单搭建一个Mysql集群,就搭建一个最简单的一主一从的主从复制吧,

    2022年8月13日
    6
  • jmeter正则提取器的使用_java正则表达式用法

    jmeter正则提取器的使用_java正则表达式用法一、正则表达式提取器各名词解(1)ApplytoMainsampleandsub-samples(作用于主节点的取样器及对应子节点的取样器)Mainsampleonly(仅作用于主节点的取样器)Sub-samplesonly(仅作用于子节点的取样器)Jmeter-VariableNametouse(作用于jmeter变量(输入框内可输入jmeter的变量名称),从指定变量值中提取需要的值)(2)要检查的响应字段1、主体:响应报文的主体,最常用2、Body(une.

    2025年10月22日
    5
  • java 取余和取整_Java取整、取余

    java 取余和取整_Java取整、取余参考链接:http://blog..net/wanlixingzhe/article/details/7359809参考链接:http://bbs..net/topics/390677448(6楼)参考链接:http://blog.sina.com.cn/s/blog_6940cab30101hji5.html最近在做一个计算的时候用到了取整取余的计算,这里对取整、取余、取模做一下总结~~~1、取…

    2022年5月30日
    164

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号