2013年深圳杯C题优秀答案

1、食品质量安全抽检数据分析摘要本文根据文章提出的不同问题，建立相应的数学模型，利用 matlab软件进行求解，对食品质量进行评价和找规律以及合理抽检方 法。针对问题一，对深圳市这三年各主要食品领域微生物、重金属、 添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对数据 处理按季节分为17个子样本点的抽样值进行分析，得到各子样本食 品安全情况。运用层次分析法决定食品安全单位指标在综合评价中权 重。可以得到深圳市近三年来的食品质量情况应是明显提高。年份201120122013食品安全系数0.270.170.095等级III级II级I级针对问题二，我们先通过MATLAB对原始数据进行检验，对残差

2、 向量进行分析，得到了残差向量分析图，剔除其中的异常点。运用 MATLAB进行编程，得到各因素的偏回方和：xix1x2x3x4x5x6x7x8Q ( x 106)0.33840.00300.36850.10090.14291.24180.14720.1963根据Q的大小可判断各因素对食品安全系数的影响程度 ix x x x x x x x63187542结论：在食品质量影响因素中食品产地影响最大，食品加工次之， 季节影响最小，抽查地点几乎无影响。针对问题三，根据问题一、二数据结果以定量比较评估的方法 分析了各类影响食品安全的因素及其可能造成的危害性的问题。改进 后的食品抽检的办法以主要食品为准

3、则层建立了层次分析（AHP）模 型，对影响食品安检的危害性因素做出定量分析如问题二个影响因素 大小。由问题一、二结果可以通过建立抽检模型，即改进后的规准型 抽样检验模型，并以蒙特卡罗法对抽检的全过程进行模拟，得到相对 误差逐渐趋向于0。关键词：方差与回归分析 残差向量分析 评价指数 层次分析问题重述“民以食为天”，食品安全关系到千家万户的生活与健康。随着人们对生活 质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民 生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比 例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、 销售以及餐饮等每一个环

4、节皆可能影响食品的质量与安全。另一方面，食品质 量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要 科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。 请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、 鸡鸭等抽检数据的获取），并根据这些资料来讨论：如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量 等安全情况的变化趋势；从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的 关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且 不过分增加监管成

5、本（食品抽检是需要费用的），例如对于抽检结果稳 定且抽检频次过高的食品领域该作怎样的调整？二、问题分析本问题是在监督最统一，最规范，最公开的城市之一深圳，通过 对食品安全抽查数据，抽样，合格，不合格以及与近几年数据比较的 数据比例进行研究，探讨使人们追求的食品质量安全与可持续科学发 展观加快吻合的方案。对于市场食品安全的分析，我们知道，和y有关的变量有8个， 研究y与变量x，x，x，*之间的定量关系的问题为多元回归问1 2 3 8题。又因为许多多元非线性回归问题都可以化为多元线性回归问题， 所以对于本问题我们建立了多元线性回归的数学模型。第一问：由于食品安全生成这一过程涉及到诸多环节，所以我们

6、必须考虑 食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮 等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与 安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都 需要科学有效的方法。运用统计分析中的方差，分析各个影响因子对 食品安全的影响程度。同时我们确定出哪些因素是主要的，哪些是次 要的。通过表格可以看出，由于抽样的各种不同方法以及子样空间的 不同导致的子样均差以及子样方差的不同，导致了食品安全系数的标 准不同。首先，考虑抽样方法的影响。通过对抽样数据和合格数据的观察， 我们发现在不同抽样方法以及不同子样本空间下，对食品安全系数有 较大的影响。通过对整个表格

7、的分析我们把抽样数据分了四个数据梯 度,，分别为2013年的，2012年的，以及2011年和2010年的，然后 列出四组数据，再应用方差分析模型进行分析，判断出评价深圳市这 三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化 趋势。其次，我们在来考虑诸多环节对食品安全系数的影响。我们把数据分 为四个组进行分析，分别为取材环节，生产环节，流通环节，餐饮环 节。然后我们同样应用方差分析模型来对数据进行处理，来判断各个 环节对评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势。第二问：首先对附表1的数据进行处理，对残差向量进行，剔除其中的异 常点。然后我们建立了多

