for k = 1:N s = binvar(4,1); objective = objective + C{k} + W{k}; Costs = [sum(s) ==1; implies(s(1), [P_G{k}>=0, C{k} == 0.055*P_G{k}]); implies(s(2), [P_G{k}<=0, C{k} == -0.0055/3*P_G{k}]); implies(s(3), [P_G{k}==0, C{k} == 0]); implies(s(4), [P_B{k}>=0, W{k} == -0.055*2*P_B{k}])];
constraints = [constraints, x{k+1} == x{k}+P_B{k}]; constraints = [constraints, P_G{k}+P_Pv(k)-d(k)-P_B{k}==0]; constraints = [constraints, -500 <= P_B{k}<= 500, -300<= P_G{k}<=300, 150<=x{k}<=850]; end
implies(s(2), [P_G{k}<=0, C{k} == 0.0055/3*P_G{k}]);