A陶子數(shù)學(xué)模型--食品安全的抽檢問(wèn)題論文

1、食品安全的抽檢問(wèn)題摘要 食品的質(zhì)量和衛(wèi)生問(wèn)題是關(guān)系到民生的大問(wèn)題，因此，對(duì)食品的檢查顯得非常重要。本文結(jié)合實(shí)際，應(yīng)用AHP方法、分層抽樣和線性目標(biāo)規(guī)化方法，建立了集時(shí)間、費(fèi)用和效果為一體的數(shù)學(xué)模型，具體如下。 對(duì)于問(wèn)題一，我們首先將主要食品進(jìn)行分類，然后將影響食品安全的因素主要分為生物性污染、化學(xué)性污染、物理性污染三大類，并將這三類污染所造成的主要危害歸納為七類，接著采用AHP法對(duì)問(wèn)題進(jìn)行定量分析，最后通過(guò)一致性檢驗(yàn)并得出其危害性的大小，得到結(jié)果細(xì)菌危害最嚴(yán)重，食品添加劑導(dǎo)致的危害次之等。 對(duì)于問(wèn)題二，針對(duì)部分主要產(chǎn)品，我們先采用了分層抽樣的方法對(duì)不同品牌不同批次的產(chǎn)品進(jìn)行抽檢，建立了樣本分配

2、率、樣本方差、總體抽樣率、分層抽樣率等函數(shù)方程，然后對(duì)上一步所抽到的批次利用線性目標(biāo)規(guī)劃的方法，建立了集時(shí)間較短、成本費(fèi)用較低和抽樣效果較好的抽檢模型線性目標(biāo)規(guī)劃模型，并利用統(tǒng)計(jì)學(xué)原理對(duì)檢測(cè)誤差進(jìn)行分析。最后，我們根據(jù)模型針對(duì)乳制品中的酸奶進(jìn)行模擬檢驗(yàn)，檢驗(yàn)的結(jié)果誤差百分比為4.24%<5%，可靠性較高。對(duì)于問(wèn)題三，我們利用問(wèn)題二所建立的模型制訂了一種較為合理的抽檢方案（根據(jù)假設(shè)總共抽檢79個(gè)批次，每個(gè)批次抽檢2個(gè)項(xiàng)目）。然后，我們進(jìn)行了可靠性分析，抽檢的誤差百分比為1.15%<5%，可靠性較高。對(duì)于問(wèn)題四，它實(shí)際是在問(wèn)題三的基礎(chǔ)上，對(duì)面粉進(jìn)行多次跟蹤抽檢。我們對(duì)問(wèn)題二所建立的模型

3、進(jìn)行了改進(jìn)，引入新的變量建立函數(shù)關(guān)系，并運(yùn)用MATLAB優(yōu)化工具箱進(jìn)行求解，得出了最佳的抽檢策略和抽檢數(shù)量（結(jié)果為跟蹤抽檢3次，共抽檢113個(gè)批次），所得結(jié)果可靠性較高、成本較低，且工時(shí)比較少，用計(jì)算機(jī)進(jìn)行模擬檢驗(yàn)時(shí)效果比較樂(lè)觀。最后，我們對(duì)模型的優(yōu)缺點(diǎn)進(jìn)行了評(píng)價(jià)，討論了其推廣應(yīng)用的價(jià)值，并主管部門寫了一份報(bào)告，提出了一些解決問(wèn)題的可行性建議，可為主管部門和市民提供一些參考。關(guān)鍵詞 ：AHP法，分層抽樣，目標(biāo)規(guī)劃、統(tǒng)計(jì)分析、可靠性 一 問(wèn)題的重述改革開放三十年來(lái)，我國(guó)人民生活水平在不斷地提高，食品安全和衛(wèi)生問(wèn)題越來(lái)越受到人們的關(guān)注。近幾年來(lái)，先后出現(xiàn)了蘇丹紅、瘦肉精、三聚氰胺等事件，以及各種不

4、利于健康的食品添加劑、強(qiáng)化劑問(wèn)題的出現(xiàn)，食品安全和衛(wèi)生的檢測(cè)已成為全社會(huì)，乃至政府有關(guān)部門重點(diǎn)關(guān)注的問(wèn)題之一。食品的質(zhì)量和衛(wèi)生問(wèn)題涉及到原材料的使用、生產(chǎn)加工、運(yùn)輸與貯存、流通與銷售等環(huán)節(jié)，在每一個(gè)環(huán)節(jié)上出現(xiàn)差錯(cuò)，都會(huì)導(dǎo)致食品出現(xiàn)安全和衛(wèi)生問(wèn)題，食品質(zhì)量和衛(wèi)生的檢測(cè)工作在實(shí)際顯得非常重要。但是，由于食品的種類、品牌和批次繁多，從生產(chǎn)加工到銷售食用中間環(huán)節(jié)復(fù)雜，質(zhì)檢部門不可能對(duì)所有食品做到全面的質(zhì)量檢測(cè)，一般只能做一定的抽檢。當(dāng)然，對(duì)食品進(jìn)行抽檢也需要一定人力、物力和財(cái)力（即成本費(fèi)用），抽檢的越多檢測(cè)效果就越好，但需要的時(shí)間就越長(zhǎng)，其成本費(fèi)用也就越高。為此，應(yīng)該如何抽檢，既能保證較好的檢測(cè)效果，

