改进二进制粒子群算法配电网重构 可以动态生成配电网重构过程目标函数为功率损耗算例为IEEE33节点系统 程序简洁明了注释详细配电网重构就像是给电力系统做一场动态手术通过调整开关状态让电流走最优路径。传统的粒子群算法处理这类离散问题容易陷入局部最优咱们今天要聊的改进版二进制粒子群算法BPSO加入了惯性权重自适应和动态变异机制配合IEEE33节点系统实战演示。改进二进制粒子群算法配电网重构 可以动态生成配电网重构过程目标函数为功率损耗算例为IEEE33节点系统 程序简洁明了注释详细先看核心代码结构。网络数据用字典存储更直观比如支路参数包含首末端节点和阻抗network { branches: [ {from:0, to:1, r:0.0922, x:0.0470, status:1}, ], loads: [ {node:1, p:100, q:60}, # 其他节点负荷数据... ] }粒子位置用二进制数组表示开关状态1代表闭合0代表断开。初始化种群时特别注意保持网络连通性def init_particle(): particle np.random.randint(0,2,33) while not is_radial(particle): # 确保辐射状结构 particle np.random.randint(0,2,33) return particle适应度计算是核心环节这里用前推回代法进行潮流计算。网损公式直接体现在目标函数中def power_loss(network, status): # 更新支路通断状态 for i in range(33): network[branches][i][status] status[i] # 执行潮流计算 loss forward_backward_sweep(network) return loss改进算法的精髓在速度更新环节。传统二进制PSO直接用sigmoid函数映射开关概率我们加入了非线性惯性权重w w_max - (w_max - w_min) * (iter/max_iter)**2 # 动态惯性权重 v w*v c1*rand()*(pbest - x) c2*rand()*(gbest - x) prob 1 / (1 np.exp(-v)) # sigmoid转换概率 x_new (prob np.random.rand()).astype(int)这里平方项让算法前期侧重全局搜索后期加强局部开发。还有个彩蛋——当连续3代最优解未更新时触发变异操作if stagnation_counter 3: idx np.random.randint(33) gbest[idx] 1 - gbest[idx] # 翻转最优解某个基因 stagnation_counter 0运行结果展示时用matplotlib绘制网损收敛曲线和拓扑变化动画。典型情况下网损从初始的202.4kW降至142.3kW降幅约30%。重构过程中会看到某些关键支路反复开合直到找到最优结构。几点编程注意使用numpy向量化操作加速计算连通性检查用并查集算法实现对于大规模系统可以采用并行计算适应度值。完整代码已开源在Github包含详细的参数注释和测试案例。