【正文】 季節(jié)因素對(duì)食品質(zhì)量的影響。城市食品的來(lái)源越來(lái)越廣泛，人們消費(fèi)加工好的食品的比例也越來(lái)越高，因此除食材的生產(chǎn)收獲外，食品的運(yùn)輸、加工、包裝、貯存、銷售以及餐飲等每一個(gè)環(huán)節(jié)皆可能影響食品的質(zhì)量與安全。本文結(jié)合實(shí)際 ， 通過(guò)建立數(shù)學(xué)模型來(lái)研究如何進(jìn)行食品安全的抽檢等問(wèn)題 。 最后 對(duì)所得結(jié)果進(jìn)行總結(jié)，預(yù)測(cè)銷售地區(qū)與食品安全之間的關(guān)系，通過(guò)觀察數(shù)據(jù)得出結(jié)論，并通過(guò)觀察圖形來(lái)預(yù)測(cè)未來(lái)的發(fā)展趨勢(shì)和檢測(cè)重點(diǎn)。我們首先對(duì) 蔬菜，肉類，水產(chǎn)品的 數(shù)據(jù)進(jìn)行處理，利用 MATLAB 對(duì)殘差向量進(jìn)行分析，剔除其中的 異常點(diǎn) 。 8） w : 表示準(zhǔn)則層對(duì)目標(biāo)層的權(quán)重向量 wi: 表示準(zhǔn)則層對(duì)目標(biāo)層的權(quán)重向量 RI: 隨機(jī)一致性指標(biāo) CI: 判斷矩陣 一致性指標(biāo) CR: 一致性比率 n : 類別個(gè)數(shù) (n=8) xi : 影響 食品安全 的 銷售地區(qū) （ i=1,2， 通常情況下，由實(shí)際得到的判斷矩陣不一定是一致的，即不一定滿足傳遞性和一致性．實(shí)際中，也 不必要求一致性絕對(duì)成立，但要求大體上是一致的，即不一致的程度應(yīng)在容許的范圍內(nèi)． 判斷矩陣一致性指標(biāo) CI(Consistency Index)： 5 1nCI n??? ? 一致性指標(biāo) CI 的值越大，表明判斷矩陣偏離完全一致性的程度越大 ， CI 的值越小，表明判斷矩陣越接近于完全一致性。 2對(duì)于蔬菜及其制品、水產(chǎn)品、肉制品三大類產(chǎn)品來(lái)說(shuō)：我們通過(guò)求得三因素對(duì)于這三大類的權(quán)重，并比較這些權(quán)重，我們得出的結(jié)論是重金屬對(duì)于 這三類產(chǎn)品影響最小，而添加劑對(duì)其影響最大。在當(dāng)今這個(gè)時(shí)代，絕大部分人都喜好吃即食食品，而該類食品的質(zhì)量安全又是岌岌可危的。通過(guò)兩者的比較，我們可以看出哪類商品可 能會(huì)造成三因素所占權(quán) 10 重的不穩(wěn)定。由于深圳市有 7個(gè)區(qū)，因此，將所有的銷售地點(diǎn)以區(qū)來(lái)劃分從而分析食品安全與銷售地點(diǎn)之間的關(guān)系。 在水產(chǎn)品的討論中，數(shù)據(jù)的局限性導(dǎo)致模型無(wú)法建立，無(wú)法判斷水產(chǎn)品的安全問(wèn)題與銷售地區(qū)之間的線性關(guān)系。因此，我們大膽猜測(cè)，蔬菜是由于氣候的原因而導(dǎo)致的不合格。③對(duì)于一些污染性比較嚴(yán)重的地區(qū)，不管是產(chǎn)地或者是銷地，適當(dāng)增加抽檢頻率，讓其危害降到最低。 七、模型的改進(jìn)與推廣 模型的改進(jìn) 本文在建立層次分析法的 判斷矩陣時(shí)， 由于數(shù)據(jù)的分散性和匱乏 ，導(dǎo)致了在各因素的相對(duì)重要性的判斷上的主觀性太強(qiáng)，因此對(duì)此缺點(diǎn)作如下改進(jìn)：在收集和整理 大量數(shù)據(jù)的基礎(chǔ)上建立一套完整的重要性評(píng)估體系，以便于進(jìn)行科學(xué)的、準(zhǔn)確的判斷。 2 1 1 1/2 1 1/4 1/2 1/5。 [x,y]=eig(A2) w4=x(:,1)/sum(x(:,1)) w4 = 30 (由 MATLAB 可得到 λ 4=，權(quán)重向量 w4=（ ))； 2020 年： B2=[1 1/4 1/4。5 1 2。3 3 1]。 [x,y]=eig(B7) w20=x(:,1)/sum(x(:,1)) w20 = 40 (由 MATLAB 可得到 λ 20=，權(quán)重向量 w20=（ ))； 2020 年： C7=[1 1/3 1/3。4 1/2 1]。2 1 5。 [x,y]=eig(B3) w8=x(:,1)/sum(x(:,1)) w8= (由 MATLAB 可得到 λ 8=，權(quán)重向量 w8=（ ))； 2020 年： C3=[1 1/5 1/4。4 5 1]。 8) (Ai， Bi， Ci分別表示 2020年， 2020 年， 2020年各分類中重金屬，微生物，添加劑的比重） 2020 年： 蔬菜 水產(chǎn) 品 肉制 品 糧食制品 飲品 休閑食品 調(diào)味 品 第八 類 w 蔬菜 1 1/2 1/2 1/3 1/2 1/5 1/3 1/6 水產(chǎn)品 2 1 1 1/2 1 1/4 1/2 1/5 肉制品 2 1 1 1/2 1 1/4 1/2 1/5 糧食制品 3 2 2 1 2 1/3 1 1/4 飲品 2 1 1 1/2 1 1/4 1/2 1/5 休閑食品 5 4 4 3 4 1 3 1/2 調(diào)味品 3 2 2 1 2 1/3 1 1/4 第八類 6 5 5 4 5 2 4 1 蔬菜 重金 屬 微生 物 添加 劑 w1 重金屬 1 1/2 1/5 微生物 2 1 1/4 添加劑 5 4 1 水產(chǎn) 品 重金屬 微生物 添加劑 w4 重金屬 1 1/4 1/3 微生物 4 1 2 添加劑 3 1/2 1 23 2020 年： 肉制 品 重金屬 微生物 添加劑 w7 重金屬 1 1/5 1/3 微生物 5 1 2 添加劑 3 1/2 1 糧食制品 重金屬 微生物 添加劑 w10 重金屬 1 1/5 1/2 微生物 5 1 2 添加劑 2 1/2 1 飲品 重金 屬 微生 物 添加 劑 w13 重金屬 1 1/5 1/2 微生物 5 1 3 添加劑 2 1/3 1 休閑食品 重金屬 微生物 添加劑 w16 重金屬 1 1/5 1/4 微生物 5 1 2 添加劑 4 1/2 1 調(diào)味品 重金屬 微生物 添加劑 w19 重金屬 1 1/5 1/4 微生物 5 1 2 添加劑 4 1/2 1 第八類 重金屬 微生物 添加劑 w22 重金屬 1 1/4 1/3 微生物 4 1 2 添加劑 3 1/2 1 蔬菜 重金 屬 微生 物 添加 劑 w2 重金屬 1 2 1/3 微生物 1/2 1 1/5 添加劑 3 5 1 水產(chǎn)品 重金屬 微生物 添加劑 w5 重金屬 1 1/4 1/4 微生物 4 1 1/2 添加劑 4 2 1 24 2020 年： 肉制 品 重金屬 微生物 添加劑 w8 重金屬 1 1/5 1/4 微生物 5 1 2 添加劑 4 1/2 1 糧食制品 重金屬 微生物 添加劑 w11 重金屬 1 5 2 微生物 1/5 1 1/2 添加劑 1/2 2 1 飲品 重金 屬 微生 物 添加 劑 w14 重金屬 1 1/2 1/2 微生物 2 1 5 添加劑 2 1/5 1 休閑食品 重金屬 微生物 添加劑 w17 重金屬 1 1/4 1/3 微生物 4 1 2 添加劑 3 1/2 1 調(diào)味品 重金屬 微生物 添加劑 w20 重金屬 1 1/4 1/3 微生物 4 1 1/2 添加劑 3 2 1 第八類 重金屬 微生物 添加劑 w23 重金屬 1 1/3 1/2 微生物 3 1 1/2 添加劑 2 2 1 蔬菜 重金 屬 微生 物 添加 劑 w3 重金屬 1 2 1/4 微生物 1/2 1 1/5 添加劑 4 5 1 水產(chǎn)品 重金屬 微生物 添加劑 w6 重金屬 1 1/3 1/4 微生物 3 1 1/2 添加劑 4 2 1 25 權(quán)重 MATLAB 源代碼及其運(yùn)行結(jié)果： A= [1 1/2 1/2 1/3 1/2 1/5 1/3 1/6。較符合實(shí)際情況，模型可靠，并且模型相對(duì)簡(jiǎn)單，利于操作；該方法不僅適用于本題，也適用于其他方面的數(shù)據(jù)預(yù)測(cè)，有實(shí)際背景，可運(yùn)用于實(shí)踐，具有廣泛適用性。 2020 年全國(guó)爆發(fā)的“瘦肉精” 事件，導(dǎo)致那年的肉制品抽檢頻率急速上升，從而肉制品的合格率也呈指數(shù)關(guān)系飛速上漲，因此對(duì)于現(xiàn)在肉制品的抽檢結(jié)果來(lái)說(shuō)，其不合格率相比之前有所下降。從所建立的模型明顯看到 7,8,9,10 月份合格率較低，符合深圳市的氣候。用最小二乘的方法給出了變量間相關(guān)關(guān)系的回歸方程，針對(duì)各因素對(duì)食品安全的影響我們與偏回歸平方和聯(lián)系起來(lái)，并將各因素的影響程度進(jìn)行了排序。設(shè) 01 8, , ,b b b??? 分別是參數(shù) 0? ， 1? ，?，7? 的最小二乘估計(jì)，則回歸方程為 7722110 xbxbxbby ????? ? （ 4） 由最小二乘法知道， 710 , bbb ? 應(yīng)使得全部觀察值 ?y 與回歸值 ?y? 的偏差平方和Q 達(dá)到最小，即使 ^ 2()Q y y???? ? ?? 最 小 (5) 所以Ｑ是 710 , bbb ? 的非負(fù)二次式，最小值一定存在。 通過(guò)三種因素的變化趨勢(shì)分析，我們可以很大膽的推測(cè)出 2020 年這三大因素對(duì)食品安全的影響權(quán)重以及變化情況，前 提如果相關(guān)部門沒(méi)有進(jìn)行有效的方案控制。 英格蘭和威爾士每年有 300人死于由微生物引起的胃腸道感染性疾病， 35000人因此住院治療。對(duì)任意兩個(gè)因素 Ci和 Cj，用 aij表示 Ci和 Cj對(duì) O的影響程度之比，按 1— 9的比例標(biāo)度來(lái)度量 aij(i,j＝ 1， 2，?， n)。 7） S2 : 平方和 S 回 : 回歸方程和 β i : 待估計(jì)系數(shù)參數(shù) （ i=1,2， 假設(shè)每年的相同季節(jié)里氣候狀況接近相同 。我們不僅可以分析出其全部產(chǎn)品中三因素的變化趨勢(shì)，還可以得到部分產(chǎn)品質(zhì)量安全的影響因素的變化趨勢(shì)，同時(shí)預(yù)測(cè)出未來(lái)一年的食品安全情況。本文用最小二乘的方法給出了變量間相關(guān)關(guān)系的回歸方程，針對(duì)各因素對(duì) 食品安全 的影響我們與偏回歸平方和聯(lián)系起來(lái)，并將各因素的影響程度進(jìn) 行了排序。另一方面，我們著重分析了蔬菜及其制品、水產(chǎn)品、肉制品 三大類質(zhì)量安全的影響因素及其所占權(quán)重。我們通過(guò)數(shù)據(jù)的整理，發(fā)現(xiàn)一些非可食性的產(chǎn)品占了很大的一部分比重，為此我們將其歸為第八類，其中主要包括保鮮膜，一次性餐具，紙巾等產(chǎn)品。 假設(shè)導(dǎo)致食品不合格的因素只包括微生物，重金屬，添加劑含量三者，其他的因素可忽略不計(jì)。 7） Qi: 變量 xj的偏回歸平方和 （ i=1,2，比較時(shí)采用相對(duì)尺度標(biāo)準(zhǔn)度量，盡可能地避免不同性質(zhì)的因素之間相互比較的困難。 對(duì)于本題我們建立了如下的層次分析法模型：目標(biāo)層為 食品質(zhì)量安全；準(zhǔn)則層為八大類產(chǎn)品，分別為： 蔬菜類及其制品（包括一些菌類），水產(chǎn)品及其制品，肉制品，糧食制品，飲品（包括茶），休閑食品，調(diào)味品，第八類；方案層為三類影響食品質(zhì)量安全的因素：微生物，重金屬，添加劑含量。在微生物污染中，細(xì)菌性污染是涉及面最廣、影響最大、問(wèn)題最多的一種污染。 由于我們所考慮的食品中大部分產(chǎn)品屬于加工性的產(chǎn)品，而非初級(jí)產(chǎn)品，因此在食 品保鮮上許多食品都加入了大量的添加劑使之保質(zhì)期延長(zhǎng)。( 721 ???? xxxy α =1,2，? ,7 （ 1） 那么這一組數(shù)據(jù)可以假設(shè)有如下的結(jié)構(gòu)式： ?????????????????????????,. . .. . . . .. . . . . . . . . .. . . . . . . . . .. . . . . . . . . .. . . . . . . . . .. . . . . . . . . .. . . . . . . . . .,. . .,. . .777772271107227722221