【正文】 disp(39。)。CI=39。t=m(i)。m=zeros(1,100)。6． 模型的評(píng)價(jià)與優(yōu)化優(yōu)點(diǎn)： 1 .該模型運(yùn)層次分析方法，模糊數(shù)學(xué)理論，內(nèi)容翔實(shí)，結(jié)構(gòu)簡(jiǎn)單明了，使得各指標(biāo)的相對(duì)重要性一目了然，同時(shí)按照分解、比較、判斷、綜合的思維方式進(jìn)行決策，而每一層的權(quán)重設(shè)置最后都會(huì)直接或間接影響到結(jié)果，而且在每個(gè)層次中的每個(gè)因素對(duì)結(jié)果的影響程度都是量化的，非常清晰、明確。6）、依據(jù)在S（a）中取最大的原理。我們現(xiàn)在根據(jù)各年的數(shù)據(jù)列的農(nóng)藥和化肥隸屬度函數(shù)圖形如下： 通過(guò)處理各指標(biāo)的測(cè)度數(shù)據(jù)，可以得到每年各指標(biāo)與上一級(jí)指標(biāo)間的模糊關(guān)系矩陣R。 符號(hào)說(shuō)明a食品安全指標(biāo)a1食品原料風(fēng)險(xiǎn)a2農(nóng)藥殘留風(fēng)險(xiǎn)a3經(jīng)營(yíng)消費(fèi)風(fēng)險(xiǎn)a11生產(chǎn)加工風(fēng)險(xiǎn)a12種植過(guò)程風(fēng)險(xiǎn)a131農(nóng)藥使用量a132化肥使用量a21加工環(huán)境安全a22輔料安全度a23加工食物合格率a24加工工具安全a31場(chǎng)所衛(wèi)生度a32設(shè)備餐具安全a33工作人員健康安全a34消費(fèi)食品合格a211環(huán)境衛(wèi)生a212工作人員健康安全a221食物防腐調(diào)味風(fēng)險(xiǎn)a222加工用水安全風(fēng)險(xiǎn)aij數(shù)量化的相對(duì)權(quán)重w權(quán)重向量CI一致性檢驗(yàn)指標(biāo) RI平均隨機(jī)指標(biāo)CR 一致性比例 最大特征值R模糊關(guān)系矩陣 S安全指數(shù) 移動(dòng)平均n移動(dòng)平均期數(shù) 模型的建立與求解食品安全指標(biāo)a食品原料風(fēng)險(xiǎn)
a1生產(chǎn)加工風(fēng)險(xiǎn)
a2經(jīng)營(yíng)消費(fèi)風(fēng)險(xiǎn)
a3農(nóng)藥殘留風(fēng)險(xiǎn)
a11種植過(guò)程風(fēng)險(xiǎn)a12加工環(huán)境安全a21輔料安全度
a22加工食物合格率a23加工工具安全
a24場(chǎng)所衛(wèi)生度
a31設(shè)備餐具安全a32工作人員健康安全a33消費(fèi)食品合格
a34農(nóng)藥使用量a131化肥使用量a132環(huán)境衛(wèi)生a211工作人員健康安全a212工作人員健康安全a212加工用水安全風(fēng)險(xiǎn)a222 目前我國(guó)還沒(méi)有成熟的食品安全風(fēng)險(xiǎn)體系，相關(guān)的統(tǒng)計(jì)信息還不完善。為了能客觀、定量、通俗、概括性地反映某地區(qū)的食品安全現(xiàn)狀，及時(shí)發(fā)現(xiàn)食品安全領(lǐng)域中存在的問(wèn)題，幫助居民合理選擇安全衛(wèi)生的食品，保障公眾安全健康，有關(guān)部門(mén)計(jì)劃定期發(fā)布“食品安全指數(shù)”，使它成為和消費(fèi)者物價(jià)指數(shù)（CPI ），空氣污染指數(shù)那樣有較大影響和指導(dǎo)意義的“民生指數(shù)”。如有違反競(jìng)賽規(guī)則的行為，我們將受到嚴(yán)肅處理。 對(duì)于問(wèn)題二，我們利用線性回歸的方法及matlab編程作出由問(wèn)題一得出的S(a)即安全指數(shù)隨年數(shù)的變化圖形，結(jié)合移動(dòng)平均法來(lái)預(yù)測(cè)未來(lái)幾年的變化趨勢(shì)，并算出2013年的食品安全指數(shù)。以近幾年國(guó)家統(tǒng)計(jì)局發(fā)布的相關(guān)統(tǒng)計(jì)信息為依據(jù)，計(jì)算出模糊關(guān)系矩陣及其相應(yīng)的安全指數(shù)。對(duì) 9, , 1 L = n ，Saaty給出了RI 的值，如表2 所示。 表7 食品安全指數(shù)與食品安全風(fēng)險(xiǎn)對(duì)應(yīng)關(guān)系 風(fēng)險(xiǎn)等級(jí)高風(fēng)險(xiǎn) 中等風(fēng)險(xiǎn) 低風(fēng)險(xiǎn)無(wú)風(fēng)險(xiǎn)對(duì)應(yīng)安全指數(shù) ~~ 結(jié)果分析 根據(jù)最后結(jié)果可得出如下結(jié)論： 1）、S（a）的值穩(wěn)定在中等風(fēng)險(xiǎn)范圍，從食品供應(yīng)鏈整體范圍看，我國(guó)食品安全狀態(tài)處于中等風(fēng)險(xiǎn)范圍，整體水平還有待提高。 對(duì)于問(wèn)題二，我們可以先使政府部門(mén)充分意識(shí)到時(shí)食品安全指數(shù)是怎樣產(chǎn)生的，它與哪些環(huán)節(jié)有關(guān)，從而為了提高該指數(shù)而采取一些相應(yīng)措施，如上述模型中提到的減少農(nóng)藥和化肥的施用，然后通過(guò)政府權(quán)威的力量對(duì)食品安全指數(shù)進(jìn)行宣傳，讓人們也意識(shí)到目前我們國(guó)家食品安全的狀況，讓人們對(duì)那些指數(shù)偏低的方面產(chǎn)生危機(jī)感，從而在政府的幫助下，一起監(jiān)督維護(hù)食品安全。)。while kp i=i+1。 %以下是一致性檢驗(yàn)CI=(tn)/(n1)。)。)。 0]s = 。B=input(39。else disp(39。)。i=2。九、參考文獻(xiàn) [1] MATLAB層次分析法教程 百度文庫(kù) 百度百科 食品論壇[2]（王琬蕪/張浩森/馬士良/何美珍/）的食品安全指數(shù)(FSI)數(shù)學(xué)建模與評(píng)價(jià)[3] [4] 參考農(nóng)業(yè)部統(tǒng)計(jì)公報(bào)、農(nóng)業(yè)統(tǒng)計(jì)年鑒、衛(wèi)生統(tǒng)計(jì)2005～2012年的統(tǒng)計(jì)數(shù)據(jù)[5] 參考李國(guó)平老師發(fā)的PPT課件中回歸分析[6] [7] [8][9][10] 附錄層次分析模型程序disp(39。移動(dòng)平均預(yù)測(cè)。根據(jù)圖2和圖3農(nóng)藥和化肥施用風(fēng)險(xiǎn)隸屬度函數(shù)，可以讀出20052012年每年農(nóng)藥和化肥施用風(fēng)險(xiǎn)的隸屬度，由隸屬度構(gòu)成模糊關(guān)系矩陣R12。