8、元线性回归的数学模型，并采用了最小二乘 法来估计参数。把模型写成矩阵的形式，化简整理得其正规方程组， 通过对正规方程组的求解，最后得到回归方程。对于各因素对食品安全的影响程度，由于利用偏回归平方和Q可i 以衡量每个变量在回归中所起的作用大小（即影响程度），我们对每 个变量气的偏回归平方和Q进行了计算，最后把影响程度由深到浅的 各因素进行了排序。回归方程的显著性检验：事先我们并不能断定随机变量y与一般 变量x，x，x，x之间是否确有线性关系。在求线性回归方程前1238线性回归模型只是一种假设，所以在求出线性回归方程之后，我们需 要对其进行统计检验。将总的平方和分解为回归平方和s回和剩余 平方和S

9、剩，运用F检验法来判定线性回归方程的显著性。回回归系数的显著性检验：由于回归方程显著并不意味着每个自变 量x，x，x，x对因变量y的影响都是重要的。而我们要找出影1238响食品安全的主要因素，即从回归方程中剔除那些次要的、可有可无 的变量，这就需要我们对每个变量进行考察。显然，如果某个变量对 y的作用不显著，那么在多元线性回归模型中，它前面的系数P .就可 以取值为零。因此，检验因子x.是否显著等价于检验7假设 H： =0。最后再运用T检验法来辨别模型中哪些因子是显著的 及从这些数据中能否找出某些规律性的东西：如食品产地与食品质量 的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素

10、等等。第三问：由于回归系数之间存在相关性，当从原回归方程中剔除一个变量 时，其他变量，特别是与它密切相关的一些变量的回归系数就会受到 影响，剔除一个变量后，这个变量对y的影响很大部分转加到另一个 变量对y的影响上。所以，我们对回归系数进行一次检验后，只能剔 除所有不显著因子中，值最小的，然后重新建立新的回归方程，再对 新的回归系数逐个进行检验，直到余下的回归系数都显著为止。三、问题假设1、土壤对取材生长状况基本相同，从而对生成食品的安全系数 的影响相同。2、假设周围自然环境因素除去土壤外，其他的对食品安全系数 的生成没有影响。3、在食品生产过程中，材料都已经转化成了食品。也就是说产生是完全的，

11、没有发生材料的剩余。4、同类食品是相同的，也就是说抽查的食品对生成食品的安全系数影响相同。5、餐饮服务等条件始终是不变的。6、抽查是在随机下进行，不存在任何的人为干扰。四、参数及符号说明表符号表示的意义单位备注Li第i个子样点与样本空间的差值个i = 1.7vi第i个空间的大小个N ij第i个子样本第j种食品的合格%N .：第i个样本空间第j子样本经管理后下一 年的合格率%w.第i个样本空间中第j子样本总量个V ij第i个子样本第j个食品合格率%T管理系数1/月0.1 t 0.5ti (i+1)样本空间与子样本之间差距个PI食品安全质量综合评价指数1Ai准则层项目因子i = 1.4Bi项目层项

12、目因子i = 1.4C.I.一致性指标R.I.平均随机一次性指标C.R.一致性比例指标Ri为第i年需要处理的不合格食品个i=1.10Hi为预测的i年的食品不合格数个i=1.10表二符号 说明X影响食品质量安全的因素。(i 1,2,.,8)y O各食品安全系数(a -1,2,.,38)_i相互独立且服从同一正态分布的随机变量(，-1,2,.,38)Q i变量X,的偏回归平方和S 总总平方和S-回回归平方和S 剩剩余平方和Pi待估计系数参数(i =1,2,.,8)yy的回归值b_i参数。的最小二乘估计(i-1,2,.,8)食材生长土壤，食材的生产收获，食品的运输、加工、包装、贮存、 销售以及餐饮八

13、个影响因素五、模型的分析、建立与求解设K取尸个水平，瓦,在水平H下总体x,服从正态分布/ = r ,这里未知，,4可以互不相同，但假定x,有相同的方差。又设在每个水平H下作了 0次独立试验，即从中抽取容量为勺的样本，记作X ，/=, X服从q2),，= 1,尸,/ = 1,0且相互独立。将这些数据列成表1 (单因素试 验数据表)的形式。表1单因素试验数据表4邑112砒a2%X22f4命&将第，行称为第，组数据。判断刀的尸个水平对指标有无显著影响，相当于要作以 下的假设检验H。:巧=/弓=“,.；: pH，,.不全相等由于的取值既受不同水平其的影响，乂受其固定下随机因素的影响，所以将它 分解为