5、又能節(jié)省時(shí)間和成本費(fèi)用？請(qǐng)根據(jù)實(shí)際情況，建立數(shù)學(xué)模型分析研究下列問(wèn)題：（1）根據(jù)主要食品的分類，試分析影響各類食品安全的因素，對(duì)其可能的危害性做出定量的比較評(píng)估分析。（2）針對(duì)部分主要食品，結(jié)合實(shí)際建立合理的抽檢模型，給出檢測(cè)誤差（即檢測(cè)的可靠性）分析的方法，并對(duì)模型進(jìn)行模擬檢驗(yàn)。（3）面粉是我國(guó)中北部地區(qū)主要的主食原材料，不妨就已推廣食用的“營(yíng)養(yǎng)強(qiáng)化面粉”抽檢問(wèn)題進(jìn)行討論，“營(yíng)養(yǎng)強(qiáng)化面粉”的配方標(biāo)準(zhǔn)如表1所示。假設(shè)某地區(qū)現(xiàn)有12個(gè)品牌的營(yíng)養(yǎng)強(qiáng)化面粉產(chǎn)品，每個(gè)品牌每月將有不少于60個(gè)批次（即同一企業(yè)、同一條生產(chǎn)線、同一批投料、同一班次生產(chǎn)的產(chǎn)品為1個(gè)批次）的產(chǎn)品在市場(chǎng)上銷售，質(zhì)檢部門要做一次全

6、面的質(zhì)量檢查，請(qǐng)你幫助制訂一種合理的抽檢方案，并分析其檢測(cè)的可靠性。所需要檢驗(yàn)項(xiàng)目、標(biāo)準(zhǔn)、成本和工時(shí)如表2所示。（4）針對(duì)問(wèn)題（3），如果質(zhì)檢部門需要連續(xù)進(jìn)行多次跟蹤抽檢，請(qǐng)你給出相應(yīng)的抽檢策略和最佳的抽檢數(shù)量，使其檢測(cè)可靠性盡量高、成本盡量低、工時(shí)盡量少，并用計(jì)算機(jī)進(jìn)行模擬檢驗(yàn)。（5）請(qǐng)根據(jù)你們的研究，深入分析食品安全存在的隱患和根源，并提出有效可行的解決問(wèn)題辦法和建議，可供主管部門和市民參考。二 模型的假設(shè)1、假設(shè)模型求解過(guò)程中所用的數(shù)據(jù)都是合理的；2、假設(shè)同一類食品每次抽檢的項(xiàng)目數(shù)相同；3、假設(shè)不考慮抽檢人數(shù)的限制；4、假設(shè)每天抽檢的最大時(shí)間為法定工作時(shí)間8小時(shí)，一個(gè)月工作22天；5、假

7、設(shè)抽檢的最大費(fèi)用有限制；三 符號(hào)約定-總體均值的估計(jì)量（st表示分層）；-層權(quán)（）,i=1,2N;L-抽樣系統(tǒng)中總的生產(chǎn)批次；-第i層生產(chǎn)批次數(shù)；-第i層的樣本均值；-第i層抽取的樣本量；-第i層的方差；-第i層的抽樣比（/）；-為L(zhǎng)在1-置信水平下的絕對(duì)誤差限；-標(biāo)準(zhǔn)整天分布的雙側(cè)分位數(shù)；-第i層每個(gè)樣本單位的平均費(fèi)用；-第i層每個(gè)樣本單位的平均檢驗(yàn)工時(shí)；-固定費(fèi)用；-每個(gè)抽檢批次的檢查項(xiàng)目；-抽檢的最長(zhǎng)時(shí)間；-抽檢的最大費(fèi)用；-抽檢的總項(xiàng)目數(shù)；-每次每項(xiàng)抽檢的平均時(shí)間；-每次每項(xiàng)抽檢的平均費(fèi)用；四 問(wèn)題的分析與模型的建立4.1對(duì)問(wèn)題一的分析與模型建立4.1.1問(wèn)題一的分析為了分析與比較影響

