概率论基础知识和其在matlab中实现.docVIP

第三章 概率论基础知识及其在matlab中的实现 概率论与应用数理统计十就是研究随机现象统计规律性的一门科学，本章学习随机 及其发生的概率，多维随机变量的分布规律，参数估计与假设检验，方差分析与回归分析等概率统计的基本方法以及它们在MATLAB中实现方法。 3．1 随机时间及其概率 3．1．1 古典概率及其模型 由古典概率的定义知，古典概率基于这样两个原则：(1) 所有可能发生的结果只有有限个；(2)每一种可能出现的结果机会是相同的。在MATLAB中提供了一个在[0,1]区间上均匀分布的随机函数rand,其命令格式为： 命令格式1： rand(N) 功能： 返回一个的随机矩阵 命令格式2： rand(N,M) 功能： 返回一个的随机矩阵 命令格式2： rand(P1,P2,…,Pn) 功能： 返回一个的随机矩阵 可以用计算机模拟掷硬币这一过程，为了模拟硬币出现正面或反面，规定随机数小于0.5时为反面，否则为正面，可以用round（）函数将其变成0—1矩阵，然后将整个矩阵的各元素值加起来再除以总的原始个数即为出现正面的概率。round（）函数的命令格式为： 命令格式： round(x) 功能：对向量或矩阵x的每个分量四舍五入取整。 现以联系掷10000次硬币为例，重复做100次试验模拟出现正面的概率。在matlab中的程序如下： for i =1：100 a(i) = sum(sum(round(rand(100))))/10000； end mx = max(a)； mn = min(a)； ma = mean(a)； a, mx, mn, ma 在该程序输出的四项中，a为实验100次中每次出现正面的频率，mx和mn分别为100次实验中出现正面频率的最大值和最小值，ma为100次试验出现正面的平均频率。运行结果如下： 这里要输入结果 下面介绍MATLAB中取整的几个函数命令： (1) 命令格式： fix(x) 功能： 对x朝零方向取整 (2) 命令格式： floor(x) 功能：对x朝负无穷大方向取整 (3) 命令格式： ceil(x) 功能：对x朝正无穷大方向取整 3．1．2 统计概率及其模型 由于古典概率是建立在事件发生的等可能基础上的概率，而现实生活中许多现象的出现并不是等可能的。例如某品种的玉米，当种下一粒后，其发芽与不发芽的机会并不相同。那么，这个概率就不能建立在等可能基础上，即不能使用古典概率的定义。 而统计概率的定义是建立在频率基础上的，就是说某事件出现 频率如果稳定在某数值附近，则称数值为该事件出现的概率。由于统计概型中的概率是一个理论上的数值，实际问题中根本无法直接得到该数值，因而通常在试验次数充分多时，利用频率值近似代替概率值。在掷硬币的试验中，在试验次数充分多的情况下，掷硬币出现正面和反面的频率均在0.5左右，故出现正面和反面的概率均为 0.5。下面来看掷硬币时，当样本容量分别为 n = 10，100，1000，10000，100000，1 000 000 时频率的变化。在MATLAB中实现时程序代码如下： for I = i：6 a(i) = sum(round(rand(1,10-i)))/10^I； end 运算结果为： 这里输入运算结果 从上面运行的结果中可以看出，当样本容量不够大时，其频率的波动范围很大，即频率不够稳定，即使有时达到 0.5，但最大时已达到0.9，然而随着样本容量的增加，频率的波动范围越来越小，相差仅有左右。 3．1．3 条件概率、全概率公式与伯努利概率 若事件B的发生会影响到事件A的发生，则事件B发生的条件下，事件A发生的概率称为条件概率，条件的概率的计算公式为： 若事件A的发与事件B的发生与否没有关系，即事件B发生与否不会影响到事件A的发生，反之亦然，则称事件A与事件B是相互独立的，这时有： 例1 袋中有10只球，其中白球7只，黑球3只。分有放回和无放回两种情况，分三次取球，每次取一个，分别求：(1) 第三次摸到了黑球的概率，(2) 第三次才摸到黑球的概率，(3) 三次都摸到了黑球的概率。 解 当有放回地摸球时，由于三次摸球互不影响，因此三次摸球相互独立，从理论上可以求得：(1) 第三次摸到黑球的概率为；(2) 第三次才摸到黑球的概率为；(3) 三次都摸到黑球的概率为。 在MATLAB中模拟这一过程时，可在[0，1]区间上产生三次随机数来模拟三次摸球，当随机数小于0.7时可认为摸到了白球，否则认为摸到了黑球。重复次分别求上述三种情况出现的概率。程序如下： a = round(rand(1000000,3) – 0.2)； for i= 1：6