14、rr夕= 一a，h =，a,=出一i = l,/即=1=1产是总均值，名是水平H对指标的效应由（1）、（2）模型可表为 孔=片+名+与 7=1%N（O,cH）, z = 1,y = 1,.,勺 原假设为（以后略去各选假设）H0 :% = OL之=叫_ =。 1.2统计分析 记 、1 ! E，=二%，元=二/与 勺 J=1n j =l J=1a是第i组数据的组平均值，元是总平均值。考察全体数据对元的偏差平方和=1 7=1 经分解可徂 记，=支勺（一元）27=1 3=1 J=1 则 S S（9） S是各组均值对总方差的偏差平方和，称为组间平方和；3反是各组内的数据对均值偏 差平方和的总和。3反映4

15、不同水平间的差异，3丘则表示在同一水平下随机误差的 大小。注意到文（孔-xj2是总体NgS）的样本方差的呵-1倍，于是有 7=1工（-打侦zW-D 7=1 由;分布的可加性知 TOC o 1-5 h z f r、I ：=i7即s/b2 /(I)且有毛=(77-*(10)对s作进一步分析可得ESa = (r-lK + 叫;(11)J=1当H成立时=(1顽(12)可知若反。成立，5只反映随机波动，而若H。不成立，那它就还反映了 乂的不同水平的效应。单从数值上看，当H。成立时，由(10)、(12)对于一次试验应有SE/(n-r)而当H。不成立时这个比值将远大于1。当丑。成立时，该比值服从白由度昭=-

16、1,n2 = (w- r)的F分布,即F =尸(井_1 /_尸)(13)为检验H,给定显著性水平】，记尸分布的分位数为Fi侦一为一心),检验 规则为F 5M-L(i)时接受丑。，否则拒绝。以上对SA,SE,ST的分析相当于对组间、组内等方差的分析，所以这种假设检验方法称 方必分析。1.3方差分析去将试验数据按上述分析、计算的结果排成表2的形式，称为单因素方差分析表 (Matlab中给出的方差分析表)表2单因素方力分析表方差来源平方和自由度均方1 pv分位数概率因素4r-1r-1Fr-l.n-r)Pv误差岛n rnr总和ST即一1最后 列给出大/值的概率Ff相方差分析一般用的显著性水平是：取1

17、= 0.01,拒绝丑0,称因素，的影响（或， 各水平的差异）非常显著；取1 = 0.01,不拒绝H。，但取a = 0.05，拒绝称因 素刀的影响显著；取& = 0.05,不拒绝H。，称因素4无显著影响。1.4 Matlab 实现Matlab统计工具箱中单因素方差分析的命令是anovaL若各组数据个数相等，称为均衡数据。若各组数据个数不等，称非均衡数据。（1）均衡数据处理均衡数据的用法为：p=anoval（x）返回值p是一个概率，当p时接受H。，X为阳X尸的数据矩阵，X的每一列是一个 水平的数据（这里各个水平上的样本容量nt = m ）o另外，还输出一个方差表和一个 Box 图。问题一的模型问题

18、一为：如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势问题的分析与模型的建立对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等 安全情况的变化趋势做出定量的综合评价，首先对数据处理按季 节分为17个子样本点的抽检值进行分析，得到各子样本食品质量 安全情况。运用层次分析法决定食品安全单位指标在综合评价中 权重。1）构造各季节食品安全质量综合评价的阶梯层次结构在构造深圳市食品领域综合评价模型时，考虑地理位置，人员结 构主要考虑老中青人员比例，经济结构等4项食品安全指标的综合效 果。构造其综合评价的递阶层次结构如下图：图：食品综合评价递阶层次结构2）建立V-A层

19、判断矩阵建立V-A层判断矩阵的依据是：考虑A层的四个因素对食品质量 安全综合评价的重要性时，认为年均指标A1最重要，且A1比微生 物、重金属、添加剂含量等值重要。而在微生物、重金属、添加剂含 量等值中，而各类食品都用到添加剂，所以，添加剂比微生物、重金 属等重要。根据判断矩阵标度内容及其含义，可得V-A层判断矩阵A：A1 A2A3A4-1233 一A11/213/23/2A2A =1/32/311A31/32/311 JA43）A-B层判断矩阵该层的矩阵元素Z确定如下：将每条准则AA2、A3、a4下各指标抽查值用分指数公式（1.1）标准化，其中C是在A.准则下第j个指标抽查值的均值（或最大 值

20、），C是第jJ个抽查指标的最大允许不合格值。对于每条准则下，令 3kj = max3 I j = 1，2,3,4i = 1,2,3,4 计算3 = min3 I j = 1，2,3,4. = 1,2,3,4 气二（3广3）/9ljij（1.2）rr =st3 -3 .-sA （ s t =A.（1.3）对于任意的3当 I r 1 1当, -1与3相比， 取3二1st 取广戏 取 3：二：1/ rt其判断矩阵元素3由以下确定st外 b = 1/ b=按式（1.1）可得相应于表1.1中各指标抽查值的标准化值由式子（1.2）（1.3）（1.4）得A-B层的判断矩阵4）用方根法计算出判断矩阵A , B

21、（）的归一化特征相量：VA、VB（n）5）计算判断矩阵B（i）,B（2）,B（3）的最大特征值人 并进行一致性检验max一致性指标人=w啤Lmax nwi=1iC.I. = y 其中，n为判断矩阵的阶数n -1CI计算一致性比例C.R.=工其中R.I.为平均一致性指标，当C.R. 0.1 R.I.时判断矩阵具有好的一直性。6）食品质量的综合评价指数的计算PI =岌 4 V . V G）j=1 i=1模型的求解及结果以2013年第一至四期食品安全抽样检验产品抽查样本的32种食 品的抽查值为例，求该抽查样本的食品质量安全综合评价指数。根据 题目给出的32种食品抽样数据值，统计得到下面的平均值表：表

22、1.12013年第一至四期食品安全抽样检验产品指标统计表表项目2010201120122013平均值8.25610.109242.43210.18286微生物平均值8.20750.111883.0250.14688重金属平均值8.20.10492.10.0675添加剂平均值8.38120.106411.41250.3125按式（1.1）可得相应于表1.1匚尸各指标抽查值的标准化值如表1.2所示。表1.2各指标抽查值的标准化值因子2010201120122013A10.688010.163850.121610.068571A20.683960.167810.151250.055077A30.68

23、3330.157350.1050.025312A40.698440.159610.0706250.11719式子（1.2）（1.3）（1.4）得A-B层的判断矩阵分别为:B（1）：B1B2B3B417.61578.22949B10.13131111.3843B20.12152111B30.111110.7223711B4B(2)：B1B2B3B417.38667.62379B10.13538111.6134B20.13117111.3763B30.111110.619820.726571B4B(3):B1B2B3B417.19417.91019B10.139111.8059B20.126421

24、11.0899B30.111110.553740.91751B4B(4):B1B2B3B417.724498.3325B10.1294611.27561B20.111110.7839311B30.12001111B44)用方根法计算出判断矩阵A ,B(i),B(2),B(3)的归一化特征相量:V =(V , V , V , V ) = (0.461539 0.230769 0.153846 a A 今 少 少0.153846);V二(V，V, 0.0802843) 12V (2) = ( V (2)，V (2)V二(V，V, 0.0802843) 12V (2) = ( V (2)，V (2)

25、，B、BiB20.0724441)V(2)，V(2) ) = (0.72337 0.104756 0.0994301B3B4V=(V，B、 BV=(V，B、 B10.0749515)V=(B0.0814922V，VBB0.08826822)VV=(B0.0814922V，VBB0.08826822)V，V)=0.734622 B3B40.0956177其中，的4个分量分别表示A层的4个因子(A1平均值、A2微生 物平均值、A3、重金属平均值、A4、添加剂平均值)在综合评价中的各自的权重。VB(1)的4个分量分别表示只考虑年均值A1的条件下4 个指标各自应占的权重；VB(2)、VB(3)、VB(

26、4)类推。5)计算判断矩阵B(i),B(2),B(3)的最大特征值人并进行一致性检max验人 =况 _L人i= 6.0i2,人2 = 6.08,人3 = 6.09 一致性指标maxnWmaxmaxmaxi=iiC.I. = 其中，n为判断矩阵的阶数n -1CI计算一致性比例C.R.=史土其中R.I.为平均一致性指标，查表R.I.R.I. =0.9，故C.R.=0 0.1C.R.（2）= 0 0.1C.R.（3）= 0 0.1当C.R.0.1时判断矩阵具有好的一致性。C.R.（4）= 0 0.1的综合评价指数的计算广.2610552.266）市场食品质量安匕 PI =岌 4 / . VB ）的综

