熵权法(matlab实现)

clear clc x 1 65 1 88 30 40 41 10 5 40 1 55 16 61 250 16 8 04 133 98 1 25 62 90 27 81 76 60 118 00 21 67 2 97 2 19 39 90 9 10 1 69 2 28 41 20 32 30 7 20 1 20 20 18 265 54 1 45 88 94 0 44 80 60 27 01 78 40 126 10 34 84 1 92 3 67 35 90 12 80 1 73 1 17 37 10 25 60 5 40 1 66 20 01 259 60 3 18 52 30 1 71 62 30 29 00 74 30 115 80 42 92 2 27 4 74 32 40 12 90 2 08 2 17 57 10 31 90 7 30 1 55 24 24 201 89 5 07 138 75 4 63 56 40 20 31 69 20 103 40 34 75 2 05 3 3 5 33 70 9 10 a min x b max x n m size x k 1 log n for i 1 n for j 1 m x i j x i j a j b j a j 100 高优指标处理 end end for i 1 n for j 1 m x i j b j x i j b j a j 100 低优指标处理 end end he sum x for i 1 n for j 1 m p i j x i j he j end end 指标归一化 for i 1 n for j 1 m if p i j 0 z i j 0 else z i j log p i j end end end e zeros 1 m for i 1 1n for j 1 m e j e j p i j z i j k end end he sum e for i 1 m g i 1 e i m he end for i 1 m w i g i sum g 计算权重 end s zeros 1 n for i 1 n 计算综合得分 for j 1 m s i s i w j p i j end end t plot t s b data a input A B C D cards 2 2 1 0 3 1 3 2 proc transpose data a out aa run proc stdize data aa out aaa method range var col1 col2 run proc transpose data aaa out result drop name var col1 col2 id name run function weights EntropyWeight 1 2 3 4 5 6 熵权法求指标权重 R 为输入矩阵 返回权重向量 weights rows cols size 2 3 输入矩阵的大小 rows 为对象个数 cols 为指标个数 k 1 log rows 求 k f zeros rows cols 初始化 fij sumBycols sum R 1 输入矩阵的每一列之和 结果为一个 1 cols 的行向量 计算 fij for i 1 rows for j 1 cols f i j R i j sumBycols 1 j end end lnfij zeros rows cols 初始化 lnfij 计算 lnfij for i 1 rows for j 1 cols if f i j 0 lnfij i j 0 else lnfij i j log f i j end end end Hj k sum f lnfij 1 计算熵值 Hj weights 1 Hj cols sum Hj end A 2 2 1 0 3 1 3 0 A 为所要计算的权重矩阵 行数为指标数 列数为专家数 m size A 1 m 存放 A 的列数 n size A 2 n 存放 A 的行数 B ones m 1 C ones m 1 b 1 C 0 For j 1 m a 1 for i 1 n a a A j i 计算第 j 行权值的乘积 End b a 1 n B j b C c b End B for j 1 m 以下为归一化处理 计算权重值 C j B j c End C C 为所求的权重值矩阵 R 1 65 1 88 30 40 41 10 5 40 1 55 16 61 250 16 8 04 133 98 1 25 62 90 27 81 76 60 118 00 21 67 2 97 2 19 39 90 9 10 1 69 2 28 41 20 32 30 7 20 1 20 20 18 265 54 1 45 88 94 0 44 80 60 27 01 78 40 126 10 34 84 1 92 3 67 35 90 12 80 1 73 1 17 37 10 25 60 5 40 1 66 20 01 259 60 3 18 52 30 1 71 62 30 29 00 74 30 115 80 42 92 2 27 4 74 32 40 12 90 2 08 2 17 57 10 31 90 7 30 1 55 24 24 201 89 5 07 138 75 4 63 56 40 20 31 69 20 103 40 34 75 2 05 3 3 5 33 70 9 10 rows cols size R k 1 log rows f zeros rows cols sumBycols sum R 1 for i 1 rows for j 1 cols f i j R i j sumBycols 1 j end end lnfij zeros rows col