8、食品安全因素所產(chǎn)生的危害程度，首先我們將主要食品分為：肉制品、乳制品、釀造食品和蛋白質(zhì)含量較高食品等；接著分析確定影響食品安全的因素，我們從微生物污染、化學(xué)污染和物理污染三方面進(jìn)行分析與確定；然后對(duì)于各因素影響食品產(chǎn)生的危害，其一是由微生物污染導(dǎo)致的細(xì)菌危害、病毒危害和寄生蟲危害；其二是由化學(xué)污染引起的食品添加劑（如非食用添加劑）導(dǎo)致的危害、環(huán)境污染（如農(nóng)藥，包括滅鼠藥和獸醫(yī)用藥）導(dǎo)致的化學(xué)危害和天然存在的化學(xué)危害；其三是由物理性污染引起的重金屬中毒等危害。因此，問(wèn)題一就可以歸結(jié)為各因素影響食品導(dǎo)致的危害對(duì)食品質(zhì)量的影響權(quán)重問(wèn)題。解決這類問(wèn)題首先要統(tǒng)計(jì)分析各類影響指標(biāo)的數(shù)值特征。然后再對(duì)其進(jìn)行

9、歸一化處理，并利用層次分析法對(duì)其權(quán)重進(jìn)行賦值。最終得到各因素對(duì)食品安全的影響排名。4.1.2問(wèn)題一的模型建立（1）根據(jù)影響食品安全的各因素與危害間的關(guān)系，我們建立的不完全層次分析算法的基本結(jié)構(gòu)圖如下： （2）根據(jù)對(duì)問(wèn)題一的分析，我們建立了如下成對(duì)比較矩陣假設(shè)要比較某一層的個(gè)因素對(duì)上一個(gè)因素O的影響，取兩個(gè)因素和，用表示和對(duì)的影響之比，全部比較的結(jié)果可用成對(duì)比較矩陣：，表示，并稱以上矩陣為正互反矩陣。（3）一致性檢驗(yàn)及權(quán)重向量的確定成對(duì)比較矩陣通常不是一致陣，為了能用它的最大特征根的特征向量作為被比較因素的權(quán)向量，其不一致程度需要控制在一定的范圍內(nèi)。計(jì)算方法如下：定義：一致性指標(biāo)；一致性比率為;

10、本文中，對(duì)應(yīng)的.當(dāng)時(shí)，認(rèn)為矩陣的不一致程度在容許的范圍內(nèi)。其中最大特征根的特征向量即為相應(yīng)的權(quán)重向量。4.2問(wèn)題二的分析與模型建立4.2.1問(wèn)題二的分析從問(wèn)題二的題設(shè)我們知道，為了建立合理的抽檢模型，我們首先需要對(duì)現(xiàn)有的統(tǒng)計(jì)信息進(jìn)行定量分析。由于問(wèn)題二的求解幾乎需要考慮文中提供的所有已知信息，所以為了簡(jiǎn)化數(shù)學(xué)模型的計(jì)算復(fù)雜度，我們可以在某些直觀問(wèn)題上進(jìn)行簡(jiǎn)單的定性分析。例如，原文附錄中給出了主要食品和主要抽檢項(xiàng)目，為此我們可以主觀上認(rèn)為我們要考慮的部分食品就是附錄1中的部分主要食品，要抽檢的項(xiàng)目就是附錄2中給的8個(gè)項(xiàng)目。在基本確定主要食品和抽檢項(xiàng)目后，首先我們利用分層抽樣模型確定這些主要食品中

11、每種品牌所要抽檢的批次數(shù)的最優(yōu)值，然后，我們利用法定的工作時(shí)間確定最大抽檢時(shí)間和費(fèi)用。最后計(jì)算出每種品牌中要抽的批次數(shù)及每個(gè)批次的藥抽檢的項(xiàng)目數(shù)，對(duì)其優(yōu)劣進(jìn)行綜合評(píng)價(jià)。4.2.2問(wèn)題二的模型（一）分層抽樣通過(guò)上述對(duì)問(wèn)題二的分析，我們利用分層抽樣法建立下列函數(shù)關(guān)系： 設(shè)整個(gè)抽樣為一個(gè)系統(tǒng)L，系統(tǒng)中含有種品牌，每個(gè)品牌的生產(chǎn)批次數(shù)為，則總體均值的估計(jì)量： （1）第i層的方差： （2）由誤差知識(shí)可得： （3）聯(lián)立上述各式可得： （4）調(diào)查費(fèi)用函數(shù)： （5）各層樣本量的最優(yōu)分配： （6）由于每個(gè)樣本的調(diào)查費(fèi)用相同，故上式可以化為 （7）聯(lián)立（4）和（7）可得，總體抽樣率為 （8）也可得到分層抽樣率為

