% 电力系统参数
num_generators = 5; % 发电机数量
load_demand = 1000; % 负荷需求(MW)
% 发电机参数
generator_capacity = [200, 300, 250, 150, 200]; % 发电机容量(MW)
generator_cost = [10, 12, 11, 9, 10]; % 发电机成本($/MWh)
% 粒子群算法参数
num_particles = 50; % 粒子数量
max_iter = 100; % 最大迭代次数
w = 0.5; % 惯性权重
c1 = 1; % 个体学习因子
c2 = 1; % 社会学习因子
% 初始化粒子群
positions = zeros(num_particles, num_generators);
velocities = zeros(num_particles, num_generators);
pbest_positions = positions;
pbest_costs = inf(1, num_particles);
gbest_position = zeros(1, num_generators);
gbest_cost = inf;
% 开始迭代
for iter = 1:max_iter
% 更新速度和位置
for i = 1:num_particles
% 更新速度
velocities(i,