上机练习六 概率统计回归分析模型doc.doc

实验六　概率统计、回归分析模型 2. 直观了解统计基本内容. 3. 掌握用数学软件包Matlab求解概率、统计问题. 4. 综合提高用Matlab软件解决实际问题的能力. 学会分析问题、建立问题的数学表达式并加以求解和推广,体会数学建模的整个过程. 例1 上海市区社会商品零售总额和全民所有制职工工资总额的数据如下 年份 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 商品零售总额(亿元) 44.1 51.8 61.7 67.9 68.7 77.5 95.9 137.4 155.0 175.0 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] 2data的数据保存在文件data1中：save data1 data 3,先用命令： load data1 data1中的数据,再用以下命令分别将矩阵data的第一、二、三行的数据赋给变量t、x、y： t=data(1,:) x=data(2,:) y=data(3,:) 若要调用矩阵data的第j列的数据,可用命令： data(:,j) 二、基本统计量 对随机变量x,计算其基本统计量的命令如下： 均值：mean(x) ：median(x) ：std(x) ：var(x) ：skewness(x) kurtosis(x) 例2 对例1中的职工工资总额x,可计算上述基本统计量. load data1 x=data(2,:) mean=mean(x) median=median(x) std=std(x) var=var(x) skewness=skewness(x) kurtosis=kurtosis(x) 三、常见概率分布的函数 常见的几种分布的命令字符为： ：norm ：exp ：poiss ：beta ：weib ：chi2 t：t F：F Matlab,其命令字符为： ：pdf ：cdf ：inv ：stat rnd (当需要一种分布的某一类函数时,将以上所列的分布命令字符与函数命令字符接起来,并输入自变量（可以是标量、数组或矩阵）和参数即可.） 例mu、标准差为sigma的正态分布,举例如下： 1、密度函数：p=normpdf(x,mu,sigma) (mu=0,sigma=1时可缺省) 例3 画出正态分布和的概率密度函数图形. 在Matlab中输入以下命令： x=-6:0.01:6; y=normpdf(x); z=normpdf(x,0,2); plot(x,y,x,z) 2：P=normcdf(x,mu,sigma) ,求. 命令为：P=normcdf(1)-normcdf(-1) 结果为：P =0.6827 3、逆概率分布：x=norminv(P,mu,sigma). x ,使得P{Xx}=P.此命令可用来求分位数. 例5 取,求 的含义是：,= 命令为： u=norminv(0.975) 时,P=0.975, norminv(0.975)=1.96 4、均值与方差：[m,v]=normstat(mu,sigma) N(3,52)的均值与方差. 命令为：[m,v]=normstat(3,5) ：m=3,v=25 5：normrnd (mu,sigma,m,n). 阶的正态分布随机数矩阵. 例7 命令：M=normrnd([1 2 3;4 5 6],0.1,2,3) M=0.9567 2.0125 2.8854 　3.8334 5.0288 6.1191 此命令产生了的正态分布随机数矩阵,各数分别服从N(1,0.12), N(2,22), N(3, 32), N(4,0.12), N(5, 22),N(6, 32) 四、频数直方图的描绘 1、给出数组data的频数表的命令为： [N,X]=hist(data,k) 此命令将区间[min(data),max(data)]分为k个小区间（缺省为10）,返回数组data落在每一个小区间的频数N和每一个小区间的中点X. 2、描绘数组da