实验八概率与频率.ppt

《实验八概率与频率.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验八概率与频率1概率，又称几率，或然率，是反映某种事件发生的可能性大小的一种数量指标，它介于0与1之间。概率论是研究随机现象统计规律的一门数学分支学科，希望通过本实验的学习，能加深对频率和概率等概念的理解和认识，并掌握一些概率统计的基本原理。概率与频率随机现象中出现的某个可能结果问题背景和实验目的2内容提要基本知识Matlab相关函数介绍六个随机试验3基本知识随机试验：满足下列三个条件试验可以在相同的情况下重复进行；试验的所有可能结果是明确可知的，且不止一个；每次试验的结果无法预知，但有且只有一个结果。概率与频率概率是

2、指某个随机事件发生可能性的一个度量，是该随机事件本身的属性。频率是指某随机事件在随机试验中实际出现的次数与随机试验进行次数的比值。频率概率随机试验进行次数4Matlab随机函数rand(m,n)满足均匀分布的随机矩阵，每个元素都在(0,1)之间randn(m,n)满足正态分布的随机矩阵randi(p)生成一个1到p之间的随机整数(均匀分布)unidrnd(p)生成一个1到p之间的随机整数(均匀分布)randperm(n)生成一个由1到n组成的随机排列perms(1:n)生成由1到n组成的全排列，共n!个随机数/随机排列的

3、生成randi(p,m,n)unidrnd(p,m,n)5Matlab取整函数fix(x)截尾取整，直接将小数部分舍去floor(x)不超过x的最大整数ceil(x)不小于x的最小整数round(x)四舍五入取整6x1=fix(3.9);x2=fix(-3.9);x3=floor(3.9);x4=floor(-3.2);x5=ceil(3.1);x6=ceil(-3.9);x7=round(3.9);x8=round(-3.2);x9=round(-3.5);x1=3x2=-3x3=3x4=-4x5=4x6=-3x7=4

4、x8=-3x9=-4取整函数举例7unique(x)合并x中相同的项，并按从小到大排序若x是矩阵，则输出为一个列向量prod(X)如果X是向量，则返回其所有元素的乘积如果X是矩阵，则作用在矩阵的每一列上其它相关函数a=[129323];b=unique(a)a=[129;323];b=unique(a)factorial(n)%计算n!8根据表达式的不同取值，分别执行不同的语句switchexprcasecase1statements1casecase2statements2......casecasemstatemen

5、tsmotherwisestatementsendswitch选择语句9这里我们主要用rand函数和randi/unidrnd函数来模拟满足均匀分布的随机试验试验方法先设定进行试验的总次数采用循环结构，进行多次试验，并统计指定事件发生的次数计算指定事件发生次数与试验总次数的比值，即频率试验方法10随机投掷均匀硬币，验证国徽朝上与朝下的概率是否都是1/2n=10000;%设定试验总次数m=0;fori=1:ny=unidrnd(2);%随机产生1或2（均匀分布）ify==1%1表示国徽朝上，2表示国徽朝下m=m+1;en

6、dendfprintf('国徽朝上的频率为：%f',m/n);试验一：投掷硬币11随机投掷骰子，验证各点出现的概率是否为1/6n=10000;m1=0;m2=0;m3=0;m4=0;m5=0;m6=0;fori=1:ny=unidrnd(6);switchycase1,m1=m1+1;case2,m2=m2+1;case3,m3=m3+1;case4,m4=m4+1;case5,m5=m5+1;otherwise,m6=m6+1;endend...%输出结果试验二：投掷骰子12用蒙特卡罗(MonteCarlo)投点法

7、计算的值n=100000;a=2;m=0;fori=1:nx=rand(1)*a/2;y=rand(1)*a/2;if(x^2+y^2<=(a/2)^2)m=m+1;endendfprintf('计算出来的pi为：%f',4*m/n);试验三：蒙特卡罗法13投点区域的选择MonteCarlo注意点完全覆盖所求区域可以扩大投点区域怎样判断点落入所求区域内？14在画有许多间距为d的等距平行线的白纸上，随机投掷一根长为l(ld)的均匀直针，求针与平行线相交的概率，并计算的值试验四：蒲丰投针实验15n=100000;L

8、=0.5;d=1;m=0;fori=1:nalpha=rand(1)*pi;y=rand(1)*d/2;if(y<=L/2*sin(alpha))m=m+1;endendfprintf('针与平行线相交的频率为：%f',m/n);fprintf('计算出来的pi为：%f',2*n*L/(m*d));试验四源程