matlab与应用多元统计分析

多元统计分析中的应用研究摘要：许多实际问题往往需要对数据进行统计分析，建立合适的统计模型，过去一般采用SAS、SPSS软件分析，本文给出Matlab软件在多元统计分析上的应用,主要介绍Matlab在聚类分析、判别分析、主成份分析上的应用,文中均给以实例,结果令人满意。关键词：Matlab软件；聚类分析；主成份分析ResearchforapplicationofMultivariateStatisticalAnalysisAbstract:ManypracticequestionsometimesneedStatisticalAnalysistodata.,andestablishappropriateStatisticalmodelSASandSPSSsoftwarewerecommonlyusedinforetime,thispapergivetheapplicationofMatlabsoftwareinMultivariateStatisticalAnalysis,mostlyintroducetheapplicationofMatlabsoftwareinpriciplecomponentanalysisandclusteranalysisanddifferentiateanalysis.Theexamplearegiveninwritingandtheresultaresatisfaction.Keywords:Matlabsoftware;clusteranalysis;priciplecomponentanalysis0引言许多实际问题往往需要对数据进行多元统计分析,建立合适的模型,在多元统计分析方面,常用的软件有SAS、SPSS、S-PLUS等。我们在这里给出Matlab在多元统计分析上的应用,在较早的版本中,统计功能不那么强大,而在Matlab6.x版本中，仅在统计工具中的功能函数就达200多个，功能已足以赶超任何其他专用的统计软件，在应用上Matlab具有其他软件不可比拟的操作简单，接口方便，扩充能力强等优势，再加上Matlab的应用范围广泛，因此可以预见其在统计应用上越来越占有极其重要的地位，下面用实例给出Matlab在聚类分析、主成份分析上的应用。1聚类分析聚类分析法是一门多元统计分类法，其目的是把分类对象按一定规则分成若干类，所分成的类是根据数据本身的特征确定的。聚类分析法根据变量（或样品或指标）的属性或特征的相似性，用数学方法把他们逐步地划类，最后得到一个能反映样品之间或指标之间亲疏关系的客观分类系统图，称为谱系聚类图。聚类分析的步骤有：数据变换，计算n个样品的两两间的距离，先分为一类，在剩下的n—1个样品计算距离，按照不同距离最小的原则，增加分类的个数，减少所需要分类的样品的个数，循环进行下去，直到类的总个数为1时止。根据类之间的距离，画出谱系聚类图。我们对杭州所辖张家港市2005年七条河流中主要污染因子(指标)即CODmn,BOD5,非离子氨,氨氮,挥发酚,石油类共6个变量(资料见表1,来源于张家港市2003年环境质量报告书),进行聚类分析。表1港2005年七条河流主要污染因子河流CODmnBOD5非离子氮氨氮挥发酚石油类类型张家港河3.14&4123.7825.794.176.473二干河5.479.5726.4823.796.425.581东横河3.14.3121.222.485.346.543横套河5.679.5410.2320.874.26.81四干河6.819.0516.1824.565.25.451华妙河6.217.0821.0531.566.15&212盐铁塘4.878.9726.5434.565.588.072我们利用Matlab6.5中的cluster命令实现,具体程序如下3.148.4123.7825.794.176.475.479.5726.4823.796.425.583.14.3121.222.485.346.545.679.5410.2320.874.26.86.819.0516.1824.565.25.456.217.0821.0531.566.158.214.878.9726.5434.565.588.07}{n,m}=size(x);Stdr=std(x);xx=x./stdr(ones(n,1),;);%标准化变换y=pdist(xx);%计算各样本间距离(这里为欧氏距离)z=linkage(y); %进行聚类(这里为最短距离法)h=dendrogram(z);%画聚类谱系图t=cluster(z,3)%将全部样本分为3类find(t==2); %找出属于第2类的样品编号执行后得到所要结果聚类谱系图见图1t={3,1,3,1,1,2,2}即全部样本分为3类。结果见表1从图1可以看出：七条河流中,二干河、横套河、四干河属于一类,污染较重,主要是CODmn、BOD5超标多;华妙河、盐铁塘属于一类,污染一般,主要是氨氮、石油类超标;张家港河、东横河属于一类,污染较轻,总的来说,各河流都存在不同程度的污染,因此全市应对各河流严格监督管理,着力实施水污染防治工作,太湖流域水污染源应限期治理达标排放,巩固水污染防治工作成果,加大投入,新建或改、扩建废水治理工程,确保达标排放。图1：聚类谱系图2主成分分析主成分分析是将多个指标化为少数几个综合指标的一种多元统计分析方法。对于实际工作中遇到的多指标系统评估问题，主成分分析可以将多个指标综合为单个指数的形式。主成分的计算步骤如下：第一步，原始数据零均值标准化。设每个指标的样本数据为X,x,Lxx,x,Lx.作如下变换，令i,1i,2i,Ni,1i,2i,Nx'ij=(xij-xi)/SiijijiipS其中k尸工k/x. >pVVQ广k?m11(i=l,2,L,P;j=l,2,L,N)第二步，计算相关矩阵R=(rij),其中第三步，求矩阵R的特征值久i与相应的标准正交化的特征向量Ai。第四步，计算第j个主成分yi的贡献率Hj=A.工.X当前q个主成分（即新的综合指标）的累积贡献率超过85%时，就提取前q个主成分作为评价指标，它们保持原始数据总信息量的85%以上。这里给出江苏省生态城市主成份分析实例城市环境生态化是城市发展的必然趁势,表现为社会、经济、环境与生态全方位的现代化水平,一个符合生态规律的生态城市应该是结构合理、功能高效和关系协调的城市生态系统所谓结构合理是指适度的人口密度,合理的土地利用,良好的环境质量,充足的绿地系统,完善的基础设施,有效的自然保护功能高效是指资源的优化配置、物力的经济投入、人力的充分发挥、物流的畅通有序、信息流的快捷关系协调是指人和自然协调、社会关系协调、城乡协调、资源利用和更新协调一个城市要实现生态城市的发展目标,关键是在市场经济的体制下逐步改善城市的生态环境质量,防止生态环境质量恶化,因此,对城市的生态环境水平调查评价很有必要。我们对江苏省十个城市的生态环境状况进行了调查,得到生态环境指标的指数值,见表2。现对生态环境水平分析和评价表2指标指数值结构功能协调度生态歼境水平權序二圾措标人口第构巧基础设施帀结构力城市盹d物质还原却配置珂生产效*^城市文明彌可持续性M无锡市078830.76330.47450,82460.87910.95380.87850.63050.8928&常州市0.73910.72870.51260.76030,873G0J2570.8M20.61870.78317锁江市0.81110.76290.8810.68880.81830.92850.85370.63130.560810张家港市0.65870.85520.89030,89770,944«0.94340.90270.74150.S4193连耳港市0.66430.75640,82880,79260.92020.91540.87290.63960.84«46摄州市0.82590.74550.785078560.92630.88710,84850,61420,76168秦州市O.84S60-780.80320.65090.91850,93570.84730.57340.82349廉州市0.68340390.88620.89020,95050,8760,9044O,S9B0.638420.84950.89180.39870.6799D磁0.95790.886«0.61860.96&44苏州市0.7B4fi0,89540.3970,98770.88730.97410.90350.73820.65141