12、（9）（二）模型的建立線性目標(biāo)規(guī)化在以上函數(shù)關(guān)系以及抽檢費(fèi)用有限制的條件下，模型以抽檢項(xiàng)目數(shù)最多為目標(biāo)函數(shù)所建立的模型如下： S.T. （10）（2）檢測(cè)誤差（即檢測(cè)的可靠性）分析的方法: （11）4.3問(wèn)題三的分析面粉是一種十分重要的原材料，因此對(duì)面份的抽檢具有代表性。根據(jù)題目的要求，質(zhì)檢部門要做一次全面的質(zhì)量調(diào)查，其抽檢方案是問(wèn)題二模型的具體應(yīng)用。所以，我們只需將問(wèn)題二建立的模型具體化，得出具體的抽檢方案，并對(duì)檢測(cè)的可靠性進(jìn)行分析即可。4.4問(wèn)題四的分析與模型的建立4.4.1問(wèn)題四的分析對(duì)于問(wèn)題四，它其實(shí)是在問(wèn)題三的基礎(chǔ)上，對(duì)面粉的質(zhì)量進(jìn)行多次跟蹤抽檢，我們對(duì)問(wèn)題二所建立的模型進(jìn)行了改進(jìn)，

13、引入新的變量時(shí)間，建立新的函數(shù)關(guān)系，并運(yùn)用MATLAB優(yōu)化工具箱進(jìn)行求解，得出了最佳的抽檢策略和抽檢數(shù)量，使檢測(cè)可靠性較高、成本較低，且工時(shí)比較少，用計(jì)算機(jī)進(jìn)行模擬檢驗(yàn)時(shí)效果比較樂(lè)觀。4.4.2問(wèn)題四的模型建立由多次跟蹤抽檢可知：抽檢批次數(shù)與時(shí)間是成反相關(guān)的，故我們可以作L與t的函數(shù)：,分別使抽檢的效果最好、費(fèi)用最低、時(shí)間最少。目標(biāo)函數(shù)為：五 模型的求解5.1模型一的求解5.1.1確定準(zhǔn)則層對(duì)目標(biāo)層的權(quán)重向量A=,=3.0385，權(quán)重向量為w=（0.6370，0.2583，0.1047）5.1.2確定方案層對(duì)準(zhǔn)則層的權(quán)重向量=,= ，=1.求解結(jié)果如下：（編程程序見(jiàn)附件一）矩陣特征值一致性指標(biāo)

14、CI隨機(jī)一致性指標(biāo)RI一致性比率CR權(quán)重向量3.03850.019250.580.033<0.1（0.6370，0.2583，0.1047）3.06490.032450.580.016225<0.1（0.6491，0.2790，0.0719）1000<0.1（1）由上表中<0.1可知矩陣、都通過(guò)了一致性檢驗(yàn)。5.1.3確定影響因素的權(quán)重向量組合權(quán)重：(i)=(0.6370*0.6370,0.2583*0.6370,0.1047*0.6370,0.6491*0.2583,0.2790*0.2583,0.0719*0.2583,1*0.1047)=(0.405769,0.1

15、645371,0.0666939, 0.16766253, 0.0720657, 0.01857177,0.1047)=（0.4058，0.1645，0.0667，0.1677，0.0721，0.0186，0.1047）CR=(0.01925*0.6370+0.03245*0.2583+0*0.1047)/(0.58*0.6370+0.58*0.2583+0*0.1047)=0.03976<0.1,通過(guò)了組合一致性檢驗(yàn)。5.1.4排名由上述組合權(quán)重可知，各因素對(duì)食品安全的影響程度由強(qiáng)到弱依次為0.40580.16770.16450.10470.07210.06670.0186細(xì)菌污染食品

16、添加劑導(dǎo)致的危害病毒危害重金屬中毒等環(huán)境污染導(dǎo)致的危害寄生蟲危害天然存在的化學(xué)危害5.2問(wèn)題二的求解（1） 層次抽樣模型:主要食品包括肉制品、乳制品、釀造食品及高蛋白食品等，針對(duì)部分主要食品進(jìn)行模擬抽檢。首先針對(duì)部分主要食品乳制品中的酸奶討論，進(jìn)行分層。第一層為：蒙牛、伊利、光明、百萬(wàn)莊園、三元、娃哈哈六個(gè)品牌的批次數(shù)，然后對(duì)每層隨機(jī)抽樣，得到樣本。分層抽樣圖如下：（2） 部分參數(shù)的假設(shè)：各品牌生產(chǎn)批次序列為：置信度取95%，根據(jù)對(duì)其他參數(shù)的討論，按照ISO質(zhì)量管理體系和AQL抽樣標(biāo)準(zhǔn)得到實(shí)際中抽檢概率在8%10%之間能夠達(dá)到抽檢效果，標(biāo)準(zhǔn)差取值范圍在0.0370.042之間，并在此區(qū)間內(nèi)等距

