节概率论与数理统计.pptVIP

5.6节 概率论与数理统计 一 各种统计分布函数 表5-6-1中列出了20种分布函数，统计工具箱提供了包括分布函数~cdf(Cumulattive Distribution Function)、概率密度函数~pdf(Probability Distribution Function)、分布函数的逆函数~inv(Inverse Cumulattive Distribution Function) 以及这些分布的理论统计特性（均值和方差）计算函数~stat。 表5-6-1 MATLAB中表示各种概率分布的前缀文字 表5-6-1的使用方法 把表中不同分布后面括号中的文字作为前缀，把所需计算的特性作为后缀，就可以组合成一个函数。例如离散类的二项式分布有binopdf, binocdf, binoinv, binostat四个函数，连续类数据的标准正态分布有normpdf, normcdf, norminv, normstat四个函数，连续类统计量的χ2分布有chi2pdf, chi2cdf, chi2inv, chi2stat四个函数，…等等，20种分布就有80个函数。由于各种分布函数的解析形式都是已知的，这些子程序的编写并不困难。 【例5-6-1】概率分布曲线的绘制 【例5-6-1】(a)求标准正态分布N(0,1)，自由度V=10的χ2分布和N=10,p=0.2的二项式分布B(10,0.2)的分布函数，并画出其概率密度曲线和分布曲线；(b) 。求出分布函数为0.05和0.95处的随机变量值；(c) 。求出这几个分布的均值和方差。 解：分别键入help normpdf, help binopdf, help chi2pdf 以了解它们的用法，特别是了解输入变元的意义，得知： f=normpdf(x,mu,sigma) 其中x为随机变量数组，mu为均值，sigma为标准差。 f=chi2pdf(x,V) 其中x为随机变量数组，V为自由度数（整数）， f=binopdf(x,N,p) 其中x为随机变量整数数组，N为次数，0≤p≤1为成功概率 概率分布参数的选择 题目中给出的变元参数应足以用来调用这些概率密度函数，只有随机变量x的取值及范围，需要事先对该分布的特性有所了解。首先要弄清它是离散量还是连续量，其次要取适当的范围。范围取小了不能显示分布的全局，取大了又可能显示不出细节。正确的取法应该使该范围内分布函数F(x)的左边界值略小于0.025，右边界值略大于0.975，这就可以基本涵盖随机变量以概率95%存在的主要区域，又不致涉及关系不大的区域。初学者往往要作几次试探才能做到。好在在计算机上改几个参数、作几次试探是很简便的事。 绘制概率分布函数的程序exn561 x1 = [-3:0.1:3]; %标准正态随机变量取值范围 f1 = normpdf(x1,0,1); % 标准正态分布的概率密度函数 F1 = normcdf(x1,0,1); % 标准正态分布的分布函数 subplot(2,2,1),plot(x1,f1,:,x1,F1,-),grid on % 绘曲线 line([-4,4],[0.025,0.025]),line([-4,4],[0.975,0.975])% 画横线 x2 = [0:0.1:20]; % 试探得出的范围 f2 = chi2pdf(x2,10); % χ2分布的概率密度函数 F2 = chi2cdf(x2,10); % χ2分布的分布函数 subplot(2,2,2),plot(x2,f2,:,x2,F2,-),grid on % 绘曲线 line([0,20],[0.025,0.025]),line([0,20],[0.975,0.975]) % 画横线); ……. 运行程序得出图5-40，其中实线是分布函数的曲线，虚线则是密度函数的曲线。上下两根横线分别为0.975和0.025， 三种分布的概率密度和分布曲线 对分布函数图形的改进 第三个子图中，随机变量是离散取值的，在两个相邻的取值点之间的概率不会变化，所以它的分布函数表现为阶梯形。密度函数是分布函数的导数，所以在阶梯突变处导数为无穷大，宽度又为无穷小，面积等于阶梯的高度，通常用一个脉冲表示。脉冲的高度表示它所包含的面积，也就等于阶梯波形的高度。而plot函数画图是把各离散点之间用直线联接，所得的曲线是不对的。应该要用stairs命令画阶梯图，用stem命令画脉冲图。改正后的程序如下： subplot(2,2,4), stairs(x3,F3,-),stem(x3,f3,:), 图中第四子图给出了改正后的结果，可见其概率密度函数是离散的脉冲。从中还可以判断，x3的取值范围不必为[0:1