1、数据分析实验报告班级： 07 信计 0 班 学号： 姓名： 实验日期 2010-3-11 实验地点： 实验楼 505 实验名称： 样本数据的特征分析 使用软件名称：MATLAB实验目的1.熟练掌握利用 Matlab 软件计算均值、方差、协方差、相关系数、标准差与变异系数、偏度与峰度，中位数、分位数、三均值、四分位极差与极差； 2.熟练掌握 jbtest 与 lillietest 关于一元数据的正态性检验；3.掌握统计作图方法；4.掌握多元数据的数字特征与相关矩阵的处理方法；实验内容安徽省 1990-2004 年万元工业 GDP 废气排放量、废水排放量、固体废物排放量以及用于污染治理的投入经费比
2、重见表 6.1.1，解决以下问题：表 6.1.1 废气、废水、固体废物排放量及污染治理的投入经费占 GDP 比重万元工业 GDP废气排放量万元工业 GDP固体物排放量万元工业 GDP废水排放量环境污染治理投资占 GDP 比重年份（立方米） （千克） （吨） （%）1990 104254.40 519.48 441.65 0.181991 94415.00 476.97 398.19 0.261992 89317.41 119.45 332.14 0.231993 63012.42 67.93 203.91 0.201994 45435.04 7.86 128.20 0.171995 46383
3、.42 12.45 113.39 0.221996 39874.19 13.24 87.12 0.151997 38412.85 37.97 76.98 0.211998 35270.79 45.36 59.68 0.111999 35200.76 34.93 60.82 0.152000 35848.97 1.82 57.35 0.192001 40348.43 1.17 53.06 0.112002 40392.96 0.16 50.96 0.122003 37237.13 0.05 43.94 0.152004 34176.27 0.06 36.90 0.131.计算各指标的均值、方差、
4、标准差、变异系数以及相关系数矩阵；2.计算各指标的偏度、峰度、三均值以及极差；3.做出各指标数据直方图并检验该数据是否服从正态分布？若不服从正态分布，利用 boxcox 变换以后给出该数据的密度函数； 4.上网查找 1990-2004 江苏省万元工业 GDP 废气排放量，安徽省与江苏省是 否服从同样的分布？程序如下：clear;clcformat bank %保留两位小数%安徽省数据%A=104254.40 519.48 441.65 0.1894415.00 476.97 398.19 0.2689317.41 119.45 332.14 0.2363012.42 67.93 203.91
5、0.2045435.04 7.86 128.20 0.1746383.42 12.45 113.39 0.2239874.19 13.24 87.12 0.1538412.85 37.97 76.98 0.2135270.79 45.36 59.68 0.1135200.76 34.93 60.82 0.1535848.97 1.82 57.35 0.1940348.43 1.17 53.06 0.1140392.96 0.16 50.96 0.1237237.13 0.05 43.94 0.1534176.27 0.06 36.90 0.13; %计算各指标的均值、方差、标准差、变异系数、偏
6、度、峰度以及极差A1=mean(A);var(A);std(A);std(A)./mean(A);skewness(A,0);kurtosis(A,0)-3;range(A)%三均值A2=1/4 1/2 1/4*prctile(A,25 50 75)%计算各指标的相关系数矩阵A3=corrcoef(A)%做出各指标数据直方图subplot(221),histfit(A(:,1),8)subplot(222),histfit(A(:,2),8)subplot(223),histfit(A(:,3),8)subplot(224),histfit(A(:,4),7)%检验该数据是否服从正态分布for
7、 i=1:4h(i),p(i),lstat(i),cv(i)=lillietest(A(:,i),0.05);endh,p%计算出前三列不服从正态分布，利用boxcox变换以后给出该数据的密度函数t1,I1=boxcox(A(:,1)t2,I2=boxcox(A(:,2)t3,I3=boxcox(A(:,3)%江苏省数据%J=5047.00 24.23 2234.005730.00 24.00 2077.005352.00 22.49 2450.005202.00 21.16 2323.005772.14 21.16 2725.007872.11 22.02 2883.007450.81 21
8、.97 2891.008164.24 22.01 2857.007487.85 20.29 2944.008353.93 20.11 2896.009078.20 20.19 3038.0013343.99 27.10 3553.0014286.46 26.30 3796.0014632.69 24.75 3893.7017818.42 26.35 4672.5320196.58 29.63 5757.37;%检验该数据是否服从正态分布for i=1:3h1(i),p1(i),lstat1(i),cv1(i)=lillietest(J(:,i),0.05);endh1,p1%计算出第一、三列不
9、服从正态分布，利用boxcox变换以后给出该数据的密度函数t11,I11=boxcox(J(:,1)t33,I33=boxcox(J(:,3)1.各指标的均值、方差、标准差、变异系数、偏度、峰度、极差和三均值如 下表：指标数字特征万元工业 GDP 废气排放量（ 3m）万元工业 GDP固体物排放（ kg）万元工业 GDP废水排放量（ t）环境污染治理投资占 GDP 比重（%）均值 51972.00 89.26 142.95 0.17方差 577053049.90 28703.47 18645.50 0.00标准差 24021.93 169.42 136.55 0.05实验结果分析变异系数 0.4
10、6 1.90 0.96 0.27偏度 1.44 2.25 1.44 0.29峰度 0.57 3.87 0.67 -0.83极差 70078.13 519.43 404.75 0.15三均值 43937.01 22.52 98.27 0.17相关系数矩阵为：165.049.62.0.49.18. 2R2.各指标数据的直方图如下：-5 0 5 10 15x 1040510-500 0 500 1000051015-500 0 500 100002468100 0.1 0.2 0.3 0.4012343.经检验可知：安徽省万元工业 GDP 废气排放量、废水排放量、固体废物排放量数据不服从正态分布，而
11、用于污染治理的投入经费数据服从正态分布。前三列数据利用 boxcox 变换以后（解得 -0.61.9,-2.3,31 ）服从正态分布，而变换公式为： /)(xy ，故密度函数为：)2(exp1)(yyf，则原始数据的密度函数为：)21(exp2)( 21xf。4.经检验可知：江苏省万元工业 GDP 废气排放量、固体废物排放量数据不服从正态分布，而废水排放量服从正态分布。附：表 1990-2005 年万元工业 GDP 废气排放量、废水排放量、固体废物排放量时间 工业废气排放量 （亿立方米） 工业废水排放量 （亿吨） 工业固体废物排放量 （万吨）1990 5047.00 24.23 2234.00
12、1991 5730.00 24.00 2077.001992 5352.00 22.49 2450.001993 5202.00 21.16 2323.001994 5772.14 21.16 2725.001995 7872.11 22.02 2883.001996 7450.81 21.97 2891.001997 8164.24 22.01 2857.001998 7487.85 20.29 2944.001999 8353.93 20.11 2896.002000 9078.20 20.19 3038.002001 13343.99 27.10 3553.002002 14286.46 26.30 3796.002003 14632.69 24.75 3893.702004 17818.42 26.35 4672.532005 20196.58 29.63 5757.37教师评语成绩日期：2010 年 3 月