我们利用Matlab6.5中的princomp命令实现,具体程序如下[0,78830,73910,81110,65870,65430.82590.84860.68340.84950,7S4fi076330.72870.76290.85520.756^0.74550.78000.94900.89180,89540.47450,5126U.88100,89030.82880.78500.80320.88620.398703970032460.76030.68880.89770.79260J8560.65090.89020.67990.98770.87910.87360.81830.94460.92020.92630.918&0.95050,86200.88730&5380.92570.92850.94340.91540.88710.9357087600.95790.97410.878&085420.85370.90270.87290.8485034730.90440.88660,9035□63050.61870.63130.74150.63980,61420.57340.8&800.61860*73820.89280.7&310.560S0.84190,84640.76160.82340.63840,96040.8514]x=x'stdr=std(x); %求各变量标准差{n,m}=size(x);Sddsta=x./stdr(ones(n,1),:);%标准化变换{p,princ,egenvalue}=princomp(sddata) %调用主成分分析程序P3=(:,1:3) %输出前三个主成分系数sc=princ(:,1:3)%输出前三个主成分得分egenvalue %输出特征根per=100*egenvalue/sum(egenvalue) %输出各个主成分贡献率执行后得到所要结果,这里是前三个主成分、主成分得分、特征根即-0*36770,37020.1364-0*36770,37020.13640.40480.14420.2313-052990.1812-0.3282-0.35350.04980』5旣P=0,3353-0.15CL0礙4-C.13i80.5273-0.02700.42360.3116-0.095&0.4815-0.0267-0.2804-0,06430.45890,5933-1.3364-0.21590.2729—1海_1矩6--205522^7540,1623O.SU43鳥沁=0J6J4-0,90021.4326-0.9266-1.66360,5837-1.8984-0.69070.54323,9332-1.5024-0.9765-1.21342.3442-0.44981.37362.7034-0,6469-0.8301 L.3H97 0J946Egenvalue=[3.8811,2.6407,1.0597]', per=[43.12,29.34,11.97]'这样,前三个主成分为Z1=-0.3677x1+0.3702x2+0.1364x3+0.4048x4+0.3355x5-0.1318x6+0.4236x7+0.4815x8-0.0643x9Z2=0.1442x1+0.2313x2-0.5299x3+0.1812x4-0.1601x5+0.5273x6+0.3116x7-0.0267x8+0.4589x9Z3=-0.3282x1-0.3535x2+0.0498x3+0.0582x4+0.5664x5-0.0270x6-0.0958x7-0.2804x8+0.5933x9第一主成分贡献率为43.12%,第二主成分贡献率为29.34%,第三主成分贡献率为11.97%,前三个主成分累计贡献率达84.24%如果按80%以上的信息量选取新因子,则可以选取前三个新因子第一新因子Z1包含的信息量最大为43.12%,它的主要代表变量为X8（城市文明）、X7（生产效率）、X4（城市绿化），其权重系数分别为0.4815、0.4236、0.4048,反映了这三个变量与生态环境水平密切相关,第二新因子Z2包含的信息量次之为29.34%,它的主要代表变量为X3（