系矩阵R。其中，农药和化肥的隶属度由图2、图3所示的函数确定，其余指标的隶属度由其测度数据直接表示。设各同级指标关于上一级指标的权重向量为W，则上一级指标的安全指数S=WR。
以此类推，通过搜集历年各基础指标的测度数据，可以逐级向上计算出各个 指标的安全指数，从而得出每年的食品安全指数，以此来定量评估每年的食品安全状态。
5.4 评估往年食品安全状态
我们从各大网页搜集到2005-2012年各个相应指标的测度值，列成表格如下：（表6） 2005 2006 2007 2008 2009 2010 2011 2012 农产品94.2 农（兽）药残留监测合格率a11（%） 农药使103.9 用量a121（万t） 93.0 95.1 97.0 94.5. 96.9 97.1 97.9 129.6 173.1 190.2 226.2 234.2 264.8 278.9 化肥使3414.8 3782.3 用量a122（万t） 89.1 加工企87.6 业卫生监督合格率a211（%） 97.3 持健康97.0 证明人数比例a212（%） 加工业食添加抽检格
4187.1 4263.2 4863.7 4521.1 4179.0 5800.0 91.2 92.1 89.9 90.1 89.3 92.2 97.6 98.2 98.0 97.2 97.6 97.5 企87.5 品剂合率90.3 93.3 94.3 94.5 93.9 94.3 94.0 9
a221（%） 饮用水89.4 的抽检合格率a222（%） 加工业食品抽样合格率a23（%） 经营/餐饮业环境卫生合格率a31（%） 餐具消毒合格率a32（%） 经营/ 餐饮业持健康证明人数比例a33（%） 经营/ 餐饮业食品抽样合格率a34（%） 89.4 81.7 85.7 88.1 86.5 90.5 91.8 96.7 90.4 91.2 92.1 91.3 91.6 91.8 93.3 88.8 87.3 91.1 93.8 92.2 93.3 93.0 92.8 79.7 82.5 83.0 84.6 83.9 84.5 83.5 84.0 89.8 96.9 97.9 98.9 98.0 97.5 98.1 97.7 86.91 91 86.9 91.2 93.3 94.9 94.3 96.3 10
80.3 87.6 81.6 83.7 85.0 91.2 91.5 加工工75.0 具设备（包装材料，容器，洗涤剂，消毒剂）的合格率a24(%) 根据图2和图3以及表6的统计数据，可计算出模糊关系矩阵R和相应的安全指数S，逐级评价，最后可得出总目标层的综合安全指数S（a）。
下面以S(a12)为例说明求解过程。
R12是表示三级指标A121农药施用风险（农药施用量）与a122化肥施用风险（化肥施用量）和它们上层的二级指标a12种植过程风险之间模糊关系的矩阵。根据图2和图3农药和化肥施用风险隶属度函数，可以读出2005-2012年每年农药和化肥施用风险的隶属度，由隶属度构成模糊关系矩阵R12。
R(12)=
?0.6537 0.5680 0.4230 0.3660 0.2460 0.2193 0.1173 0.4037 ???0.4112 0.3479 0.2781 0.2650 0.1614 0.2205 0.2795 0.0000?? ??则相应的安全指数为： S（a12）=wR12
?0.6537 0.5680 0.4230 0.3660 0.2460 0.2193 0.1173 0.4037 ????=(0.667 0.333)?0.4112 0.3479 0.2781 0.2650 0.1614 0.2205 0.2795 0.0000??
=(0.5729 0.4947 0.3747 0.3124 0.2178 0.2197 0.1713 0.2693)
同理，可得： S（a1)
?0.9420 0.9300 0.9500 0.9700 0.9450 0.9690 0.9710 0.9790?=(0.667 0.333)??0.5729 0.4947 0.3747 0.3124 0.2178 0.2197 0.1713 0.2693??
??=(0.8191 0.7850 0.7584 0.7510 0.7082 0.7195 0.7047 0.7427)
运用同样的方法，可算得如下：
S(a21)=(0.9073 0.9183 0.9333 0.9413 0.9260 0.9246 0.9206 0.9396)
S(a22)=0.8845 0.8600 0.8950 0.9120 0.9050 0.9220 0.9305 0.9535）
S(a2)=( 0.8557 0.8662 0.9009 0.8942 0.8929 0.9021 0.9208 0.9359)
S(a3) =(0.8548 0.8848 0.8843 0.9095 0.9083 0.9164 0.9118 0.9181)
11
S(a) =(0.8447 0.8451 0.8559 0.8542 0.8390 0.8502 0.8552 0.8756) 最后，按照表6的对应关系确定食品安全风险等级。
表7 食品安全指数与食品安全风险对应关系 风险等级 高风险 中等风险 低风险 无风险 0.8~0.9 0.9~0.95 >0.95 对应安全指数 <0.8
结果分析
根据最后结果可得出如下结论： 1）、S（a）的值稳定在中等风险范围，从食品供应链整体范围看，我国食品安全状态处于中等风险范围，整体水平还有待提高。 2）、食品安全综合指数估算值逐年增加，说明我国食品安全状态呈现缓缓提高趋势， 这与我国食品安全的实际情况相符。
3)、从结果看，S(a12)的值很低，风险很大，要限制农药化肥的施用才行。S（a1）在高风险上下波动，说明风险值比较大，要调节兽药和激素的使用，降低风险。 4）、S(a21) 、S(a22)、S(a2)的值较大。属于低风险范畴，说明我国还是比较注重食品加工安全，继续保持。 5）、S（a3)表示的是经营消费环节的安全风险，其值在中等风险和低风险上下波动，有待提升。 6）、依据在S（a）中取最大的原理，我国现在的安全指数为0.8756。
5.4 利用回归分析和移动平均法预测未来食品安全状态
由上面建立的模型，我们得出了粗略的各年的食品安全指数，历史是用来借鉴并可以用来推测未来趋势发展，要想给人们达到一个预警，我们采取了以下方法来预测短期性的未来，首先我们利用线性回归画出S（a）s随时间推移而变化的函数曲线，如下：（相关程序设计在附录给出）
12
图4
从图形可以看出安全指数波动不大又逐年上升的。于是我们采用了移动平均法来预测估计接下来几年的安全指数值。
移动平均法 ? 基本思想
— 假定在一个较短的时间间隔里，序列值之间的差异主要是由随机波动造成的。根据这种假定，我们可以用一定时间间隔内的平均值作为某一期的估计值。 ? 分类
— n期中心移动平均 — n期移动平均
n期中心移动平均
????????t??????????~x1(xt?n?1?xt?n?1?1???xt???xt?n?1?1?xt?n?1)，n为基数2222n1(1xt?n?xt?n?1???xt???xt?n?1?xt?n)，n为偶数2222n213