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

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

《数学建模食品质量安全抽检数据分析Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《数学建模食品质量安全抽检数据分析Word格式文档下载.docx（29页珍藏版）》请在冰豆网上搜索。

1.评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；

2.从这些数据中找出某些规律性的东西：

如食品产地与食品质量的关系；

食品销售地点（即抽检地点）与食品质量的关系；

季节因素等等；

3.改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作出的调整.

二、问题分析

针对问题一，随着中国经济的快速发展，人们的生活质量逐渐提高，因而人们对于食品质量的要求也日渐增高，因此对食品进行产地、加工地等方面可能影响食品质量的因素进行抽检。

又因为深圳是食品抽检、监督最统一、最规范、最公开的城市之一，所以分析深圳市2010—2012三年个主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势，具有典型代表意义。

从问题二可知，为了建立更合理的模型，我们需要考虑影响食品抽检不合格的其他因素，并进行相应的定性分析。

深圳三年食品抽检的数据中给出了抽检地点，主要可分为八大区，抽检季度，四个季度。

计算各地区各因素的不合格数，并可以借助SPSS软件对数据之间的关系进行分析和刻画，最后找到它们之间的规律和联系。

针对问题三，为了能改进食品抽检的办法，可以考虑对影响食品质量的各因素进行层次分析，建立一个抽检个数与类别、地区的关系模型，从而有效地降低监管的成本。

三、基本假设

1.假设主要食品仅能分为六大类，其他没有被分类的食品对食品抽检的不合格性所造成的影响忽略不计；

2.假设影响主要食物抽检不合格的因素主要有四大类，其他没有被分类的因素对食品抽检不合格性的影响忽略不计；

3.假设抽样过程中其他地理或人为因素对抽样食品不合格性的影响忽略不计；

4.假设对不同食品领域的各种指标是抽检是随机的；

5.假设抽检间隔时间相同；

6.假设检测的不同环节、不同因素的成本和工时相同。

四、定义与符号说明

符号

含义

成对比较矩阵

判断矩阵

权向量

特征根

一致性指标

随机一致性指标

自变量

因变量

参数

误差项

多重判定系数

回归平方和

偏回归平方和

第层含有不合格产品的批次

从第层内抽取的样本数

五、模型的建立与求解

5．1问题一的模型

5.1.1数据处理

为了提高数据处理过程中的效率，我们从原来的所有食品抽检数据中抽取了其中不合格产品的数据，以此作为参照数据并对它进行分类汇总和统计，然后分别得到2010年、2011年以及2012年这三年各个季度的统计数据。

表1三年来各个季度分类统计表

年份

因素

第一季度

第二季度

第三季度

第四季度

2010年

微生物

41

21

46

添加剂

28

9

重金属

4

6

3

其他

36

2011年

214

14

37

94

205

15

23

61

138

11

16

113

48

112

2012年

65

19

8

58

10

92

1

24

12

5.1.2模型的建立与求解

（1）进行图表分析

通过上面的数据处理，我们得到了2010、2011和2012这三年在四个影响因素下的不合格产品数，并以直方图的形式进行呈现。

图12010年各因素引起的不合格产品数变化的趋势图

图22011年各因素引起的不合格产品数变化的趋势图

图32012年各因素引起的不合格产品数变化的趋势图

从上面三幅图中，我们大致可以知道，微生物、添加剂、重金属以及其他因素不达标而造成的食品的不合格数三年来在不断地波动。

（2）进行层次分析

为了能更好地了解深圳市这三年来各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势，我们有必要对食品质量的因素进行深入分析。

虽然我们已经知道了影响食品质量的因素是上述四方面，但我们仍不清楚它们到底是怎样影响食品的质量的。

因此，我们还要注重研究它们在影响食品质量中所占的权重，而在此情况下，层次分析法能得到很好的应用。

层次分析法是是一种定性和定量相结合的、系统化的、层次化的分析方法，它的思路主要体现在分层上，从最底层开始分析各层对上一层的权重，一直到目标层，最后才综合得出最底层对目标层的总权重，从而能达到我们的预定目标。