17、選取六個(gè)數(shù)進(jìn)行計(jì)算，六個(gè)數(shù)分別為，d取0.01。 由標(biāo)準(zhǔn)正態(tài)分布表查表可知：=1.64。其他參數(shù)的計(jì)算：N = =440 （批） 所以由有： ； （3） 分層抽樣率的計(jì)算由（9）計(jì)算得到各層的抽樣率如下表：酸奶品牌蒙牛伊利光明百萬(wàn)莊園三元哇哈哈抽樣率0.08130.08350.08570.08790.09010.0923假設(shè)t=1.5,p=60,Pmax=6000，Tmax=176, =11項(xiàng)目數(shù)為2.56413項(xiàng)。誤差：=440/6=73.3=10.8013=1.6530誤差百分比：1.6530/39=4.24%<5%權(quán)重=（0.2175，0.2175，0.0725，0.0725，0.

18、0435，0.0435，0.0362，0.0311，0.2175，0.0242，0.0242）排序：三聚氰胺,維生素A,維生素B11，維生素B1、維生素B2、鐵、鋅、維生素B5、鈣、水分、蛋白精、硫氰酸鈉隨及抽檢三個(gè)項(xiàng)目，前三項(xiàng)抽檢概率最大，分別為三聚氰胺,維生素A,維生素B11.5.3 問(wèn)題三的求解與結(jié)果分析5.3.1 問(wèn)題三的求解部分參數(shù)的假設(shè)：各品牌生產(chǎn)批次序列為：置信度取95%，d取0.01。根據(jù)對(duì)其他參數(shù)的討論，按照ISO質(zhì)量管理體系和AQL抽樣標(biāo)準(zhǔn)得到實(shí)際中抽檢概率在8%10%之間能夠達(dá)到抽檢效果，標(biāo)準(zhǔn)差取值范圍在0.03250.0365之間，并在此區(qū)間內(nèi)選取12個(gè)數(shù)進(jìn)行計(jì)算，12

19、個(gè)數(shù)分別為（5）其他參數(shù)的計(jì)算：N = =880（批）所以由有： ； 用MATLAB對(duì)函數(shù)進(jìn)行優(yōu)化可得：項(xiàng)目數(shù)為1.01272項(xiàng)。誤差：=880/12=73.3=8.4782=0.9106誤差百分比：0.9106/79=1.15%<5%權(quán)重：=10，CI=0,CR=0,通過(guò)了一致性檢驗(yàn)（程序見(jiàn)附件6）=（0.2779，0.2779，0.0926，0.0926，0.0556，0.0556，0.0463，0.0397，0.0309，0.0309）排序：三聚氰胺,維生素A,維生素B11，維生素B1、維生素B2、鐵、鋅、維生素B5、過(guò)氧化苯甲酞、鈣、水分隨及抽檢三個(gè)項(xiàng)目，前三項(xiàng)抽檢概率最大，分別

20、為三聚氰胺,維生素A,維生素B115.3 問(wèn)題四的求解與結(jié)果分析用MATLAB（程序見(jiàn)附件7）做出函數(shù)當(dāng)n=1、n=2、n=3、n=4的圖形，如圖示：由圖形可得：當(dāng)取=4時(shí)，即跟蹤抽檢四次時(shí)，圖號(hào)效果費(fèi)用時(shí)間圖11658250330圖21135650226圖3934650186圖4854250170當(dāng)取=3時(shí)，即跟蹤抽檢三次時(shí)，圖號(hào)效果費(fèi)用時(shí)間圖11457050290圖21085400216圖3924600184圖4854250170由表中內(nèi)容可得：L與的函數(shù)為,共需跟蹤抽檢三次，共抽檢113個(gè)批次。六 模型的評(píng)價(jià)6.1 模型一的評(píng)價(jià)模型一采用的AHP方法是一種定性與定量相結(jié)合的、系統(tǒng)的、層次

21、化的分析方法，能夠解決涉及到社會(huì)、經(jīng)濟(jì)、人文等因素的決策、評(píng)價(jià)、分析、預(yù)測(cè)等問(wèn)題。但這種方法同樣具有局限性，它只能從原有方案選優(yōu)，不能生成新方案，且它的判斷結(jié)果都是比較粗糙的，不適于精度要求很高的問(wèn)題，另外模型中的成對(duì)比較矩陣人的主觀因素的作用較大，可能有時(shí)難以為眾人接受。6.2 模型二的評(píng)價(jià) 模型二采用的是分層抽樣的方法，它對(duì)不同品牌不同批次的產(chǎn)品進(jìn)行抽檢，建立了樣本分配率、樣本方差、總體抽樣率、分層抽樣率等函數(shù)方程，再次利用統(tǒng)計(jì)學(xué)原理對(duì)檢測(cè)誤差進(jìn)行分析。這種抽樣調(diào)查的方法適用性特別強(qiáng)，對(duì)任何特征分層來(lái)比較明顯的系統(tǒng)均可以采用這種方法進(jìn)行抽樣。而且這種可靠性較高。七 模型的推廣本文中所建立的

