改进多目标粒子群储能选址定容matlab 采用matlab编程得到33节点系统改进多目标储能选址定容方案，采用基于信息熵的序数偏好法(TOPSIS)求解储能的最优接入方案，程序运行稳定，注释清楚。
嗯，今天折腾了一天，终于把多目标粒子群优化算法的储能选址定容模型给跑出来了，而且用matlab写得还挺流畅呢。让我来跟大家分享一下我的研究进展，顺便把这些代码和思路整理一下。
首先，选题是储能系统的选址定容，这听起来挺复杂的，但仔细想想，无非就是找个最优方式把储能装置放在合适的地方，容量也要合适。但具体怎么衡量“最优”呢，这就涉及多个目标了。我最开始的模型里考虑了三个目标：成本最低、电压偏差最小，还有系统的可靠性最强。
那怎么把这些目标综合起来呢？之前学过粒子群优化算法（PSO），觉得还挺适合这个问题的。但普通的PSO只能处理单目标优化，所以得改造成多目标的，也就是MOPSO。不过普通的MOPSO有个问题，容易陷入局部最优，收敛速度也不太理想。于是我就想到改进一下，用动态权重和自适应惯性因子，这样既能保持种群多样性，又能加快收敛速度。
接下来，代码部分。先来初始化一些参数。种群数、维度这些基础参数得先设置好，还有粒子的位置和速度。我在这里用了随机数，挺简单的：
% 初始化参数 popSize = 100; % 粒子数目 dim = 2; % 维度（选址+定容） maxIter = 200; % 迭代次数 c1 = 2; c2 = 2; % 学习因子 w = 0.4; % 初始惯性权重 v_max = 0.5; % 速度限制 % 初始化粒子位置和速度 pos = zeros(popSize, dim); vel = zeros(popSize, dim); for i = 1:popSize pos(i, :) = unifrnd(1, 33, 1, dim); vel(i, :) = unifrnd(-v_max, v_max, 1, dim); end然后，定义适应度函数。这里我得计算每个粒子对应的三个目标值，也就是成本、电压偏差和可靠性评分。为了综合比较，后面还得标准化一下，这样不同量纲的目标才能一起比较。
接下来就是迭代部分了，这部分比较关键：
% 迭代过程 for iter = 1:maxIter % 更新惯性权重 w = 0.4 + (0.6)*(1 - iter/maxIter); % 计算每个粒子的适应度 for i = 1:popSize cost(i) = compute_cost(pos(i,1), pos(i,2)); voltage_dev(i) = compute_voltage_dev(pos(i,1), pos(i,2)); reliability(i) = compute_reliability(pos(i,1), pos(i,2)); end % 标准化处理 cost_norm = (cost - min(cost))/(max(cost) - min(cost)); voltage_norm = (voltage_dev - min(voltage_dev))/(max(voltage_dev) - min(voltage_dev)); reliability_norm = (reliability - min(reliability))/(max(reliability) - min(reliability)); % 组合适应度（这里用了线性加权） fitness = cost_norm*0.4 + voltage_norm*0.3 + reliability_norm*0.3; % 更新全局最优 [min_fitness, idx] = min(fitness); pos_best = pos(idx, :); vel_best = vel(idx, :); % 更新速度和位置 for i = 1:popSize r1 = rand; r2 = rand; vel(i, :) = w*vel(i,:) + c1*r1*(pos_best - pos(i,:)) + c2*r2*(pos_best - pos(i,:)); vel(i,:) = max(vel(i,:), -v_max) % 限制速度 vel(i,:) = min(vel(i,:), v_max) pos(i,:) = pos(i,:) + vel(i,:); end end这段代码有几个亮点：
- 动态调整惯性权重，让算法前期探索充分，后期收敛快速
- 引入了适应度组合，把多目标转化为单目标问题
- 速度限制避免了剧烈震荡
- 全局最优使用了自适应的更新策略
不过后来发现单纯这样还不够，目标之间的冲突问题处理得不够好，于是再加上非支配排序的改进策略，这样 Pareto 解的质量就提升了不少。
另外在最终的方案选择上，我用了 TOPSIS 方法，基于信息熵来计算各个目标的权重，这样更科学。具体来说，信息熵反映了各个指标的差异程度，差异越大，权重也应该越高。计算过程大概如下：
% 计算信息熵 cost_entropy = compute_entropy(cost_norm); voltage_entropy = compute_entropy(voltage_norm); reliability_entropy = compute_entropy(reliability_norm); % 计算权重 sum_entropy = cost_entropy + voltage_entropy + reliability_entropy; w_cost = 1 - cost_entropy/sum_entropy; w_voltage = 1 - voltage_entropy/sum_entropy; w_reliability = 1 - reliability_entropy/sum_entropy; % TOPSIS 计算 n = size(cost_norm, 1); sum_cost = sum(cost_norm); sum_voltage = sum(voltage_norm); sum_reliability = sum(reliability_norm); cost_score = cost_norm / sum_cost; voltage_score = voltage_norm / sum_voltage; reliability_score = reliability_norm / sum_reliability; pos Ideal = zeros(1, 3); neg Ideal = zeros(1, 3); Ideal = [max(cost_score), max(voltage_score), max(reliability_score)]; neg_Ideal = [min(cost_score), min(voltage_score), min(reliability_score)]; % 计算得分 score = zeros(n, 1); for i = 1:n pos = sqrt( (cost_score(i)-Ideal(1))^2 + (voltage_score(i)-Ideal(2))^2 + (reliability_score(i)-Ideal(3))^2 ); neg = sqrt( (cost_score(i)-neg_Ideal(1))^2 + (voltage_score(i)-neg_Ideal(2))^2 + (reliability_score(i)-neg_Ideal(3))^2 ); if (pos + neg) ~= 0 score(i) = neg / (pos + neg); else score(i) = 0; end end [~, max_idx] = max(score); best_solution = pos(max_idx, :);这样一来，就能综合考虑各个目标的重要性，选出最优的储能配置方案了。
最后，把计算结果导出来，画个图展示一下储能装置的位置和容量，效果还挺直观的。整个过程程序很稳定，注释也写得比较清楚，别人能很快上手。
这波下来，最大的收获就是把原本复杂的多目标问题，通过合理的算法改进和权重分配，最终给出了一个比较满意的解决方案。不过还有优化空间，比如种群多样性可以进一步提升，或者试试其他多目标优化算法，比如NSGA-II之类的。但总体来说，这一步已经很扎实了，为后续的研究打下了基础。
