【正文】 )== f54=1。else if x(i,5)amp。amp。x(i,5) f54=(x(i,5))/()。 else if x(i,5)amp。amp。x(i,5) f54=(x(i,5))/()。 else f54=0。 end endendif x(i,5)= f55=0。else if x(i,5)amp。amp。x(i,5) f55=(x(i,5))/()。 else f55=1。 endend %if x(i,6)= f61=1。else if x(i,6)amp。amp。x(i,6) f61=(x(i,6))/()。 else f61=0。 endendif x(i,6)== f62=1。else if x(i,6)amp。amp。x(i,6) f62=(x(i,6))/()。 else if x(i,6)amp。amp。x(i,6) f62=(x(i,6))/()。 else f62=0。 end endendif x(i,6)== f63=1。else if x(i,6)amp。amp。x(i,6) f63=(x(i,6))/()。 else if x(i,6)amp。amp。x(i,6) f63=(x(i,6))/()。 else f63=0。 end endendif x(i,6)== f64=1。else if x(i,6)amp。amp。x(i,6) f64=(x(i,6))/()。 else if x(i,6)amp。amp。x(i,6) f64=(x(i,6))/()。 else f64=0。 end endendif x(i,6)= f65=0。else if x(i,6)amp。amp。x(i,6) f65=(x(i,6))/()。 else f65=1。 endend%if x(i,7)= f71=1。else if x(i,7)amp。amp。x(i,7) f71=(x(i,7))/()。 else f71=0。 endendif x(i,7)== f72=1。else if x(i,7)amp。amp。x(i,7) f72=(x(i,7))/()。 else if x(i,7)amp。amp。x(i,7) f72=(x(i,7))/()。 else f72=0。 end endendif x(i,7)== f73=1。else if x(i,7)amp。amp。x(i,7) f73=(x(i,7))/()。 else if x(i,7)amp。amp。x(i,7) f73=(x(i,7))/()。 else f73=0。 end endendif x(i,7)== f74=1。else if x(i,7)amp。amp。x(i,7) f74=(x(i,7))/()。 else if x(i,7)amp。amp。x(i,7) f74=(x(i,7))/()。 else f74=0。 end endendif x(i,7)= f75=0。else if x(i,7)amp。amp。x(i,7) f75=(x(i,7))/()。 else f75=1。 endend%PM10一小時(shí)的隸屬度if x(i,8)= f81=1。else if x(i,8)amp。amp。x(i,8) f81=(x(i,8))/()。 else f81=0。 endendif x(i,8)== f82=1。else if x(i,8)amp。amp。x(i,8) f82=(x(i,8))/()。 else if x(i,8)amp。amp。x(i,8) f82=(x(i,8))/()。 else f82=0。 end endendif x(i,8)== f83=1。else if x(i,8)amp。amp。x(i,8) f83=(x(i,8))/()。 else if x(i,8)amp。amp。x(i,8) f83=(x(i,8))/()。 else f83=0。 end endendif x(i,8)== f84=1。else if x(i,8)amp。amp。x(i,8) f84=(x(i,8))/()。 else if x(i,8)amp。amp。x(i,8) f84=(x(i,8))/()。 else f84=0。 end endendif x(i,8)= f85=0。else if x(i,8)amp。amp。x(i,8) f85=(x(i,8))/()。 else f85=1。 endend%PM10二十四小時(shí)的隸屬度if x(i,9)= f91=1。else if x(i,9)amp。amp。x(i,9) f91=(x(i,9))/()。 else f91=0。 endendif x(i,9)== f92=1。else if x(i,9)amp。amp。x(i,9) f92=(x(i,9))/()。 else if x(i,9)amp。amp。x(i,9) f92=(x(i,9))/()。 else f92=0。 end endendif x(i,9)== f93=1。else if x(i,9)amp。amp。x(i,9) f93=(x(i,9))/()。 else if x(i,9)amp。amp。x(i,9) f93=(x(i,9))/()。 else f93=0。 end endendif x(i,9)== f94=1。else if x(i,9)amp。amp。x(i,9) f94=(x(i,9))/()。 else if x(i,9)amp。amp。x(i,9) f94=(x(i,9))/()。 else f94=0。 end endendif x(i,9)= f95=0。else if x(i,9)amp。amp。x(i,9) f95=(x(i,9))/()。 else f95=1。 endend R=[f11,f12,f13,f14,f15。 f21,f22,f23,f24,f25。 f31,f32,f33,f34,f35。 f41,f42,f43,f44,f45。 f51,f52,f53,f54,f55。 f61,f62,f63,f64,f65。 f71,f72,f73,f74,f75。 f81,f82,f83,f84,f85。 f91,f92,f93,f94,f95]。% R為模糊關(guān)系矩陣s1=1/5*(++++)。s2=1/5*(5+10+35+60+90)。s3=1/5*(++++)。s4=1/5*(++++)。s5=1/5*(++++)。s6=1/5*(++++)。s7=1/5*(++++)。s8=1/5*(++++)。s9=1/5*(++++)。k1=(x(i,1)/s1)。k2=(x(i,2)/s2)。k3=(x(i,3)/s3)。k4=(x(i,4)/s4)。k5=(x(i,5)/s5)。k6=(x(i,6)/s6)。k7=(x(i,7)/s7)。k8=(x(i,8)/s8)。k9=(x(i,9)/s9)。w0=k1+k2+k3+k4+k5+k6+k7+k8+k9。f1=k1/w0。f2=k2/w0。f3=k3/w0。f4=k4/w0。f5=k5/w0。f6=k6/w0。f7=k7/w0。f8=k8/w0。f9=k9/w0。W=[f1,f2,f3,f4,f5,f6,f7,f8,f9]。% W為權(quán)重矩陣S=[50,100,150,200,300]。% S為空氣質(zhì)量標(biāo)準(zhǔn)類別向量FCI(i)=round(W*R*S39。)。% FCI為四舍五入取整后的模糊綜合指數(shù)End 程序三：x=[1 2 3 4 5 6 7 8 9 10 11]。xi=linspace(min(x),max(x))。yi1=interp1(x,y1,xi,39。spline39。)。 yi2=interp1(x,y2,xi,39。spline39。)。 yi3=interp1(x,y3,xi,39。spline39。)plot(xi,yi1,xi,yi2,xi,yi3)。legend(39。原始數(shù)據(jù)39。,39。國(guó)標(biāo)39。,39。模型39。)。32