自考概率论与数理统计Word文件下载.docx

1、12X （ 0.1+0.7+0.2+1+1.9+1.3+1.8 ）= 2例7 4在一批产品取样 n件，发现其中有 m件次品，试用此样本求该批产品的次品 率p的矩估计。10070103针对该题提问】因为/? = 二例如抽样总数n=100，其中次品m=5.= = 0.05则. :.J.例7 5电话总机在一分钟间隔内接到呼唤次数 XP （上）。观察一分种接到呼唤次数共观察40次，结果如下接到呼唤次数345观察次数10128求未知参数的矩估计-.10070104针对该题提问】（1）v X P （二） EX=-由矩法二二儿 .T（2）计算 - （0X5+1 X10+2X12+3X8+4X3+5X2）=

2、2-. = 27.1.2极大似然估计为了叙述极大似然原理的直观想法，先看例 7 6例7 6设有外表完全相同的两个箱子， 甲箱中有99个白球和1个黑球，乙箱中有99个黑球和1个白球，现随机地抽取一箱， 并从中随机抽取一球，结果取得白球，问这球是从 哪一个箱子中取出的？10070105针对该题提问】不管是哪一个箱子，从箱子中任取一球都有两个可能的结果： A表示取出白球，B表示取出黑球，如果我们取出的是甲箱，则 A发生的概率为0.99，而如果取出的是乙箱，则A发生的概率为0.01，现在一次试验中结果 A发生了，人们的第一印象就是： 此白球（A） 最像从甲箱取出的 ”或者是说，应该认为试验条件对事件

3、A出现有利，从而可以推断这球是从甲箱中取出的，这个推断很符合人们的经验事实，这里 最像”就是极大似然”之意。本例中假设的数据很极端，一般地，我们可以这样设想，在两个箱子中各有 100个球，甲箱中白球的比例是 Pi，乙箱中白球的比例是 P2,已知Pi P2,现随机地抽取一个箱子并 从中抽取一球，假定取到的是白球，如果我们要在两个箱子中进行选择， 由于甲箱中白球的比例高于乙箱，根据极大似然原理，我们应该推断该球来自甲箱。下面分别给出离散型随机变量和连续型随机变量的极大似然估计求未知参数 的估计的步骤（1）离散型随机变量第一步，从总体 X取出样本Xi,X2,x第二步，构造似然函数L （ Xi,X2，

4、x，幻）=P （ X = Xi） P （ X = X2）P （X = xn）第三步，计算ln L （Xi,X2,x, J）并化简第四步，当=卜时ln L （Xi,X2,Xi, ）取最大值则取=-1常用方法是微积分求最值的方法。（2）连续型随机变量若Xf （ x卫）第一步 从总体X取出样本Xi,X2，,X第二步构造似然函数L （Xi,X2，X, ）= f （Xi, 4） f （X2,旳彳（Xn, 9）第三步 计算In L （xi,x2,x,）并化简第四步 当一时In L （Xi,X2,x., 0）取最大值则取常用方法是微积分求最值的方法例7 7设总体XB （i , P）即ri*若a发生I0若A不

5、发生设P （A ）=，从总体X中抽样Xi,X2，,x，问最大似然法求i0070i06针对该题提问】当 X B （i, P）时，应有 - - p （x = i）= p, p （X = 0） =i P第一步构造似然函数L （Xi,X2,，X, P）= P （ X = Xi） P （ X = X2）P （ X = Xn）（1 -左厂（1 -戸厂卜气1弋厂=左倨丹严）（- ”严）第二步 计算In L （Xi,X2,X, P）并化简 =（Xi +Xn） Inp+ （ n- （Xi +Xn） In （i-p）第三步求二丄二 H册+花+ -八斗 囘一（焉+耳）1-刀Zj + -+ 嚣一（血+耳）驻点为 匚

6、-;化简为（x什+Xn） （ 1-p ） =p n- （ Xi +Xn）（Xi +Xn） =np驻点p=-（x因为只有一个驻点丁是最大点取一-例抽样n次A发生m次，则在Xi, X2xn中有m个1,其余为0,.P=_ 上例7 8 （1）设总体X服从泊松分布p U）,求兄的极大似然估计；（2）设总体X服从指数分布E （兄）求兄的极大似然估计10070107针对该题提问】 p （ X=k ） =! 从总体X中取样本X1 , X2Xn。血）=口尸（丙丿）=口刀3 =可）匕d:旷丄, 用=1dA AA = = x驻点 :InL（A）=迟再）血几-曲Tn（召!兀孟口 dlnL（A） *解得的极大似然估计1

