指标权重计算方法源代码

Jude

权重计算的四种方法，CRITIC法、AHP层次分析法、熵权法EWM、Kantiray权重法KWM。
欢迎大神们下载，讨论，使用。

Jude

AHP是网上找到的程序，其他三个都是自己编的。

kingjames

fenxiangdaimahuxiangxuexithanks

sciren

大神   CRITIC法方法有推荐的资料看看吗 百度都好少 看不明白

sciren

大神   CRITIC法方法导入 %导入评价指标X，列表示不同指标
怎么导入啊。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

sciren

146745        6361        1126        18.86        78.86        96.39        78.2        20        0.53
150114        6953        1188        17.76        80.26        96.57        72.81        20.4        0.44
167809        7582        1322        19.63        80.42        97.03        73.73        19.6        0.47
190186        7923        1314        21.48        77.54        96.66        80.43        18.3        0.34
211179        8668        1254        22.99        72.51        96.92        63.95        16.2        0.48
这是我的数据 为啥权重得出来是
0.9824    0.0162    0.0012    0.0000    0.0000    0.0000    0.0001    0.0000    0.0000
用的是CRITIC方法

Jude

sciren 发表于 2018-12-24 19:36
大神   CRITIC法方法有推荐的资料看看吗 百度都好少 看不明白

《考虑不同水文年及 IHA 指标相关性的水文特征评估方法》我参考的是这篇文献以及https://wiki.mbalib.com/wiki/Cri ... crieria_Correlation
这个上面的资料

Jude

sciren 发表于 2018-12-24 21:11
146745        6361        1126        18.86        78.86        96.39        78.2        20        0.53
150114        6953        1188        17.76        80.26        96.57        72.81        20.4        0.44
...

我算的是这个值：
0.1076        0.0942        0.09414        0.1026        0.1284        0.0924        0.120         0.1366        0.1240
应该是程序里数据标准化出了点问题

Jude

sciren 发表于 2018-12-24 21:11
146745        6361        1126        18.86        78.86        96.39        78.2        20        0.53
150114        6953        1188        17.76        80.26        96.57        72.81        20.4        0.44
...

%CRITIC权重法
%导入评价指标X，列表示不同指标
function [W]=CRITIC(X)

%数据标准化
Y=zeros(size(X,1),size(X,2));
for j=1:size(X,2)
    for i=1:size(X,1)
        Y(i,j)=(X(i,j)-min(X(:,j)))/(max(X(:,j)-min(X(:,j))));
    end
end
X=Y;

for j=1:size(X,2)
    for i=1:size(X,2)
        z=corrcoef(X(:,j),X(:,i));
        r(i,j)=z(1,2);
    end
    C(j)=std(X(:,j),1)*sum((1-r(:,j)));
end
%计算权重
for j=1:size(X,2)
    W(j)=C(j)/sum(C);
end
end

用这个试一下吧

sciren

我和网上的一个例子对比了一下下0.4830        13.2682        0.0000        4.3646        5.1070
0.4035        13.4909        39.0131        3.6151        5.5005
0.8979        25.7776        9.0513        4.8920        7.5342
0.5927        16.0245        13.2935        4.4529        6.5913
这是他的数据 然后他各指标的权重的结果是
0.12191        0.16497        0.38126        0.13543        0.19683
但是小可爱代码跑出来
0.1285    0.1337    0.4307    0.1790    0.1281