有功优化问题

%潮流计算 + 发电调度
%输入有功出力的功率  发电成本的计算
%4台可调发电器(支路6-9、6-10、4-12、28-27)、2个无功补偿装置(节点10、24)。
%发电机的无功出力上、下限成在1序点、发电成本与总有功发电量成正比。
%负荷节点电压在允许的范围内0.95~1.05；发电机组允许电压范围0.9~1.1。
%步长0.02；补偿装置的最大可调范围为4组，节点10、24的装置电容器额定容量分别为0.05和0.01；
%4台发电器 11 12 15 36 6台发电机(节点1分别对应节点前、节点2、5、8、11、13为PV节点)
%控制变量的顺序：T * 4 V5 I C2 Pg * 5
function y = f_user(x)
mpc = case30;

mpc.branch([11 12 15 36],9) = 0.90 + ceil(x(1,1:4))*0.02;
%变压器的变比为第9列
mpc.gen(2:6,6) = x(1,5:9);%发电机的端电压为第六列
mpc.gen(2:6,2) = x(1,12:16);

mpc.bus([2 5 8 11 13],8) = x(1,5:9);
mpc.bus([10 24],4) = mpc.bus([10 24],4) - ceil(x(1,10:11)) .* [0.05 0.01];

%01 潮流计算
res = runpf(mpc);

Ploss = res.branch(:,14) + res.branch(:,16);

%02 发电调度
cost = sum(x(1,12:16)*1.5);

end