[理学]matlab课件--第5讲 概率统计实验.ppt

* Matlab 软件实习 第五讲 概率统计实验 一、古典概型 1. 阶乘 n! 的计算函数：factorial(n) 例1： 求8!。 调用 factorial(8) 计算得 40320 2. 排列 的计算函数：paily(n,r) 编辑paily.m文件： function y=paily(n,r) y=factorial(n)/factorial(n-r) 例2： 求在15个元素中取6个的排列。 程序如下： n=15;r=6; A=paily(n,r) 3. 组合 的计算函数：nchoosek(n,r) 编辑nchoosek.m文件： function y=nchoosek (n,r) y=paily(n,r)/factorial(r) 例3： 求在100个元素中任取6个的组合。 程序如下： A=nchoosek(100,6) 例4： 一个盒子中有10个产品,其中有7个正品3个次品,任取3个,求恰有1个是次品的概率。 程序如下： p=nchoosek(7,2)*nchoosek(3,1)/nchoosek(10,3) 4. 重要公式: 条件概率公式、乘法公式、全概率公式、贝叶斯公式 5. n重贝努力试验 计算函数：binopdf(k,n,p) 例5： 某人射击命中率为0.7,求其射击10次恰有4次命中的概率. 程序如下： p=binopdf(4,10,0.7) 二、随机数的产生 所有分布的随机数的产生方法都始于均匀分布的随机数.在MATLAB工具箱中提供了通用的随机数产生函数random和特定分布的随机数产生函数(以rnd结尾).可以直接调用这些函数来获得所需要的随机数. 统计工具箱中的随机产生函数及调用格式 R=frnd(V1,V2,m,n) frnd F-分布 R=chi2rnd(V,m,n) chi2rnd ?2-分布 R=trnd(V,m,n) trnd t-分布 R=lognrnd(MU,SIGMA,m,n) lognrnd 对数正态分布 R=normrnd(MU,SIGMA,m,n) normrnd 正态分布 R=exprnd(MU,m,n) exprnd 指数分布 R=unifrnd(N,m,n) unifrnd 连续均匀分布 R=poissrnd(LAMBDA,m,n) poissrnd 泊松分布 R=hygernd(M,K,N,m,n) hygernd 超几何分布 R=geornd(P,m,n) geornd 几何分布 R=binornd(N,P,m,n) binornd 二项分布 R=unidrnd(N,m,n) unidrnd 离散均匀分布 函数调用格式 函数名称 分布类型名称 1. random函数 功能：产生可选分布的随机数 调用格式：y=random(‘name’,A1,A2,m,n) 说明：random函数产生统计箱中任意分布的随机数.’name’为相应的分布的名称.A1,A2为分布参数.m,n为确定了运行结果y的行数与列数. 例6: 产生服从正态分布N(0,1)的2行4列的随机数. 程序如下： y=random(‘Normal’,0,1,2,4) 2. exprnd函数 例7: 产生4行5列的指数分布的随机数. 程序如下： y=exprnd(3,4,5) %参数?=3 三、随机变量与概率分布密度 1. 几个常用的离散型分布密度函数(…pdf ) (1)均匀分布 P(X=xn)=1/n 密度函数调用格式：y=unidpdf(X,N) 例8：求X取值为1,2,3,4,5,6,7,8时服从均匀分布的概率值. 程序如下: X=1:8,N=8; Y=unidpdf(X,N) (2)二项分布 P(X=k)= 密度函数调用格式：y=binopdf(X,N,P) 例9：求X取值为0,1,2,3,4,5,6,7,8时服从二项分布b(X;12,0.4)的概率值. 程序如下: X=0:8,N=12;P=0.4; Y=binopdf(X,N,P) (3)泊松分布 P(X=k)= 密度函数调用格式：y=poisspdf(X,LAMBDA) 例10:求X取值为0,1,2,3,4,5,6,7,8. ?=3时服从泊松分布的概率值. 程序如下: X=0:8; Y=poisspdf(X,3) (4)几何分布 P(X=k)=p(1