22、模型不僅對(duì)于解決題目中的問(wèn)題具有很強(qiáng)的實(shí)用性，而且具有很好的推廣性。例如，本文運(yùn)用的AHP方法適用于其他類似一類選優(yōu)排序問(wèn)題、合理安排住房問(wèn)題等都具有很強(qiáng)的適用性。本文針對(duì)問(wèn)題二所建立的模型具有較強(qiáng)的普適性，能夠?qū)Ω黝愔饕称愤M(jìn)行精確度比較搞的抽檢，而且對(duì)批量生產(chǎn)的各種零件、衣服等的質(zhì)量抽檢提供了一種較為簡(jiǎn)單高效的方法，具有較強(qiáng)的研究意義。八 解決問(wèn)題的可行辦法和建議近年來(lái)，我國(guó)食品工業(yè)始終保持持續(xù)快速的增長(zhǎng)趨勢(shì)。據(jù)有關(guān)資料介紹，我國(guó)食品安全狀況得到有效改善和提高。具體表現(xiàn)在國(guó)家抽檢結(jié)果顯示，2008年上半年我國(guó)糧食加工品、副食調(diào)味品、肉制品、酒類、飲料飲品、休閑食品33類3288家企業(yè)生產(chǎn)的

23、3813種食品抽樣合格率為98.4%，同比提高4個(gè)百分點(diǎn)，比2006年提高近9個(gè)百分點(diǎn)。 但是，目前我國(guó)食品安全狀況距離廣大群眾的要求還存在不小差距，主要存在五大問(wèn)題：一是菌落指數(shù)超標(biāo)和農(nóng)藥殘留超標(biāo)；二是存在違法濫用非食品用添加劑情況；三是部分企業(yè)食品安全保障管理體系不完善，檢測(cè)系統(tǒng)不健全，易發(fā)生質(zhì)量事故；四是食品流通領(lǐng)域經(jīng)營(yíng)秩序不規(guī)范，假冒偽劣食品屢禁不絕；五是消費(fèi)者食品安全知識(shí)和維權(quán)教育亟待加強(qiáng)。經(jīng)過(guò)分析，我們發(fā)現(xiàn)食品安全問(wèn)題存在的原因主要有：1、生產(chǎn)經(jīng)營(yíng)者法制意識(shí)淡薄是產(chǎn)生食品安全問(wèn)題的根源。如個(gè)別經(jīng)營(yíng)戶見(jiàn)利忘義，為一些劣質(zhì)食品提供了可乘之機(jī)。比如食用油、餅干、糕點(diǎn)等到食品國(guó)家

24、實(shí)行生產(chǎn)許可證制度，但糧油店特別是農(nóng)村流動(dòng)經(jīng)營(yíng)戶大量經(jīng)銷大桶散裝油，有的摻雜使假并游擊作戰(zhàn)，執(zhí)法部門難以跟蹤監(jiān)管，消費(fèi)者上當(dāng)受騙的情況時(shí)有發(fā)生；2、食品安全法律體系結(jié)構(gòu)不科學(xué)、有法難依。3、食品安全監(jiān)督管理體制存在弊端，導(dǎo)致對(duì)食品安全問(wèn)題監(jiān)管不力。比如一些食品經(jīng)銷店，特別是農(nóng)村小店大量購(gòu)進(jìn)小作坊“三無(wú)”產(chǎn)品，并銷售霉變過(guò)期食品；4、食品安全標(biāo)準(zhǔn)和監(jiān)測(cè)手段落后，不能適應(yīng)監(jiān)管要求。主要表現(xiàn)在食品安全標(biāo)準(zhǔn)技術(shù)指標(biāo)落后，且交叉重疊，質(zhì)監(jiān)局制定質(zhì)量標(biāo)準(zhǔn)，各行業(yè)主管部門制定行業(yè)標(biāo)準(zhǔn)，衛(wèi)生部制定衛(wèi)生標(biāo)準(zhǔn)，難免使食品標(biāo)準(zhǔn)政出多門，指標(biāo)不統(tǒng)一。食品檢測(cè)是保障食品安全的重要手段，檢驗(yàn)設(shè)備、檢測(cè)人員分散在衛(wèi)生、質(zhì)監(jiān)

