水資源短缺風險綜合評價(修復的)-資料下載頁
2025-06-26 18:29本頁面
【正文】 11 26 24 38 作用水量月缺水量對時間的圖形走勢程序：x1=[1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008]x3=[ ]x4=[ ]x2=[ ]x5=[ ]x6=[ 26 24 38 ]x7=x2x6plot(x1,x2,x1,x3,x1,x4,x1,x5,x1,x6,x1,x7)問題二中標準計算程序：x6=[ 26 24 38 ]average=mean(x6)I1=average*。I2=average*。I3=average*。求Q和e的相關程序：function c=judge(x)%輸入?yún)?shù)x為監(jiān)測值三維矩陣for k=1:6 for i=1:3 for j=1:5 c(i,j,k)=(x(i,j,k)min(x(i,:,k)))./(max(x(i,:,k))min(x(i,:,k)))。 end endendfunction w=powershan(h)%輸入為熵值矩陣，輸出為各評價指標熵權矩陣for k=1:6 for i=1:3 w(i,1,k)=(1h(i,1,k))/(3sum(h(:,k)))。 endendfunction h=shan(c)%輸入為所有監(jiān)測值判斷矩陣，輸出為各評價標準熵值for k=1:6 for i=1:3 for j=1:5 f(i,j,k)=(1+c(i,j,k))/(sum(c(i,:,k))+1)。 end endend for k=1:6 for i=1:3 h(i,1,k)=f(i,:,k)*(log(f(i,:,k)))39。./log(5)。 endendfunction Q=testexp(x,w,I2)%輸入為監(jiān)測數(shù)值矩陣和各指標熵權矩陣，輸出為綜合評測指數(shù)for k=1:6 for i=1:3 d(i,:,k)=x(i,:,k)I2(i)。 D(i,1,k)=sum(d(i,:,k))。 endendfor k=1:6 Q(1,1,k)=(w(:,1,k))39。*D(:,1,k)。End%y1=[ 。 。 ]。%19791983y2=[ 。 。 ]。%19841988y3=[ 。 。 ]。%19891993y4=[ 。 。 ]。%19941998y5=[ 。 。 ]。%19992003y6=[ 。 。 ]。%20042008a=cat(3,y1,y2,y3,y4,y5,y6)。 %監(jiān)測值矩陣，頁為分組，行為指標，列為年份I1=[,12]。I2=[,]。%評價標準c=judge(a)。 %歸一化后的判斷矩陣h=shan(c)。 %各分組中各個評價指標的熵w=powershan(h) %評價指標熵權矩陣for i=1:3 w_(i)=mean(w(i,1,:))。 %各評價指標的熵權平均值endQ=testexp(a,w,I2) %各分組綜合評測指數(shù)Qs=abs(I1I2).*5*w_39。 %衡量標準e=sum(s)/sum(abs(s)) %程度系數(shù)19 / 19
環(huán)評公示相關推薦
鄂ICP備17016276號-1