别再手动算权重了！用MATLAB+熵权法优化你的TOPSIS评价模型（附完整代码）-洪萨配资

熵权法赋能TOPSIS：用MATLAB实现科学客观的多指标决策

在数据分析与决策科学领域，我们常常面临一个经典难题：如何从多个候选方案中选出最优解？当评价指标超过三个时，人脑已经难以直观比较，这时TOPSIS（优劣解距离法）便成为一盏明灯。但传统TOPSIS有个致命弱点——它默认所有指标权重相同，这在实际业务场景中几乎不可能成立。想象一下，在供应商评估时，交货准时率的重要性怎么可能和办公环境装修成本等同？

1. 权重困境：TOPSIS模型的阿喀琉斯之踵

1.1 等权重假设的现实悖论

我曾参与过一个医疗器械供应商评估项目，当使用传统TOPSIS时，发现某供应商因包装美观度得分高而排名靠前，但其灭菌合格率却低于行业标准。这种反直觉的结果正源于等权重假设的缺陷：

技术指标（如灭菌合格率）与非技术指标（如付款周期）被同等对待
关键绩效指标（KPI）与普通指标权重没有区分
不同量纲的指标简单相加（比如将ppm浓度与评分等级直接运算）

% 传统TOPSIS权重向量示例（等权重） equal_weights = [0.25, 0.25, 0.25, 0.25]; % 四个指标权重相同

1.2 主观赋权法的局限性

常见的AHP（层次分析法）虽然能设定不同权重，但面临三大挑战：

专家依赖性强：不同专家给出的判断矩阵可能差异巨大
一致性检验复杂：CR>0.1时需要反复调整判断矩阵
动态适应性差：当指标数据分布变化时，权重无法自动调整

实践发现：当评价指标超过7个时，AHP的判断矩阵一致性很难满足CR<0.1的要求

2. 熵权法：让数据自己说话

2.1 信息熵的核心哲学

熵权法的精妙之处在于其物理学基础——信息熵理论。指标数据的离散程度越高，说明该指标对方案区分度越大，理应获得更高权重。这就像学生考试成绩：

如果全班数学成绩都在85-90分之间（低离散度），则数学科目区分力弱
如果英语成绩从40分到95分分布（高离散度），则英语科目权重应该更高

2.2 熵权法四步实现框架

步骤1：数据标准化处理

function Z = standardize(X) [n,m] = size(X); Z = X ./ repmat(sum(X.^2).^0.5, n, 1); % 向量化运算提升效率 end

步骤2：计算概率矩阵

P = Z ./ repmat(sum(Z), size(Z,1), 1); % 每列和为1的概率分布 P(P==0) = 1e-11; % 避免log(0)错误

步骤3：计算信息熵

e = -1/log(n) * sum(P .* log(P)); % 指标信息熵向量

步骤4：确定权重

d = 1 - e; % 信息效用值 weights = d / sum(d); % 归一化权重

3. MATLAB实战：水质评估案例深度解析

3.1 数据准备与预处理

以某流域20个监测点的水质数据为例：

监测点	溶解氧(ppm)	PH值	大肠菌群(个/mL)	氨氮(mg/L)
A	4.69	6.59	5111	1.94
B	2.03	7.86	196	6.46
...	...	...	...	...

load('water_quality.mat'); [normalized, weights] = entropy_weight(X); disp(['各指标权重：', num2str(weights)]);

3.2 权重敏感度分析

通过蒙特卡洛模拟测试权重稳定性：

% 进行1000次随机扰动测试 weight_distribution = zeros(1000,4); for i = 1:1000 noisy_X = X .* (1 + 0.1*randn(size(X))); % 添加10%噪声 [~, temp_weights] = entropy_weight(noisy_X); weight_distribution(i,:) = temp_weights; end

3.3 与传统方法对比

下表对比三种权重确定方法的结果差异：

方法类型	计算耗时(s)	权重波动率	业务解释性
等权重法	0.01	0%	差
AHP法	32.7	15.2%	中等
熵权法	0.83	4.8%	强

4. 工程化应用中的进阶技巧

4.1 混合权重策略

对于既需要专家经验又需要数据驱动的场景，可采用组合权重：

alpha = 0.6; % 熵权法权重占比 combined_weights = alpha*entropy_weights + (1-alpha)*ahp_weights;

4.2 动态权重调整

当数据随时间变化时，建议采用滑动窗口更新权重：

function update_weights() global historical_data current_weights new_data = acquire_latest_data(); % 获取最新监测数据 historical_data = [historical_data(end-99:end,:); new_data]; % 保留最近100期 current_weights = entropy_weight(historical_data); % 重新计算权重 end

4.3 缺失值处理方案

当数据存在缺失时，推荐采用列均值填充+权重补偿策略：

missing_cols = any(isnan(X),1); X_filled = fillmissing(X, 'constant', mean(X,'omitnan')); weights(missing_cols) = weights(missing_cols) * 0.8; % 缺失指标权重打折

5. 避坑指南：来自实战的经验结晶

量纲陷阱：熵权法对指标量纲敏感，务必先标准化

% 错误做法：直接计算原始数据的熵权 % 正确做法：先进行极差标准化或Z-score标准化

零值处理：当数据含零时，概率计算会出错
```
P(P==0) = eps; % 用极小值替代零
```
权重稀释：当指标过多时，单个权重可能过小
```
% 解决方案：先通过聚类降维，再计算权重
```
结果验证：建议用Kendall相关系数检验排序稳定性
```
corr(original_rank, new_rank, 'type', 'Kendall')
```

在最近一次城市智慧水务项目中，我们采用熵权TOPSIS对50个泵站进行优先级排序。相比专家打分法，新方法将决策时间缩短70%，且选出的重点改造泵站在后续审计中确实发现了83%的设备隐患。这种数据驱动的决策方式，正在重新定义我们的评估范式。