改进二进制粒子群算法配电网重构 可以动态生成配电网重构过程，目标函数为功率损耗，算例为IEEE33节点系统 程序简洁明了，注释详细
配电网重构就像是给电力系统做一场"动态手术"，通过调整开关状态让电流走最优路径。传统的粒子群算法处理这类离散问题容易陷入局部最优，咱们今天要聊的改进版二进制粒子群算法（BPSO）加入了惯性权重自适应和动态变异机制，配合IEEE33节点系统实战演示。
先看核心代码结构。网络数据用字典存储更直观，比如支路参数包含首末端节点和阻抗：
network = { 'branches': [ {'from':0, 'to':1, 'r':0.0922, 'x':0.0470, 'status':1}, # 其他32条支路数据... ], '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，包含详细的参数注释和测试案例。
