mathematics教程第9章概率統(tǒng)計計算課件

第九章概率統(tǒng)計計算

北京交通大學9.1概率統(tǒng)計軟件包Mathematica可以處理概率統(tǒng)計方面的計算,有關的命令都在Mathematica自帶的統(tǒng)計軟件包中,這些軟件包存放在Mathematica系統(tǒng)自己帶有程序包,存放在C:\wnmath22\Packges\Statisti目錄中,用戶可以在Mathematica的工作窗口鍵入Ctrl+O,調出Open窗口,將該窗口左下腳的文件類型選為Packages(*.m),并用鼠標雙擊文件夾packages打開其中的子文件夾,然后任意雙擊Statisti文件夾,就可以在窗口左上部分看到很多以.m為擴展名的Mathematica所有自帶的概率統(tǒng)計軟件包文件:(見圖)下一頁返回9.2Mathematica概率統(tǒng)計軟件包中最常用的命令為了使用的方便,下面寫出一些概率統(tǒng)計軟件包中最常用的內容及其調用文件名需調用Statistics`DescriptiveStatistics`軟件包才能使用的函數(shù):Mean[data]計算樣本數(shù)據(jù)data的均值Median[data]計算樣本數(shù)據(jù)data的中值Variance[data]計算樣本數(shù)據(jù)data的方差StandardDeviation[data]計算樣本數(shù)據(jù)data的標準差

注意:

data是由離散數(shù)據(jù)組成的表例1：1)已知樣本數(shù)據(jù)為dat={3.2,5.1,1,4,2},試計算dat的均值、中值、方差、標準差。2)產生[0，1]上的20個隨機實數(shù)，并計算它們的均值、中值、方差、標準差。解:In[1]:=<<Statistics`DescriptiveStatistics`*調用統(tǒng)計軟件包In[2]：=dat={3.2,5.1,1,4,2}；In[3]:=Mean[dat]Out[3]:=3.06In[4]:=Median[dat]Out[4]:=3.2In[5]:=Variance[dat]Out[5]:=2.608In[6]:=StandardDeviation[dat]Out[6]:=1.61493In[7]:=dat1=Table[Random[],{20}]Out[7]:={0.93234,0.439331,0.407442,0.469035,0.741679,0.884562,0.111029,0.696056,0.0591917,0.622276,0.825287,0.540449,0.594691,0.597846,0.490196,0.463414,0.404672,0.19069,0.105273,0.942455}In[8]:=Mean[dat1]Out[8]:=0.525896In[9]:=Median[dat1]Out[9]:=0.515323In[10]:=Variance[dat1]Out[10]:=0.0724088In[11]:=StandardDeviation[dat1]Out[11]:=0.269089需調用Statistics`DiscreteDistributions`軟件包才能使用的概率分布和函數(shù):

BernoulliDistribution[p]表示均值為p的離散伯努力分布BinomialDistribution[n,p]表示參數(shù)為n,p的二項分布b(n,p)GeometricDistribution[p]表示參數(shù)為p的幾何分布HypergeometricDistribution[n,nsucc,ntot]表示參數(shù)為n,nsucc,ntot的超幾何分布PoissonDistribution[mu]表示參數(shù)為mu的F泊松分布PDF[distribution,k]離散分布distribution的分布律P{=k}CDF[distribution,x]概率分布為distribution且隨機變量小于值x的概率P{<x}Mean[distribution]計算離散分布distribution的均值Variance[distribution]計算離散分布distribution的方差StandardDeviation[distribution]計算離散分布distribution的標準差Random[distribution]產生具有概率分布為distribution一個偽隨機數(shù)例5：假設投擲一個均勻硬幣只能出現(xiàn)正面和反面兩種情況,用Mathematica命令來驗證投擲出現(xiàn)正面的概率為0.5。解：設X表示投擲一個均勻硬幣出現(xiàn)正面和反面的隨機變量，它只取兩個值0和1,采用具有概率分布均值為0.5的離散伯努力分布BernoulliDistribution[0.5]產生的偽隨機數(shù)Random[BernoulliDistribution[0.5]]來模擬實際投擲一個均勻硬幣的情況，規(guī)定出現(xiàn)隨機數(shù)是1表示投擲硬幣出現(xiàn)正面;0表示投擲硬幣出現(xiàn)反面。命令中分別用產生的100個偽隨機數(shù)、500個偽隨機數(shù)和1000個偽隨機數(shù)出現(xiàn)數(shù)1的頻率來驗證投擲出現(xiàn)正面的概率為0.5的結論，命令為：In[1]:=<<Statistics`DiscreteDistributions`*調用統(tǒng)計軟件包In[2]：=sy[n_]:=Module[{face,s},*定義模擬函數(shù)s=BernoulliDistribution[0.5];For[face=0;i=1,i<=n,i=i+1,If[Random[s]==1,face=face+1]];N[face/n]]

