function [Result,Cost,PDE,PWind,rate] = Yalmip_Cplex(Load,Wind,mpc,tan,wind)
%% 初始定义
yalmip;
T = 96; %调度周期
% 决策变量
P_DE1 = sdpvar(1,T,'full');
P_DE2 = sdpvar(1,T,'full');
P_Wind = sdpvar(1,T,'full');
rate = sdpvar(1,1,'full'); % 风电机组的装机容量系数
%% 约束
St = [];
St = [St, 0 <= rate(1,1) <= 1];
for t = 1:T
St = [St, mpc(2,1) <= P_DE1(t) <= mpc(1,1)];
St = [St, mpc(2,2) <= P_DE2(t) <= mpc(1,2)];
St = [St, 0 <= P_Wind(t) <= rate*Wind(t)];
St = [St, -0.1 <= P_DE1(t) + P_DE2(t) + P_Wind(t) - Load(t) <= 0.1];
end
%% 目标
Object = 0; % 总成本
F = 0; Tan = 0;
WIND_yunxing = 0; WIND_qifeng = 0; % 风电
% 煤耗成本
for t = 1:T
F = F + 0.25*(mpc(6,1)*P_DE1(t)^2 + mpc(5,1)*P_DE1(t) + mpc(4,1));
F = F + 0.25*(mpc(6,2)*P_DE2(t)^2 + mpc(5,2)*P_DE2(t) + mpc(4,2));
end
Object = ((F/1000)*1.5*700)/4; % 运行成本(成本转化)
% 产碳量
for t = 1:T
Tan = Tan + P_DE1(t)*mpc(3,1)*0.25*tan;
Tan = Tan + P_DE2(t)*mpc(3,2)*0.25*tan;
end
Object = Object + Tan;
% 风电成本
for t = 1:T
WIND_yunxing = WIND_yunxing + (rate*Wind(t)*1000) * wind(1)*0.25;
WIND_qifeng = WIND_qifeng + ((rate*Wind(t)-P_Wind(t))*1000) * wind(2)*0.25;
end
Object = Object + WIND_yunxing; % 风电运行
Object = Object + WIND_qifeng; % 弃风成本
%% 求解
Option = sdpsettings('solver','cplex','debug',0);
tic;
Result = optimize(St,Object,Option);
fprintf('模型求解时间：')
toc;
%% 表达
fprintf('全天运行费用：');
disp(value((F/1000)*1.5*700+Tan+WIND_yunxing+WIND_qifeng));
PDE = [value(P_DE1);value(P_DE2)];
PWind = value(P_Wind);
rate = value(rate);
Cost = [value((F/1000)*1.5*700),value(Tan),value(WIND_yunxing),value(WIND_qifeng)];
end
