%% 目标函数
objective=0;
%发电成本
% for t=1:24
% objective=objective+(k_f*Pf(t)+b_f+one*Pf(t)+const);
% end
for t=1:24
objective=objective+(two*Pf(t)^2+b_f+one*Pf(t)+const);
end
%投资运维
for t=1:24
objective=objective+r_est*((1+r_est)^T)*(1+lamda_e)*(E0*100+P0*(Pdis(t)+Pchr(t)))/(((1+r_est)^T)-1) ;
end
%弃风惩罚
for t=1:24
objective=objective+(Pwind(t)-Pw(t))*fa_wind ;
end
%
%调峰不足惩罚成本
for t=1:24
objective=objective+Ploss(t)*fa_peak ;
end
%
% %延缓火电投资
% for t=1:24
% objective=objective-Pdis(t)*r_ther*((1+r_ther)^T_fire)*(1+lamda_th)/(24*k*((1+r_ther)^T_fire)-1)-Pdis(t)*gas*mei+Pdis(t)*dirt ;