In[3]={sy[100],sy[500],sy[1000]}Out[3]={0.53,0.514,0.472}

從模擬試驗結果可以看到投擲出現(xiàn)正面的概率在0.5附近波動。需調用Statistics`ContinuousDistributions`軟件包才能使用的概率分布和函數(shù)BetaDistribution[,]表示參數(shù)為和的Beta連續(xù)分布CauchyDistribution[,]表示參數(shù)和的柯西連續(xù)分布ChiSquareDistribution[n]表示有n個自由度的2連續(xù)分布ExponentialDistribution[lambda]表示參數(shù)為的指數(shù)連續(xù)分布"FRatioDistribution[n1,n2]表示分子參數(shù)為n1和分母參數(shù)為n2的F連續(xù)分布NormalDistribution[,]表示均值為標準差為的正態(tài)分布N(,2)RayleighDistribution[]表示參數(shù)為的瑞利連續(xù)分布"StudentTDistribution[n]表示有n個自由度的t連續(xù)分布UniformDistribution[min,max]表示[min,max]區(qū)間上的均勻分布PDF[distribution,x]概率分布為distribution的分布密度函數(shù)f(x)CDF[distribution,x]概率分布為distribution且隨機變量小于值x的概率P{<x}Mean[distribution]計算概率分布為distribution均值Variance[distribution]計算概率分布為distribution方差StandardDeviation[distribution]計算概率分布為distribution標準差Random[distribution]產生具有概率分布為distribution一個偽隨機數(shù)例3：設隨機變量服從正態(tài)分布N(0,32),（１）求出對應的分布密度函數(shù),并畫出對應的分布密度函數(shù)圖形（２）求隨機變量<2的概率解:Mathematica命令為:In[1]:=<<Statistics`ContinuousDistributions`In[2]:=dis=NormalDistribution[0,3]Out[2]:=NormalDistribution[0,3]In[3]:=PDF[dis,x]1Out[3]=-------------------x2/183ESqrt[2Pi]

In[4]:=Plot[PDF[dis,x],{x,-10,10},PlotRange->All]Out[4]:=-Graphics-In[5]:=CDF[dis,2]*求隨機變量<2的概率Out[5]=0.747507實驗3生成自由度為12的t分布的連續(xù)型隨機變量及其概率密度函數(shù)，分布函數(shù)，并用圖形顯示。Mathematica命令In[1]:=<<Statistics`ContinuousDistributions`In[2]:=rv=StudentTDistribution[12];In[3]:=f=PDF[rv,x]Out[3]:=(*t(12)的概率密度函數(shù)*)

In[4]:=Plot[f,{x,-5,5}]In[5]:=g=CDF[rv,x];Out[5]:=(*t(12)的分布函數(shù)*)In[6]:=Plot[g,{x,-4,4}]實驗4某地區(qū)18歲女青年的血壓（收縮壓，以mm-Hg計）服從N(110,122)。在該地區(qū)任選一個18歲的女青年，測量她的血壓X。求P(X≤105)和P(100<X≤120)，畫出血壓X概率密度函數(shù)的圖像。Mathematica命令In[1]:=<<Statistics`ContinuousDistributions` rv=NormalDistribution[110，12]; f[x_]:=PDF[rv，x]; F[x_]:=CDF[rv，x]; N[F[105]]Out[5]=0.338461(*P(X105)=0.338461*)In[6]:=N[F[120]-F[100]]Out[6]=0.595343(*P(100<X100)=0.595343*)In[7]:=Plot[f[x],{x,110-12*3,110+12*3}]實驗5設隨機變量X-b(20-0.4),計算（1）P{X=0}（2）P{X=1}（3）P{X<2}（4）P{X≤6}（5）P{X>10}（6）P{X≥15}Mathematica命令

