理学matlab概率统计实验.ppt

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