27、合评价指数的计算广.2610552.26j=1 i=1根据深圳市食品安检局标准处推荐的评价标准（GB3838-2002 食品安检质量标准以及食品质量分级）为参考依据，见下表 1.3表1.3深圳市食品质量系数（PI）分级质量安全值级别食品质量评价0.80劣V严重劣食品这样就得到了 2013年第一至四期食品安全抽样检验产品抽查样本的 食品安全的优劣状况，并做出定量的综合评价。其它年份的食品安全 情况同理可求。求解程序见附录1。表1.4为求得的各个抽查样本的食品质量安 全状况和定量综合评价。其中C.R.为一致性指标（当C.R.0.1时有好 的一致性）详细的计算结果见附录2。2013年第一至四期食品安

28、全抽样检验产品抽查样本的食品安全情况 的综合评价可由加权平均法来作出判断。表1.4 2013年各抽查产品合格评价表抽查产品食品综合评 价指数C.R.（食品质量安全系数）超市速食面制 品0.2651150.00000000 000.000000000 00.00000008830.0000000000预包装调味0.2659770.00000000 000.00000008830.00000008830.0000000000预包装即食小 食品0.2846380.00000008830.000000000 00.00000008830.0000000000一次性纸杯塑 料杯0.2486960.000

29、00000 000.00000008830.00000000 000.0000000883饼干及膨化食 品0.258630.00000000 000.000000000 00.00000000 000.0000000000速冻面制品0.2478460.00000008830.00000008830.00000000 000.0000000000巧克力及其制 品0.2610550.00000000 000.000000000 00.00000000 000.0000000000乳制品及含乳 制品0.3574770.00000000 000.000000000 00.00000000 000.00

