食品安全问题数学建模论文

葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅芀蒅薄螄羀芇蒀螄肂蒃莆袃膅芆蚄袂袄蒁薀袁羇芄蒆袀腿葿蒂衿芁莂螁袈羁膅蚇袈肃莁薃袇膆膃葿羆袅荿莅羅羈膂蚃羄肀莇蕿羃节膀薅羂羂蒅蒁羂肄芈螀羁膆蒄蚆羀艿芇薂聿羈蒂蒈蚅肁芅莄蚅膃蒀蚃蚄羃芃虿蚃肅蕿薅蚂膇莁蒁蚁芀膄蝿蚀罿莀蚅蚀肂膃薁蝿膄莈蒇螈袄膁莃螇肆莆螂螆膈艿蚈螅 食品安全模型承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们授权数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. xd 2. pjp 3. lck 指导教师或指导教师组负责人 (打印并签名)： 日期： 2013 年 8 月 7日赛区评阅编号（由赛区组委会评阅前进行编号）：编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：食品安全指数摘要 食品安全问题近年来渐受全社会关注。提高食品安全程度，让人民群众吃得放心，已成为当前主要的民生问题之一。本文研究了食品安全指数的建立及其深度利用方法。 对于问题一，我们针对我国食品产业链的现状，将食品供应链划分为供应源头、食品加工和经营消费三个环节，建立了“从生产到消费”的评价指标体系。然后，用层次分析法计算出各同级指标之间的权重，并通过一致性指标进行验证。接着，应用模糊数学的理论处理 20052012年各项指标的数据，计算出各同级指标与其上一级指标之间的模糊矩阵，根据各个指标间的权重，计算出各指标的安全指数，以此对食品安全问题做定量评估。然后通过各种媒体进行发布宣传。 对于问题二，我们利用线性回归的方法及matlab编程作出由问题一得出的S(a)即安全指数随年数的变化图形，结合移动平均法来预测未来几年的变化趋势，并算出2013年的食品安全指数。最后查阅相关资料，得到关于政府部门如何提高食品安全的具体措施。 对于问题三，我们简单介绍了食品安全指数，向群众普及了这一概念，并叙述了如何运用该指数提高食品安全。关键词：食品安全指数 层次分析法 模糊数学理论 线性回归 移动平均法 未来预测20 一、问题重述 食品安全问题近年来渐受全社会关注。“民以食为天”，食品安全关系人民群众生活，关系社会稳定和谐。食品安全重大事件甚至可能损害国家形象、影响对外交往。提高食品安全程度，让人民群众吃得放心，已成为当前主要的民生问题之一。为了能客观、定量、通俗、概括性地反映某地区的食品安全现状，及时发现食品安全领域中存在的问题，帮助居民合理选择安全卫生的食品，保障公众安全健康，有关部门计划定期发布“食品安全指数”，使它成为和消费者物价指数（CPI ），空气污染指数那样有较大影响和指导意义的“民生指数”。作为此项工作的基础，有关部门常年对主要食品的质量和所含有毒有害物质进行日常性抽检，并可根据需要获取其他与食品安全有关的已有资料和数据，以及开展相关的调查工作。 1请你建立完整的食品安全指数数学模型，并确定信息发布形式，使这一指数能发挥预期的作用。 2 如何对食品安全指数进行深度利用，为政府部门加强和改善食品安全工作提供决策参考，提出改进食品安全的政策措施，指导居民改善膳食结构，加强对有毒有害物质的预警和食品安全重大事件的预防。 3 试用不超过两页的篇幅向公众介绍这一新生事物，以及如何利用这一指数控制食品风险、提高生活质量。 2、 问题分析 我们的生活中，存在着各种各样的食物。在提倡绿色食品、健康食品的今天，还是有一些食品存在严重安全隐患，而食品的产生具体通过原料，原料加工，食品销售这三个主要流程，从而食品的安全问题也是在其中产生，基于此，我们想到用层次分析法来进行逐次分层分析，得到这些方面对食品安全的权重（即对上层的影响值），计算出各要素相对于其上一层某要素的权重向量，然后根据搜集到的近几年的统计信息，通过运用模糊数学理论，用模糊矩阵来表示指标与其上一级指标之间的模糊关系。接着，评估往年食品安全状态。引入安全指数，并应用模糊综合运算进行结果评价。以近几年国家统计局发布的相关统计信息为依据，计算出模糊关系矩阵及其相应的安全指数。最后通过线性回归与移动平均法来进行未来预测，从而达到预警的效果。3、 问题假设假设一：我国的食品安全程度主要受模型中所研究的19种指标的影响。 假设二：模型中所研究的各同级指标之间相互独立，忽略他们之间的相互影响。 假设三：我国每年的食品生产量保持在一个稳定的水平。 假设四：历年的食品安全指数具有可预测性。 假设五：居民的饮食习惯在一段时间内不发生突变。4、 符号说明a食品安全指标a1食品原料风险a2农药残留风险a3经营消费风险a11生产加工风险a12种植过程风险a131农药使用量a132化肥使用量a21加工环境安全a22辅料安全度a23加工食物合格率a24加工工具安全a31场所卫生度a32设备餐具安全a33工作人员健康安全a34消费食品合格a211环境卫生a212工作人员健康安全a221食物防腐调味风险a222加工用水安全风险aij数量化的相对权重w权重向量CI一致性检验指标 RI平均随机指标CR 一致性比例 最大特征值R模糊关系矩阵 S安全指数 移动平均n移动平均期数 5、 模型的建立与求解5.1.建立基于食品产业链的评价指标体系食品安全指标a食品原料风险a1生产加工风险a2经营消费风险a3农药残留风险a11种植过程风险a12加工环境安全a21辅料安全度a22加工食物合格率a23加工工具安全a24场所卫生度a31设备餐具安全a32工作人员健康安全a33消费食品合格a34农药使用量a131化肥使用量a132环境卫生a211工作人员健康安全a212工作人员健康安全a212加工用水安全风险a222 目前我国还没有成熟的食品安全风险体系，相关的统计信息还不完善。为使食品安全风险评价指标具有可操作性，将食品供应链简化成三个环节，即：农产品供应源头、食品加工环节、经营消费环节，这里的“经营消费环节”包括除食品加工之外的其他食品经营服务环节，例如食品的配送、物流、销售、餐饮业以及食堂等集体用餐配送环节，因而体现了“从生产到消费”的食品供应链全过程特点。 根据系统性及层次性原则，构建基于整体监控的食品安全评价指标体系的层次结构模型，层次结构图如下：图1. 一级指标层二级指标层三级指标层 5.2 以层次分析法为基础的对各项指标的权重求取 a.采取对因子进行两两比较建立成对比较矩阵的办法。即每次取两个因子xi 和xj，以aij表示xi和xj对Z的影响大小之比，全部比较结果用矩阵A=(aij )n*n 表示，称A为Z X 之间的成对比较判断矩阵（简称判断矩阵）。容易看出，xj与xi对Z的影响之比应为aji =aij/1.若矩阵A =(aij )n n 满足 （i ）aij 0 ，（ii ）aji = aij/1 （i,j=1,2.n ） 则称之为正互反矩阵(易见aii =1，i=1.n ) 。 关于如何确定aij的值,Saaty 等建议引用数字 19 及其倒数作为标度。表 1 列出了19 标度的含义： 表1 标度的含义 1表示两个因素相比，具有相同重要性3表示两个因素相比，前者比后者稍重要5表示两个因素相比，前者比后者明显重要 7表示两个因素相比，前者比后者强烈重要9表示两个因素相比，前者比后者极端重要2，4，6，8 表示上述相邻判断的中间值 倒数若因素i 与因素j 的重要性之比为aij,那么因素j 与因素i 重要性之比为aij=aji/1。 b.计算权重向量及最大特征根 根据构造的判断矩阵用matlab编程计算出各层次对上层的权重向量及特征根。c.进行一致性检验（i ）计算一致性指标CI （ii ）查找相应的平均随机一致性指标 RI 。对 9, , 1 L = n ，Saaty给出了RI 的值，如表2 所示。 表2 RI 的值 n 1 2 3 4 5 6 7 8 9 RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 RI的值是这样得到的，用随机方法构造 500 个样本矩阵：随机地从 19 及其倒数中抽取数字构造正互反矩阵，求得最大特征根的平均值，并定义 。 （）计算一致性比例CR CR=当 10 . 0 CR 时，认为判断矩阵的一致性是可以接受的，否则应对判断矩阵作适当修正。 按照上面步骤，我们列得求得如下表格中的结果。 第一指标层（1）判断矩阵及权重向量和特征根（表3）AA1A2A3wA111/220.2973.01A22130.540A31/21/310.163 第二指标层（2）判断矩阵及权重向量和特征根（表4）A1A11A12wA1111/20.6672A12210.333A2A21A22A23A24wA2113/42/31/20.174.1596A224/31210.315A233/21/213/20.254A24212/310.261A3A31A32A33A34wA3111/23/24/50.2134.0679A32214/37/60.324A332/33/411/20.173A345/46/7210.29第三指标层（3）判断矩阵及权重向量和特征根（表5）A12A121A122 wA13A131A132WA121120.6672A131120.6672A1221/210.333A1321/210.333A21A211A212 wA22A221A222WA211120.6672A221110.52A2121/210.333A222110.5A24A241A242WA241110.52A242110.55.3 应用模糊数学的理论处理评价信息得到指标 根据模糊数学理论，用模糊隶属度来表示某一指标对其上一级指标而言的评价表现。 我们查出相关指标测度值，除了农药和化肥施用量要建立隶属度函数，其余都是在0-1之间的百分数，可以直接用作隶属度。然后根据图形中的隶属关系结合各层相应的权重，求出安全指数。我们现在根据各年的数据列的农药和化肥隶属度函数图形如下： 图2.农药施用量隶属度函数图3.化肥施用量隶属度函数通过处理各指标的测度数据，可以得到每年各指标与上一级指标间的模糊关系矩阵R。其中，农药和化肥的隶属度由图2、图3所示的函数确定，其余指标的隶属度由其测度数据直接表示。设各同级指标关于上一级指标的权重向量为W，则上一级指标的安全指数S=WR。 以此类推，通过搜集历年各基础指标的测度数据，可以逐级向上计算出各个指标的安全指数，从而得出每年的食品安全指数，以此来定量评估每年的食品安全状态。5.4 评估往年食品安全状态 我们从各大网页搜集到2005-2012年各个相应指标的测度值，列成表格如下：（表6）20052006200720082009201020112012农产品农（兽）药残留监测合格率a11（%）94.293.095.197.094.5.96.997.197.9农药使用量a121（万t）103.9129.6173.1190.2226.2234.2264.8278.9化肥使用量a122（万t）3414.83782.34187.14263.24863.74521.14179.05800.0加工企业卫生监督合格率a211（%）87.689.191.292.189.990.189.392.2持健康证明人数比例a212（%）97.097.397.698.298.097.297.697.5加工企业食品添加剂抽检合格率a221（%）87.590.393.394.394.593.994.394.0饮用水的抽检合格率a222（%）89.481.785.788.186.590.591.896.7加工业食品抽样合格率a23（%）89.490.491.292.191.391.691.893.3经营/餐饮业环境卫生合格率a31（%）88.887.391.193.892.293.393.092.8餐具消毒合格率a32（%）79.782.583.084.683.984.583.584.0经营/ 餐饮业持健康证明人数比例a33（%）89.896.997.998.998.097.598.197.7经营/ 餐饮业食品抽样合格率a34（%）86.919186.991.293.394.994.396.3加工工具设备（包装材料，容器，洗涤剂，消毒剂）的合格率a24(%)75.080.387.681.683.785.091.291.5 根据图2和图3以及表6的统计数据，可计算出模糊关系矩阵R和相应的安全指数S，逐级评价，最后可得出总目标层的综合安全指数S（a）。下面以S(a12)为例说明求解过程。 R12是表示三级指标A121农药施用风险（农药施用量）与a122化肥施用风险（化肥施用量）和它们上层的二级指标a12种植过程风险之间模糊关系的矩阵。根据图2和图3农药和化肥施用风险隶属度函数，可以读出2005-2012年每年农药和化肥施用风险的隶属度，由隶属度构成模糊关系矩阵R12。R(12)=则相应的安全指数为：S（a12）=wR12=(0.667 0.333)=(0.5729 0.4947 0.3747 0.3124 0.2178 0.2197 0.1713 0.2693)同理，可得：S（a1)=(0.667 0.333)=(0.8191 0.7850 0.7584 0.7510 0.7082 0.7195 0.7047 0.7427)运用同样的方法，可算得如下：S(a21)=(0.9073 0.9183 0.9333 0.9413 0.9260 0.9246 0.9206 0.9396)S(a22)=0.8845 0.8600 0.8950 0.9120 0.9050 0.9220 0.9305 0.9535）S(a2)=( 0.8557 0.8662 0.9009 0.8942 0.8929 0.9021 0.9208 0.9359)S(a3) =(0.8548 0.8848 0.8843 0.9095 0.9083 0.9164 0.9118 0.9181)S(a) =(0.8447 0.8451 0.8559 0.8542 0.8390 0.8502 0.8552 0.8756) 最后，按照表6的对应关系确定食品安全风险等级。 表7 食品安全指数与食品安全风险对应关系 风险等级高风险 中等风险 低风险无风险对应安全指数 0.95 结果分析 根据最后结果可得出如下结论： 1）、S（a）的值稳定在中等风险范围，从食品供应链整体范围看，我国食品安全状态处于中等风险范围，整体水平还有待提高。 2）、食品安全综合指数估算值逐年增加，说明我国食品安全状态呈现缓缓提高趋势， 这与我国食品安全的实际情况相符。 3)、从结果看，S(a12)的值很低，风险很大，要限制农药化肥的施用才行。S（a1）在高风险上下波动，说明风险值比较大，要调节兽药和激素的使用，降低风险。4）、S(a21) 、S(a22)、S(a2)的值较大。属于低风险范畴，说明我国还是比较注重食品加工安全，继续保持。5）、S（a3)表示的是经营消费环节的安全风险，其值在中等风险和低风险上下波动，有待提升。6）、依据在S（a）中取最大的原理，我国现在的安全指数为0.8756。5.4 利用回归分析和移动平均法预测未来食品安全状态 由上面建立的模型，我们得出了粗略的各年的食品安全指数，历史是用来借鉴并可以用来推测未来趋势发展，要想给人们达到一个预警，我们采取了以下方法来预测短期性的未来，首先我们利用线性回归画出S（a）s随时间推移而变化的函数曲线，如下：（相关程序设计在附录给出）图4从图形可以看出安全指数波动不大又逐年上升的。于是我们采用了移动平均法来预测估计接下来几年的安全指数值。移动平均法l 基本思想 假定在一个较短的时间间隔里，序列值之间的差异主要是由随机波动造成的。根据这种假定，我们可以用一定时间间隔内的平均值作为某一期的估计值。l 分类 n期中心移动平均 n期移动平均n期中心移动平均如5期中心移动平均 移动平均期数确定的原则l 事情的发展有无周期性 - 以周其长度作为移动平均的间隔长度，以消除周期效应 的影响。l 对趋势平滑的要求 - 移动平均的期数越多，拟合趋势越平滑对趋势为反应近期变化的敏感程度，要求移动的平均期数越少，拟合趋势越敏感。移动平均预测运用该方法我们求得在2013年的安全指数值为0.8535。问题解答： 对于问题一，前面我们已经建立了相应了模型来解决如何确定食品安全指数的问题，关于信息的发布于宣传，我们可以通过各个传播媒介如电视、网络、广播、报纸等来进行，让群众们认识到这一指数，相关部门也可利用此指数来提高人们的食品安全意识，增强食品安全的防范，从而让政府有关部门对食品生产各个流程的严厉审查。 对于问题二，我们可以先使政府部门充分意识到时食品安全指数是怎样产生的，它与哪些环节有关，从而为了提高该指数而采取一些相应措施，如上述模型中提到的减少农药和化肥的施用，然后通过政府权威的力量对食品安全指数进行宣传，让人们也意识到目前我们国家食品安全的状况，让人们对那些指数偏低的方面产生危机感，从而在政府的帮助下，一起监督维护食品安全。从模型中我们也提出了未来安全指数的预测，这可以给人们起到一个预警的作用，在宣传方面，要做到尽量利用权威媒体，这样人们才可以信服，加大宣传力度，使这一概念深入人心，从而逐渐为人们所接受。政府相关政策如下：（1） 建立分工明确、协调一致的食品安全管理体制；（2） 建立中国食品安全需要的科技支撑体系；（3） 完善中国食品安全标准体系；（4） 整合现有检验检测监测机构；（5） 制定食品安全基本法，增强执法能力；（6） 建立健全食品突发事件应急反应机系；（7） 建立统一、规范的食品认证认可体系；（8） 建立健全食品安全信息交流体系；（9） 完善食品供应组织体系；（10） 实行规范的市场准入制度，完善食品溯源管理；（11） 充分发挥行业协会、媒体的作用；（12） 加强教育的培训；（13） 加强国际交流与合作。 对于问题三,我们自撰了一篇文章如下： 长期以来，食品安全事故频繁发生，有关劣质食品的报道层出不穷，比如说关于劣质奶粉、劣质面粉、劣质大米、劣质豆制品、染白粉丝、注水肉、苏丹红，福寿螺、毒火腿、瘦肉精、毛发酱油、黑心月饼及2008年9月份曝光的“三鹿婴幼儿奶粉事件”，一起起怵目惊心，轻者致病，重者导致死亡。各级监管部门针对于此的执法检查，始终没有停止过，并在每年的元旦、春节、国庆、中秋等重大节日期间加大执法检查的力度，同时也查处了不少违法单位。但食品安全事件时有发生，严重威胁着人民的生命健康，时时令我们提心吊胆。 基于此，我们引入了食品安全指数这一概念，并为此建模得出其标度值，首先我们把食品可能发生污染的环节给列出来，主要是源头，加工过程与销售过程，接着列出各个环节的子环节，其中包括农药化肥的使用，加工过程的环境与辅料的安全度，各环节食品合格水平，甚至包括从业人员的健康状况等，这些指标很大程度上的反应了我们所吃的食品是否安全，我们查阅了大量的相关资料，然后通过数学的方法将其量化处理，得到食品安全指数，它的产生不是无缘无故的，而是比较合理的，我们可以通过其值而让人们意识到目前我国食品安全的现状，从而采取相应措施，人们也会配合共同监督食品安全，努力提高其指数。6 模型的评价与优化优点： 1 .该模型运层次分析方法，模糊数学理论，内容翔实，结构简单明了，使得各指标的相对重要性一目了然，同时按照分解、比较、判断、综合的思维方式进行决策，而每一层的权重设置最后都会直接或间接影响到结果，而且在每个层次中的每个因素对结果的影响程度都是量化的，非常清晰、明确。 2.指标的相对重要性一目了然模型兼具评估及预测功能，可以鉴定过往绩效，又能为来年的安检工作提供参考价值。 3.移动平均法方法简单，只需要因变量的历史资料，减少了计算量。缺点1.人为确定权重的主观性强 。2.所查数据可能准确度不高，影响模型结果。3. 移动平均法对于未来的预测不够精确。模型的优化：在层次分析中可以多加几项指标，考虑更多对食品安全威胁的因素，从而使指数的确定更加准确而全面，对未来的预测中可以采取其它更高级的方法如指数平滑法，使预测更加准确。九、参考文献 1 MATLAB层次分析法教程 百度文库 百度百科 食品论坛2（王琬芜/张浩森/马士良/何美珍/）的食品安全指数(FSI)数学建模与评价3 中层次分析模型4 参考农业部统计公报、农业统计年鉴、卫生统计20052012年的统计数据5 参考李国平老师发的PPT课件中回归分析6 7 891010、 附录层次分析模型程序disp(请输入判断矩阵A(n阶);A=input(A=);n,n=size(A);x=ones(n,100);y=ones(n,100);m=zeros(1,100);m(1)=max(x(:,1);y(:,1)=x(:,1);x(:,2)=A*y(:,1);m(2)=max(x(:,2);y(:,2)=x(:,2)/m(2);p=0.0001;i=2;k=abs(m(2)-m(1);while kp i=i+1; x(:,i)=A*y(:,i-1); m(i)=max(x(:,i); y(:,i)=x(:,i)/m(i); k=abs(m(i)-m(i-1);enda=sum(y(:,i);w=y(:,i)/a;t=m(i);disp(权向量);disp(w);disp(最大特征值);disp(t); %以下是一致性检验CI=(t-n)/(n-1);RI=0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59;CR=CI/RI(n);if CR0.10 disp(此矩阵的一致性可以接受!); disp(CI=);disp(CI); disp(CR=);disp(CR);else disp(此矩阵的一致性不可以接受!);end 回归分析程序及图像x=2005 2006 2007 2008 2009 2010 2011 2012;X=ones(8,1),x;Y=0.8447 0.8451 0.8559 0.8542 0.8390 0.8502 0.8522 0.8756;b,bint,r,rint,stats=regress(Y,X);rcoplot(r,rint)层次分析求权重举例：disp(请输入矩阵A(n阶);A=input(A=);disp(请输入判断矩阵B(n阶);B=input(B=);s=A*B例：isp(请输入矩阵A(n阶);A=input(A=);disp(请输入矩阵B(n阶);B=input(B=);s=A*B请输入矩阵A(n阶)A=0.667 0.333请输入矩阵B(n阶)B=0.6537 0.5680 0.4230 0.366 0.246 0.2193 0.1173 0.4037;0.4112 0.3479 0.2781 0.2650 0.1614 0.2205 0.2693 0s = 0.5729 0.4947 0.3747 0.3324 0.2178 0.2197 0.1679 0.2693 膇芁蒁蚁袇膄莇蚀罿莀芃蚀肂膃薁虿螁莈蒇蚈袄膁莃螇羆莆艿螆肈腿薈螅螈羂薄螄羀芇蒀螄