基于MATLAB概率统计数值实验.pptVIP

基于MATLAB的概率统计数值实验 一、古典概型 主讲教师：董庆宽 副教授 研究方向：密码学与信息安全 电子邮件：qkdong@ 个人主页：/qkdong/ 内容介绍 一、古典概型 MATLAB常用的及与随机数产生相关的函数 实验1：计算超几何分布 实验2：频率稳定性实验 实验3：利用频率估计自然对数底e 实验4：蒲丰投针实验，利用频率估计圆周率? 实验5：生日悖论实验 */21 一、古典概型 利用MATLAB 软件的图形可视功能将概率统计的内容用图形表示出来，以加深对概率的理解 MATLAB常用的及与随机数产生相关的函数 factorial(n) ：阶乘，n!，可通过阶乘来计算排列组合数 1.rand(m,n)：生成m×n的随机矩阵，每个元素都在(0,1)间，生成方式为均匀分布。 2.randn(m,n)：生成m×n的随机矩阵，每个元素都在(0,1)间，生成方式为正态分布 3.randperm(m)：生成一个1~m的随机整数排列 4.perms(1:n)：生成一个1~n的全排列，共n!个 5.取整函数系列：（1）fix(x)：截尾法取整；（2）floor(x)：退一法取整（不超过x的最大整数）；向负方向舍入（3）ceil(x)：进一法取整（= floor(x)+1）； 向正方向舍入（4）round(x)：四舍五入法取整。 6.unique(a)：合并a中相同的项 7.prod(x)：向量x的所有分量元素的积 */21 一、古典概型 示例： rand(1) %生成一个(0,1)间的随机数 ans = 0.8147 rand(2,2) %生成一个2×2阶(0,1)间的随机数矩阵ans = 0.9134 0.0975 0.6324 0.2785 randperm(5) %生成一个1~5的随机整数排列ans = 4 1 5 2 3 a=[1 2 4 2 3 3 2];unique(a)ans = 1 2 3 4 */21 实验 1：计算超几何分布的结果 设有N件产品，其中D件次品，今从中任取n件 问其中恰有k(k?D)件次品的概率是多少？ （令N=10，D=3，n=4，k=2） 解：编辑组合函数zuhe.m文件 function y=Com(n,r) y=factorial(n)/(factorial(r)*factorial(n-r)) 计算如下： N=10; D=3; n=4; k=2; p=Com(3,2)*Com(10-3,4-2)/Com(10,4)=0.3 */21 实验2 频率稳定性实验 随机投掷均匀硬币，观察国徽朝上与国徽朝下的频率 解 n= 3000~100000000;m=0; for i=1:n t=randperm(2); %生成一个1~2的随机整数排列 x=t-1; %生成一个0~1的随机整数排列 y=x(1); %取x排列的第一个值 if y==0; m=m+1; end end p1=m/n p2=1-p1 */21 试验次数n 3000 5000 1万 2万 3万 国徽朝上频率 0.5040 0.5006 0.4879 0.4999 0.5046 国徽朝下频率 0.4960 0.4994 0.5121 0.5001 0.4954 试验次数n 5万 10万 100万 100万 1亿 国徽朝上频率 0.5021 0.4999 0.4999 0.5001 0.5000 国徽朝下频率 0.4979 0.5001 0.5001 0.4999 0.5000 可见当 时， */21 实验3 用频率估计自然对数e 某班有n个人，每人各有一支枪，这些枪外形一样。某次夜间紧急集合，若每人随机地取走一支枪，求没有一个人拿到自己枪的概率？ 解：记事件Ai为第i个人拿到自已枪，事件 为第i个人没拿到自己枪，易知： ； ， 又记 p0为没有一个人拿到自己枪的概率。 由乘法公式可知 … … */21 于是 所以  特别地，当n较大时， 。 因此，可随机模拟出没有人拿到自己枪的频率，根据频率的稳定性，近似当做概率，然后去估计自然对数e。并考虑估计精度与人数是否有关系，为什么。算法如下： */21 1、产生n个随机数的随机序列； 2、检验随机列与自然列是否至少有一个配对； 3、对没有一个配对的序列进行累积 t； 4、