请解释这段程序：%%%%机组组合%%%%; %%%%线性化MILP模型%%%; %%%考虑风电场景、不同电价场景、电动汽车充放电%%% %%%电动汽车数量按照10倍压缩%%%% clear clc %%%%%机组组合数据%%%%%%%%%%%; Ji=10;%机组数量; Time=24;%时间尺度; SS=20;%场景数量; Pmax=[455,455,130,130,162,80,85,55,55,55];%机组最大出力; Pmin=[150,150,20,20,25,20,25,10,10,10];%机组最小出力; a=[1000,970,700,680,450,370,480,660,665,670]; b=[16.19,17.26,16.60,16.50,19.7,22.26,27.74,25.92,27.27,27.79]; c=[0.00048,0.00031,0.002,0.0021,0.00398,0.00712,0.00079,0.00413,0.00222,0.00173]; Ton=[8,8,5,5,6,3,3,1,1,1];%最小开机时间; Toff=[8,8,5,5,6,3,3,1,1,1];%最小停机时间; Tcs=[5,5,4,4,4,2,2,0,0,0];%冷启动时间; Sh=[4500,5000,550,560,900,170,260,30,30,30];%热启动费用; Sc=[9000,10000,1100,1120,1800,340,520,60,60,60];%冷启动费用; T=[8,8,-5,-5,-6,-3,-3,-1,-1,-1];%初始运行状态; Xbefore=zeros(8,10); for t=1:8 for j=1:10 if T(j)+t<=0 Xbefore(t,j)=0; else Xbefore(t,j)=1; end end end Xf=Xbefore(1,:);%初始序列; Xbefore=[zeros(1,10);Xbefore]; PL=[700,750,850,950,1000,1100,1150,1200,1300,1400,1450,1500,1400,1300,1200,1050,1000,1100,1200,1400,1300,1100,900,800];%日负荷; delta_hot=[130,130,60,60,90,40,40,40,40,40];%爬坡速率 delta_cold=[150,150,20,20,25,20,25,10,10,10];%开停机爬坡速率 R=0.1*PL;%备用容量，这里取10%PL; Pf=[89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84,89.84;99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76,99.76;100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100,100;86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52,86.52;93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54,93.54;85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14,85.14;76.14,76.14,76.14,76.14,76.14,68.36,68.36,68.36,61.72,61.72,61.72,61.72,61.72,61.72,58.32,58.32,58.32,58.76,58.76,58.76;57.72,57.72,57.72,57.72,57.72,42.48,42
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
