数据的统计分析与描述-matlab.pptVIP

三、常见概率分布的函数**Matlab工具箱对每一种分布都提供五类函数，其命令字符为：概率密度：pdf概率分布：cdf逆概率分布：inv均值与方差：stat随机数生成：rnd（当需要一种分布的某一类函数时，将以上所列的分布命令字符与函数命令字符接起来，并输入自变量（可以是标量、数组或矩阵）和参数即可.）12在Matlab中输入以下命令：x=-6:0.01:6;y=normpdf(x);z=normpdf(x,0,2);plot(x,y,x,z)1、密度函数：p=normpdf(x,mu,sigma)(当mu=0,sigma=1时可缺省)如对均值为mu、标准差为sigma的正态分布，举例如下：逆概率分布：x=norminv(P,mu,sigma).即求出x，使得P{Xx}=P.此命令可用来求分位数.2、概率分布：P=normcdf(x,mu,sigma)4、均值与方差：[m,v]=normstat(mu,sigma)例5求正态分布N(3,52)的均值与方差.命令为：[m,v]=normstat(3,5)结果为：m=3,v=255、随机数生成：normrnd(mu,sigma,m,n).产生m?n阶的正态分布随机数矩阵.例6命令：M=normrnd([123;456],0.1,2,3)结果为：M=0.95672.01252.88543.83345.02886.1191此命令产生了2?3的正态分布随机数矩阵，各数分别服从N(1,0.12),N(2,22),N(3,32),N(4,0.12),N(5,22),N(6,32)四、频数直方图的描绘**给出数组data的频数表的命令为：[N,X]=hist(data,k)此命令将区间[min(data),max(data)]分为k个小区间（缺省为10），返回数组data落在每一个小区间的频数N和每一个小区间的中点X.描绘数组data的频数直方图的命令为：hist(data,k)Matlab相关命令介绍load从matlab数据文件中载入数据S=load(数据文件名)hist绘制给定数据的直方图hist(x,m)table=tabulate(x)绘制频数表，返回值table中，第一列为x的值，第二列为该值出现的次数，最后一列包含每个值的百分比。对于给定的数据集，假设它们满足以上十种分布之一，如何确定属于哪种分布？01例1：某次笔试的分数见data1.txt，试画出频数直方图03从图形上看，笔试成绩较为接近正态分布05x=load(data1.txt);x=x(:);hist(x)02绘制频数直方图，或列出频数表04频数直方图或频数表频数直方图或频数表x=load(data2.txt);x=x(:);hist(x)例2：某次上机考试的分数见data2.txt，试画出频数直方图从图形上看，上机考试成绩较为接近离散均匀分布x=load(data3.txt);x=x(:);hist(x)例3：上海1998年来的月降雨量的数据见data3.txt，

试画出频数直方图从图形上看，月降雨量较为接近?2分布频数直方图或频数表x=load(data4.txt);

x=x(:);tabulate(x)hist(x)在重复数据较多的情况下，我们也可以利用Matlab自带的tabulate函数生成频数表，并以频数表的形式来发掘数据分布的规律。例4：给出数据data4.txt，试画出其直方图，并生成频数表频数直方图或频数表x=load(data5.txt);x=x(:);hist(x)fiugrehistfit(x)%加入较接近的正态分布密度曲线例5：现累积有100次刀具故障记录，当故障出现时该批刀具完成的零件数见data5.txt，试画出其直方图。从图形上看，较为接近正态分布当我们可以基本确定数据集X符合某种分布后，我们还需要确定这个分布的参数。01由于正态分布情况发生的比较多，故我们主要考虑正态分布的情形。02对于未知参数的估计，可分两种情况：点估计区间估计03参数估计参数估计：点估计构造样本X与某个统计量有关的一个函数，作为该统计量的一个估计，称为点估计。Matlab统计工