基于微量元素的主成分分析对茶叶的评价毕业论文.doc

基于微量元素的主成分分析对茶叶的评价摘要：以微量元素含量为指标，应用软件和软件对微量元素进行基本统计分析和主成分分析，筛选出四个主成分，找出微量元素的相关关系，并对不同的茶叶进行了排名，为研究茶叶的品质和功效提供依据。关键词：微量元素；基本统计；主成分分析；茶叶评价引言现代医学已经证明，微量元素与人体的健康密切相关,已经发现了铁、镁、锌、铜、硒、锰等微量元素。这些微量元素都参与了生命体内的许多重要生理过程，与生命体内的许多疾病的发生密切相关。因此，微量元素的研究极为重要 曾自林,陈启霞,居明乔微量元素与中药功效相关性分析j时珍国医国药,2001,12(7)：658661。茶是我们生活的必需品，具有润肺止渴，安神定性的功效。为了更好地研究茶的特性功效，通过基本统计分析和主成分分析对茶叶中的微量元素进行研究，对不同品种的茶进行综合排名。一、 主成分分析概况1.1 基本思想主成分分析是一种数学降维的方法，找出几个综合变量来代替原来众多的变量，使这些综合变量能尽可能地代表原来变量的信息量，而且彼此之间互不相关，这种将把多个变量化为少数几个相互无关的综合变量的统计分析方法叫做主成分分析或主分量分析。它对原变量进行转换，使少数几个变量是原变量的线性组合，同时，这些变量又尽可能多的表征原变量的数据结构特征而不失信息 许禄,邵学广.化学计量学方法m .北京:科学出版社.2006:130-133.。1.2 数学模型对于一个样本资料，观测个变量，个样品的数据资料阵为：其中主成分分析就是将个观测变量综合成为个新的变量（综合变量），即简写为其中互不相关，的方差大于的方差大于的方差，以此类推，为第一主成分，为第二主成分，以此类推，有第个主成分。称为主成分系数。1.3 计算步骤第一步：对原始数据进行标准化处理。，其中，。第二步：计算样本相关系数矩阵。假定原始数据标准化后仍为表示，则经标准化处理后的数据的相关系数为：第三步：用雅克比方法求相关系数矩阵的特征值和相应的特征向量。特征值为，特征向量第四步：选择重要的主成分，并写出出成分表达式。贡献率=贡献率越大，说明该主成分所包含的原始变量的信息越强。主成分个数k的选取，主要根据主成分的累积贡献率来决定，即一般要求累积贡献率达到85%以上，这样才能保证综合变量能包括原始变量的绝大多数信息 张文霖.主成分分析在spss中的操作应用j.市场研究,2005(12)：3134.。第五步：计算主成分得分根据标准化的原始数据，按照各个样品，分别代入主成分表达式，就可以得到各主成分下的各个样品的新数据，即为主成分得分。第六步：依据主成分得分的依据，进行进一步的统计分析，如综合评价。对主成分进行加权综合。将主成分的权数根据它们的方差贡献率为确定，因为方差贡献率反映了各个主成分的信息含量多少。综合评价得分：二、数据处理过程2.1 数据来源数据来自文献不同产地茶叶矿物质元素含量的调查分析 苏冰霞,郑亚军,吴学进等不同产地茶叶矿物质元素含量的调查分析j.微量元素与健康研究,2012,29(1)：29322.2 数据基本统计分析与预处理 首先对原始数据进行描述性统计，包括最值、均值和方差的统计；其次由于原始数据差别较大，因此对数据进行标准化处理，以消除由于原始数据存在的量纲不同和数量级不同而对分析结果造成的影响。标准化后的数据具备可比性，并遵从正态分布规律 曹红翠,吴启勋.中药甘草中微量元素的主成分分析和系统聚类分析j.甘肃联合大学学报:自然科学版,2010,24 (2 )：61 63。三、数据处理结果3.1 基本统计分析结果对原始数据进行描述性统计后，得到各元素在茶叶中的含量分析。数据中各微量元素的最大值与最小值相差很大，最小相差约4倍，最大达20余倍。表 1 数据描述性统计表最小值最大值标准差ba11.59352.90103.10co0.160.510.10cr0.531.580.32li0.080.610.17zn20.7074.1019.84mn251.503048.00760.68fe60.30844.00231.07ca2540.0018200.004710.73na11.40213.0060.58cu6.3135.208.28s1580.004530.00841.02k14000.0020100.001927.12la0.0015.705.21sn23.0033.704.02mg2040.004540.00755.56sc0.308.412.44ni2.9311.302.52在对数据描述性统计后，数量级别差别便暴露了出来，的数量级相比其他的数据大。若不进行标准化处理，在进行主成分分析时，数据的级别差将会影响主成分的提取。进行标准化处理，使得数据的均值为0，方差为1，可以满足正态分布检验。这也为主成分的提取做好了准备。3.2 主成分的提取通过软件的数据分析，我们将原始数据多变量进行降维处理，得到数据的特征根和方差贡献率。从下表特征值及贡献率可以看出，各成分的贡献率依次递减，前4个成分的累积贡献率达87.461%，大于85%，它表示前4个成分已经可以解释87.461%的数据信息，因此可以萃取前4个为主成分。明显地，第一个主成分的特征根与第二个主成分的特征根相差很大，但仅取一个成分只能解释50.784%的原始数据，信息遗漏过大。表 2 特征根和方差贡献率表序号特征根方差贡献率（%）累积贡献率（%）18.63350.78450.78423.25419.14369.92631.92711.33381.25941.0546.20287.46150.9395.52592.98760.5433.19696.18270.3241.90698.08880.2361.38799.47590.0890.525100100010011001001200100130010014001001500100160010017001003.3 主成分载荷矩阵主成分的载荷矩阵表示某一成分中各变量与成分的相关关系，正值表示正相关，负值表示负相关。表 3 因子载荷阵zscore(ba)0.9750.0920.070.134zscore(co)-0.681-0.080.1850.674zscore(cr)0.1160.757-0.565-0.004zscore(li)0.866-0.008-0.2450.154zscore(zn)-0.6010.7290.1510.06zscore(fe)-0.0990.344-0.7580.121zscore(mn)-0.344-0.7320.2550.264zscore(cu)-0.6280.4540.1280.529zscore(na)0.9180.1460.0990.256zscore(ca)0.9860.0230.0360.142zscore(mg)0.987-0.030.030.115zscore(k)-0.609-0.582-0.356-0.027zscore(s)-0.7940.3760.253-0.046zscore(la)0.7710.4220.2770.044zscore(sc)0.9910.0020.0250.106zscore(sn)-0.4120.7010.06-0.081zscore(ni)0.1810.3220.738-0.261根据表2，可以看出主成分覆盖了反映茶叶微量元素价值程度的各个维度 王鹏，李贵全.不同大豆杂交后代新品系产量性状的主成分分析j.山西农业大学学报（自然科学版）,2012,32（5）:408-411.，且各主成分相对独立。由第一主成分可以看出具有较大的载荷，说明与这几种元素相关性较大。同时，我们发现的载荷的绝对值也较大，说明与第一主成分负相关关联性大。第二主成分中，具有较高的载荷，说明与这三种正相关大，而载荷的绝对值大，说明与第二主成分的负相关关联性大。第三主成分中，具有较高的正载荷，有较高的负载荷。第四主成分中，具有较高的正载荷，表明这两种元素与第四主成分关联性大。3.4 特征向量矩阵根据因子载荷阵与特征根，求出特征向量矩阵。表 4 特征向量矩阵z1z2z3z40.3320.0510.0510.131-0.232-0.0450.1330.6570.0400.420-0.407-0.0040.295-0.004-0.1760.150-0.2050.4040.1090.059-0.0340.190-0.5460.118-0.117-0.4060.1840.257-0.2140.2510.0920.5160.3130.0810.0710.2500.3360.0130.0260.1390.336-0.0170.0210.112-0.207-0.323-0.257-0.026-0.2700.2080.182-0.0450.2630.2340.1990.0430.3370.0010.0180.103-0.1400.3880.043-0.0790.0610.1780.532-0.255通过特征向量矩阵的计算，便于我们得到主成分的表达式：3.5 综合得分及排名根据，计算综合得分并对茶叶的价值给予排名。表 5 茶叶微量元素主成分对茶叶价值的综合得分与排名茶叶排名海南山苦茶（野生）8.035-0.0700.9420.4974.811南海白茶绿茶-1.519-1.256-0.4180.297-1.1910南海佳人（有机乌龙茶）-0.6151.331-0.262-0.423-0.133巴陵春春绿茶-0.3091.2470.770-0.7030.142乌龙茶-0.731-1.615-1.5870.473-0.959福建铁观音0.819-1.901-1.819-0.875-0.244武功绿茶-1.397-0.500-0.4530.845-0.928井冈翠绿-1.169-0.3320.542-0.915-0.757庐山云雾-1.3091.2061.504-1.287-0.396婺源毛尖-1.8071.8900.7812.091-0.395根据四个主成分的表达式与主成分的权数我们对各种茶叶计算得分 徐敏,万德光.川黄柏中微量元素的主成分分析和聚类分析j.时珍国医国药,2008 ,19 (1) 151153。海南山苦茶最好；其次为巴陵春春绿茶，南海佳人，福建铁观音，婺源毛尖，庐山云雾，井冈翠绿，武功绿茶，乌龙茶，得分最低的为南海白茶绿茶。根据得分不难发现海南山苦茶与其余品种茶叶得分相差很大，这也说明野生茶叶在微量元素含量上最高，价值更大。这也符合客观事实，野生的茶叶价值大于人工培养的茶叶，在第一主成分得分上，海南山苦茶具有明显的优势，它所含的主要微量元素符合人体所需的矿物质元素，可以说越是对人体有益的其价值的凸显越明显。四、总结及建议现代医学表明钙、镁、钠、锰、锌、铁、