找最合适的管教方案，只不过这次我们把管教权交给了粒子群算法)
基于PSO粒子群PID控制器参数整定粒子群PID/psopid/。 基于粒子群算法整定PID控制器，实现PID控制器参数的自整定(PSO-PID)。 matlab编写，源码注释详细具体如图，评价指标详细，直接运行即可。
先看段核心代码——粒子群初始化部分：
% 粒子群参数初始化 n_particles = 20; % 粒子侦察队数量 max_iter = 50; % 最大搜捕次数 w = 0.6; % 惯性权重(决定粒子有多恋旧) c1 = 1.7; % 自我认知系数 c2 = 1.7; % 社会认知系数 % PID参数取值范围 [Kp_min, Ki_min, Kd_min; Kp_max, Ki_max, Kd_max] param_range = [0 0 0; 20 20 20]; % 初始化粒子位置(参数组合)和速度 particles = rand(n_particles,3) .* (param_range(2,:)-param_range(1,:)) + param_range(1,:); velocities = zeros(n_particles,3);这里有几个关键点：粒子侦察队数量不宜过多（20个刚刚好），惯性权重相当于给粒子们保留30%的"老经验"。参数范围设置得像孙悟空给唐僧画的圈——Kp/Ki/Kd最大不超过20，防止参数暴走。
适应度函数是整个算法的灵魂，我们选用ITAE(时间乘绝对误差积分)作为考核指标：
function fitness = evaluatePID(parameters) % 解包PID三参数 Kp = parameters(1); Ki = parameters(2); Kd = parameters(3); % 搭建PID控制器 controller = pid(Kp,Ki,Kd); sys = feedback(controller * plant, 1); % 仿真获取系统响应 [y,t] = step(sys, sim_time); % 计算ITAE指标 error = ref - y; itae = sum(t .* abs(error) .* (t(:)'-t(end)+1)); % 时间加权误差积分 % 惩罚项：超调量超过5%要扣分 overshoot = max(y) - ref; if overshoot > ref*0.05 itae = itae * 1.2; % 超调惩罚系数 end fitness = itae; end这个评价函数暗藏玄机：时间加权让系统更关注稳态误差，超调量惩罚项就像交通罚单，防止参数组合飙车失控。注意那个(t(:)'-t(end)+1)的写法，其实是个时间归一化的小技巧，避免仿真时间过长导致ITAE爆炸。
基于PSO粒子群PID控制器参数整定粒子群PID/psopid/。 基于粒子群算法整定PID控制器，实现PID控制器参数的自整定(PSO-PID)。 matlab编写，源码注释详细具体如图，评价指标详细，直接运行即可。
粒子更新环节是算法最激动人心的部分，看着参数们像赶集的蚂蚁一样寻找最优路径：
for iter = 1:max_iter for i = 1:n_particles % 速度更新公式 velocities(i,:) = w * velocities(i,:) + ... c1*rand(1,3).*(pbest_pos(i,:) - particles(i,:)) + ... c2*rand(1,3).*(gbest_pos - particles(i,:)); % 位置越界处理 particles(i,:) = particles(i,:) + velocities(i,:); particles(i,:) = max(particles(i,:), param_range(1,:)); particles(i,:) = min(particles(i,:), param_range(2,:)); % 更新个体最优 current_fitness = evaluatePID(particles(i,:)); if current_fitness < pbest_val(i) pbest_val(i) = current_fitness; pbest_pos(i,:) = particles(i,:); end end % 更新全局最优 [min_val, idx] = min(pbest_val); if min_val < gbest_val gbest_val = min_val; gbest_pos = pbest_pos(idx,:); end % 动态调整惯性权重 w = w * 0.98; % 每次迭代减少2%的惯性 end速度更新公式里的rand(1,3)给算法注入了必要的随机性，就像炒菜时撒的不定量盐。动态惯性权重调整让算法前期大胆探索，后期精细调整，这个0.98的衰减系数是调参老手们用无数杯咖啡换来的经验值。
最终我们得到的可能是这样的参数组合：
最优PID参数: Kp = 8.742, Ki = 15.326, Kd = 2.114 ITAE指标: 1.257 超调量: 3.8% 调节时间: 2.3秒对比传统Ziegler-Nichols法4.2秒的调节时间和12%的超调量，PSO整定的参数就像给控制系统做了套精准按摩——既放松了肌肉(快速响应)又避免了抽筋(抑制超调)。
代码包里还藏着几个小惊喜：
- 实时收敛曲线可视化，看着粒子们逐渐聚拢的过程特别治愈
- 三维参数搜索轨迹图，像夜空中的萤火虫群
- 不同惯性权重策略的对比测试模块
- 支持自定义被控对象，把plant换成你的实际系统模型就能直接开整
最后说句实在话：这算法就像智能版的"瞎猫碰死耗子"，只不过碰得很有策略。下次调PID别手动较劲了，放粒子们出去遛弯，说不定就有惊喜。
)