7、 、=一2%7易知上的矩估计亦为： （2）v X E （二）第一步，从中取样本值 X1 , x2xn,应有X1 0, X2 0X n 0 似然函数L（ X1 , X2X n） = f（X1）f（X2 ）f（ Xn）=几巳跖）W巳空）（肮 第二步 计算: ： :丁第三步求討、二一厂 A= L驻点 叫亠- 是最大点S=i 取.：在例7-2中用矩法估计也是同样结果八J 3 X f（x） = & 例7 - 9设即 其它从中取样Xi , X2Xn,试用最大似然法求210070108针对该题提问】因为样本Xi , X2xn已经取出。所以应有OX W1,0三2,0rX所以。的取值范围为 m做 X. xa T

8、 - - Xfl） 0，很明显，似然函数- 丁-二: 是J的单调减函数，因此当 J最小时，似然函数二一.最大，由条件二知J的最小值为 -mass（忑,X2所以x 时-二二 J最大。取一心:厂-二、: 这一结果与用矩法估计（例 7-3）的结果不同。例7 10若-丨i -从中抽样x1, x2Xn,试用最大似然估计法求： ，-10070109针对该题提问】解: X的似然函数- i - 二 ： I （鬲-qHi-1=（加,戸exp 哙評TInZO- =-乔郭召-莎-守2 yln（2r）将丄.一宀|分别关于两个分量求偏导并令其为 0即得到似然方程组,（1）nL（t dr2） 1 & 卫 用 门 - .-

9、,（2）- 菲吃W解此方程组，由（1）可得驻点-，的极大似然估计为 ： 一将之代入（2）给出匚的极大似然估计i1 =7.2点估计的评价标准我们已经看到，点估计有各种不同的求法，为了在不同的点估计间进行比较选择， 就必须对各种点估计的好坏给出评价标准。数理统计中给出了众多的估计量评价标准， 对同一估计量使用不同的评价标准可能会得到完全不同的结论，因此，在评价某一个估计好坏时首先要说明是在哪一个标准下， 否则所论好坏毫无意义。但在诸多标准中，有一个基本标准是所有的估计都应该满足的， 它是衡量一个估计是否可行的必要条件，这就是估计的相合性，我们就从相合性开始介绍。7.2.1 相合性我们知道，点估计是

10、一个统计量，因此它是一个随机变量，在样本量一定的条件下，我 们不可能要求完全等同于参数的真实取值， 但如果我们有足够的观测值，根据格里纹科定理,随着样本量的不断增大，经验分布函数逼近真实分布函数， 因此完全可以要求估计量随着样本量的不断增大而逼近参数真值，这就是相合性，严格定义如下，定义7- 2 设J - ?为未知参数是J的一个估计量，n是样本容量，右对任何一个 0,有I J （ 7.2.1）则称一为参数J的相合估计相合性被认为是对估计的一个最基本要求， 如果一个估计量，在样本量不断增大时， 它都不能把被估参数估计到任意指定的精度， 那么这个估计是很值得怀疑的， 通常，不满足相合性要求的估计一

11、般不予考虑，证明古计的相合性一般可应用大数定律或直接由定义来证。例11用大数定律证明丁是的相合估计10070110针对该题提问】是；的相合估计为了避免用定义判断相合性的困难，下面介绍一个判断相合性很有用的定理:定量：设-J - 是的估计量若（ 1）鳥（2） 则一是J的相合估计。例12证明丁 是J的相合估计10070111针对该题提问】 证：在前面我们已经证明（1）- 7 n（s1）= z?y（x-x）（2）D迟（召-z） -0（ -co）iT 二是J的相合估计7.2.2 无偏性相合性是大样本下估计量的评价标准， 对小样本而言，需要一些其他的评价标准， 无偏性便是一个常用的评价标准。设是的一个估

12、计，的参数空间为口，若对任意的 口，有门了 - 则称丁是J的无偏估计，否则称为有偏估计。例7 -13对任一总体而言，样本均值是总体均值的无偏估计，当总体 k阶矩存在时,样本k阶原点矩s是总体k阶原点矩二的无偏估计，但对 k阶中心矩则不一样，例如，二阶样本中心矩就不是总体方差 /的无偏估计，事实上,10070112针对该题提问】对此，有如下两点说明（1）当样本量趋于无究时，有 我们称为的渐近无偏估计，这表明当样本量较大时，可近似看作J的无偏估计它比訂更常用，这是则是总体方差的无偏估计，这种简章的修正方法在一些场合常被采用,因为在n2时，且二 因此用上估计F有偏小的倾向，特别在小样本场合要使用 0

