基于主成分分析质量安全评价模型构建.docVIP

基于主成分分析质量安全评价模型构建 　　摘要：本文针对分析深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势的问题，首先确定评价食品安全的检测指标，建立评价指标体系，本文利用已有数据得出三年主要食品领域如蔬菜、肉制品以及水产品领域，微生物、重金属、添加剂含量以及其他影响食品安全因素检查的合格率；再利用主成分分析方法建立的数据矩阵进行求解，得出影响食品质量安全的主成分因子，以及相应的方差贡献率（又称信息贡献率）作为主成分指标变量的权重，进而定义食品安全综合指数FSCI（Food Safety Comprehensive Index ，FSCI），通过食品安全综合指数的比较，对三年来深圳市的食品安全水平进行综合评价并得出结论：在蔬菜、肉类以及水产品三个主要食品领域中，食品质量安全水平2012最优，2011在三年来相对较差，2010年居中。 　　关键词：主成分分析 信息贡献率 特征向量 综合评价指数 　　中图分类号：R155.5 文献标识码：A 文章编号：1672-5336（2014）06-0064-03 　　针对评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势的问题，我们首先确立指标体系，并对相关数据进行标准化处理，建立如下基于主成分分析的综合评价模型，并利用Matlab编程进行求解。 　　假设进行主成分分析的指标变量有个[1]，分别为，共有个评价对象，第个待评对象的第个指标值为，得到观测数据矩阵：。 　　将各指标值转换成标准化的指标值，有： 　　其中： 　　其中即和为第个指标的样本均值和样本标准差。对应地，称为标准化的指标变量。 　　根据公式概率统计以及线性代数的相关知识[2]可以知道：即为标准化的指标变量和的相关系数矩阵，相关系数矩阵中的元素 　　其中是第个指标与第个指标的相关系数，显然易见相关系数矩阵为是对称矩阵。 　　由上述相关系数矩阵的性质知道：相关系数矩阵为是实对称矩阵，因此相关系数矩阵进行特征分解得到[3]：，其中，是由相关系数矩阵的特征值组成的对角阵，是由的标准正交化的特征向量按列并排组成的正交矩阵，称为主成分载荷矩阵，它是用主成分指标变量表示中心标准化的原指标变量时的系数矩阵，即用矩阵可表示为：。由于主成分载荷矩阵是正交矩阵，所以满足，故又有，即可知：，可见的转置矩阵是用原变量表示主成分变量时的系数矩阵。 　　另外还可以知道矩阵为主成分得分阵，每行的数据就是对应各次观测对应的主成分变量的取值。因而主成分的样本协方差矩阵为 　　，由此可见主成分指标变量的样本方差就是特征值，且不同主成分之间的样本协方差都等于零。 　　为计算特征值的信息贡献率和累计贡献率，定义为主成分指标变量的信息贡献率；同时，有为主成分指标变量的累计贡献率。当接近于1（一般取）时，则选取前个指标变量作为主成分指标变量，代替原来的个指标变量，从而可以对个主成分指标变量进行综合评价[4]。 　　具体到本问题中时，我们定义了食品安全综合指数FSCI； 　　其中：为第个主成分指标变量的信息贡献率，根据食品安全综合指数FSCI，我们就可以实现对2010、2011以及2012年三年来深圳市主要食品领域食品质量安全的变化趋势。 　　为综合评价深圳市这2010、2011以及2012三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势，我们确定了由微生物检测合格率、重金属检测合格率、添加剂含量检测合格率以及其他因素检测合格率四项所组成的评价指标体系，2010、2011、2012年食品质量安全水平组成的评价对象。为综合考虑三年的食品质量安全水平，我们主要以蔬菜、肉制品以及水产品三个食品领域作为指标变量数据获取的样本空间[5]，得到如下的表1，表1中各项数据为三个主要食品领域中不同年份四项检测指标的合格率情况。 　　注：此处仅考虑蔬菜、肉制品以及水产品三个主要食品领域的各项指标。 　　因此，指标变量主要食品领域微生物检测指标的合格率、主要食品领域重金属检测指标的合格率、主要食品领域添加剂检测指标的合格率、主要食品领域其他检测指标的合格率，得到相应的数据矩阵[6]为 　　对相关系数矩阵进行特征值和相应的特征向量求解，并计算各个特征值的贡献率和累计贡献率，如下表2： 　　由表2得到的主成分分析结果，可以利用主成分的贡献率作为综合评价指标的权重系数。由表中数据可以得知前两个主成分指标变量的累计贡献率就达到了100%，主成分分析效果较好。 　　主成分载荷矩阵，是由的标准正交化的特征向量按列并排组成的正交矩阵，即；进而可以得到表3――标准化的4个主成分指标变量对应的特征向量[7]。 　　由此可得2个主成分变量分别为： 　　据主成分系数可以看出第一主成分主要反映了前两个指