30、00000000酒店奶制品0.2901320.00000008830.00000008830.00000000 000.0000000000大米0.3788740.00000000 000.00000008830.00000008830.0000000000食用油0.2380920.00000000 000.00000008830.00000017660.0000000883广式烧卤熟肉 制品0.3214150.00000000 000.000000000 00.00000000 000.0000000000烟熏肉制品0.2946080.00000008830.000000000 00.000

31、00000 000.0000000883预包装茶叶0.2751090.00000000 000.000000000 00.00000008830.0000000883大型企事业单 位食堂0.8441380.00000000 000.00000008830.00000000 000.0000000000速冻汤圆食品0.2857430.00000008830.000000000 00.00000000 000.0000000000幼儿园食堂食 品0.2705070000.0000000883其中权值七可以反映不同食品质量对整个深圳市食品质量的 影响的不同作用：我们采用目标应用最多的因子实测值与标准

32、值为双重判定依据的赋 权方法一超标指数法或指数赋权法，即根据各食品的质量分级指数 C来确定权重。其公式为：L(i=1,2, n)乎Iii=1L(i=1,2, n)乎Iii=1I = ioi式中PI为各食品质量安全分级的指数，C为PI对应分级指数 的最大允许值。将我们计算出的PI值代入式1.o5计算得i47PIPI = t a j = 0.317 x x x x x x x 63187542结论：在食品质量影响因素中食品产地影响最大，食品加工次之， 季节影响最小，抽查地点几乎无影响。问题三：规准型抽样检验“民以食为天，食以安为先”，食品安全和卫生问题一直是关系国计民生 的重大问题。根据问题一、二

33、数据结果以定量比较评估的方法分析了各类影响 食品安全的因素及其可能造成的危害性的问题。改进后的食品抽检的办法以主 要食品为准则层建立了层次分析（AHP）模型，对影响食品安检的危害性因素 做出定量分析如问题二个影响因素大小。由问题一、二结果可以通过建立抽检 模型，即改进后的规准型抽样检验模型，并以蒙特卡罗法对抽检的全过程进行 模拟，得到相对误差逐渐趋向于0。就目前“营养强化.（查看更多）面粉”如 何抽检问题，本题在将各检测项目与各类食品安全影响因素一一对应的基础上， 建立逐步优化模型，最终得出一次全面检验的最小总成本。依据问题一二结果 通过参照计量统计检验程序给出了一套完备的连续多次跟踪抽检策略

34、，建立更 加切合实际的多目标规划的抽检策略。最终得到当检测数量达到100以上时， 该部门进行抽检工作时可达到检测可靠性高、成本低、工时少的目的。故针对 至少720个总批次的食品本题最终建议抽检125批次。综上，本文较好地解决 了当前食品安全问题中存在的矛盾关系。不仅提高了消费者的满意度，也兼顾 了生产者的利益，具有较高的社会经济价值。6模型的评价6.1模型的优点本文对于各种因素对于食品质量安全的影响建立了多元线性回 归模型，全面综合考虑了各个方面的因素，避免了单一因素分析的不 准确性，得出了合理的数学模型。并且通过各因素的显著性分析，找 到了影响食品质量安全的主要因素，较符合实际情况，模型可靠

