博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
高斯 到 正态分布 的前世今生
阅读量:4611 次
发布时间:2019-06-09

本文共 2160 字,大约阅读时间需要 7 分钟。

一、正态分布,熟悉的陌生人

学过基础统计学的同学大都对正态分布非常熟悉。这个钟型的分布曲线不但形状优雅,其密度函数写成数学表达式

 也非常具有数学的美感。其标准化后的概率密度函数

更加的简洁漂亮,两个最重要的数学常量π,e都出现在了公式之中。在我个人的审美之中,它也属于top-N的最美丽的数学公式之一,如果有人问我数理统计领域哪个公式最能让人感觉到上帝的存在,那我一定投正态分布的票。因为这个分布戴着神秘的面纱,在自然界中无处不在,让你在纷繁芜杂的数据背后看到隐隐的秩序。

 

【正态分布曲线】

正态分布又通常被称为高斯分布,在科学领域,冠名权那是一个很高的荣誉。早年去过德国的兄弟们还会发现,德国的钢镚和10马克的纸币上都留有高斯的头像和正态密度曲线。正态分布被冠名高斯分布,我们也容易认为是高斯发现了正态分布,其实不然,不过高斯对于正态分布的历史地位的确立是起到了决定性的作用。

 

【德国马克上的高斯头像和正态分布曲线】

正态曲线虽然看上去很美,却不是一拍脑袋就能想到的。我们在本科学习数理统计的时候,课本一上来介绍正态分布就给出密度分布函数,却从来不说明这个分布函数是通过什么原理推导出来的。所以我一直搞不明白数学家当年是怎么找到这个概率分布曲线的,又是怎么发现随机误差服从这个奇妙的分布的。我们在实践中大量的使用正态分布,却对这个分布的来龙去脉知之甚少,正态分布真是让人感觉既熟悉又陌生。直到我读研究生的时候,我的导师给我介绍了陈希儒院士的《数理统计学简史》这本书,看了之后才了解了正态分布曲线从发现到被人们重视进而广泛应用,也是经过了几百年的历史。

正态分布的这段历史是很精彩的,我们通过讲一系列的故事来揭开她的神秘面纱。

 

二、邂逅,正态曲线的首次发现

第一个故事和概率论的发展密切相关,主角是棣莫弗(De Moivre)和拉普拉斯(Laplace)。拉普拉斯是个大科学家,被称为法国的牛顿;棣莫弗名气可能不算很大,不过大家应该都熟悉这个名字,因为我们在高中数学学复数的时候我们都学过棣莫弗定理

 

(cosθ+isinθ)n=cos(nθ)+isin(nθ)(cos⁡θ+isin⁡θ)n=cos⁡(nθ)+isin⁡(nθ).

古典概率论发源于赌博,惠更斯、帕斯卡、费马、贝努利都是古典概率的奠基人,他们那会研究的概率问题大都来自赌桌上,最早的概率论问题是赌徒梅累在1654年向帕斯卡提出的如何分赌金的问题。统计学中的总体均值之所以被称为期望(Expectation),就是源自惠更斯、帕斯卡这些人研究平均情况下一个赌徒在赌桌上可以期望自己赢得多少钱。

 

棣莫弗(De Moivre) 拉普拉斯 (Laplace)

有一天一个哥们,也许是个赌徒,向棣莫弗提了一个和赌博相关的一个问题:A,B两人在赌场里赌博,A,B各自的获胜概率是pp和q=1pq=1−p,赌nn局,若A赢的局数X>npX>np,则A付给赌场XnpX−np元,否则B付给赌场npXnp−X元。问赌场挣钱的期望值是多少?

问题并不复杂,本质上是一个二项分布,最后求出的理论结果是

2npqb(n,p,np),2npqb(n,p,np),

其中

b(n,p,i)=(ni)piqnib(n,p,i)=(ni)piqn−i

是常见的二项概率。但是对具体的nn,要把这个理论结果实际计算出数值结果可不是件容易的事,因为其中的二项公式中有组合数.这就驱动棣莫弗寻找近似计算的方法。

与此相关联的另一个问题,是遵从二项分布的随机变量XB(n,p)X∼B(n,p),问XX落在二项分布中心点一定范围的概率Pd=P(|Xnp|d)Pd=P(|X−np|≤d)是多少?

对于p=12p=12的情形,棣莫弗做了一些计算并得到了一些近似结果,但是还不够漂亮,幸运的是棣莫弗和斯特林(Stirling)处在同一个时代,而且二人之间有联系,斯特林公式是在数学分析中必学的一个重要公式:(事实上斯特林公式的形式其实是棣莫弗最先发现的,但是斯特林改进了这个公式,改进的结果为棣莫弗所用)

n!2πn−−−√(ne)nn!∼2πn(ne)n

1733年,棣莫弗很快利用斯特林公式进行计算并取得了重要的进展。考虑nn是偶数的情形,令二项概率

b(i)=b(n,12,i)=(ni)(12)nb(i)=b(n,12,i)=(ni)(12)n

通过斯特林公式做一些简单的计算容易得到,

b(n2)2πn−−−√  b(n2+d)b(n2)exp(2d2n).b(n2)∼2πn ⇒ b(n2+d)b(n2)∼exp(−2d2n).

于是有

b(n2+d)22πn−−−√exp(2d2n).b(n2+d)∼22πnexp(−2d2n).

使用上式的结果,并在二项概率累加求和的过程中近似的使用定积分代替求和,很容易就能得到

 

  P(∣∣∣Xn12∣∣∣cn−−√)=   =   cn√icn√b(n2+i)cn√icn√22πn−−−√exp(2i2n)                                        (1)2c2in√2c12π−−√exp(12(2in−−√)2)2n−−√

转载于:https://www.cnblogs.com/nucdy/p/6343617.html

你可能感兴趣的文章
BZOJ-3289 Mato的文件管理
查看>>
自旋锁和互斥锁的区别
查看>>
react混合开发APP,资源分享
查看>>
入门篇
查看>>
【洛谷1829】 [国家集训队] Crash的数字表格(重拾莫比乌斯反演)
查看>>
[转]免费api大全
查看>>
git 认证问题之一的解决 : http ssh 互换
查看>>
sql where 1=1作用
查看>>
搜索算法----二分查找
查看>>
Python语言编程
查看>>
[poj 1469]Courses
查看>>
Xcode8出现AQDefaultDevice(173):Skipping input stram 0 0 0x0
查看>>
数据结构(二十四)二叉树的链式存储结构(二叉链表)
查看>>
Material Design Lite,简洁惊艳的前端工具箱 之 布局组件。
查看>>
关于bootstrap Modal弹窗 滚动条的问题
查看>>
Django----------路由控制
查看>>
将数字转化为字符串的快捷方式
查看>>
java23种设计模式
查看>>
冲刺周期一--站立会议04
查看>>
支持IE6以上阴影效果纯CSS
查看>>