In[1]:=<<Statistics`DiscreteDistributions`brv=BinomialDistribution[20,0.4];f[x_]:=PDF[brv,x];df[x_]:=CDF[brv,x];f[0]Out[5]=0.0000365616(*得P(=0)=0.0000365616*)In[6]:=f[1]Out[6]=0.000487488(*得P(=1)=0.000487488*)In[7]:=df[2]-f[2]Out[7]=0.000524049(*得P(<2)=0.000524049*)In[8]:=df[6]Out[8]=0.250011(*得P(

6)=0.250011*)In[9]:=1-df[10]Out[9]=0.127521(*得P(>10)=0.127521*)In[10]:=1-df[15]+f[15]Out[10]=0.00161152(*得P(

15)=0.00161152*)實驗7給出20個服從均值為0、標準差為3的正態(tài)分布N(0,32)隨機數(shù)組成的表Mathematica命令

In[1]:=<<Statistics`ContinuousDistributions`In[2]:=rv=NormalDistribution[0,3];RandomArray[rv,20]Out[3]={0.636589,-4.25557,2.04924,1.58478,0.0244065,0.371864,-0.933664,3.54688,-0.888601,-0.650029,-2.49356,-3.07764,-2.44536,-0.512286,-1.68181,3.8912,-4.28302,-2.01939,-0.294215,2.13797}實驗8n個人每人攜帶一件禮物參加聯(lián)歡會。聯(lián)歡會開始后，先把所有的禮物編號，然后每人任意抽取一個號碼，按號碼領取禮物。請分別就參加聯(lián)歡會的人數(shù)n=1到20人求所有人都得到別人贈送禮物的概率，并從這些概率值推斷隨著參加聯(lián)歡會的人數(shù)增加是否會出現(xiàn)所有人都得到別人贈送禮物的概率會不斷變小的情況？Mathematica命令

In[1]:=p[n_]:=Sum[(-1)^k*1/k!,{k,2,n}]In[2]:=Table[N[p[k],18],{k,1,20}]Out[2]={0,0.500000000000000000,0.333333333333333333,0.375000000000000000,0.366666666666666667,0.368055555555555556,0.367857142857142857,0.367881944444444444,0.367879188712522046,0.367879464285714286,0.367879439233605900,0.367879441321281599,0.367879441160691161,0.367879441172161906,0.367879441171397190,0.367879441171444985,0.367879441171442173,0.367879441171442329,0.367879441171442321,0.367879441171442322}從計算結果可以看到，隨著參會人數(shù)的增加，所有人都得到別人贈送禮物的概率不會不斷變小，而是會收斂到一個約為0.367879,也就是e-1。實驗9在某紡織廠中，一個工人要照顧800個紗錠。每個紗錠旋轉時，由于偶然的原因，紗會被扯斷。假設在某一段時間內，每個紗錠的紗被扯斷的概率為0.005，求在這段時間內，紗被扯斷次數(shù)不大于10的概率。分析：相當于進行800次獨立試驗，用X表示紗被扯斷次數(shù)，則有X服從b(800,0.005)的二項分布，而所求概率為P{X≤10}可以用求b(800,0.005)的分布函數(shù)得到。Mathematica命令

In[1]:=<<Statistics`DiscreteDistributions`In[2]:=rvb=BinomialDistribution[800,0.005];f[k_]:=CDF[rvb,k]f[10]Out[4]=0.997239所以在這段時間內，紗被扯斷次數(shù)不大于10的概率為0.997239。實驗10設樣本數(shù)據(jù)為{110.1,25.2,50.5,50.5,55.7,30.2,35.4,30.2,4.9,32.3,50.5,30.5,32.3,74.2,60.8}

