概率论基础知识及其在matlab中的实现.doc
《概率论基础知识及其在matlab中的实现.doc》由会员分享,可在线阅读,更多相关《概率论基础知识及其在matlab中的实现.doc(11页珍藏版)》请在咨信网上搜索。
1、_ 第三章 概率论基础知识及其在matlab中的实现 概率论与应用数理统计十就是研究随机现象统计规律性的一门科学,本章学习随机 及其发生的概率,多维随机变量的分布规律,参数估计与假设检验,方差分析与回归分析等概率统计的基本方法以及它们在MATLAB中实现方法。31 随机时间及其概率311 古典概率及其模型由古典概率的定义知,古典概率基于这样两个原则:(1) 所有可能发生的结果只有有限个;(2)每一种可能出现的结果机会是相同的。在MATLAB中提供了一个在0,1区间上均匀分布的随机函数rand,其命令格式为:命令格式1: rand(N)功能: 返回一个的随机矩阵命令格式2: rand(N,M)功
2、能: 返回一个的随机矩阵命令格式2: rand(P1,P2,Pn)功能: 返回一个的随机矩阵可以用计算机模拟掷硬币这一过程,为了模拟硬币出现正面或反面,规定随机数小于0.5时为反面,否则为正面,可以用round()函数将其变成01矩阵,然后将整个矩阵的各元素值加起来再除以总的原始个数即为出现正面的概率。round()函数的命令格式为:命令格式: round(x)功能:对向量或矩阵x的每个分量四舍五入取整。现以联系掷10000次硬币为例,重复做100次试验模拟出现正面的概率。在matlab中的程序如下:for i =1:100 a(i) = sum(sum(round(rand(100)/100
3、00;endmx = 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朝正无穷大方向取整312 统计概率及其模型由于古典概率是建立在事件发生的等可能基础上的概率,而现
4、实生活中许多现象的出现并不是等可能的。例如某品种的玉米,当种下一粒后,其发芽与不发芽的机会并不相同。那么,这个概率就不能建立在等可能基础上,即不能使用古典概率的定义。而统计概率的定义是建立在频率基础上的,就是说某事件出现 频率如果稳定在某数值附近,则称数值为该事件出现的概率。由于统计概型中的概率是一个理论上的数值,实际问题中根本无法直接得到该数值,因而通常在试验次数充分多时,利用频率值近似代替概率值。在掷硬币的试验中,在试验次数充分多的情况下,掷硬币出现正面和反面的频率均在0.5左右,故出现正面和反面的概率均为 0.5。下面来看掷硬币时,当样本容量分别为 n = 10,100,1000,100
5、00,100000,1 000 000时频率的变化。在MATLAB中实现时程序代码如下:for I = i:6 a(i) = sum(round(rand(1,10-i)/10I;end运算结果为:这里输入运算结果从上面运行的结果中可以看出,当样本容量不够大时,其频率的波动范围很大,即频率不够稳定,即使有时达到 0.5,但最大时已达到0.9,然而随着样本容量的增加,频率的波动范围越来越小,相差仅有左右。313 条件概率、全概率公式与伯努利概率若事件B的发生会影响到事件A的发生,则事件B发生的条件下,事件A发生的概率称为条件概率,条件的概率的计算公式为: 若事件A的发与事件B的发生与否没有关系,
6、即事件B发生与否不会影响到事件A的发生,反之亦然,则称事件A与事件B是相互独立的,这时有: 例1 袋中有10只球,其中白球7只,黑球3只。分有放回和无放回两种情况,分三次取球,每次取一个,分别求:(1) 第三次摸到了黑球的概率,(2) 第三次才摸到黑球的概率,(3) 三次都摸到了黑球的概率。解 当有放回地摸球时,由于三次摸球互不影响,因此三次摸球相互独立,从理论上可以求得:(1) 第三次摸到黑球的概率为;(2) 第三次才摸到黑球的概率为;(3) 三次都摸到黑球的概率为。在MATLAB中模拟这一过程时,可在0,1区间上产生三次随机数来模拟三次摸球,当随机数小于0.7时可认为摸到了白球,否则认为摸
7、到了黑球。重复次分别求上述三种情况出现的概率。程序如下: a = round(rand(1000000,3) 0.2); for i= 1:6 b = a(1:10i,3; c(i) =sum(b)/(10i); end c for i =1:6 b = (a(1:10i,1)&(a(1:10i,2)&a(1:10i,3); d(i) = sum(b)/(10i); end d for I = 1:6 b=a(1:10i,1)&a(1:10i,2)&a(1:10i,3); e(i) = sum(b)/(10i); end e运行结果为:这里加入运行结果执行结果中c为第三次摸到黑球的概率,d为第
8、三次才摸到黑球的概率,e为三次都摸到黑球的概率.可以看到,随着试验次数的增加,其频率都会逐渐稳定在理论值附近.当无放回地摸球时,由于第二次摸球会受到第一次的影响,而第三次摸球又会受到前两次的影响,因而三次摸球相互影响,并不独立.从理论上可求得:(1) 第三次摸到黑球的概率为,(2) 第三次才摸到黑球的概率为(3) 三次都摸到了黑球的概率为 用计算机模拟该过程时,在0,1区间模拟第一次摸球,当值小于0.7时认为摸到了白球,否则认为摸到了黑球;第二次摸球时由于少了一个球,故可在区间长度为0.9的区间上模拟,若第一次摸到白球,可将区间设为0.1,1,否则区间设为0,0.9;第三次摸球可依次类推,其模
9、拟程序如下:a = rand(1000000,3);a(:,1) = round(a(:,1) 0.2);a(:,2) =round(a(:,2)*0.9 0.2 0.1*(a:,1) 1);a(:,3) =round(a(:,3)*0.8 0.2 0.1*(a(:,1) 1) 0.1*(a(:,2) 1); for i=1:6 b =a(1:10i,3); c(i) = sum(b)/(10i); end c for i =1:6 b = (a(1:10i,1) & (a(1:10i,2) & a(1:10i,3); d(i) =sum(b)/(10i); end d for i =1:6
10、b =a(1:10i,1) & a(1:10i,2) & a(1:10i,3); e(i) =sum(b)/(10i); end e运行结果为上面在理论上计算第三次摸到黑球的概率时,用到了全概率公式:若构成一个完备事件组,且事件的发生总是伴随着事件中的某一个发生而发生,则下面将用到伯努利概型,所谓伯努利概型是指:在相同条件下,进行次独立重复试验,每次试验只有事件A发生或不发生两种结果,且 这里第三次摸到黑球的四种情况分别是:白,白,黑,白,黑,黑,黑,白,黑,黑,黑,黑.这四种情况构成了完备事件组.现考虑下面问题: (1) 当不放回时,已知第三次摸到了黑球,问前两次是黑球的概率为多少? (2)
11、 若有放回地连续摸10次,则恰有三次摸到黑球的概率是多少?第一问是一逆概率问题,由逆概率公式即贝叶斯公式得到其概率应为第二问则属伯努利概型,这里A为摸到的是黑球,故 , .于是由二项概率公式有,10次有放回摸球中,恰有三次摸到黑球的概率为在MATLAB中实现这两个过程的程序如下: a = rand(100000,3); a(:,1) = round(a(:,1) 0.2); a(:,2) = round(a(:,2)*0.9 0.2 0.1*a(:,1) 1); a(:,3) = round(a(:,3)*0.8 0.2 0.1*(a(:,1) 1) 0.1*(a(:,2) 1); for i
12、=1:6 b =a(1:10I,3); c(i) =sum(b); b =a(1:10i,1) & a(1:10i,2) & a(1:10i,3); d(i) = sum(b); e(i) =d(i)./c(i); end 这里添加运行结果 a = round(rand(1000000,10) 0.2); for I =1:6 b =sum(a(1:10i, :),2) 3; c(i) =sum(-b)/910I); end c3.2 随机变量的分布及其数字特征随机变量的统计行为完全决定于其概率分布,按随机变量的取值不同,通常可将其分为离散型连续型和奇异型三大类.由于奇异型在实际应用中很少遇到
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 概率论 基础知识 及其 matlab 中的 实现
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。