食品质量安全抽检数据分析模型

1、.食品质量安全抽检数据分析模型摘要本文通过对所给食品抽检数据进行分析，根据问题建立不同的模型， 利用相应数据处理和建模软件进行求解，评价了食品安全情况， 寻找相应规律， 得到改进的食品抽检模型。问题一 , 评价各主要食品领域的安全情况，分析所给数据，将食品主要分为八大类，将安全指标分为四大类。对数据进行统计、筛选，用多项式拟合得出各主要食品三年的安全情况的走势。利用改进的层次分析法的模糊综合评价模型（ AHPFCE）定量的得出各类食品的综合评价指标。 综合评价八类食品中冷食饮品类、粮食类、蔬果类、水产类、调味品类等食品领域的食品安全质量对比为：Q2012Q2010Q2011；肉类、甜点小吃的食

2、品安全质量对比为：Q2010Q2012Q2011；整合三年总的食品综合评价指数可以得出三年内食品的安全质量： Q2012Q2010Q2011。问题二，根据提供的数据，统计出按照季节、销售地、产地分别抽查的样本，用层次分析法给出影响食品安全系数的四个指标的权重值，构建食品安全系数的公式 RW * P ，用相关性分析分别给出季节、销售地、产地与食品安全系数的相关性：10.8866, 20.3839,30.4468，由此结果可知，季节与食品安全系数有很强的相关性， 销售地点与食品安全系数的相关性弱，产地与食品安全系数的相关性强。问题三，为了更科学更有效的反映食品质量情况，结合问题一与问题二得出的相关

3、结论， 在经费一定的条件下， 采用多层分层抽样的方法，减少食品质量稳定的样本的抽查次数， 相应增加食品质量不稳定的抽查次数，适当采用抽查的记忆性来抽查食品质量。关键字： 拟合层次分析相关性分析分层抽样.一、问题重述城市食品的来源越来越广泛， 人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。问题要求：1、如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2、从这些数据中能

4、否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3、能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的） ，例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？二、问题分析对于问题一， 由于调查中食物种类复杂， 为了方便问题的讨论， 我们将食品分为八大类，即冷食饮料类、粮食及其制品类、 肉及其制品类、水产类，蔬果类，蛋及乳制品类、 调味品类和甜点小吃类， 对数据进行筛选。 在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属、添加剂和其

5、他（包括食品变质、主成份不合格等） 。从附录中得到深圳市2010，,2011,2012 三年抽检样本数据，并用 EXCEL制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。 根据相应的数据制作出每个因素安全情况三年的变化趋势。再用改进的层次分析法的模糊综合评价模型 （AHP FCE），从定量的角度求出三年的八大类食品及总体的变化趋势。对于问题二，第一问中涉及的影响因素是食品是否安全的直接因素， 但这些因素又是由其它相关因素决定。 分析数据，我们将影响食品安全的这些因素归结为销售地（检测地点）、食品产地、抽查季节三个要素。用层次分析法给出影响食品安全系数的四个指标的权重值，构建食品安

6、全系数的公式RW* P，用相.关性分析分别给出季节、销售地、产地与食品安全系数的相关性。对于问题三， 要改进原有抽检方法， 找到一个新方法， 使其在成本相同的情况下，测得数据结果更准确， ，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据，得到不同因素的影响程度是不同的。采用多层分层抽样的方法， 减少稳定数据的测量量， 加大不稳定数据的测量量， 减少测量环节，加重必要环节。三、模型假设1、抽查是在随机下进行，不存在任何的人为干扰;2、假设影响食品安全性的因素能且仅能分为四大类，其他没被分类的因素对食品安全性所造成的影响忽略不计;3、食品可以分为八大类，每一种食品都有明确的分