13、估计E3无偏性不具有不变性。即若 &是&的无偏估计，一般而言， g （）不是g （4）的无偏估计，除非g 3）是J的线性函数，例如，是J的无偏估计，但s不是h的无偏估计 例14证明近住+是口的无偏估计二也- -。其中6是X的样本10070113针对该题提问】证：匕=碍码 +务+敢.=.：.:- 一 - 行=:.-=-1工M - ; - - - +特别情形* - 丁是J的无偏估计 例15证明丁 是匚的无偏估计10070114针对该题提问】证二】；一二厂证 一 ：童（兀-xf = 乞彳-/）必；-恥 =左（士 +/）-超（3+（遲力，）z a（超-1好敲十+用/）-用（丄/ +站7.2.3 有效性

14、参数的无偏估计可以有很多， 那么如何在无偏估计中进行选择？直观的想法是希望该估 计围绕参数真值的波动越小越好， 波动的大小可以用方差来衡量， 因此人们常用无偏估计的方差的大小作为度量无偏估计优劣的标准，这就是有效性。定义7-4 设二，宀是j的两个无偏估计，如果对任意的 有则称囚比r有效例16设九呦是取自某总体的样本,记总体均值为.:，总体方差为/,则-. 丁都是二的无偏估计，但 显然，只要n 1, c比效，这表明，用全部数据的平均估计总体均值要比只使用部分数据更有效。亠 1 2 1 2 1 2 应再=&包曲十&再）=y目罚+ E愁-i- p = p 与门都是的无偏估计卫1 = Q（亍兀+憨j

15、= Q疋J + D（亏用）1 n 1 2 1 . 2 2 1 2=I - 丨 一7 1 7。卫a=询+亍花）=西）+。（彳花）0（WCO）-2-=-x 二是。7.3参数的区间估计用点估计去估计总体的参数， 即使是无偏且有效的， 也会由于样本的随机性， 使得从一个样本X1, X2, X3,，Xn算得的估计值不一定是被估计的参数的真实值，而且估计值的可 靠性并不知道，这是一个重大的问题，因此， 必须解决根据估计量的分布， 在一定可靠性的程度下指出被估计的总体参数的取值范围，这正是本节要介绍的参数的区间估计问题。7.3.1 置信区间概念为了引入置信区间的概念，请看下面的引例。引例 设某种绝缘子抗扭强

16、度 X服从正态分布 ；，其中未知，J已知（口 =45 公斤 米），试对总体均值作区间估计。10070201针对该题提问】对于区间估计，要选择一个合适的统计量，若在该总体取一个容量为 n的样本X1, X2,X3,，Xn,样本均值为-的点估计即：，然而我们要给出 门的一个区间估计，以体现出估丁叫）计的误差，我们知道 。在区间估计问题中，要选取一个合适的估计函数。这时，u_c可取 J人，它是：的标准化随机变量，且具备下面两个特点:（1） u中包含所要估计的未知参数 匸（其中匚已知）；（2） u的分布为N（0, 1），它与未知参数 启无关。因为uN （0, 1）,因而有Puu * = cj （0 cs

17、 1） 刃 ，根据uN （0, 1）的概率密度I的对称性（见下图）在引例中，若a=160，壬=40， n=16。则有- 40a-1.96 - = 160- 1.96x = 140.4 16+1.96 - = 160+1.96x41 = 179.6 血 716、（140.4 179.6） = 0.95说明该绝缘子抗扭强度 X的期望在（140.4, 179.6）内的可靠度为 0.95。下面，引出置信区间的概念。定义7-5设已为总体的未知参数: i _ j |是由样本 二定出的两个统计量，若对于给定的概率 1-a （0V aV 1）,有4 4 長- _匚 匚 二I为置信下限，已称为置信上限。者说，随

