1MATLAB软件在基本统计中的命令.pptVIP

MATLAB软件在基本统计中的命令 一、数据的录入、保存和调用 统计中的数据往往较大，在交互环境输入后，往往需要保存为数据文件，以便调用来作各种统计分析。下面通过一个实例说明。 例1 上海市区社会商品零售总额和全民所有制职工工资总额的数据如下： 1、年份数据是以1为增量的，可用产生向量的方法来输入： 命令格式：x=a:h:b,即从a到b，以h为增量的行向量，当h缺省时，认为增量为1. 如年份：t=78：87 2、分别以x、y表示职工工资总额和商品零售总额，以向量的形式输入： X=[23.8,27.6,31.6,32.4,33.7,34.9,43.2,52.8,63.8,73.4]; Y=[41.4,51.8,61.7,67.9,68.7,77.5,95.9,137.4,155.0,175.0] 3、将变量t,x,y的数据保存在文件data中，用以下命令：sava data t x y 4、进行统计分析时，用以下命令调用数据文件data中的数据：load data 注：也可用矩阵的形式输入和保存数据。 （1）输入矩阵 data=[78,79,80,81,82,83,84,85,86,87;23.8,27.6,31.6,32.4,33.7,34.9,43.2,52.8,63.8,73.4;[41.4,51.8,61.7,67.9,68.7,77.5,95.9,137.4,155.0,175.0] (2)将矩阵data的数据保存在文件data1中：save data1 data (3)进行统计分析时，先用命令：load data1 在用命令将矩阵data中第1、2、3行的数据赋给变量t,x,y:t=data(1,:);x=data(2,:);y=data(3,:) 即若要调用矩阵data的第j列的数据，可用命令：data(:,j) 二、基本统计量 均值： mean(x) 中位数：medaian(x) 标准差：std(x) 方差： var(x) 偏度： skewness(x) 峰度： kurtosis(x) 三、常见概率分布的函数 1、正态分布：norm 2、指数分布：exp 3、泊松分布：poiss 4、β分布：beta 5、威布尔分布：weib 6、X2 分布：chi2 7、t分布：t 8、F分布：F MATLAB工具箱对每一种分布提供五类函数，其命令字符为： 概率密度：pdf 概率分布：cdf 逆概率分布：inv 均值与方差：stat 随机生成：rnd 当需要一种分布的某一类函数时，将以上所列的分布命令字符与函数命令字符接起来，并输入自变量（可以是标量、数组或矩阵）和参数即可。 1、密度函数：p=normpdf(x,mu,sigma) (当mu=0,sigma=1时可缺省） 例2 画出正态分布N(0,1),N(0,22)的概率密度函数图形。 在MATLAB中输入命令： X=-6:0.1:6; y=normpdf(x); z=normodf(x,0,2); Plot(x,y,x,z) 2、概率分布：p=normcdf(x,mu,sigma) 例3 （1）设X～N(0,1)，求P{-1X1} 命令： p=normcdf(1)-normcdf(-1) 结果：p = 0.6827 （2）设 X～N(2,52)，求P{1X3} p=normcdf(3,2,5)-normcdf(1,2,5) p=0.1585 （3）设X～λ=2的指数分布，求p{1x3}. 解：∵mu=1/2,sigma=1/4 ∴命令：p=expcdf(3,1/2,1/4) 结果：0.1329 3、逆概率分布：x=norminv(p,mu,sigma) 即求出x，使得P{Xx}=p，此命令可用来求分位数。 例3 (1)取α=0.05,求u1-α/2. 即使得p{X u1-α/2}=1-α/2=0.975 u=norminv(0.975) u =1.9600 (2) X～λ=2的指数分布，p=0.5 命令： x=expinv(0.5,1/2,1/4) 结果： x=0.3466 4、均值与方差：[m,v]=normstat(mu,sigma) 例 设X～N(1,32)，求其均值与方差。 [m,v]=normstat(1,3) m =1,v =9 四、频数直方图 1、给出数组data的频数表的命令： [N,x]=hist(data,k) 此命令将区间[min(data),max(data)]分为k个小区间（缺省时为10），返回数组在每个小区间的频数N和小区间的中点x。 2、描绘数组data的频数直方图 Hiat(data,k) 例4 某校60名学生的一次考试成绩如下： 93 75