35、，并 且模型相对简单，利于操作；该方法不仅适用于本题，也适用于其他 方面的数据预测，有实际背景，可运用于实践，具有广泛适用性。6.2模型的缺点本文忽略了除了所给因素之外的因素对食品安全的影响，与实际 问题存在偏差。同时是在假设各因素相互独立的情况下对食品质量的 影响进行分析，可能会导致误差7模型的改进与推广应用统计推断先判断食品的安全性，又有中心极限定律假设食品的合格服从正 态分布使其食品安全量化，并且量化成数据。比如超过量化后的数据为不合格， 没有超过为合格。各个国家的食品安全等级不一样，反映深圳市也是如此。运 用2010年、2011年和2012年深圳市的食品抽检数据统计量化后的数据分别作为

36、 各类产品安全等级。生活食品质量标准是指按特定生产方式生产，并经国家有 关的专门机构认定，或许使用绿色食品标志的无污染、无公害、安全、优质、 营养型的食品。在许多国家，绿色食品又有著许多相似的名称和叫法，诸如“生 态食品”、“自然食品”、“蓝色天使食品”、“健康食品”、“有机农业食品”等。影 响因素有原料生长土壤环境及食品的生产收获，此外食品的运输，加工，包装， 储存，销售，以及餐饮等每个环节皆可能影响食品的质量与安全。运用统计分 析中的方差，分析各个影响因子对食品安全运用统计分析中的方差，分析各个 影响因子对食品安全的影响程度的影响程度。最终使人们心中无污染、无公害、 安全、优质、营养型的食

37、品与市场食品安全等级数据差极限为零。相反如果人 们认识的食品安全数据比市场标准化低，或没有很好认识安全食品概念及本质。 我们可以通过媒体，广告，食品安全讲座等。让人们认识和提高对安全食品实 质。如果人们心目中的标准已经是市场安全等级，若是的话，就可以得出合理 的结论：人们追求更加持续性，更加科学化的生活质量标准。正好和可持续科 学发展观相吻合。模型中得到最优回归方程的方法是从包含全部变量的回归方程 中逐次剔除不显著因子，这种方法是在不显著因子不多时采用，当不 显著因子较多时，则工作量将会相当大，因为每剔除一个变量就得重 新计算回归系数。鉴于以上问题，我们引入了逐步回归分析的方法，它的基本思想

38、是将因子一个个引入，引入因子的条件是，该因子的偏回归平方和经 检验时显著的。同时，每引入一个新因子后，要对老因子逐个检验， 将偏回归平方和变为不显著的因子剔除。这种方法不需要计算偏相关 系数，计算较简便，并且由于每步都作检验，因而保证了最后所得的 方程中所有因子都是显著的。若回归方程是拟合好的，就可以进一步利用它来进行预报和控 制。预报问题，用统计数学的语言来说就是一个区间估计问题。在建 立气象预报、地震预报、自动控制等数学模型时，都可以用到本文的 模型。8参考文献马新民，王逸迅.概率与数理统计M.北京：机械工业出版社， 2010.刘卫国.MATLAB程序设计与应用M.北京：高等教育出版社，2

39、009.茹诗松.回归分析及其试验设计M.上海：华东师范大学出版社，1986.陈巧玲，沈晓斌，长江水质污染预测的灰色模型，泉州师范学院学报(自然科学)，2005.11孙靖南，邹志红，任广平.模糊综合评价在天然水体水质评价中的应用研究,环境污染治理技术与设备,2005 ,6 (2) :45248.李小燕灰色神经网络预测模型的优化研究, HYPERLINK /p-67589 /p-67589378.html邓聚龙，灰理论基础M，武汉:华中科技大学出版社,2002.傅立，灰色系统理论及其运用M，北京:科学技术文献出版社,1992.刘荣珍，赵军，模糊评价模型在长江水质评价中的应用，兰州交通大学学报,20

40、07.1212C.Henry Edwards &David E.penney微分方程及边值问题，北京：清华大 学出版社，2007附录：一数据统计名称抽取 合格 不合格2012.26即食蔬 菜学校学 生餐 娱乐场 所食品 水发水 产潮式打 冷食品 自制米 面油炸食 品即食水 果快餐食 品城中村 快餐食 品中央厨 房半成 品餐饮火 锅汤底 学校食 堂面制 品西式菜 肴餐前小 吃食品 餐饮领 域盐焗 类食品 工业区 餐饮 餐饮领 域乳及 其乳制2002000300300030255343406054634340606001001000150150015015006160183830100871320