25、、農(nóng)業(yè)等部門，難以滿足現(xiàn)有的食品安全監(jiān)督綜合執(zhí)法的需要，如流通領(lǐng)域是食品安全監(jiān)管的重要環(huán)節(jié)，執(zhí)法部門檢查食品，識(shí)別假冒偽劣商品，基本上是借助經(jīng)驗(yàn)，依靠眼看手摸，看是否是“三無(wú)”產(chǎn)品，很難發(fā)現(xiàn)食品內(nèi)在質(zhì)量上的問(wèn)題。近年來(lái)雖配備了食品快速檢測(cè)設(shè)備，但檢測(cè)的食品種類有限，不能適應(yīng)食品市場(chǎng)監(jiān)管要求，嚴(yán)重影響了工作效率。通過(guò)查閱相關(guān)文獻(xiàn)和資料，我們提出一些解決食品安全問(wèn)題的對(duì)策和建議，希望可以為有關(guān)主管部門對(duì)食品的抽檢和廣大市民對(duì)食品安全問(wèn)題的了解以及抵制存在安全問(wèn)題的食品提供一些參考和幫助。針對(duì)食品安全存在的問(wèn)題，我們應(yīng)做好以下幾方面工作：1、加強(qiáng)宣傳教育，提高全民素質(zhì)；2、建立健全和完善食品安全管理

26、法規(guī)體系和標(biāo)準(zhǔn)體系；3、加大監(jiān)督力度，堅(jiān)決打擊制假售假等違法行為。對(duì)一些嚴(yán)重違法違規(guī)生產(chǎn)、經(jīng)營(yíng)不合格食品的企業(yè)及責(zé)任人要從嚴(yán)從重予以查處，形成一種高壓態(tài)勢(shì)和震懾力量；4、加強(qiáng)自律意識(shí)，充分發(fā)揮行為協(xié)會(huì)的作用；5、加強(qiáng)部門協(xié)調(diào)配合，共筑食品安全防線。食品安全是關(guān)系到千家萬(wàn)戶身體健康的基本民生問(wèn)題，涉及多個(gè)監(jiān)管部門，任務(wù)繁重，需用各部門配合協(xié)作，齊抓共管，為此，要樹立部門協(xié)作意識(shí)，加強(qiáng)與各有關(guān)部門的溝通協(xié)調(diào)，通過(guò)制度與機(jī)制上的創(chuàng)新，達(dá)到長(zhǎng)效監(jiān)管的目標(biāo)，共同推動(dòng)工作的開展。因?yàn)橹挥羞@樣，才能真正有效的避免食品安全問(wèn)題對(duì)社會(huì)造成嚴(yán)重危害。九 參考文獻(xiàn)【1】劉衛(wèi)國(guó) MATLAB程序設(shè)計(jì)教程，中國(guó)水利水電

27、出版社，2006年3月第三版【2】韓中庚 數(shù)學(xué)建模方法及其應(yīng)用，高等教育出版社，2005年6月第一版【3】食品安全體系的抽樣理論研究，李兵,陳國(guó)華 ,楊滌塵 ,朱寧，湖南人文科技學(xué)院數(shù)學(xué)與應(yīng)用數(shù)學(xué)系,湖南婁底 417000;2.桂林電子科技大學(xué)數(shù)學(xué)與計(jì)算科學(xué)學(xué)院 ,廣西桂林 541004【4】居民出行調(diào)查抽樣率模型，石飛，陸建，王煒，葛宏偉 （東南大學(xué) 交通學(xué)院，江蘇 南京210096）【5】朱德通，優(yōu)化模型與實(shí)驗(yàn)，統(tǒng)計(jì)大學(xué)出版社，2003年4月第一版【6】【7】【8】【9】【10】附錄附件1：>> A=1 3 5;1/3 1 3;1/5 1/3 1A = 1.0000 3.000

