浅谈几个著名的大数定律及应用李 蕊(青海大学成教学院,青海 西宁,810001)摘均后得到的随机变量 ,当n充分大时, 以概率1趋于 要 大数定律以严格的数学形式表达了随机现象最根本的性质——平均结果的稳定性,是随机现象统计规律性的具体 。这样,切比雪夫大数定律给出了平均值稳定性科学的表现,本文介绍了几种常用的大数定律,并给出一些简单应用。描述。[1]关键词 大数定律 随机变量 数学期望 概率推论1:设随机变量ξ,ξ,…ξ相互独立,它们具有相12n2同的数学期望和方差,Eξ=a,Dξ=σ(i=1,2,…)则对于任意1 引言ii“大数定律”本来是一个数学概念,又叫做“平均法给定的正数ε,有则”。在随机事件的大量重复出现中,往往呈现几乎必然的规此推论表明:n个相互独立的具有相同数学期望和方差的律,这个规律就是大数定律,通俗地说,这个定理就是在试验随机变量,当n很大时,它们的算术平均值几乎是一个常数,这不变的条件下,重复试验多次,随机事件的频率以概率为稳定个常数就是它们的数学期望。[2]值。比如,我们向上抛一枚硬币,硬币落下后哪一面朝上本来定理2(贝努力大数定律)是偶然的,但当我们向上抛硬币的次数足够多时,达到上万次设是u次独立试验中事件A发生的次数,p是事件A在每次试n甚至几十万几百万次之后,我们就会发现,硬币向上的次数约验中发生的概率,则对任意给定的正数ε,有占总次数的二分之一。偶然中包含着必然。此定理表明:当n很大时,n重贝努力试验中事件A发生的从概率的统计定义中可以看出:一个事件发生的频率具有频率几乎等于事件A在每次试验中发生的概率,这个定理以严格稳定性,即随着试验次数的增多,事件的频率逐渐稳定在某个的数学形式刻画了频率的稳定性,因此,在实际应用中,当试常数附近.人们在实践中观察其他一些随机现象时,也常常会验次数很大时,便可以用事件发生的频率来代替事件的概率。发现大量随机个体的平均效果的稳定性。这就是说,无论个别贝努力大数定律是切贝晓夫大数定律的特例,在它们的证明随机个体以及它们在试验进行过程中的个别特征如何,大量随中:都是以切贝晓夫不等式为基础的,所以要求随机变量具有机个体的平均效果与每一个体的特征无关,且不再是随机的深方差,但是进一步的研究表明,方差这个条件并不是必要的。[2]入考虑后,人们会提出这样的问题:稳定性的确切含义是什么?定理3(辛钦大数定律)在什么条件下具有稳定性?这就是大数要研究的问题。设ξ,ξ,…ξ是相互独立的随机变量,而且有相同的分12n布,具有有限的数学期望Eξ(i=1,2,…),则对任意给定的正2 几个大数定律i在介绍大数定律之前,先介绍几个相关定义。数ε,有 ,其中 a=Eξi[1]定义1设ξ(n=1,2,……)为概率空间(Ω,F,P)上定义的前面说过贝努力大数定律表明了当n很大时,事件发生的n随机变量序列(简称随机序列),若存在随机变数ξ,使对任频率会“靠近”概率,而这里的辛钦大数定律表明,当n很大意ε>0,恒有: 则称随机序时,随机变量在n次观察中的算术平均值会靠近它的期望。[1]列 依概率收敛于随机变量ξ(ξ也可以是一个常数),并用定理4 (泊松大数定律)下面的符号表示:设ξ,ξ,…ξ是相互独立的随机变量,12n[2]定义2设 为一随机序列,数学期望E(ξ)存在,令 (其中p=1-q),则 服从大数nnn ,若 ,则称随机序列 服从大数定律。定律,或者说大数法则成立。泊松大数定律是贝努力大数定律的推广,贝努利大数定切比雪夫不等式律证明了事件在完全相同条件下重复进行的随机试验中频率的设随机变量X的数学期望E(X)与方差D(X)存在,则对稳定性,而泊松定理表明,当独立进行的随机试验的条件变化于任意正数ε,不等式时,频率仍然具有稳定性.随着n的无限增大,在n次独立试验 都成立。不等式(1)和(2)称为切比雪夫不等式。中,事件A的频率在各次试验中事件A出现概率的算术平均值处切比雪夫不等式给出了在随机变量X的分布未知的情况取得稳定值。[1]下,只利用J的数学期望和方差即可对J的概率分布进行估值的定理5(马尔可夫大数定律)方法,这就是切比雪夫不等式的重要性所在。设{ξ}是随机变量序列,若 ,则{ξ}ii大数定律形式很多,我们仅介绍几种最常用的大数定律。服从大数定律。[1]定理1 (切比雪夫大数定律)3 应用设随机变量ξ,ξ,…ξ相互独立,它们的数学期望依次 误差领域的应用12n222为a,a,…a方差依次为σ,σ,…σ而且存在正常数k,使例1 仪器测量已知量A时,设n次独立得到的测量数据为 12n12n2得对一切i=1,2,…,有σ<k则对任意给定的数ε,恒有x,x…x如果仪器无系统误差,问:当n充分大时,是否可取 i12n 作为仪器测量误差的方差的值?切比雪夫大数定律说明:独立随机变量序列ξ,ξ,…ξ解:把x视作n个独立同分布的随机变量(i=1,2…n)的观12ni2的数学期望与方差都存在,且方差一致有上界,则经过算术平察值,则E(X)=u,D(X)=σ(i=1,2,…,n),仪器第i次测量的 4
2误差X-A的数学期望E(X-A)=u-A,方差D(X-A)=σ,设y=(X-似值。iiiii2A),i=1,2…n,则y,也相互独立服从于统一分布。在仪器无例3[3].假设 求其极限。i系统误差时E(X-A)=0,既有u=A解:假设随机变量ξ(n=1,2,…n),在[0,1]上均匀分布,in222E(y)=E[(X-A)]=E[(X-E)],D(X)=σ,i=1,2,…n而且相互独立,有iiiXii由切比雪夫定律,可得: 即 。从而,当n→∞时,随机变量 2依概率收敛于σ,即当n充分大时,可以取 作为仪 器测量误差的方差。 根据大数定律,对于随机误差ξ,ξ,…ξ,应有 12n这说明当测量次数较多时,实测数据的平均值 和预测真值的差以很大概率趋于0,因此,用求样本数据平均值的方法来22进行测量是可行的。由ξ,ξ,…ξ独立同分布,可见ξξ独立同分布,12n1…2…根据辛钦大数定律知: 在分布型未知的情况下估计数学期望E(ξ)及方差D(ξ) 若ξ及{ξ}都是随机变量,则有:k[3]例4.在贝努力试验中,事件A出现的概率率为p,令x=1,若在第n次及第n+1次试验中A都出现;x=0其他证明{x}nnn服从大数定律。证:{x}为同分布随机变量序列,其共同分布为n22p(x=0)=1-p,p(x=1)=pnn2222且E(x)=E(x)=p,从而Var(x)=p(1-p)<1,又当│i-j│nnn>2 时,x与x独立所以ij 在数学分析中的应用例2.计算定积分 的近似值。为了解这种近似计算的依据,先进行如下分析:若令 为均匀分布的概率密度函数,即若令 为均匀分布即马尔可夫条件成立,故服从大数定律。的概率密度函数,即 ,则 而函参考文献数g(x)的数学期望 根据大数定律应用可对该[1]黄清龙,阮宏顺.概率论与数理统计[M].北京:北数学期望值进行估计,即京大学出版社,2005. 故可用 [2]杨亚非.概率与数理统计基础[M].北京:北京工业这种近似计算的具体过程如下:出版社,2003.欲计算定积分 的近似值,则应先取样本数列[3]峁诗松,程依明,濮晓龙.概率论与数理统计教程习{X}→求函数序列(x)→求出 ,即作为J的近题与解答[M].北京:高等教育出版社,(上接第63页)低的条件,在这三个限制条件下选择最佳的承台安全经济高度。(~)、(~)、(~) 、(~)从以上的承台材料费用计算表中的计算结果可知,该桩基时的每立方米混凝土含筋率为~%(采用HRB335热实例承台的安全经济高度在 ~ 左右。大于此高度时,轧钢筋),(当采用HPB235热轧钢筋时每立方米混凝土含筋率为虽然降低了承台配筋数量,但承台抗冲切满足系数超过,抗~%),桩基承台材料费用最低,此时相应的承台高剪切满足系数k>,k>,这说明承台高度设计,十分保度为安全经济高度。B1B2守,增大了混凝土材料用量,致使材料总费用偏高;当设计的笔者认为,轴心受压矩形承台的安全经济高度可由公式 承台高度小于此高度时,虽然抗冲切系数及抗剪切系数都大于 初步计算确定(F为竖向力设计值,b、h为柱截cc1,满足规范要求,但没有合理地利用承台高度的抗弯能力,使面尺寸);偏心受压矩形承台的安全经济高度一般可用公式得承台配筋量大为增加,也导致承台材料费用偏高。初步计算确定;三角形三桩承台的安全经济高度根据笔者研究,当采用C15、C20、C25、C30、C35强度等可由公式初步计算确定。级的混凝土,承台的抗冲切满足系数在分别为(~)、