1、数学建模校内竞赛论文 论文题目：食品质量安全抽检数据分析组号：284#成员：王林炜 陈小红 葛敏雪 选题：A 题 姓名 学院 年级 专业 学号 联系电话 数学分析 高等代数 高等数学 线性代数 概率统计 数学实验 数学模型 CET4 CET6 王林炜 自动化 2011 物流工程 20115172 15922768472 / / 86 86 / 85 / 567 / 陈小红 自动化 2011 物流工程 20115223 13637985570 / / 95 86 / 92 / / / 葛敏雪 自动化 2011 物流工程 20115174 13512354364 / / 91 90 82 / /
2、519 / 2013-05-211食品质量安全抽检数据分析摘要“民以食为天” ，食品安全关系千家万户，本文根据深圳市给出 2010年、2011年、2012 年的食品质量安全抽检数据，针对深圳市的食品质量安全做出了分析。通过对数据的研究，我们对数据进行了处理，将主要食品分为六大类，分别是：蔬果类、肉蛋类、水产类（、淀粉类、调味类、豆乳类。在建立数学模型时，应用了数据筛选和均值统计法、层次分析法、回归分析法、数据拟合法、概率统计分析。针对问题一，我们将影响食品安全的因素分为微生物、重金属、添加剂、其他这四大类。先用数据筛选和均值法对数据进行处理，从而得到 2010年、2011年和 2012年每个季
3、度食品因微生物、重金属、添加剂、其他而导致的不合格率，并用图表表示出每一类三年的变化趋势。然后用层次分析法，定量地求出三年食品安全系数。三年内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。年份 2010 2011 2012食品安全系数0.1932 0.5676 0.2238等级 针对问题二，由于微生物、重金属、添加剂、其他这四方面的不合格导致了食品质量的不过关，但是这四类直接因素又由其它因素影响，通过对数据的分析，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。通过建立回归模型，运用 MATLAB软件，求出六大要素与食
4、品质量间的线性关系，定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。针对问题三，根据第一问和第二问，我们可以得到不同因素的影响程度是不同的。为了科学有效地反映食品安全情况，合理简化抽检方法，我们减少稳定数据的测量量，增加不稳定数据的测量量，使测量环节有效减少，重要环节的监测力度加大。通过使用 MATLAB软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于 0，新方法具有科学性。最终得到，新的抽检方法为减
5、少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数关键词：统计 层次分析 回归分析 拟合 2一、 问题重述“民以食为天” ，食品安全关系到千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市
6、之一。请下载 2010年、2011 年和 2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取） ，并根据这些资料来讨论：1. 如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2. 从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3. 能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的） ，例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？二、问题分析对于问题一，由于调查中食物种类复杂，为了方便问题的讨论，我们
7、将食品分为六大类，即淀粉类，水产类，蔬果类，肉蛋类豆乳类和调味品，对数据进行筛选。在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属、添加剂和其他。从附录中得到深圳市2010，,2011,2012 三年抽检样本数据，并用 EXCEL制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。根据相应的数据制作出每个因素安全情况三年的变化趋势。再用层次分析法，从定量的角度求出三年的总体变化趋势。对于问题二，第一问中涉及的影响因素是食品是否安全的直接因素，但这些因素又是由其它相关因素决定。分析数据，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽
8、查季节、生产环节、流通环节、餐饮环节六个3要素。这一问要求的是这些要素与食品质量的规律性关系，即线性关系。通过建立回归模型，来定量求出各因素在影响食品安全中所占影响度大小和影响规律。对于各因素对食品安全的影响程度，利用偏回归平方和 衡量每个变量在iQ回归中所起的作用大小（即影响程度） ，对每个变量 的偏回归平方和 进行ixi了计算，最后把影响程度由深到浅的各因素进行了排序。最后进行回归方程的显著性检验，将总的平方和 分解为回归平方和 和剩余平方和 ，运用 F总S回S剩S检验法来判定线性回归方程的显著性。对于问题三，要改进原有抽检方法，找到一个新方法，使其在成本相同的情况下，测得数据结果更准确，
9、 ，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据，得到不同因素的影响程度是不同的。减少稳定数据的测量量，加大不稳定数据的测量量，减少测量环节，加重必要环节。由于我们获取的数据均是在原抽检方法下得到的数据，而没有真正标准真实的数据，所以无论我们用什么样的新方法，只能在原抽检下的数据（或有关预测、推测数据）进行验证。故新方法在测量不稳定数据的结果与原方法不具可比性。但稳定数据的结果两种方法下的测量数据基本一致。新方法得到后，检验新方法的科学性、合理性。拟合曲线确定其未来检测值，采用蒙特卡洛法抽检，带入回归方程进行验证新方法是否合理。三、问题假设1、 假设模型求解过程中所用的
10、数据都是合理的2、假设不考虑抽检时抽检人员的人为影响3、假设每次抽检时遵循公平性原则4、假设所有食品取主要食品即六大类，其余食品对质量安全分析无影响5、假设抽检时所有环节、所有因素的成本和工时相同4、参数及符号说明4符号 表示的意义 备注Ai 目标层指标 i=1Bi 准则层指标 i=1,2,3,4Ci 方案层指标 i=1,2,3Wi第 i 层权重 i=1,2,3i0标准化后的第 i 层权重 i=1,2,3入 一致性指标过度数C.I 一致性指标R.I. 平均随机一次性指标C.R. 一致性比例指标bij 每个因子的评分 i=1,2,3j=1,2,3,4Ci 综合评价指数 i=1,2,3PI 食品安
11、全质量综合评价指数ix影响食品质量安全的因素。 i=1,2,.6y各食品安全系数 a=1,2.6i相互独立且服从同一正态分布的随机变量 1,2.38iiQ变量 的偏回归平方和ixi=1,2,.6总S回归平方和回剩余平方和i待估计系数参数 1,2.8iy的回归值yib参数 的最小二乘估计i,.i5五、 模型的分析、建立与求解（1）问题一的分析与模型的建立1、统计、计算、绘制图表来表示出六大类食品，在2010、2011、2012三年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势；通过数据筛选和统计得出附表一，根据附表一，计算求解出各个季度下 微生物、重金属、添加剂、其它等因素不
12、合格的概率年份 2010 2011 2012季度 二 三 四 一 二 三 四 一 二 三 四微生物 0.033 0.051 0.141 0.105 0.02 0.187 0.342 0.067 0.059 0.036 0.074重金属 0.039 0.03 0 0.389 0.149 0.178 0.162 0 0.122 0.077 0.134添加剂 0.038 0.03 0 0.117 0.038 0.178 0.119 0.098 0.086 0.037 0.092其它 0.05 0.048 0.141 0.041 0.015 0 0.027 0.102 0.067 0.023 0.02
13、86后根据表格，绘制图表来表示三年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势：从图表得出结论：三年来，因微生物 、重金属、添加剂导致的不合格率变化趋势有波动，食品安全整体有下降，从总体上说，食品安全情况没有改善。2、 根据三年中六类食品的食品安全情况，采用层次分析法，定量结合定性的评价出由于微生物、重金属、添加剂超标而导致食品不安全的年度食品安全情况。对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况各因素不合格的概率年平均 2010 2011 2012微生物 0.075 0.164 0.059重金属 0.023 0.220 0.083添加剂 0.023
14、0.113 0.078其它 0.080 0.021 0.0557的变化趋势做出定量的综合评价，首先对3年的数据处理按季节分为12个子样本点的抽检值进行分析，得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。建立层次结构模型（1）构造各季节食品安全质量综合评价的阶梯层次结构根据层次分析法的基本步骤,建立食品质量安全评价递阶层次模型,如下图所示。图：食品综合评价递阶层次结构（2）建立两两比较的 A-B层判断矩阵并进行一致性检验。建立 A-B层判断矩阵的依据是：考虑 B层的四个因素对食品质量安全综合评价的重要性时，根据在查找数据知年均指标微生物最重要，且 B1比重金属
15、、添加剂、其它含量等值重要。而在重金属、添加剂、其它含量等值中，而各类食品都用到添加剂，所以，添加剂比重金属、其它重要。而其它比重金属重要。利用 19 标度法(见下表)进行成对比较,同时参考专家意见,确定各因素之间的相对重要性并赋以相应的分值,构造出各层次中的所有判断矩阵,并计算权向量和一致性检验。表 2 19 标度的意义标度 aij 意义1 Ci 与 Cj 的影响相同3 Ci 比 Cj 的影响稍强5 Ci 比 Cj 的影响强7 Ci 比 Cj 的影响明显地强9 Ci 比 Cj 的影响绝对地强2 ,4 ,6 为上述两判断级的中食品安全综合评价A1 AA微生物 B1 重金属 B2 添加剂 B3
16、其它 B42010 年食品安全情况 C1 2012 年食品安全情况C32011 年食品安全情况C2目标层 A准则层 B方案层 C8,8 间值1 ,1/ 2 , ,1/ 9 Ci 较 Cj 的影响之比与上述说明相反1、 根据判断矩阵标度内容及其含义，可得 A-B层判断矩阵 B： B1 B2 B3 B41 4 2 3 B1 1/4 1 1/2 3/4 B2B= 1/2 2 1 3/2 B31/3 4/3 2/3 1 B4 其中 B矩阵中的每一个元素 Bij根据调查相关数据定性判断出：2、用方根法计算出判断矩阵 B的权重值并进行归一化：方根法（几何平均法）：Wi= i=1,2,3,.,n归一化： =
17、计算步骤：1、B 中的元素按行相乘得一新向量；2、将新向量的每一个分向量开 n次方；3、将所有向量归一化即为权重向量。3、计算判断矩阵 B的最大特征值 并进行一致性检验.maxax1()niiBWw计算一致性指标 C.I 其中，n 为判断矩阵的阶数max.1CI计算一致性比例 其中 为平均一致性指标，当 时判.CIR.I .01CR断矩阵具有好的一致性。找查相应的平均随机一致性指标 R.I。平均随机一致性指标n 1 2 3 4 5 6 7 8 9 10njbi)(1inji1i0njji19经计算得图表：(3)建立两两比较的B-C层判断矩阵并进行一致性检验。根据途径一的表格中数据 （表 1-1
18、）可完成表格：方法和计算 B判断矩阵一致。R.I 0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49A B1 B2 B3 B4 Wi 入 miB1 1 4 2 3 2.21 0.48 4.00B2 1/4 1 1/2 3/4 0.55 0.12 4.02B3 1/2 2 1 3/2 1.11 0.24 3.99B4 1/3 4/3 2/3 1 0.74 0.16 3.99=4C.I.=0R.I.=0.89C.R.0.1B1 C1 C2 C3 Wii0入 miC1 1 0.457 1.271 0.834 0.25 3.00C2 2.187 1 2.780 1.825 0.55 3.00C3 0.787 0.360 1 0.657 0.20 3.00=3C.I.=0R.I.=0.52C.R.0.1 B2 C1 C2 C3 ii0入 miC1 1 0.105 0.277 0.308 0.07 3.00C2 9.565 1 2.651 2.938 0.67 3.00C3 3.609 0.377 1 1.108 0.25 3.00=3C.I.=0R.I.=0.52C.R.0.1入 axi0入 ax入 ax