28、0 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000>> x,y=eig(A)x = 0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304iy = 3.0385 0 0 0 -0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i>> w=x(:,1)/sum(x(:,1)w = 0.6370 0.25830.1047（=3.0385，權(quán)重向

29、量為w=（0.6370，0.2583，0.1047）=(3.0385-3)/(3-1)=0.01925RI=0.58=0.01925/0.58=0.033<0.1,表明A通過(guò)了一致性驗(yàn)證。>> B1=1 3 5;1/3 1 3;1/5 1/3 1B1= 1.0000 3.0000 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000>> x,y=eig(B1)x = 0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753

30、 - 0.1304i -0.0753 + 0.1304iy = 3.0385 0 0 0 -0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i>> w=x(:,1)/sum(x(:,1)w = 0.6370 0.25830.1047（1=3.0385，權(quán)重向量w1=（0.6370，0.2583，0.1047）>> B2=1 3 7;1/3 1 5;1/7 1/5 1B2 = 1.0000 3.0000 7.0000 0.3333 1.0000 5.0000 0.1429 0.2000 1.0000>> x,y=eig(B2)x

31、= 0.9140 0.9140 0.9140 0.3928 -0.1964 + 0.3402i -0.1964 - 0.3402i 0.1013 -0.0506 - 0.0877i -0.0506 + 0.0877iy = 3.0649 0 0 0 -0.0324 + 0.4448i 0 0 0 -0.0324 - 0.4448i>> w2=x(:,1)/sum(x(:,1)w2 = 0.6491 0.27900.0719（2=3.0649，權(quán)重向量w2=（0.6491，0.2790，0.0719）附件2： >> f=0.037*(9/44*0.037+2/11*0.0

32、38+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0813>> f=0.038*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7

33、/44*0.0422)+440*(0.01/1.64)2)f = 0.0835>> f=0.039*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0857>> f=0.04*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+

34、7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.0879>> f=0.041*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.09

35、01>> f=0.042*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/(9/44*0.0372+2/11*0.0382+3/22*0.0392+15/88*0.042+13/88*0.0412+7/44*0.0422)+440*(0.01/1.64)2)f = 0.092390，80，60，75，65，707.317，6.68，5.142，6.5925，5.8565，6.461=38.327n=8+7+5+7+6+6=39 假設(shè)t=1.5,p=60,Pmax=6000，Tmax=176,

36、=11附件3：c=-39; A=39;39;1; b=117.3;100;11; Aeq=; beq=; vlb=0;0;0; vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)x = 2.5641fval = -100.0000>> w=1 1 3 3 5 5 6 7 1 9 9;1 1 3 3 5 5 6 7 1 9 9;1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3;1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3;1/5 1/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5;1/5 1

37、/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5;1/6 1/6 1/2 1/2 5/6 5/6 1 7/6 1/6 3/2 3/2;1/7 1/7 3/7 3/7 5/7 5/7 6/7 1 1/7 9/7 9/7;1 1 3 3 5 5 6 7 1 9 9;1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1;1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1;w = 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0000 9.0000 9.000

38、0 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0000 9.0000 9.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 0.3333 3.0000 3.0000 0.3333 0.3333 1.0000 1.0000 1.6667 1.6667 2.0000 2.3333 0.3333 3.0000 3.0000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 0.2000 1.80

39、00 1.8000 0.2000 0.2000 0.6000 0.6000 1.0000 1.0000 1.2000 1.4000 0.2000 1.8000 1.8000 0.1667 0.1667 0.5000 0.5000 0.8333 0.8333 1.0000 1.1667 0.1667 1.5000 1.5000 0.1429 0.1429 0.4286 0.4286 0.7143 0.7143 0.8571 1.0000 0.1429 1.2857 1.2857 1.0000 1.0000 3.0000 3.0000 5.0000 5.0000 6.0000 7.0000 1.0

40、000 9.0000 9.0000 0.1111 0.1111 0.3333 0.3333 0.5556 0.5556 0.6667 0.7778 0.1111 1.0000 1.0000 0.1111 0.1111 0.3333 0.3333 0.5556 0.5556 0.6667 0.7778 0.1111 1.0000 1.0000>> x,y=eig(w)x = Columns 1 through 7 -0.9883 -0.5443 0.7667 0.7667 -0.9887 0.9780 0.9841 0.0988 -0.5443 -0.2898 - 0.4162i -

41、0.2898 + 0.4162i 0.1232 -0.0768 0.0071 0.0329 -0.1814 0.1120 + 0.0340i 0.1120 - 0.0340i 0.0352 -0.0552 -0.0970 0.0329 -0.1814 0.1120 + 0.0340i 0.1120 - 0.0340i 0.0352 -0.0552 -0.0970 0.0198 -0.1089 -0.1255 + 0.0472i -0.1255 - 0.0472i -0.0104 -0.0576 -0.0638 0.0198 -0.1089 -0.1255 + 0.0472i -0.1255 -

42、 0.0472i -0.0104 -0.0576 -0.0638 0.0165 -0.0907 0.0560 + 0.0170i 0.0560 - 0.0170i 0.0176 -0.0276 -0.0485 0.0141 -0.0778 0.0471 + 0.0559i 0.0471 - 0.0559i 0.0203 -0.0342 0.0111 0.0988 -0.5443 -0.2684 - 0.0868i -0.2684 + 0.0868i -0.0443 0.1501 0.0331 0.0110 -0.0605 -0.0162 - 0.0370i -0.0162 + 0.0370i 0.0309 0.0145 0.0227 0.0110 -0.0605 -0.0162 -