41、019822001928880100100027270品69654555322001991797366564165623一次性 餐饮具 腌腊肉 制品 餐饮领 域茶叶 2012.25 饼干干果及 坚果食 品预包装 肉制品动物性水产干制 品、藻类 干制品 交通枢纽 商铺预包 装食品8019870189109糕点面包1009010预包装龟27270苓膏茶叶1201200焙烤咖啡27270及其制品饮料酒47470乳制品及1461460含乳制品火锅原料20016238食品罐头50482调味品100964塑料奶瓶橡胶16160奶嘴大米1501500食用植44440物油火锅底料产品27252腌腊肉25250

42、制品餐饮自2012.242002000. 制饮料餐饮大1001000米餐饮食 用油100982餐饮西式菜肴100973餐饮生食水产1361333品企事业单位食堂分装1001000菜肴.饭盒罗湖区.龙岗区18162餐饮肉丸食用植2012.2380773物油腐竹及淀粉制55532品批发商批发环50491节粮油士多店预包装40400食品乳制品及含乳53530制品批发商批发环节调味100919品集体食2012.223002991. 堂食品快餐食 品2002000使用冰100982粉丝腐100991竹即食果-H-71710疏2012.212012.(20)2012.19餐前小 吃 凉茶 即食蛋 制品 即

43、食熟 肉制品 即食水 产品 方便面 凉茶 小麦粉 茶叶 乳制品 及含乳 制品 生产领 域月饼 流通领 域月饼 餐饮领 域月饼 瓶装饮 用水 固体饮 料 饮料食 品 饮料酒 淀粉及 淀粉制 品 冷冻饮 品 速冻米 面制品 乳制品 及含乳 制品 餐具洗 涤剂流通领 域畅销 大米1321257484801010035323262333838061601393901009825050013212936058211811803533250491807826663379772109127270494904040077770食品和食品用2012.1875573817相关产品预包装2012.1750491罐头

44、谷物粉类制成80791品食糖60600蛋制品50482食用植86860物油乳制品及含乳50500制品福田区学校周37370边食品餐饮茶2012.16叶2001991凉拌菜1531530现配奶茶2001964米面制 品917912沙县小吃1691690商业中心餐饮48471原料2012.15市场4生产领2012.1411110域粽子流通领域粽子16151餐饮领域粽子88826餐饮酒2012.13类97943生食水1261242主口 产品沙津2001991寿司1001000餐饮原料饭盒120501205000煎炸用食物油100991保鲜膜学校食50446堂米油南山区商业中4004000心餐饮原料r

45、e 预包装11011002012.12豆制品预包装26224饼干食 品80737膨化食 品80782糖果40391湿纸巾白酒洋20200酒学校周37343边小卖铺食品63612葡萄酒乳制品27261及含乳 制品盐焗食505002012.11品45414调味品100964油炸食 品89890小麦粉福田区商业区1001000专项执法抽查494902012.（10）蜂产品46442调味品796811炒货38371餐巾纸35350十货45441乳制品 及含乳 制品乳制品505002012.9及含乳制品40400巧克力及制品50500预包装果冻50500预包装水产品40355餐饮乳2012.8制品及含乳

46、制40400品大米食用植物油100200100199012012.7散装水残品66633加油站便利店28315268蜜饯473215清洗产 品50500罐头食 品交通枢20018614纽便利店散装和200181192012.6预包装食用油1531503腌腊肉原料79709菜肴98620966一次性用5预包装 肉制品15014010糖果10090102012.42012.32012.22012.12011.192011.18炒货食 品 乳制品 大米 油炸食 品 鲜榨果 蔬汁 即食果 蔬 佐餐饮 料 速冻米 面食品 水饺和 汤圆 预包装 乳制品 及含乳 制品 米面制 品

47、培炒咖 啡及其 制品 龟苓膏 食用油 工业食 堂原料 商业中 心餐饮 单位原 料饮料单 位酒类 油条 集中式 消毒餐 具 预包装 乳制品 及含乳 制品 果冻96942299299010710618281110099130030001009911131130121121015014916059150500257251628127652982944200197324186549510393939073649巧克力及其制57543品方便面68671食用植艾2092063物油即食豆_r制品31301腌品肉80755十货坚79736果调味品15513124预包装熟肉制81801品2011.172011.

48、16食品1559151148乳制品40400豆制品80800烧卤熟18217111肉一次性包装材604911料预包装2011.15乳制品80及含乳80制品大米78762茶叶80800腐竹粉77前m丝776710饼干806812十货食68662品蜂产品806218膨化食如,80755品小麦粉48462预包装糕点面80746包糖果40400葡萄酒20200熟制鸭脖.生食水102594314 柠60591产品寿司60600西餐菜60600肴冰淇淋44431餐前小80719吃小麦粉47470大米50500即食现配奶茶50500沙拉44440凉拌菜63630茶叶80791凉茶53530预包装13乳制品1

