XXXX年深圳杯数学建模竞赛A题食品质量安全抽检数据分析

..答卷编号：论文题目：A题：食品质量安全抽检数据分析组别：本科生参赛队员信息(必填)：姓名学院、班级学号联系电话参赛队员1教育实验学院实验1101班参赛队员2教育实验学院实验1101班参赛队员3教育实验学院实验1101班参赛学院：教育实验学院1A题：食品质量安全抽检数据分析摘要“民以食为天”，食品安全问题越来越引起社会各界的重视，因此食品的抽检对了解食品安全情况就起到了非常重要的作用，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节都可能影响食品的质量与安全。本文主要对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；食品产地与食品质量的关系，食品销售地点（即抽检地点）与食品质量的关系，季节因素与食品质量的关系；以及如何改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的）等问题进行了分析研究，建立了相应的数学模型，运用了SPSS及MATLAB软件工具对模型进行了求解。对于问题一，首先将三年的不合格数据进行统计分类，由相关标准将不合格食品按主要食品领域分为十类，将不合格的项目分为微生物、重金属、添加剂、食物固有成分四类。然后对三年不合格主要食品按照此类别进行筛选，计算出每年各主要食品领域中每类不合格项目在总的不合格项目中所占比例，并根据此比例对年份做出折线图，由此得到食品安全情况的变化趋势。对于问题二，首先本文运用统计学的方法把三年来食品的产地、抽检地点、季节因素进行了分类并统计。然后运用归一化原理分别计算出了每年各个食品产地、抽检地点、季节因素占总不合格数的比例。再对这些比值进行K-均值聚类分析，聚为三类，由此把这三个因素对食品质量的影响分为良好、一般、严重三个等级，以表示食品产地、抽检地点、季节因素与食品质量的关系。对于问题三，首先将所有食品进行分类，然后运用了统计学的方法统计出了每年在各主要食品领域抽检的总数目以及其中的合格数、不合格数，并计算出各主要食品领域的不合格率，再配合问题一中所统计出的各不合格项目在该食品领域所占的比例，得到了各主要食品领域不合格项目的不合格率，再以此不合格率为基础建立基于实际数据的层次分析法来确定各主要食品领域和不合格项目的权重，最后基于此权重来调整食品的抽检方法。关键词：食品的分类统计学归一化折线图K-均值聚类分析基于实际数据的层次分析法权重21问题重述国以民为本，民以食为天，食以安为先。食物是人类赖以生存和发展的基本物质条件，也是国家安定，社会发展的根本要素。在任何一个国家，食品质量及其安全性是百姓共同关注的一个永恒主题。食品安全不仅是社会关注的热点，也是政府民生的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品比例也越来越高，因此食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。根据题意以及对附件表中的数据分析，本文需要解决的问题有：（1）如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；（2）从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；（3）能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？2模型假设（1）假设食品能且仅能分为十大类，其他没有被分类的食品对食品安全质量所造成的影响忽略不计。（2）假设影响食品质量的因素能且仅能分为四大类，其他没有被分类的因素对食品质量所造成的影响忽略不计。（3）食品抽检数据上传时间的季节即为抽检季节。（4）对于食品产地，只有来自于深圳八个区的数据才为有效数据。（5）假设所有食品生产商的信誉度均相同。（6）假设抽检不受国家相关政策的影响。3符号说明A………………目标层；C………………准则层；P………………措施层；ija………………准则层两个因素ic和jc对目标层的影响程度之比；..CI……………判断矩阵一致性指标；ijb………………方案层两个因素ip和jp对准则层的影响程度之比；..CR……………随机一致性比率；max……………最大特征值；3..RI……………平均随机一致性指标；4模型的建立与求解4.1问题一首先，按照相关标准把主要食品领域分为了:[1]粮食及其制品、肉类及其制品、食用油、豆制品、水产品、调味品、蔬果及其制品、冷食饮料酒水乳制品、糕点、糖果蜜饯这十类，把不合格项目分为了：微生物、重金属、添加剂及食物固有成分四类。然后分别把每年的不合格数据进行汇总，并按照食品类别进行筛选，统计出每类食品总的不合格数目，再进一步计算出该食品类别中各不合格项目的个数及其在总不合格数目中所占比例，然后根据此比例和年份作出了深圳市这三年来各主要食品领域微生物、重金属、添加剂，食物固有成分安全情况的变化趋势图。本文在下面列举出了豆制品、粮食及其制品的图表。例：三年中豆制品不合格项目的数量及比例，如表1：表1：三年中豆制品不合格项目的数量及比例201020112012豆制品类别个数比例类别个数比例类别个数比例微生物2450%微生物7459.2%微生物433.3%重金属00重金属00重金属00添加剂2143.75%添加剂3931．2%添加剂541.67%食物固有成分36.25%食物固有成分129.6%食物固有成分325%总数48100%总数125100%总数12100%根据表1豆制品中微生物、重金属、添加剂、食物固有成分所占的比例用[2]MATLAB做出折线图，如图1：图1：豆制品中不合格项目的变化趋势图2010201120120102030405060比例年份微生物添加剂食物固有成分4由图1分析可知这三年豆制品中微生物对食品安全情况的影响是先增加后减少；添加剂对食品安全情况的影响是先减少后增加；食物固有成分对食品安全情况的影响是一直增加；豆制品无重金属问题出现。例：三年中粮食及其制品中不合格的数量及比例如表2：表2：三年中粮食及其制品中不合格的数量及比例201020112012粮食及其制品类别个数比例类别个数比例类别个数比例微生物836.36%微生物178.17%微生物1317.81%重金属313.64%重金属12861.54%重金属2230.14添加剂731.82%添加剂5928.37%添加剂3852.05%食物固有成分418.18%食物固有成分41.92%食物固有成分00%总数22100%总数208100%总数73100%根据表3粮食及其制品中微生物、重金属、添加剂、食物固有成分所占的比例做出折线图，如表4：图2：粮食及其制品中不合格项目的变化趋势图201020112012010203040506070比例年份微生物重金属添加剂食物固成分有由图2分析可知这三年粮食及其制品中微生物对食品安全情况的影响是先减少后增加；重金属对食品安全情况的影响是先增加后减少；添加剂对食品安全情况的影响是先减少后增加；食物固有成分对食品安全情况的影响是一直减少。剩余图表见附录1。本文通过对附录1中深圳市这三年来剩余的八种主要食品领域得到的图和表分析可得如下结论：（1）这三年糖果蜜饯中微生物对食品安全情况的影响是先增加后减少；重金属对食品安全情况的影响是一直减少；添加剂对食品安全情况的影响是基本上没有发生变化；食物固有成分对食品安全情况的影响是一直增加。5（2）这三年糕点中微生物对食品安全情况的影响是一直增加；重金属对食品安全情况的影响是一直减少；添加剂对食品安全情况的影响是先减少后增加；食物固有成分对食品安全情况的影响先增加后减少。（3）这三年冷食饮料酒水乳制品中微生物对食品安全情况的影响是一直减少；重金属对食品安全情况的影响是先增加后减少；添加剂对食品安全情况的影响是一直增加；食物固有成分对食品安全情况的影响是先减少后增加。（4）这三年调味品中微生物对食品安全情况的影响是一直减少；添加剂对食品安全情况的影响是一直减少；食物固有成分对食品安全情况的影响是一直增加。（5）这三年水产品中微生物对食品安全情况的影响是先减少后不变；重金属对食品安全情况的影响是先增加后不变；添加剂对食品安全情况的影响是先增加后不变。（6）这三年蔬果及其制品中微生物对食品安全情况的影响是先增加后减少；重金属对食品安全情况的影响是先增加后减少；添加剂对食品安全情况的影响是先不变后增加；食物固有成分对食品安全情况的影响是一直减少。（7）这三年食用油中微生物对食品安全情况的影响是一直增加；添加剂对食品安全情况的影响是先增加后减少；食物固有成分对食品安全情况的影响是一直减少。（8）这三年肉类及其制品中微生物对食品安全情况的影响是一直减少；重金属对食品安全情况的影响是先增加后减少；添加剂对食品安全情况的影响是一直增加；食物固有成分对食品安全情况的影响是一直增加。4.2问题二4.2.1食品产地与食品质量的关系假设只有产地为深圳宝安、福田、龙岗、罗湖、南山、光明新区、坪山新区这八个区的数据有效。先统计出各年不合格食品的产地，再对数据进行筛选分类，分别得到这八个区的不合格数目及其在总数目中所占比例。分别见表3、表4：表3：产地不合格数目统计表产地不合格数目统计表201020112012深圳宝安5111846深圳福田195719深圳光明新区8286深圳龙岗5112041深圳罗湖136822深圳南山154511深圳坪山新区074深圳盐田6228总和163465157表4：产地不合格比例表（%）产地不合格比例表（%）201020112012深圳宝安31.2925.3829.29深圳福田11.6612.2612.10深圳光明新区4.916.023.82深圳龙岗31.2925.8026.11深圳罗湖7.9714.6214.02深圳南山9.209.677.016深圳坪山新区0.001.512.55深圳盐田3.684.745.10总比例100.00100.00100.00对表4数据用[3]SPSS进行K-均值聚类分析，分析结果如下：IterationHistoryaIterationChangeinClusterCenters12313.0973.1843.9322.000.000.000FinalClusterCentersCluster123地点28114NumberofCasesineachClusterCluster16.00028.000310.000Valid24.000Missing.000得到分类结果见附录2表10，经过数据分析得如下结果：InitialClusterCentersCluster123地点311507表5：分析结果良好一般严重坪山新区罗湖宝安盐田福田龙岗光明新区南山产地为宝安和龙岗的食品质量问题最为严重，罗湖、福田和南山的食品质量问题一般，而坪山新区、盐田和光明新区的食品质量问题良好。4.2.2食品销售地点（抽检地点）与食品质量的关系将抽检地点具体到深圳的八个区：宝安、福田、龙岗、罗湖、南山、盐田、坪山新区和光明新区。再将数据筛选统计，分别得到这八个区的不合格数目及其在总数目中所占比例。分别见表6、表7：表6：销售地点不合格数目统计表销售地点不合格数目统计表宝安福田光明新区龙岗罗湖南山坪山新区盐田总和2010693897022294122532011273100402551041012540938201211832128853391214368表7：销售地点不合格比例表（%）销售地点不合格比例表（%）宝安福田光明新区龙岗罗湖南山坪山新区盐田总和201027.2715.023.5627.678.7011.461.584.74100.00201129.1010.664.2627.1911.0910.772.674.26100.00201232.078.703.2623.9114.4010.603.263.80100.00对表7数据用SPSS进行K-均值聚类分析，分