现在根据深圳市的实际情况，将主要食品分为六大类，通过层次分析模型（AHP），建立了三个层：

食品安全综合评价、六类食品、四个影响因素，具体的分层结构图如下所示：

图4食品安全综合评价层次结构图

构造层次分析模型的建立具体应该包括以下几个过程。

Ⅰ.构造判断矩阵。

判断矩阵表示针对上一层次某因素而言，本层次与之有关的各因素之间的相对重要性。

其形式如下：

，

其中表示对上一层A层而言，该B层中因素对的相对重要程度。

在这里，我们通常使用1-9尺度可以方便地表示如表1.

表21-9尺度的含义

尺度

Cｉ与Cｊ的影响相同

Cｉ比Cｊ的影响稍强

5

Cｉ比Cｊ的影响强

7

Cｉ比Cｊ的影响明显的强

Cｉ比Cｊ的影响绝对的强

2,4,6,8

Cｉ与Cｊ的影响之比在上述两个相邻等级之间

1,1/2，…,1/9

Cｉ与Cｊ的影响之比为上面aｉｊ的互反数

针对本文中的问题一，通过以上的步骤建立模型之后，本文用成对比较法和1-9比较尺度对层次结构模型中的准则层对于目标层建立的6*6成对比矩阵为

。

Ⅱ.进行层次单排序和一致性检测。

采用和法计算准则层的因素相对于目标层的层次单排序。

首先，对的每一列向量归一化，得到对按行求和得：

其次，将归一化得

最后，计算矩阵的最大特征根

其中表示向量的第个元素。

此外，还需要计算一致性指标：

，其中。

RI为平均随机一致性指标，当CR<

0.1，认为判断矩阵的一致性是可以接受的，否则应对判断矩阵作适当调整。

Ⅲ.进行层次总排序。

需要将计算出来的层次单排序再次进行适当计算。

假设层次结构模型是由目标层（A）、准则层（B）和方案层（C）所组成，准则层有6个因素，方案层有4个因素。

已知B层对A层的层次单排序为：

以为列向量构成矩阵，则C层对B层的准则Bk的层次单排序为：

同时，C层各方案对A层的层次总排序的方法为：

更一般地，若共有s层，则第k层对第1层（设只有1个因素）的组合权向量，即层次总排序满足其中

。

其中是以第k层对第k-1的权向量为列向量组成的矩阵，于是最下层（第s层）对最上层的层次总排序为

再根据上面的步骤进行逐步计算，最后可以得到2010年各主要食品受各种重要因素影响比重表为：

表32010年各主要食品受各种重要因素影响比重表

矩阵

权重向量

CI

RI

CR

A-B

6.1140

0.0228

1.24

0.0184

B1-P

4.0329

0.0110

0.90

0.0122

B2-P

4.0276

0.0092

0.0102

B3-P

4.0435

0.0145

0.0161

B4-P

4.0734

0.0245

0.0272

B5-P

4.0104

0.0035

0.0038

B6-P

4.0042

0.0014

0.0015

其中准则层相对于目标层的排序为

方案层相对于目标层的排序向量为

（其中权重结果的排序依次是微生物，添加剂，重金属，其它。

）

Ⅳ.组合一致性检验。

我们可以逐层进行组合一致性检验，若第p层的一致性的指标为（n是第p-1层因数的数目），随机一致性指标为。

定义

，，

则第p层的组合一致性比率为

同时，第p层通过组合一致性检验的条件为

从而最下层（第s层）对第1层的组合一致性比率为

当且仅当适当的小时，才认为整个层次的比较判断通过一致性检验。

我们对上述结果进行组合一致性收敛，得到第2层的组合一致性比率为

而在第3层中，有

从而有,均通过检验。

同时最下层（第3层）对第1层的组合一致性比率为

组合一致性检验通过，前面得到的组合权向量W可作为最终决策影响安全因素权重依据。

同理，按照上述的方法分析得出2011年各主要食品受各种重要因素影响比重表。

表42011年各主要食品受各种重要因素影响比重表

6.1367

0.0273

0.0244

4.0078

0.0026

0.0029

3.9320

-0.0227

-0.0252

4.0697

0.0232

0.0258