数学建模食品质量安全抽检数据分析Word格式文档下载.docx

1、1. 评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2. 从这些数据中找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3. 改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定 且抽检频次过高的食品领域该作出的调整.二、 问题分析针对问题一，随着中国经济的快速发展，人们的生活质量逐渐提高，因而人们对于食品质量的要求也日渐增高，因此对食品进行产地、加工地等方面可能影响食品质量的因素进行抽检。又因为深圳是食品抽检、监督最统一、最规范、最公开

2、的城市之一，所以分析深圳市20102012三年个主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势，具有典型代表意义。从问题二可知，为了建立更合理的模型，我们需要考虑影响食品抽检不合格的其他因素，并进行相应的定性分析。深圳三年食品抽检的数据中给出了抽检地点，主要可分为八大区，抽检季度，四个季度。计算各地区各因素的不合格数，并可以借助SPSS软件对数据之间的关系进行分析和刻画，最后找到它们之间的规律和联系。针对问题三，为了能改进食品抽检的办法，可以考虑对影响食品质量的各因素进行层次分析，建立一个抽检个数与类别、地区的关系模型，从而有效地降低监管的成本。三、 基本假设1.假设主要食品仅能分

3、为六大类，其他没有被分类的食品对食品抽检的不合格性所造成的影响忽略不计；2.假设影响主要食物抽检不合格的因素主要有四大类，其他没有被分类的因素对食品抽检不合格性的影响忽略不计；3.假设抽样过程中其他地理或人为因素对抽样食品不合格性的影响忽略不计；4.假设对不同食品领域的各种指标是抽检是随机的；5.假设抽检间隔时间相同；6.假设检测的不同环节、不同因素的成本和工时相同。四、 定义与符号说明 符号含义成对比较矩阵判断矩阵权向量特征根一致性指标随机一致性指标自变量因变量参数误差项多重判定系数回归平方和偏回归平方和第层含有不合格产品的批次从第层内抽取的样本数五、 模型的建立与求解51问题一的模型5.1

4、.1数据处理为了提高数据处理过程中的效率，我们从原来的所有食品抽检数据中抽取了其中不合格产品的数据，以此作为参照数据并对它进行分类汇总和统计，然后分别得到2010年、2011年以及2012年这三年各个季度的统计数据。表1 三年来各个季度分类统计表年份因素第一季度第二季度第三季度第四季度2010年微生物412146添加剂289重金属463其他362011年2141437942051523611381116113481122012年65198581092124125.1.2模型的建立与求解 （1）进行图表分析通过上面的数据处理，我们得到了2010、2011和2012这三年在四个影响因素下的不合格产

5、品数，并以直方图的形式进行呈现。图1 2010年各因素引起的不合格产品数变化的趋势图图2 2011年各因素引起的不合格产品数变化的趋势图图3 2012年各因素引起的不合格产品数变化的趋势图从上面三幅图中，我们大致可以知道，微生物、添加剂、重金属以及其他因素不达标而造成的食品的不合格数三年来在不断地波动。（2）进行层次分析为了能更好地了解深圳市这三年来各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势，我们有必要对食品质量的因素进行深入分析。虽然我们已经知道了影响食品质量的因素是上述四方面，但我们仍不清楚它们到底是怎样影响食品的质量的。因此，我们还要注重研究它们在影响食品质量中所占的权

6、重，而在此情况下，层次分析法能得到很好的应用。层次分析法是是一种定性和定量相结合的、系统化的、层次化的分析方法，它的思路主要体现在分层上，从最底层开始分析各层对上一层的权重，一直到目标层，最后才综合得出最底层对目标层的总权重，从而能达到我们的预定目标。现在根据深圳市的实际情况，将主要食品分为六大类，通过层次分析模型（AHP），建立了三个层：食品安全综合评价、六类食品、四个影响因素，具体的分层结构图如下所示：图4 食品安全综合评价层次结构图构造层次分析模型的建立具体应该包括以下几个过程。.构造判断矩阵。判断矩阵表示针对上一层次某因素而言，本层次与之有关的各因素之间的相对重要性。其形式如下：，其中

7、表示对上一层A层而言，该B层中因素对的相对重要程度。在这里，我们通常使用1-9尺度可以方便地表示如表1.表2 1-9尺度的含义尺度C与C的影响相同C比C的影响稍强5C比C的影响强7C比C的影响明显的强C比C的影响绝对的强2,4,6,8C与C的影响之比在上述两个相邻等级之间1,1/2，,1/9C与C的影响之比为上面a的互反数针对本文中的问题一，通过以上的步骤建立模型之后，本文用成对比较法和1-9比较尺度对层次结构模型中的准则层对于目标层建立的6*6成对比矩阵为。.进行层次单排序和一致性检测。采用和法计算准则层的因素相对于目标层的层次单排序。首先，对的每一列向量归一化，得到对按行求和得：其次，将归

8、一化得最后，计算矩阵的最大特征根其中表示向量的第个元素。此外，还需要计算一致性指标：，其中。RI为平均随机一致性指标，当CR0.1，认为判断矩阵的一致性是可以接受的，否则应对判断矩阵作适当调整。.进行层次总排序。需要将计算出来的层次单排序再次进行适当计算。假设层次结构模型是由目标层（A）、准则层（B）和方案层（C）所组成，准则层有6个因素，方案层有4个因素。已知B层对A层的层次单排序为：以为列向量构成矩阵，则C层对B层的准则Bk的层次单排序为：同时，C层各方案对A层的层次总排序的方法为：更一般地，若共有s层，则第k层对第1层（设只有1个因素）的组合权向量，即层次总排序满足其中 。其中是以第k层

9、对第k-1的权向量为列向量组成的矩阵，于是最下层（第s层）对最上层的层次总排序为再根据上面的步骤进行逐步计算，最后可以得到2010年各主要食品受各种重要因素影响比重表为：表3 2010年各主要食品受各种重要因素影响比重表矩阵权重向量CIRICRA-B6.11400.02281.240.0184B1-P4.03290.01100.900.0122B2-P4.02760.00920.0102B3-P4.04350.01450.0161B4-P4.07340.02450.0272B5-P4.01040.00350.0038B6-P4.00420.00140.0015其中准则层相对于目标层的排序为方案

10、层相对于目标层的排序向量为（其中权重结果的排序依次是微生物，添加剂，重金属，其它。）.组合一致性检验。我们可以逐层进行组合一致性检验，若第p层的一致性的指标为（n是第p-1层因数的数目），随机一致性指标为。定义，则第p层的组合一致性比率为同时，第p层通过组合一致性检验的条件为从而最下层（第s层）对第1层的组合一致性比率为当且仅当适当的小时，才认为整个层次的比较判断通过一致性检验。我们对上述结果进行组合一致性收敛，得到第2层的组合一致性比率为而在第3层中，有,从而有,均通过检验。同时最下层（第3层）对第1层的组合一致性比率为组合一致性检验通过，前面得到的组合权向量W可作为最终决策影响安全因素权重依据。 同理，按照上述的方法分析得出2011年各主要食品受各种重要因素影响比重表。表4 2011年各主要食品受各种重要因素影响比重表6.13670.02730.02444.00780.00260.00293.9320-0.0227-0.02524.06970.02320.0258