《新编MATLAB&Simulink自学一本通》课件第16章 概率分布与随机数.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《《新编MATLAB&Simulink自学一本通》课件第16章 概率分布与随机数.ppt》由用户(momomo)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 新编MATLAB&Simulink自学一本通 新编MATLAB&Simulink自学一本通课件第16章 概率分布与随机数 新编 MATLAB Simulink 自学 课件 16 概率 分布 随机数
- 资源描述:
-
1、2023-5-5 概率分布概率分布 生成一元分布随机数生成一元分布随机数 生成多元分布随机数生成多元分布随机数主要内容:主要内容:2023-5-5第一节第一节 概率分布概率分布2023-5-5一、随机变量的分布函数一、随机变量的分布函数设设X是一随机变量,对任意的实数,称是一随机变量,对任意的实数,称()()F xP Xx为随机变量为随机变量X的的分布函数分布函数。1.1.定义定义2.2.性质性质 单调性:单调非降单调性:单调非降 有界性:有界性:右连续性右连续性0()1F xlim()0,lim()1xxF xF x2023-5-5二、离散随机变量的概率函数(或分布列,分布律)二、离散随机变
2、量的概率函数(或分布列,分布律)设设X是一随机变量,称是一随机变量,称X取可能值的概率取可能值的概率()(),1,2,iiipp xP Xxin为随机变量为随机变量X的的概率函数(或分布列,分布律)概率函数(或分布列,分布律)。列。列表如下:表如下:1.1.定义定义2023-5-52.2.分布律性质分布律性质 非负性:非负性:()0,1,2,ip xi 正则性:正则性:1()1iip x3.3.分布律与分布函数的关系分布律与分布函数的关系()()iixxF xp x【例例16.1-001】由随机变量的分布律求分布函数。由随机变量的分布律求分布函数。1230.250.50.25XP0,10.25
3、,12()0.75,231,3xxF xxx ()F xxo123112023-5-5三、连续随机变量的密度函数三、连续随机变量的密度函数 设随机变量设随机变量X的分布函数为的分布函数为F(x),若存在非负可积函数,若存在非负可积函数f(x),使得对任意实数,使得对任意实数x有有()()xF xf x dx则称则称X为为连续随机变量连续随机变量,f(x)为为X的的密度函数密度函数。1.1.连续随机变量及密度函数定义连续随机变量及密度函数定义2023-5-52.2.密度函数性质密度函数性质 非负性:非负性:正则性:正则性:()0f x()1f x dx3.3.利用分布函数或密度函数求概率(连续随
4、机变量)利用分布函数或密度函数求概率(连续随机变量)()()()()P aXbP aXbP aXbP aXb()()()baF bF af x dx2023-5-5【例例16.1-002】已知随机变量的密度为已知随机变量的密度为34,01()0,xxf x其他3400,0()4,011,1xxF xx dxxxx求(求(1)概率)概率P(0 x 0.5);(2)分布函数)分布函数F(x)。解:解:0.530(00.5)4116PXx dx2023-5-5四、数学期望的定义四、数学期望的定义1.1.离散随机变量情形离散随机变量情形 设离散随机变量设离散随机变量X的分布律为的分布律为p(xi),i
5、=1,2,,则,则X的数的数学期望为学期望为()()iiiE Xx p x2.2.连续随机变量情形连续随机变量情形 设连续随机变量设连续随机变量X的密度函数为的密度函数为f(x),则,则X的数学期望为的数学期望为()()E Xxf x dx2023-5-5【例例16.1-003】航班每次飞行坠机概率为十万分之一,每位航班每次飞行坠机概率为十万分之一,每位乘客保费为乘客保费为20元,死亡赔付金额为元,死亡赔付金额为40万。问保险公司从每位万。问保险公司从每位顾客手中平均获取多大利润顾客手中平均获取多大利润解:令解:令 ,用,用Y 表示保险公司从一位顾客手中表示保险公司从一位顾客手中获取的利润。则
6、获取的利润。则X 和和Y 的分布律为的分布律为 0,1,X平安坠机 01 0.99999 0.00001 2020400000XPY从而可得从而可得Y 的期望为:的期望为:()20 0.99999(20400000)0.0000116E Y 2023-5-5五、方差的定义五、方差的定义212()(),var()()(),iiixE Xp xXxE Xf x dx离散情形连续情形2.2.离散和连续情形离散和连续情形1.1.定义定义222var()()()()XEXE XE XE X方差用来描述随机变量取值的波动(集中与分散)程度方差用来描述随机变量取值的波动(集中与分散)程度2023-5-5【例
7、例16.1-004】在在M电子公司生产的简易二极管中,按质量电子公司生产的简易二极管中,按质量等级可分为等级可分为5级,其中级,其中1级最差,级最差,5级最好。现统计了今年级最好。现统计了今年1月月份生产的二极管质量各等级所占比率,如下表所列。求平均份生产的二极管质量各等级所占比率,如下表所列。求平均质量等级和质量等级的方差。质量等级和质量等级的方差。解:平均质量等级解:平均质量等级质量等级的方差质量等级的方差()1 0.12 0.23 0.34 0.35 0.13.1E X X12345P0.10.20.30.30.12var()1 0.14 0.29 0.3 16 0.325 0.1 3.
8、11.29X 2023-5-5六、六、0-10-1分布分布1.1.定义定义 抛一枚硬币一次,用抛一枚硬币一次,用X表示出现正面的次数,则表示出现正面的次数,则X所服所服从的分布就是从的分布就是0-1分布分布(两点分布两点分布)。)。两点分布的分布律为两点分布的分布律为X01P1-pp其中其中0 p 1.2023-5-5七、二项分布七、二项分布1.1.定义定义 抛一枚硬币抛一枚硬币n次,用次,用X表示出现正面的次数,则表示出现正面的次数,则X所服从所服从的分布就是的分布就是二项分布二项分布。2023-5-52.2.实例实例【实例实例1】一袋中有一袋中有N个大小形状相同的球,其中有个大小形状相同的
9、球,其中有M个个白球,从中白球,从中有放回抽取有放回抽取n个球,记个球,记X为取到的白球数,为取到的白球数,X服从的分布即为二项分布服从的分布即为二项分布B(n,M/N)。【实例实例2】在在n次独立试验中,若每次只有次独立试验中,若每次只有“成功成功”和和“失败失败”两种结果,且每次成功概率均为两种结果,且每次成功概率均为p,则,则n次试验次试验中成功次数服中成功次数服X从二项分布从二项分布B(n,p)。2023-5-5八、泊松分布八、泊松分布1.1.定义定义 泊松分布是常见的,例如中午时分,每分钟进入肯德基泊松分布是常见的,例如中午时分,每分钟进入肯德基的顾客数;一定时间内接错电话的次数;一
10、个铸件上的缺的顾客数;一定时间内接错电话的次数;一个铸件上的缺陷数;一平方米玻璃上的气泡数;一页书上的错字书等陷数;一平方米玻璃上的气泡数;一页书上的错字书等.2023-5-52.2.实例实例【实例实例1】二十世纪初卢瑟福和盖克两位科学家在观察二十世纪初卢瑟福和盖克两位科学家在观察与分析放射性物质放出的与分析放射性物质放出的 粒子个数的情况时粒子个数的情况时,他们做了他们做了2608次观察次观察(每次时间为每次时间为7.5秒秒)发现放射性物质在规定的发现放射性物质在规定的一段时间内一段时间内,其放射的粒子数其放射的粒子数X 服从泊松分布服从泊松分布.2023-5-5【实例实例2】von Bor
11、tkiewicz记录了记录了1875-1894年间普鲁士年间普鲁士骑兵军团被马踢死的士兵数。这些数据和骑兵军团被马踢死的士兵数。这些数据和 的泊松的泊松分布的对比如下表:分布的对比如下表:0.612023-5-5九、超几何分布九、超几何分布1.1.定义定义2023-5-52.2.实例实例【实例实例】一袋中有一袋中有N个大小形状相同的球,其中有个大小形状相同的球,其中有M个个白球,从中白球,从中不放回抽取不放回抽取n个球,记个球,记X为取到的白球数,为取到的白球数,X服从的分布即为超几何分布。服从的分布即为超几何分布。2023-5-5十、几何分布十、几何分布1.1.定义定义 若独立试验中仍只有若
12、独立试验中仍只有“成功成功”和和“失败失败”两种结果,且两种结果,且每次成功概率均为每次成功概率均为p,则直到首次出现,则直到首次出现“成功成功”为止所进为止所进行的试验次数服从行的试验次数服从几何分布几何分布.2023-5-52.分布律图形分布律图形3.3.实例实例【实例实例】掷一枚骰子,直到掷一枚骰子,直到1点朝上,记录投掷的次数点朝上,记录投掷的次数X,则则X服从参数服从参数 p=1/6 的几何分布。的几何分布。2023-5-5十一、负二项分布十一、负二项分布1.1.分布律分布律 若独立试验中仍只有若独立试验中仍只有“成功成功”和和“失败失败”两种结果,且两种结果,且每次成功概率均为每次
13、成功概率均为p,则直到出现,则直到出现r次次“成功成功”为止所进行为止所进行的试验次数服从的试验次数服从负二项分布负二项分布.11()(1),1,rrx rxp xCppxr r 2.分布律图形分布律图形2023-5-5十二、连续均匀分布十二、连续均匀分布1.1.定义定义2023-5-52.2.密度函数图密度函数图3.3.分布函数图分布函数图2023-5-54.4.实例实例【实例实例1】用用X表示四舍五入取整的误差,则表示四舍五入取整的误差,则X 服从服从-0.5,0.5上的均匀分布。上的均匀分布。【实例实例2】如图所示,轮盘赌指针停留位置与水平正向如图所示,轮盘赌指针停留位置与水平正向的夹角
14、记为的夹角记为X,则,则X 服从服从0,2p p上的均匀分布。上的均匀分布。X2023-5-5十三、指数分布十三、指数分布1.1.定义定义 某些元件或设备的寿命服从指数分布某些元件或设备的寿命服从指数分布.例如无线电元例如无线电元件的寿命件的寿命、电力设备的寿命、动物的寿命等都服从指数、电力设备的寿命、动物的寿命等都服从指数分布分布.2023-5-52.2.密度函数图密度函数图3.3.分布函数图分布函数图2023-5-5十四、一元正态分布十四、一元正态分布一个游戏:一个游戏:高尔顿钉板游戏高尔顿钉板游戏考察某一学科考察某一学科考试成绩的分布考试成绩的分布考察人类考察人类身高的分布身高的分布情况
15、情况思考:以上分布具有什么样的特点?思考:以上分布具有什么样的特点?2023-5-51 1、一元正态分布的定义、一元正态分布的定义则称则称 x 服从参数为服从参数为 ,2 的的正态分布正态分布,记作记作 x N(,2)定义定义 若随机变量若随机变量 X 的密度函数为的密度函数为xexfx222)(21)(p,为常数,为常数,0其中其中 亦称高斯亦称高斯(Gauss)分布分布2023-5-52 2、标准正态分布、标准正态分布=0,=1 的正态分布称为的正态分布称为标准正态分布标准正态分布,记作记作 x N(0,1)xexx2221)(p密度函数记为密度函数记为3 3、标准正态分布与一般正态分布之
16、间的关系、标准正态分布与一般正态分布之间的关系记记 u N(0,1),则则 x=+u N(,2)2023-5-5记记 u N(0,1),对于给定的,对于给定的0 x=0:10;%产生一个向量产生一个向量 Y=normpdf(x,1.2345,6)%求密度函数值求密度函数值 P=normcdf(x,1.2345,6)%求分布函数值求分布函数值2023-5-52023-5-5 u=norminv(1-0.05,0,1)t=tinv(1-0.05,50)chi2=chi2inv(1-0.025,8)f1=finv(1-0.01,7,13)f2=finv(1-0.99,13,7)2023-5-5第二节
17、第二节 生成一元分布随机数生成一元分布随机数2023-5-5一、均匀分布随机数和标准正态分布随机数一、均匀分布随机数和标准正态分布随机数调用格式:调用格式:Y=randY=rand(n)Y=rand(m,n)Y=rand(m n)Y=rand(m,n,p,)Y=rand(m n p)Y=rand(size(A)1.rand函数函数2023-5-5在在MATLAB7.7以前的版本中,以前的版本中,rand函数还可以这样调用:函数还可以这样调用:rand(method,s)s=rand(method)其中其中method是字符串变量,它的可能取值如下表所列:是字符串变量,它的可能取值如下表所列:2
18、023-5-5调用格式:调用格式:与与rand函数类似函数类似2.randn函数函数2023-5-5 x=rand(10)y=x(:);hist(y)xlabel(0,1上均匀分布随机数上均匀分布随机数);ylabel(频数频数);【例例16.2-1】调用调用rand函数生成函数生成1010的随机数矩阵,并将矩阵的随机数矩阵,并将矩阵按列拉长,然后调用按列拉长,然后调用hist函数画出频数直方图。函数画出频数直方图。2023-5-5%设置随机数生成器的算法为设置随机数生成器的算法为Mersenne Twister算法,初始种算法,初始种子为子为1 rand(twister,1);%生成生成2行
19、行6列的随机数矩阵,其元素服从列的随机数矩阵,其元素服从0,1上均匀分布上均匀分布 x1=rand(2,6)【例例16.2-1续续】设置随机数生成器的算法为设置随机数生成器的算法为Mersenne Twister算算法,生成均匀分布随机数矩阵法,生成均匀分布随机数矩阵2023-5-5二、常见一元分布随机数二、常见一元分布随机数MATLAB统计工具箱中函数名以统计工具箱中函数名以rnd三个字符结尾的三个字符结尾的函数用来生成常见分布的随机数。例如:函数用来生成常见分布的随机数。例如:betarndBeta分布分布exprnd指数分布指数分布gamrndGamma分布分布lognrnd对数正态分布
20、对数正态分布normrnd正态分布正态分布poissrnd泊松分布泊松分布randsample从有限总体中随机抽样从有限总体中随机抽样random指定分布指定分布2023-5-5%调用调用normrnd函数生成函数生成1000行行3列的随机数矩阵列的随机数矩阵x,其元素服从均值为,其元素服从均值为75,标准差为标准差为8的正态分布的正态分布 x=normrnd(75,8,1000,3);hist(x)%绘制矩阵绘制矩阵x每列的频数直方图每列的频数直方图 xlabel(正态分布随机数(正态分布随机数(mu=75,sigma=8));%为为X轴加标签轴加标签 ylabel(频数频数);%为为Y轴加
21、标签轴加标签 legend(第一列第一列,第二列第二列,第三列第三列)%为图形加标注框为图形加标注框【例例16.2-2】调用调用normrnd函数生成函数生成10003的正态分布随机数矩的正态分布随机数矩阵,其中均值阵,其中均值 为为75,标准差为,标准差为8,并作出各列的频数直方图,并作出各列的频数直方图405060708090100110050100150200250300正态分布随机数(=75,=8)频数 第一列第二列第三列2023-5-5%调用调用normrnd函数生成函数生成1000行行3列的随机数矩阵列的随机数矩阵x,其各列元素分别服从不同,其各列元素分别服从不同的正态分布的正态分
22、布 x=normrnd(repmat(0 15 40,1000,1),repmat(1 2 3,1000,1),1000,3);hist(x,50)%绘制矩阵绘制矩阵x每列的频数直方图每列的频数直方图 xlabel(正态分布随机数正态分布随机数);%为为X轴加标签轴加标签 ylabel(频数频数);%为为Y轴加标签轴加标签%为图形加标注框为图形加标注框 legend(mu=0,sigma=1,mu=15,sigma=2,mu=40,sigma=3)【例例16.2-3】调用调用normrnd函数生成函数生成10003的正态分布随机数矩的正态分布随机数矩阵,其中第各列均值分别为阵,其中第各列均值分
23、别为0,15,40,标准差分别为,标准差分别为1,2,3,并作出各列的频数直方图并作出各列的频数直方图2023-5-5%调用调用random函数生成函数生成10000行行1列的随机数向量列的随机数向量x,其元素服从二项分布,其元素服从二项分布B(10,0.3)x=random(bino,10,0.3,10000,1);fp,xp=ecdf(x);%计算经验累积概率分布函数值计算经验累积概率分布函数值 ecdfhist(fp,xp,50);%绘制频率直方图绘制频率直方图 xlabel(二项分布(二项分布(n=10,p=0.3)随机数)随机数);%为为X轴加标签轴加标签 ylabel(f(x);%
24、为为Y轴加标签轴加标签【例例16.2-4】调用调用random函数生成函数生成100001的二项分布随机数向的二项分布随机数向量,然后作出频率直方图。其中二项分布的参数为量,然后作出频率直方图。其中二项分布的参数为 n=10,p=0.3 01234567891000.511.5二项分布(n=10,p=0.3)随机数f(x)2023-5-5 x=random(chi2,10,10000,1);fp,xp=ecdf(x);%计算经验累积概率分布函数值计算经验累积概率分布函数值 ecdfhist(fp,xp,50);%绘制频率直方图绘制频率直方图 hold on t=linspace(0,max(x
25、),100);y=chi2pdf(t,10);plot(t,y,r,linewidth,3)xlabel(x (chi2(10);%为为X轴加标签轴加标签 ylabel(f(x);%为为Y轴加标签轴加标签 legend(频率直方图频率直方图,密度函数曲线密度函数曲线)%为图形加标注框为图形加标注框【例例16.2-5】调用调用random函数生成函数生成100001的卡方分布随机数向的卡方分布随机数向量,然后作出频率直方图,并与自由度为量,然后作出频率直方图,并与自由度为10的卡方分布的密度的卡方分布的密度函数曲线作比较。其中卡方分布的参数(自由度)为函数曲线作比较。其中卡方分布的参数(自由度)
展开阅读全文