7、类，即我们所列举出的食品必属于食品分类中的某一类型，不存在模糊分类的情况;4、假设食品生产、加工、运输各方面的公司和工厂在食品运营过程中均为正常运作，不存在某个环节脱节情况;5、检测过程中的系统误差在所难免，也难以定量计算，忽略检测过程中因仪器等原因而引起的系统误差有利于考虑问题；四、符号说明B比较判断矩阵RI随机一致性指标CI一致性指标CR一致性指标比率Z食品安全系数Si编号向量Wi层次分析单排序的权重向量W层次分析总排序的权重向量P四大安全指标的不合格率的最小最优值max特征值的最大值.五、模型建立、求解5.1 问题一5.1.1 如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量

8、等安全情况的变化趋势。5.1.2 根据对提供的资料的分析，本文将食品主要分为八大类：冷食饮料酒水、粮食及其制品 (含豆类 )、肉类及其制品、蔬果、水产、调味品、乳制品、甜点小吃，主要研究其中主要的三大类；安全指标分为四大类：微生物、添加剂、重金属、其他。通过具体的统计筛选得到如下数据（只提供部分）：冷食饮料酒水微生物重金属添加剂其他2010 上0.006593000.0021982010 下0.03030300.0079740.0271132011 上0.0071650.0020470.0122820.0030712011 下0.0063170.0034010.0068030.00340120

9、12 上0.00875700.00175102012 下0.0035310.0014120.0035310.002825根据已有数据， 用最小二乘法进行多项式拟合。经多次实验发现四次多项式较为贴切的反映数据的趋势。 应用 matlab 进行求解作图（部分图，其他见附录）如下 (具体程序见附录 )：绘制粮食及制品数据图形.通过定性的分析可得到：（ 1） 粮食类制品其重金属的不合格率、添加剂的不合格率的波动较大。2010、2012 年的这两项指标的不合格率较低，而 2011 年的这两项指标的不合格率明显较高。（ 2） 粮食类制品其微生物的不合格率有明显的下降趋势。由图可以得出一些危害物的定性的变化

10、趋势，但不能显示出三年中粮食类食品的总的变化趋势，也不能定量的给出食品安全的综合评价标准。下面我们用改进的层次分析法的模糊综合评价模型，定量的对三年的各类食品的安全情况进行评价分析；5.1.3基于改进的层次分析法的模糊综合评价模型（AHPFCE）5.1.3.1 AHP FCE模型建立（1）根据所研究评价系统的实际情况，从代表性、系统性和适用性等的角度，建立模糊综合评价的评价指标体系，由于个评价指标的样本数据建立单评价指标.的相对隶属度的模糊评价矩阵。设有n 个 评 价 指标 组成 对全 体m 个方 案的 评 价指 标样 本集 数 据x(i , j ) i1 : n, j1 : m ，各项指标

11、x(i , j ) 均为非负值。为确定单个评价指标的相对隶属度的模糊评价矩阵，消除各指标的量纲效应， 使建模具有通用性， 需对样本数据集x(i , j ) 进行标准化处理。 为了尽可能保持各评价指标值的变化信息，对越大约优型指标的标准化处理可取为r (i , j ) x(i , j ) / xmax (i ) xmin (i)(1)对于越小约优型指标的标准化处理公式可取为r (i , j ) xmax (i ) xmin (i ) x(i, j ) / xmax (i)xmin (i) (2)对于越中越优型指标的标准化处理公式可取为x(i, j ) / xmid (i)xmin (i)， xm