18、机区间以100 （1-a） %的概率包含旧。粗略地说，当 a =0.05时，在100次的抽样中，大致有 95次旧包含在L J中，而其余5次可能不在该区间中。a常取的数值为0.05, 0.01，此时置信度1-a分别为0.95, 0.99。置信区间的长度可视为区间估计的精度，下面分析置信度与精度的关系。n固定时，置信区间长度增大，即区间估计精度减低 置信区间长度减小，即区间估计精度提高。n增大时，置信区间减小（如引例中，置信区间长度% 丁为：），区间估计精度提高。7.3.2单个正态总体参数的置信区间正态总体-1 是最常见的分布，本小节中我们讨论它的两个参数的置信区间。1已知时的置信区间设总体X服从

19、正态分布3小,其中/已知，而已未知，求的置信度1-a的置信区间。 这一问题实际上已在引例中的讨论中解决，得到斗-盂“盂盂十肚。莘= 1迄2 U为 2 7就JP所以八的置信度1-a的置信区间为例1某车间生产滚珠，从长期实践知道，滚珠直径X服从正态分布。从某天产品里随 机抽取6个，测得直径为（单位：毫米）：14.6，15.1，14.9，14.8，15.2，15.1。若总体方差C72=0.06，求总体均值P的置信区间（a =0.05 a =0.01。10070202针对该题提问】解- 1，a =0.05时，置信度为95%的置信区间为14 75J5 15a =0.01时，置信度为99%的置信区间为从此

20、例知，在样本容量 n固定时，当置信度1- a较大时，置信区间长度较大；当置信度1-a较小时，置信区间较小。 _例2用天平称量某物体的质量 9次，得平均值为1=15.4 （g）,已知天平称量结果为正态分布，其标准差为 0.1g，试求该物体质量的 0.95置信区间。10070203针对该题提问】解 此处1-a =0.95 a =0.05查表知U0.025=1.96,于是该物体质量 的0.95的置信区间为早=15 4 1 96x1 = 15.4 0 0653从而该物体质量的 0.95置信区间为15.3347 , 15.4653。例3设总体为正态分布 上：打亠为得到匸的置信水平为0.95的置信区间长度

21、不超过 1.2,样本容量应为多大？10070204针对该题提问】解由题设条件知的0.95置信区间为量至少为11时才能使得,的置信水平为0.95的置信区间长度不超过 1.2。2匸未知时匸的置信区间解得P （-1）仝徑3-曲=1 -os7 2戸（一 3 1）且上企=一 E Pa 一 J 仗 一 1）莘兰 “ X+4 （料-1）卓）= l-fl宀-4（再-分2其中 1.-. 是.一的无偏估计。例4假设轮胎的寿命服从正态分布。为估计某种轮胎的平均寿命，现随机地抽 12只轮胎试用，测得它们的寿命（单位：万千米）如下：4.68 4.85 4.32 4.85 4.61 5.02 5.20 4.60 4.58

22、 4.72 4.38 4.70试求平均寿命的0.95置信区间。10070205针对该题提问】 _解此处正态总体标准差未知，可使用 t分布求均值的置信区间。本例中经计算有 A=4.7092 , s2=0.0615。取 a =0.05 查表知 to.025 （11 ） =2.2010,于是平均寿命的 0.95 置信区间 为（单位：万千米）的置信区间此时虽然也可以就是否已知分两种情况讨论 J的置信区间，但在实际问题中未知时 匸已知的情况是极为罕见的，所以我们只在 门未知的条件下讨论 J的置信区间。-厂中包含未知参数又它的分布与-J无关，以上作为估计函数，可用于-J的区间估 计。由于一厂分布是偏态分布

23、，寻找平均长度最短区间很难实现，一般都改为寻找等尾置信a区间：把a平分为两部分，在一r分布两侧各截面积为1的部分，即采用一r的的两个分位数 处帥1）和迁口-1）2 1_2尸C?、佥-1） = 1） = 1-它们满足 ： : - O（见下图）F1）-11_7 23-1）兰F莖丘心一D） = 1亠比 F 75孑守習十G.Pts 兰 y 0.0325=0.26，这里 a =0.95 查表知Zafl（8） = 2 1797,xt（8）=门 5345二 . 代入公式可得的0.95置信区间为从而F的0.95置信区间为0.1218 , 0.3454。以上关于正态总体参数的区间估计的讨论列表如表 7-1所示。表正态邑体掺魏的区间估计未条件佶计参数置信区闾戶【1-財卡7丰斗卡】Sf 0-皆 13 片本章小结本章考核要求为（1）点估计（1）知道点估计的概念（2）会用矩法求总体参数的矩估计值，主要依据是fl - X（3）