食品质量安全抽检数据分析论文

数学建模校内竞赛论文论文题目：食品质量安全抽检数据分析组号：284#成员：王林炜陈小红葛敏雪选题：A题姓名学院年级专业学号联系电话数学分析高等代数高等数学线性代数概率统计数学实验数学模型CET4CET6王林炜自动化2011物流工程2011517215922768472//8686/85/567/陈小红自动化2011物流工程2011522313637985570//9586/92///葛敏雪自动化2011物流工程2011517413512354364//919082//519/2013-05-211食品质量安全抽检数据分析摘要“民以食为天”，食品安全关系千家万户，本文根据深圳市给出2010年、2011年、2012年的食品质量安全抽检数据，针对深圳市的食品质量安全做出了分析。通过对数据的研究，我们对数据进行了处理，将主要食品分为六大类，分别是：蔬果类、肉蛋类、水产类（、淀粉类、调味类、豆乳类。在建立数学模型时，应用了数据筛选和均值统计法、层次分析法、回归分析法、数据拟合法、概率统计分析。针对问题一，我们将影响食品安全的因素分为微生物、重金属、添加剂、其他这四大类。先用数据筛选和均值法对数据进行处理，从而得到2010年、2011年和2012年每个季度食品因微生物、重金属、添加剂、其他而导致的不合格率，并用图表表示出每一类三年的变化趋势。然后用层次分析法，定量地求出三年食品安全系数。三年内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。年份201020112012食品安全系数0.19320.56760.2238等级ⅠⅤⅡ针对问题二，由于微生物、重金属、添加剂、其他这四方面的不合格导致了食品质量的不过关，但是这四类直接因素又由其它因素影响，通过对数据的分析，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。通过建立回归模型，运用MATLAB软件，求出六大要素与食品质量间的线性关系，定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。针对问题三，根据第一问和第二问，我们可以得到不同因素的影响程度是不同的。为了科学有效地反映食品安全情况，合理简化抽检方法，我们减少稳定数据的测量量，增加不稳定数据的测量量，使测量环节有效减少，重要环节的监测力度加大。通过使用MATLAB软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于0，新方法具有科学性。最终得到，新的抽检方法为减少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数关键词：统计层次分析回归分析拟合2一、问题重述“民以食为天”，食品安全关系到千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取），并根据这些资料来讨论：1.如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2.从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3.能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？二、问题分析对于问题一，由于调查中食物种类复杂，为了方便问题的讨论，我们将食品分为六大类，即淀粉类，水产类，蔬果类，肉蛋类豆乳类和调味品，对数据进行筛选。在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属、添加剂和其他。从附录中得到深圳市2010，,2011,2012三年抽检样本数据，并用EXCEL制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。根据相应的数据制作出每个因素安全情况三年的变化趋势。再用层次分析法，从定量的角度求出三年的总体变化趋势。对于问题二，第一问中涉及的影响因素是食品是否安全的直接因素，但这些因素又是由其它相关因素决定。分析数据，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。这一问要求的是这些要素与食品质量的规律性关系，即线性关系。通过建立回归模型，来定量求出各因素在影响食品安全中所占影响度大小和影响规律。对于各3因素对食品安全的影响程度，利用偏回归平方和iQ衡量每个变量在回归中所起的作用大小（即影响程度），对每个变量ix的偏回归平方和iQ进行了计算，最后把影响程度由深到浅的各因素进行了排序。最后进行回归方程的显著性检验，将总的平方和总S分解为回归平方和回S和剩余平方和剩S，运用F检验法来判定线性回归方程的显著性。对于问题三，要改进原有抽检方法，找到一个新方法，使其在成本相同的情况下，测得数据结果更准确，，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据，得到不同因素的影响程度是不同的。减少稳定数据的测量量，加大不稳定数据的测量量，减少测量环节，加重必要环节。由于我们获取的数据均是在原抽检方法下得到的数据，而没有真正标准真实的数据，所以无论我们用什么样的新方法，只能在原抽检下的数据（或有关预测、推测数据）进行验证。故新方法在测量不稳定数据的结果与原方法不具可比性。但稳定数据的结果两种方法下的测量数据基本一致。新方法得到后，检验新方法的科学性、合理性。拟合曲线确定其未来检测值，采用蒙特卡洛法抽检，带入回归方程进行验证新方法是否合理。三、问题假设1、假设模型求解过程中所用的数据都是合理的2、假设不考虑抽检时抽检人员的人为影响3、假设每次抽检时遵循公平性原则4、假设所有食品取主要食品即六大类，其余食品对质量安全分析无影响5、假设抽检时所有环节、所有因素的成本和工时相同四、参数及符号说明符号表示的意义备注4Ai目标层指标i=1Bi准则层指标i=1,2,3,4Ci方案层指标i=1,2,3Wi第i层权重i=1,2,3Wi0标准化后的第i层权重i=1,2,3入一致性指标过度数C.I一致性指标R.I.平均随机一次性指标C.R.一致性比例指标bij每个因子的评分i=1,2,3j=1,2,3,4Ci综合评价指数i=1,2,3PI食品安全质量综合评价指数ix影响食品质量安全的因素。i=1,2,..6y各食品安全系数a=1,2...6i相互独立且服从同一正态分布的随机变量1,2,...,38iiQ变量ix的偏回归平方和i=1,2,..6总S回归平方和回S剩余平方和i待估计系数参数1,2,...,8iˆyy的回归值ib参数i的最小二乘估计1,2,...,8i5五、模型的分析、建立与求解（1）问题一的分析与模型的建立1、统计、计算、绘制图表来表示出六大类食品，在2010、2011、2012三年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势；通过数据筛选和统计得出附表一，根据附表一，计算求解出各个季度下微生物、重金属、添加剂、其它等因素不合格的概率年份201020112012季度二三四一二三四一二三四微生物0.0330.0510.1410.1050.020.1870.3420.0670.0590.0360.074重金属0.0390.0300.3890.1490.1780.16200.1220.0770.134添加剂0.0380.0300.1170.0380.1780.1190.0980.0860.0370.092其它0.050.0480.1410.0410.01500.0270.1020.0670.0230.0286后根据表格，绘制图表来表示三年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势：从图表得出结论：三年来，因微生物、重金属、添加剂导致的不合格率变化趋势有波动，食品安全整体有下降，从总体上说，食品安全情况没有改善。2、根据三年中六类食品的食品安全情况，采用层次分析法，定量结合定性的评价出由于微生物、重金属、添加剂超标而导致食品不安全的年度食品安全情况。对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对3年的数据处理按季节分为12个子样本点的抽检值进行分析，得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。各因素不合格的概率年平均201020112012微生物0.0750.1640.059重金属0.0230.2200.083添加剂0.0230.1130.078其它0.0800.0210.0557建立层次结构模型（1）构造各季节食品安全质量综合评价的阶梯层次结构根据层次分析法的基本步骤,建立食品质量安全评价递阶层次模型,如下图所示。图：食品综合评价递阶层次结构（2）建立两两比较的A-B层判断矩阵并进行一致性检验。建立A-B层判断矩阵的依据是：考虑B层的四个因素对食品质量安全综合评价的重要性时，根据在查找数据知年均指标微生物最重要，且B1比重金属、添加剂、其它含量等值重要。而在重金属、添加剂、其它含量等值中，而各类食品都用到添加剂，所以，添加剂比重金属、其它重要。而其它比重金属重要。利用1～9标度法(见下表)进行成对比较,同时参考专家意见,确定各因素之间的相对重要性并赋以相应的分值,构造出各层次中的所有判断矩阵,并计算权向量和一致性检验。表21～9标度的意义标度aij意义1Ci与Cj的影响相同3Ci比Cj的影响稍强5Ci比Cj的影响强7Ci比Cj的影响明显地强9Ci比Cj的影响绝对地强2,4,6,8为上述两判断级的中间值1,1/2,⋯,1/9Ci较Cj的影响之比与上述说明相反食品安全综合评价A1AA微生物B1重金属B2添加剂B3其它B42010年食品安全情况C12012年食品安全情况C32011年食品安全情况C2目标层A准则层B方案层C81、根据判断矩阵标度内容及其含义，可得A-B层判断矩阵B：B1B2B3B41423B11/411/23/4B2B=1/2213/2B31/34/32/31B4其中B矩阵中的每一个元素Bij根据调查相关数据定性判断出：2、用方根法计算出判断矩阵B的权重值并进行归一化：方根法（几何平均法）：Wi=i=1,2,3,.....,n归一化：=计算步骤：1、B中的元素按行相乘得一新向量；2、将新向量的每一个分向量开n次方；3、将所有向量归一化即为权重向量。3、计算判断矩阵B的最大特征值max并进行一致性检验.max1()niiiBWnw计算一致性指标C.Imax..1nCIn其中，n为判断矩阵的阶数计算一致性比例......CICRRI其中..RI为平均一致性指标，当..0.1CR时判断矩阵具有好的一致性。找查相应的平均随机一致性指标R.I。经计算得图表：平均随机一致性指标n12345678910R.I000.520.891.121.261.361.411.461.49nnjbij)(1ninnjbij11Wi0njjiWW19(3)建立两两比较的B-C层判断矩阵并进行一致性检验。根据途径一的表格中数据（表1-1）可完成表格：方法和计算B判断矩阵一致。AB1B2B3B4Wi入mi=4C.I.=0R.I.=0.89C.R.0.1B114232.210.484.00B21/411/23/40.550.124.02B