【正文】 } 附錄 3 深圳食品安全三年的評(píng)價(jià) 2021 年 △ i 0. 00444444 0. 00119586 0. 00259151 0. 00143554 求判斷矩陣 *************************************** 1 1. 28571 0. 12963 1. 28571 0. 777778 1 0. 111111 1 7. 71429 9 1 9 0. 777778 1 0. 111111 1 *************************************** 1 1. 95652 0. 150327 2. 34783 0. 511111 1 0. 116162 1 6. 65217 8. 6087 1 9 0. 425926 1 0. 111111 1 *************************************** 1 3. 375 0. 177778 2. 25 0. 296296 1 0. 111111 0. 888889 5. 625 9 1 7. 875 0. 444444 1. 125 0. 126984 1 *************************************** 1 7. 97143 0. 972222 7. 45714 0. 125448 1 0. 111111 1 1. 02857 9 1 8. 48571 0. 1341 1 0. 117845 1 VA 0. 25 0. 25 0. 25 0. 25 VB 0. 100641 0. 0801792 0. 739 0. 0801792 26 0. 139288 0. 0745828 0. 714942 0. 0 0. 170342 0. 0654139 0. 686532 0. 0 0. 428627 0. 0534368 0. 462789 0. 0 VBW 0. 402565 0. 320717 2. 956 0. 320717 0. 557152 0. 298331 2. 85977 0. 284747 0. 681366 0. 261656 2. 74613 0. 310849 1. 71451 0. 213747 1. 85116 0. 220588 PI=0. 0127681 2021 年 △ i 0. 00261041 0. 00267164 0. 00259233 0. 00160037 求判斷矩陣 *************************************** 1 2. 32653 4. 71428 9 0. 429825 1 2. 38775 6. 67347 0. 212121 0. 418803 1 4. 28571 0. 111111 0. 149847 0. 233333 1 *************************************** 1 4. 06849 0. 202778 1. 60274 0. 245791 1 0. 111111 0. 405556 4. 93151 9 1 6. 53425 0. 623932 2. 46575 0. 15304 1 *************************************** 1 9 1. 59231 5. 95385 0. 111111 1 0. 134995 0. 328283 0. 628019 7. 40769 1 4. 36154 0. 167959 3. 04615 0. 229277 1 *************************************** 1 8. 38842 1 9 0. 119212 1 0. 129707 。j++) { sum+=va[i]*vb[i][j]*a[j][i]。i++) for(j=0。 for(i=0。 } void last() { float sum=0。\t39。j4。i++) { 25 coutendl。 for(i=0。 } coutendl。j++) coutvb[i][j]39。 for(j=0。i4。 } ///////////////////////////////////////////////////// //out coutendlendlVB。j4。j++) sum+=vbw[i][j]。 for(j=0。i4。k++) vbw[i][j]+=t[j][k]。j++) for(k=0。 } for(j=0。k4。k++) sum+=e[i][k][j]。 for(k=0。j4。i4。 float sum。 } void forth() { float t[4][4]。\t39。i4。 ///////////////////////////////////////////////////// //out coutendlendlVAendl。i4。i++) sum+=va[i]。 for(i=0。j++) va[i]+=t[i][j]。i++) for(j=0。 } for(i=0。j4。j++) sum+=aa[j][i]。 for(j=0。i4。 float sum。 } } } void third() { float t[4][4]。\t39。k4。j4。i++) { cout***************************************endl。 for(i=0。 continue。 continue。amp。k++) { t=e[i][j][k]。j++) for(k=0。i++) for(j=0。 for(i=0。k++) e[i][j][k]=(a[i][j]a[i][k])/d[i]。j++) for(k=0。i++) { for(j=0。 for(i=0。 //求出 △ i //////////////////////////////////////////// } void second() { int i, j, k。i++) coutd[i] 。 23 for(i=0。 } } d[i]=(maxmin)/9。 } if(a[i][j]max) { max=a[i][j]。j++) { if(a[i][j]min) { min=a[i][j]。 for(j=1。i++) { min=a[i][0]。 for(i=0。 } void first() { int i, j。 forth()。 second()。 //求 VB, VBW void last()。 void third()。 //C. R. void first()。 //B的最大化特征向量 float ci[4]={0}。 //A歸一特征向量 float vb[4][4]={0}。 //△ i float e[4][4][4]={0}。*///2021年 /*float a[N][M]={0. 001445087, 0, 0. 00433526, 0. 00433526, 0. 017286915, 0, 0. 020218067, 0. 002881152, 0. 005790917, 0. 004571777, 0. 011581835, 0, 0. 010071415, 0. 004028566, 0. 024903864, 0}。 //在 2021年用此數(shù)據(jù) /*float aa[4][4]={1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1}。 %指定模型的結(jié)構(gòu) [coeffX,errorsX,LLFX] = garchfit(spec2,da) %擬合參數(shù) [sigmaForecast,w_Forecast] = garchpred(coeffX,da,4) %計(jì)算 4步預(yù)報(bào)值 x_pred=a(end)+cumsum(w_Forecast) %計(jì)算原始數(shù)據(jù)的 4 步預(yù)測(cè)值 附錄 2 食品安全的評(píng)價(jià)程序 includeiostream. h define N 4 //季度 define M 4 //主要的食品領(lǐng)域 float aa[4][4]={1, 1, 1, 0. 5, 1, 1, 1, 0. 5, 1, 1, 1, 0. 5, 2, 2, 2, 1}。off39。Display39。M39。R39。)。m=input(39。輸入階數(shù) R＝ 39。,i,j,aic,bic)。 fprintf(39。 %擬合參數(shù) num=garchcount(coeffX)。)。,39。,j,39。,i,39。 %計(jì)算差分后的數(shù)據(jù)個(gè)數(shù) for i=0:3 for j=0:3 spec= garchset(39。 %按照原來數(shù)據(jù)的順序去掉零元素 r11=autocorr(a) %計(jì)算自相關(guān)函數(shù) r12=parcorr(a) %計(jì)算偏相關(guān)函數(shù) da=diff(a)。%第一列為微生物的不合格率，二列為重金屬，三列為添加劑，四 列為其它 a=nonzeros(a39。 所聚類的 結(jié)果見附 錄 5. 最終統(tǒng)計(jì)為表 14. 表 14 食品種類數(shù)據(jù)統(tǒng)計(jì) 類別 合格總數(shù) 不合格總數(shù) 合格率 不合格率 1 744 144 0. 837838 0. 162162 2 6640 57 0. 991489 0. 008511 3 6917 797 0. 896681 0. 103319 4 2889 224 0. 928044 0. 071956 5 1052 157 0. 870141 0. 129859 6 10011 293 0. 971564 0. 028436 再由 卡方統(tǒng)計(jì)量判別 , 結(jié)果見 表 15. 表 15 食品種類與食品合格率之間的卡方統(tǒng)計(jì)量 ?? 1 2 3 4 5 6 1 0 717. 5154 28. 3018 67. 3158 4. 3455 389. 6461 2 0 577. 9853 307. 4469 572. 6331 79. 9131 3 0 25. 5426 7. 7106 435. 2385 4 0 36. 3214 122. 2230 5 0 296. 3570 6 0 由表 15 知 食品種類的變化是很顯著 , 說明食品種類的不同決定著有不同程度的不合格率 . 且由表 14 可知第一類 (即豆制品 )的不合格率是最高的 , 為 0. 162162；第二類的不合格率是最低的 , 為 0. 008511. 由高到低依次為第一類