12、in (i )x(i, j ) xmid (i )r (i , j )x(i, j ) / xmax (i)xmid (i ), x(i) midx(i, j )x(i) max xmax (i) xmid (i )(3)式中： xmin (i) 、 x(i ) max 、 xmid (i) 分别为方案中第 i 个指标的最小值、最大值和中间最适值； r (i, j ) 为标准化后的评价指标值，也就是第 j 个方案第 i 个评价指标从属于优的相对隶属度值，i =1n, j = 1 m 。以这些 r (i , j ) 值为元素可组成单评价指标的模糊评价矩阵R = (r (i, j)nm 。（2）根

13、据模糊评价矩阵 R = ( r (i, j )nm 构造用于确定各评价指标权重的判断矩阵 B = (bij )nm 。模糊综合评价的实质是一种优选过程， 从综合评价的角度看，若评价指标 i1的样本系列 r (i2, j ) j = 1 m 的变化程度比评价指标 i 2的样本系列 r (i 2, j ) j = 1 m 的变化程度大，则评价指标 i1传递的综合评价信息比评价指标 i2 传递的综合评价信息多。基于此，可用各评价指标的样本.m0.5轾2s(i)=犏?%标准差犏反映各评价指标对综合评价指标的影响臌j= 1m程度，并用于构造判断矩阵 B 。其中 r%=? r (i, j ) / m 为个

14、评价指标下样本系j = 1列的均值， i = 1 n。于是可得判断矩阵s(i )s( j ) b1,s(i) s( j)smaxsminm(4)bijs( j ) b1/ s(i )1 , s(i) s( j )smaxmsmin式中： smin 、 smax 分别为 s(i) i = 1 n 的最大值和最小值； 相对重要性程度参数值 bm = smax / smin 。（3）判断矩阵 B 的一致性检验及其权重 w (i = 1 n) 的计算。i1、判断矩阵的特征值和特征向量的求解（ 1）先求判断矩阵 B 中每行元素之积 M in, 有 Mia(i 1,2,L, n)ijj 1（ 2）再求 M

15、 i 的 n 次方根 , 得 Win M i (i1,2,L,n)(3)然 后对向量W(W1 ,W2 ,L ,Wn )T进 行 归一化 ,得nLTWW /W12n.(i 1,2, , n), 从而得到特征向量iiiW(w ,w Lw )i 1（ 4）计算判断矩阵的最大特征根（值）n ( BW )imaxi1nWi2、进行一致性检验的步骤如下(1) 计算CImaxnn1.nRI.(2) 计算根据 RI(random index)随机一致性指标 , 可查表确定，如表 1所示。表1随机一致性指标 RI 值1、 2345678910111200.580.91.121.261.361.411.461.4

16、91.521.54CRCIRI若判断矩阵不满足一致性条件（CR0.1 ），则需要修改 （ 4）把各评价指标的权重值 w 与各方案相应评价指标的相对隶属度r (i , j ) 相乘i并累加，可得到模糊评价的综合指标值z( j )mz( j )= ? wi r (i , j ),( j = 1 m)（ 8）i= 1综合指标值 z( j ) 越大说明第 j个方案越优，据此可进行科学决策。5.1.3.2 对于问题一的求解对于该题以粮食及其制品类为例给出食品安全评价体系。粮食类食品安全系数微生物重金属添加剂其他2010 年2011 年2012 年粮食类食品安全评价体系（ 1） 以四类危害物三年的不合格率

17、作为样本数据，构造样本矩阵X ( x(i, j )43 :.轾0.0007350.0161650.01396犏犏0.0073820.0359890.005878X=犏犏0.0033840.0013250.02307犏犏0.0141340.0039750.0053臌再利用（ 2）式r (i , j) xmax (i ) xmin (i )x(i, j ) / xmax (i ) xmin (i) 对数据进行标准化处理得到模糊评价矩阵R =(r (i , j )43 ：轾0.51700.91330.0867犏犏0.02000.85570.9800R=犏犏0.37990.62010.5655犏犏0.

18、63470.57080.3653臌m0.5轾)2/ m（ 2） 求出各评价指标的样本标准差s(i)= 犏? (r (i , j) - ri得：犏%臌j= 1S = 0.3376 0.42630.10280.1149(3) 根据（ 4）式可得判断矩阵 B(b(i, j )44 :轾1.00000.46794.00943.8542犏犏2.13721.00005.14664.9914B = 犏犏0.24940.19431.00000.8657犏犏0.25950.20031.15521.0000臌.(4) 判断矩阵 B 的一致性检验及其权重 wi (i = 1 n) 的计算。由 MATLAB可求出判断

19、矩阵的最大特征值为：max4.0336对应的特征向量为： W = 0.3089 0.5138 0.0847 0.0926T所以：CI =max -n4.0336-40.0112n - 1=4 - 1=取RI = 0.9,则：CR= CI =0.0112 = 0.0124 0.1RI0.9判断矩阵符合一致性检验。（ 5）把各评价指标的权重值wi 与各方案相应评价指标的相对隶属度r (i , j ) 相乘并累加，可得到模糊评价的综合指标值z( j )mz( j )= ? wi r (i , j ),( j = 1 m)i = 1即：Z = 0.61200.26100.8271同理可求出各类食品的综

20、合评价指数和总食品的评价指数，如图表2 所示。图表 2 各类食品的综合评价指数年份冷食饮粮食类肉类蔬果类水产类调味品蛋乳制甜点小品品吃20100.55020.61200.88150.37330.53460.58450.54930.906120110.39500.26100.33850.53820.48660.22420.73300.149620120.75450.82710.61200.92140.68870.62640.95040.2539总的食品综合评价指数：.年份201020112012综合评价指数0.67730.27050.7284Qi 代表三年中食品的安全质量，i =20102011

21、2012；结果分析：(1) 八类食品中冷食饮品类、粮食类、蔬果类、水产类、调味品类等食品领域的食品安全质量对比为： Q2012 Q2010 Q2011；(2) 八类食品中肉类、甜点小吃类食品的食品安全质量对比为：Q2010Q2012 Q2011；(3) 整合三年总的食品综合评价指数可以得出三年内食品的安全质量：Q2012Q2010 Q20115.2 问题二5.2.1 层次分析 +相关性分析模型在进行食品安全抽检之前， 为了提高抽检工作的效率以及准确性，首先分析并找出影响食品安全的主要因素。 本文首先建立层次分析模型的主要食品添加剂、重金属、微生物等因素的权重( W ) ，找出关键性因素，这有利

22、于把握抽检工作的关键，使抽检结果更科学化。5.2.1.1理论基础由于食品种类的不同 , 影响食品安全因素也很多，因此首先确定影响食品安全问题的主要因素 , 在通过层次分析法 （AnalyticHierarchyProcess, 简称 AHP）分析影响食品安全各个因素的相对重要性。这一方法的特点是对复杂问题的本质、影响因素以及内在关系等进行深入分析之后，构建一个层次结构模型递阶层次结构，然后利用有限的定性与定量信息，把决策的思维过程数学化。因此，应用 AHP方法可以科学地、系统地确立各决策权重。5.2.1.2确立递阶层次结构在深人分析实际问题的基础上， 将有关的各个因素按照不同属性自上而下地分解

23、成若干层次。 同一层的诸因素从属于上一层的因素，同时支配下一层的因素或受到下层因素的作用， 而同一层的各个因素之间尽量相互独立。最上层为目标.层，最下层通常为方案层，中间为准则层。5.2.1.3构造判别矩阵从层次结构模型的第二层开始，对同一层的各个因素， 用成对比较法和比较尺度构造成对比较矩阵， 直到最下层。 以上一层某个因素作为比较准则时，可用标度“表示对下一层中第i 元素与第j 元素的相对重要性，由构成的矩阵称A (aij ) 为比较判别矩阵。比较判别矩阵 A 是正互反矩阵。5.2.1.4计算权向量并做一致性检验计算比较判别矩阵 A 的最大特征值所对应的特征向量并归一化， 就是所求的权重向

24、量。对成对比较进行一致性检验， 在一致性检验。 在一致性检验中需要引出 Satty 的随机一致性指标 RI ; 对于不同的 n ， RI 的数值不同。当时 n 1,2 ，RI 0 ，是因为 1，2 阶的正互反矩阵总是一致阵。对于n 3 的成对比较矩阵 A,将它的一致性指标 CImaxn 与同阶的随机一致性指标CR, 当 CRCI0.1n1RI时认为 A 的不一致程度在容许的范围之内，可用其特征向量作为权向量。5.2.1.4食品安全系数 R通过所给数据，可以构建食品安全系数RW * P ，其中 P 为各抽样样本中微生物、重金属、添加剂、其他不合格率的最小最优值。最小最优公式为：Pimax( P)

25、min( P) Pi ，其中 P 是各样本中微生物、重金属、添加剂、其他的max( P)min( P)概率向量。5.2.2 用层次分析模型分析影响各类食品安全因素权重5.2.2.1建立方案评价的递阶层次结构模型.主要食品安全指标权重冷食饮品粮食类肉类微生调味品甜点小 .微生物重金属添加物其他安全指标权重层次结构图该模型的最高层为总目标层：A确定影响各主要食品安全指标权重.第二层为方案评价的准则层，分别为：B1 ：冷食饮品； B2 ：粮食及制品；B3 ：肉及制品； B4 ：蔬果； B5 ：水产品；B6 ：调味剂； B7 ：乳制品； B8 ：甜点小吃最底层为方案层，分别为：P：微生物含量； P：重

26、金属含量； P ：添加剂； P ：其他12345.2.2.2构造比较判断矩阵利用 1 9标度法 ( 见下表 ) ，以A为比较准则， B层次各因素的两两比较判断矩阵为 AB，如表下表所示。 同理，以 每一个 Bi 为比较矩阵， P 层次各因素的两两比较判断矩阵为 Bi P(i 1,2,.8) ，进行成对比较 , 同时参考专家意见 , 确定各因素之间的相对重要性并赋以相应的分值, 构造出各层次中的所有判断矩阵,并计算权向量和一致性检验。表三 1 9 标度的意义aij135792、4、6、 1 ,1/2 ,.8?,1/9aij 的物ci 和cjci比 cjci比c jci比cjci比c jci比c

27、jci比c j 的影响理意义影响相影响影响影响影响绝为上述之比与上述说同稍强强明显对地强判断的明相反地强中间级A- B之间构成的判断矩阵 :113264572214375681111423534211215346A B231111111267453211113124452311112113563421111111178562435.2.2.3层次单排序以及其一致性检验对于各比较矩阵， 求出其最大特征值及其对应的特征向量， 将特征向量经过归一化后，即可得到相应的层次单排序的相对重要性权重向量， 以及一致性指标 CI 和一致性比例 CR, 具体 MATLAB程序详见附录，其中特征值最大值n(BW)

28、 i；结果如下表所示maxnwii 1矩阵层次单排序的权重向量 WimaxCIRICRA-B(0.2304,0.3307,0.1059,0.1572,0.0336,0.04541.410.0710,0.0477,0.0235)8.31800.0322B1P(0.4673,0.0954,0.2772,0.1601)4.03100.01030.900.0114B2P(0.2772,0.1601,0.4673,0.0954)4.03100.01030.900.0114.B3P(0.4673,0.1601,0.2772,0.0954)4.03100.01030.900.0114B4P(0.1601,0

29、.0954,0.2772,0.4673)4.03100.01030.900.0114B5P(0.2772,0.1601,0.4673,0.0954)4.03100.01030.900.0114B6P(0.1601,0.0954,0.2772,0.4673)4.03100.01030.900.0114B7P(0.0954,0.4673,0.1601,0.2772)4.03100.01030.900.0114B8P(0.2772,0.0954,0.4673,0.1601)4.03100.01030.900.01145.2.2.4计算层次总排序权重向量并做一致性检验由上表可得 B 层对 A 层的权重

30、向量为：W 2(0.2304,0.3307,0.1059,0.1572,0.0336,0.0710,0.0477,0.0235)P 层对 A 层的组合一致性指标为：CI 3(0.0103,0.0103,0.0103,0.0103,0.0103,0.0103,0.0103,0.0103) gW 20.0103组合随机一致性指标为：RI 3(0.90,0.90,0.90,0.90,0.90,0.90,0.90,0.90) gW 20.90组合一致性比率指标为：CR3CR2CI 30.0436RI 3可知 CR30.10 ，所以整个层次的判断矩阵通过一致性检验，有组合权重向量为W( 0.3057,0

31、.1436,0.3453,0.2054)即微生物、重金属、添加剂、其他所占的权重分别为：0.3057 、0.1436、0.34530.2054，5.2.2.5计算各样本的食品安全系数通过数据的整理分析，得到有关季节、销售地点、产地的样本数据，按照R的公式，用 MATLAB编程（程序见附录）得到计算结果为：四个季节的食品安全系数为：R1( 0.3247 ，0.4154，0.5997，0.5513).销售地点的食品安全系数为：R2 = （0.5279 ，0.7297，0.5119，0.4436 ，0.5413 ）产地的食品安全系数为：R3 =()5.2.4 对季节、销售地点、产地与食品安全系数进行

32、相关性分析要进行相关性分析，本文首先对季节、销售地、产地编号量化，规则如下表：季节第一季度第二季度第三季度第四季度编号 S11234销售地点商场及各类各类食堂及加油站便利餐饮类场所交通枢纽便专营店学校周边小店利店店编号 S212345产东莞佛山福建广东湖南山东上海深圳四川浙江重庆地S31234567891011对各编号 Si 进行归一化处理，公式为 hiSi, i 1,2,3，其中 hi 为归一max( Si )化后的数值，由于篇幅限制，结果不在详细贴出。在各类型的相关分析中， 只有两个变量的线性相关关系的分析是最简单的。两个变量之间线性相关程度可以用简单线性相关系数去度量，这种相关系数是最常

33、用的，也简称为相关系数。 两个相互联系的变量的相关系数称为总体相关系数，通常用表示，计算公式如下：Cov(Si , Ri )iVar (Si )Var ( Ri ).其中 Var ( Si ) 表示各样本编号的方差，Var ( Ri ) 表示各样本食品安全系数的方差，Cov (Si , Ri ) 表示二者的协方差。通过 MATLAB编程（见附录），分别求出各样本的相关系数为：10.8866, 20.3839,30.44685.2.5 模型结果分析本文设相关系数的等级比较，认为时0.85 ，两者具有很强的相关性，当0.4 时，认为相关性弱。根据此标准，观察本模型的结果可知，季节与食品安全系数有很

34、强的相关性， 销售地点与食品安全系数的相关性弱，产地与食品安全系数的相关性强。5.3 问题三5.3.1 分层抽样方法分层抽样也称类型抽样或者分类抽样，通俗地说就是异质性较强的总体分成若干个同质性较强的子总体，再从不同的子总体中抽取样本分别代表该子总体，进而将所有各总体中抽取的样本合并成一个样本代表总体。它的优点是每一层上都可以采用不同的抽样方法，并且提高了参数估计的精度。即先将包含 N 个单元的总体分成各包含 N1, N 2,. N L 个单元的子总体。这些子总体互不重复合起来就是整个的总体， 因此有 N1 N2 . N L N 。这些子总体就成为层。层被确定后，就从每一层抽一个样本，抽样是在

35、各层独立地进行的，各层内的样本含量分别用n1, n2 ,.nl 表示这就是分层抽样5.3.2 食品质量安全的分层抽样设在某个时间短，对主要食品抽检的批次为N，每个批次的抽检次数相同，每次的抽检费用相同，抽检时间相同。本文按照食品类型进行分层， 分为八层：冷食饮料酒水、粮食及其制品 ( 含豆类 ) 、肉类及其制品、蔬果、水产、调味品、乳制品、甜点小吃，每层的权重有第二问的层次分析模型给出，分别为：（ 0.2304,0.3307,0.1059,0.1572,0.0336,0.0710,0.0477,0.0235），则每层的抽检批次为.N*（0.2304,0.3307,0.1059,0.1572,0

36、.0336,0.0710,0.0477,0.0235）。对每层食品按照季节、产地、四大食品安全指标的比例权重分批进行抽样检验。抽检工作并不能就此结束， 我们应该从长远考虑， 对某个品牌进行多次跟踪抽检，成本费及时间比一次抽检明显是增大。 在这里将采用放宽加严策略， 即每次抽检都具有记忆性， 前几次抽检结果都将对以后抽检抽取样本数有影响。 如果前一次的检测效果比较好， 后一次检测就采用放宽政策， 减少产品的抽检批次和抽检项目，这样就使得成本及检测时间大大减小； 如果前一次的检测效果并不理想，则后一次就采用加严政策， 增加下一次的抽检力度， 厂商就会重视产品的质量，等下一次再抽检时，产品的抽检效果

37、会提高，当达到一定要求后，就又可进行放宽政策， 减少产品的抽检批次及抽检项目， 这样，产品的抽检成本及时问也会大大减小。六、模型的优缺点及改进方向优点：本文对于各种因素对于食品质量安全的影响建立了改进的层次分析法的模糊综合评价模型、 层次分析模型和食品安全系数， 全面综合考虑了各个方面的因素，避免了单一因素分析的不准确性， 得出了合理的数学模型。 并且通过各因素的显著性分析， 找到了影响食品质量安全的主要因素， 较符合实际情况， 模型可靠，并且模型相对简单，利于操作；该方法不仅适用于本题，也适用于其他方面的数据预测，有实际背景，可运用于实践，具有广泛适用性。缺点：本文在做第一个问题的时候， 数

38、据处理存在漏洞， 对结果的定量分析有一定的误差 ; 在做第二问的时候食品安全系数问题的时候没有食品产地、食品加工环节、食品抽检环节因素对食品安全的权重， 并且默认了各个环节都符合标准，但实际上各个环节有待检验， 第三问抽检的数量有限， 没做具体的数据处理。改进方向：基于问题一的改进的层次分析法的模糊综合评价模型、 问题二的食品安全度模型以及问题三的分层抽样模型的缺点以及优点， 需要对对抽检模型进行修正，通过 MATLAB 软件进行仿真，做出最佳抽检测策略模型有利于对问题的定性分析和定量分析。六、参考文献1 马新民，王逸迅 . 概率与数理统计 M. 北京：机械工业出版社， 2010.2 刘卫国

39、. MATLAB程序设计与应用 M. 北京：高等教育出版社， 2009.3 茆诗松 . 回归分析及其试验设计 M. 上海：华东师范大学出版社，1986.4 姜启源，数学模型 M, 北京，高等教育出纳版社， 20065 刘荣珍 , 赵军 , 模糊评价模型在长江水质评价中的应用 , 兰州交通大学学报 ,2007.12附录问题一.改进的层次分析法的模糊综合评价模型（AHPFCE）求解程序：.'ht','sheet1'表格截图function y=zuixiaozuiyou(x)a=xlsread('ht','sheet1','J

40、83:J85');b=xlsread('ht','sheet1','K83:K85');c=xlsread('ht','sheet1','L83:L85');d=xlsread('ht','sheet1','M83:M85');R=a' b' c' d'RMax=max(R')'RMin=min(R')'for i=1:4for j=1:3R(i,j)=(RMax(i)+RMin(i)-R(i,j)/(RMax(i)+RMin(i); endendRfor i=1:4g=0;sp=sum(R(i,:)/3;for j=1:3k(j)=(R(i,j)-sp)2;g=g+k(j);endS(i)=(g/3)0.5;.endSsmin=min(S);smax=max(S)