49、551550及含乳制品调味品60591预包装即食豆987810制品蛋制品49472食糖60600凉茶60573谷物及淀粉制60600品罐头食51501品饮料80791瓶装饮48453用水固体饮50491料冷冻饮16151品即使水503812产品散装糕4030102011.2011.2011.122011.112011(10)2011.92011.82011.7点面包 保鲜膜 月饼 预包装 乳制品 及含乳 制品 儿童食 品 肉制品 十货 烘烤食 品 餐饮单 位厨房 一次性 用品 餐饮单 位使用 十货 餐饮单 位制作 油炸食 品 餐饮经 营单位 使用乳 制品 餐饮经 营单位 使用调 味品 市民日

50、常消费 的各种 原辅料 保鲜膜 鲜榨果 蔬汁 使用油 预包装 乳制品 及含乳 制品31310280278239390100982797904847150464494545047350428353505049111951150453030010010002011983404002011.62011.52011.4蜂蜜 餐饮服 务单位 的即食 熟肉制 品 餐饮经 营单位 的食品 原料腌 腊肉 餐饮经 营单位 的食品 原料使 用油 餐饮经 营单位 使用粮 食和调 味品 餐饮经 营单位 使用水 产十制 品 预包装 茶叶 预包装 饮料酒 预包装 糕点面 包饼干 预包装 乳制品 及含乳 制品 十货一次性

51、餐具 餐饮服 务单位 使用纸 品52484115833220164201914040010098291910969249585107877197841383785200175252011.32011.22011.12011.元旦前2010 豆制品2010 小麦粉2010 商业中心餐饮经 营单位 制作糕 点八、 餐饮经 营单位 制作面 类食品 加油站 便利店 鸡精燕 窝饮品 米面制 品 十果坚 果 休闲食 品婴幼儿 用奶瓶 奶嘴 预包装 冷冻食 品生产环 节抽查 流通环 节抽查 餐饮消 费环节 抽查 生产环 节抽查 流通环 节抽查 餐饮消 费环节 抽查豆制品小麦粉 商业中 心餐饮29817812

52、01981663239535936343133823691350473649522274646020219844194109126212085421802085954694541551549916615588278573125050026825810餐饮食 品原料 2010 双 节前2011 双节全2010.32010.22010.32012.2单位原 料 应节食 品 生产环 节抽查 流通环 节抽查 餐饮服 务环节 应节食 品 生产环 节抽查 流通环 节抽查 餐饮服 务环节 大米 生产企 业 市场畅 销 餐饮单 位 中秋佳 节将全 月饼 月饼生 产企业 餐饮单 位制作 市场畅 销 酱油和 耗油

53、 乳饮料 和植物 蛋白饮 料 啤酒 预包装 沙拉酱 预包装159315326119519148458202555352132143613756119519148458202539636432333332178771747401811810377376112812711981980515102220217170111101414040382固体饮 料 食饮具 瓶装饮 用水 洗涤消 毒剂 保鲜膜 西餐菜 肴 凉菜 罐头食 品 现冻动 物性水 产品 预包装 粗粮谷 物 预包装 饮料 鲜榨果 蔬汁饮 料 水食水 产品 预包装 水产品 预包装 肉制品 预包装 冷冻饮 品 预包装 淀粉及 谷物制 品 洗涤

54、剂 产品 预包装 凉茶 保鲜膜 一次性 饭盒固体饮 料 食饮具 瓶装饮 用水 洗涤消 毒剂 保鲜膜 西餐菜 肴 凉菜 罐头食 品 现冻动 物性水 产品 预包装 粗粮谷 物 预包装 饮料 鲜榨果 蔬汁饮 料 水食水 产品 预包装 水产品 预包装 肉制品 预包装 冷冻饮 品 预包装 淀粉及 谷物制 品 洗涤剂 产品 预包装 凉茶 保鲜膜 一次性 饭盒1331312312921816232302251246518218116059116116015050080800200198260582342861391345505001651650474704848029290917615一次性筷子16315211餐饮使 用油2482426餐饮使用茶叶32275生产环 节粽子18162流通环 节粽子44440流通环 节咸蛋 餐饮服40400务环节49481粽子 生产环 节儿童1121075食品 流通环 节儿童食品2192127生产环节食品1483141568流通环 节食品3022948餐饮消费环节79790生产环节应品流通环节应节食品1391390餐饮消 费环节 食端午前2010 六一将全2010 春节