数学建模——熵权法步骤及程序详解
权重的求解一直都是数学建模的重点关注对象，所以学好建模论文的重要一步就是如何确定权重，今天是来介绍一种客观确定几个指标各自所占的权重的方法——熵权法。之前的数学建模实战里有提到用熵权法确定了每个指标各自的权重，这里展开详细的写一下。
文章目录
- 数学建模——熵权法步骤及程序详解
- 前言
- 一、熵权法的介绍
- 1、熵权法的应用场景
- 2、熵权法的基本思想
- 3、熵权法的算法步骤
- 二、代码程序
- 总结
前言
按照信息论基本原理的解释，信息是系统有序程度的一个度量，熵是系统无序程度的一个度量；根据信息熵的定义，对于某项指标，可以用熵值来判断某个指标的离散程度，其信息熵值越小，指标的离散程度越大， 该指标对综合评价的影响（即权重）就越大，如果某项指标的值全部相等，则该指标在综合评价中不起作用。因此，可利用信息熵这个工具，计算出各个指标的权重，为多指标综合评价提供依据。
一、熵权法的介绍
我们分析某个因素时，主观地列出了一些有影响力的指标，这些指标已经与因变量存在某种的相关性，不管是正相关还是负相关。指标的值变化会直接影响因素的变化，变化量越大，说明指标对于因素的变化作用也应该是越明显的。
简单的讲，就是在评价对象时，往往每个对象会有几个指标。那这几个指标哪个指标所占的权重最大呢？我们自己捏造一个好像又不太合理，所以我们瞎分析一下，好像一些变化不是那么大，整体相对平滑的数据，举一个极端的例子：对于所有的样本而言，这个指标都是相同的数值，那么我们可认为这个指标的权值为0，即这个指标对于我们的评价起不到任何帮助。依据这种思想我们就诞生了熵权法。
1、熵权法的应用场景
熵权法既然是客观赋权法，那么应用有一个前提条件，那就是必须得有数据，没有数据支撑的赋权方法是主观赋权法。
这种方法虽然简单，但是用于企业的实际问题当中是非常有用的，可以用于确定初步权重，后续可通过其他方法进行优化。
2、熵权法的基本思想
我们对于一件事情的普遍看法是，越有可能发生的事情，信息量越少；越不可能发生的事情，信息量就越多。而熵的定义我们在高中化学中已经接触过了，有兴趣的同学可以回顾一下，熵的定义是系统的混乱程度，熵越高，系统的状态越混乱，这个定义就很符合我们对事情的普遍看法的定义，所以诞生了名词——信息熵。
一个随机变量的信息熵越高，则他所带来的信息量越高，那么相对应的，其他变量的信息量会比较低（因为总体信息量(权重)是固定的，为1）
啊，说到这里，大家不是很理解也很正常滴，主要写论文会用，会抄就可以了。下面进入重点，
3、熵权法的算法步骤
第一步为数据标准化。
首先需要正向化指标（后续需要用到），当然，如果数据都是非负的，那么可以跳过正向化步骤，公式如下。
就是指标的值减去最小值，除以指标的振幅。这样子就可以得到正向化后的矩阵Z。
第二步为进行归一化处理
处理公式如下。
y i j = z i j / ∑ i = 1 n x i j y_{ij}=z_{ij}/\sum_{i=1}^{n}{x_{ij}} yij=zij/i=1∑nxij
此外还可以采用最大最小法、Z-score等方法进行标准化处理。由此可得到行列的标准化矩阵(共有n行m列，一行代表一个样本，一列代表一个指标)
Y = ( y i j ) n m Y=(y_{ij})_{nm} Y=(yij)nm
第三步计算每个指标的信息熵，并计算信息效用值，并归一化得到每个指标
公式如下：
e j = k ∑ i = 1 n y i j ln y i j e_j=k\sum_{i=1}^n{y_{ij}\ln y_{ij}} ej=ki=1∑nyijlnyij
这里为什么要除以一个常数 k k k呢？ k k k的取值是为了保证 e j e_j ej可以在区间[0,1]上，当指标的样本值完全一样的时候，k= 1 ln n \frac{1}{\ln\text{\ }n} ln n1刚好可以令信息熵位于[0,1]之间，所以我们 k k k值一般取 1 ln n \frac{1}{\ln\text{\ }n} ln n1
熵值越小，信息量越大,有兴趣大家可以试一下画他的函数图像。
第四步计算信息效用值（信息熵冗余度）
计算方法如下：
d j = 1 − e j d_j=1-e_j dj=1−ej
信息效用越大，信息量越大。
第五步计算指标的权重系数
将信息效用值进行归一化，就可以得到每个指标的权重：
w j = d j / ∑ i = 1 m d i w_j=d_j/\sum_{i=1}^m{d_i} wj=dj/i=1∑mdi
二、代码程序
Matlab代码如下：
clc,clear %行为样本，列为指标 X=[39414 2823 34877 44562 2036 603 322 934936 929914 1492 29811 54934 1911 52242 35262 3862 908 396 1075563 1030664 1780 29811 96442 2743 88737 303221 4307 1596 694 1104835 1010146 1936 32678 107079 3036 98513 478883 3956 2530 1089 909220 862077 2160 36063 124359 3326 116897 378318 4102 2669 1179 1117851 1123109 2349 38951 140167 3900 130355 261203 4180 3538 1991 1116429 1100510 2446 40324 161523 3989 153722 444755 4309 3727 1593 878466 880226 2637 43211 177681 4669 167161 422267 4630 6629 1867 1048053 1003952 2904 47116 124969 4416 111415 286399 3829 5665 2591 1142395 1112661 3092 49406 146015 3200 129997 228695 5308 4911 2506 1202365 1112475 3252 51119]; [n,m]=size(X); %% 正向化 z=find(X