求該樣本的均值、方差、標準差、中位數(shù)、眾數(shù)。Mathematica命令

In[1]:=<<Statistics`DescriptiveStatistics`In[2]:=d1={110.1,25.2,50.5,50.5,55.7,30.2,35.4,30.2,4.9,32.3,50.5,30.5,32.3,74.2,60.8};In[3]:=Mean[d1]Out[3]=44.8867(*均值為44.8867*)In[4]:=var=Variance[d1]Out[4]=614.89(*方差為614.89*)In[5]:=Sqrt[var]Out[5]=24.797(*標準差為24.797*)In[6]:=Median[d1]Out[6]:=35.4(*中位數(shù)為35.4*)In[7]:=Mode[d1]Out[7]:=50.5(*眾數(shù)為50.5*)實驗11設樣本數(shù)據(jù)為{16.5,13.8,16.6,15.7,16.0,16.4,15.3}，求該樣本的均值、幾何均值和調和均值。Mathematica命令

In[1]:=<<Statistics`DescriptiveStatistics`In[2]:=d2={16.5,13.8,16.6,15.7,16.0,16.4,15.3};In[3]:=Mean[d2]Out[3]=15.7571(*均值為15.7571*)In[4]:=GeometricMean[d2]Out[4]=15.7296(*幾何均值為15.7296*)In[5]:=HarmonicMean[d2]Out[5]=15.7007(*調和均值為15.7007*)實驗12設樣本數(shù)據(jù)為{6.5,3.8,6.6,5.7,6.0,6.4,5.3}，畫出該樣本的條形圖和餅形圖。Mathematica命令

In[1]:=<<Graphics`Graphics`In[2]:=d3={6.5,3.8,6.6,5.7,6.0,6.4,5.3};In[3]:=BarChart[d3](*畫樣本條形圖10-1*)In[4]:=PieChat[d3](*畫樣本餅形圖10-2*)實驗14設有如下30個樣本數(shù)據(jù)

0.192,-1.382,0.508,-0.813,0.531,-0.536,0.826,1.404,-1.372,-0.349,1.054,1.372,1.624,0.709,1.034,1.670,-0.205,-0.017,-0.204,0.056,-1.179,-0.645,1.201,0.453,0.304,-1.832,0.058,1.870,0.912,-1.769

（1）畫出具有10個等距子區(qū)間的直方圖；（2）畫出具有16個等距子區(qū)間的直方圖。（1）Mathematica命令

In[1]:=d={0.192,-1.382,0.508,-0.813,0.531,-0.536,0.826,1.404,-1.372,-0.349,1.054,1.372,1.624,0.709,1.034,1.670,-0.205,-0.017,-0.204,0.056,-1.179,-0.645,1.201,0.453,0.304,-1.832,0.058,1.870,0.912,-1.769}In[2]:=h1=-2;(*左端點值為-2*)h=(2-(-2))/10;(*小區(qū)間的長度為0.4*)num=10;(*小區(qū)間的個數(shù)為10*)dnum=30；(*樣本容量為30*)fpn=Table[m=0;Do[If[d[[i]]>=h1&&d[[i]]<h1+h,m=m+1],{i,1,dnum}];h1=h1+h;m,{k,1,num}]Out[2]={2,2,2,2,4,4,4,4,3,3}In[3]:=<<Graphics`Graphics`In[4]:=h1=-2;GeneralizedBarChart[Table[{h1+h*i,fpn[[i+1]],h},{i,0,num-1}],PlotRange->{0,5}]（2）Mathematica命令In[1]:=d={0.192,-1.382,0.508,-0.813,0.531,-0.536,0.826,1.404,-1.372,-0.349,1.054,1.372,1.624,0.709,1.034,1.670,-0.205,-0.017,-0.204,0.056,-1.179,-0.645,1.201,0.453,0.304,-1.832,0.058,1.870,0.912,-1.769}

In[2]:=h1=-2;(*左端點值為-2*)

h=(2-(-2))/16;(*小區(qū)間的長度為0.25*)num=16;(*小區(qū)間的個數(shù)為16*)

dnum=30；