clear
clc
yalmip;
Cplex;
%%系統(tǒng)參數(shù)
%所有參數(shù)均用有名值表示
paragen=xlsread（‘excel2017‘‘機(jī)組參數(shù)‘）;
loadcurve=xlsread（‘excel2017‘‘負(fù)荷曲線‘）;
netpara=xlsread（‘excel2017‘‘網(wǎng)絡(luò)參數(shù)‘）;
branch_num=size（netpara）;%網(wǎng)絡(luò)中的支路
branch_num=branch_num（11）;
PL_max=netpara（:6）;%線路最大負(fù)荷
PL_min=netpara（:7）;%線路最小負(fù)荷
limit=paragen（:3:4）;%機(jī)組出力上下限//limit（:1）表示上限，limit（:2）表示下限
para=paragen（:5:7）;%成本系數(shù)//para（:1）表示系數(shù)apara（:2）表示系數(shù)bpara（:3）表示系數(shù)c。
price=100;
para=price*para;%價(jià)格換算
lasttime=paragen（:9）;%持續(xù)時(shí)間
Rud=paragen（:8）;%上下爬坡速率//因題中簡(jiǎn)化上坡下坡速度相同
H=paragen（:10）;%啟動(dòng)成本
J=paragen（:11）;%關(guān)停成本
u0=[1?1?1?1?1?1];%初始狀態(tài)
%%?規(guī)模變量
%機(jī)組數(shù)
gennum=size（paragen）;
gennum=gennum（11）;
%節(jié)點(diǎn)數(shù)
numnodes=size（loadcurve）;
numnodes=numnodes（11）-1;
%時(shí)間范圍
T=size（loadcurve）;
T=T（12）-1;
%線性化分段數(shù)（按需要更改）
m=4;
%各時(shí)刻節(jié)點(diǎn)總負(fù)荷
PL=loadcurve（numnodes+12:T+1）;
%%
%決策變量
u=binvar（gennumT‘full‘）;%狀態(tài)變量
p=sdpvar（gennumT‘full‘）;%即各機(jī)組實(shí)時(shí)功率p（it）
Ps=sdpvar（gennumTm‘full‘）;%分段出力
costH=sdpvar（gennumT‘full‘）;%啟動(dòng)成本
costJ=sdpvar（gennumT‘full‘）;%關(guān)停成本
sum_PowerGSDF=sdpvar（Tbranch_numnumnodes‘full‘）;%發(fā)電機(jī)的輸出功率轉(zhuǎn)移總和
%%?目標(biāo)函數(shù)線性化
MaxPs=zeros（gennumTm）;%這里表示分段出力的上限
st=[];%st約束初始化
for?i=1:gennum???%目標(biāo)函數(shù)線性化后分段出力的不等式約束
???for?t=1:T
?????for?s=1:m
	MaxPs（its）=（limit（i1）-limit（i2））/m;
????st=st+[Ps（its）>=0Ps（its）<=MaxPs（its）];
?????end
???end
end
K=zeros（gennumm）;%煤耗函數(shù)的斜率值
for?i=1:gennum
for?s=1:m
K（is）=2*para（i1）*（2*s-1）*MaxPs（i11）+para（i2）;%推導(dǎo)簡(jiǎn)化后的煤耗斜率
end
end
?%目標(biāo)函數(shù)線性化后分段出力的等式約束
for?i=1:gennum?
????for?t=1:T
st=st+[p（it）==（sum（Ps（it:）3）+u（it）*limit（i2））];
????end
end
%%?目標(biāo)函數(shù)
totalcost=0;%機(jī)組費(fèi)用成本最小
%線性化的最優(yōu)成本目標(biāo)
for?i=1:gennum
for?t=1:T
for?s=1:m
????totalcost=totalcost+K（is）*Ps（its）;%線性化煤耗成本
end
????totalcost=totalcost+u（it）*（para（i2）*limit（i2）+para（i1）*limit（i2）^2+para（i3））;%加上表示機(jī)組開機(jī)并以最小出力?運(yùn)行產(chǎn)生的煤耗
????totalcost=totalcost+costH（it）+costJ（it）;%加上機(jī)組啟停產(chǎn)生的開停機(jī)成本
end
end
%原二次函數(shù)式的最優(yōu)成本目標(biāo)
%?for?i=1:gennum
%?????for?t=1:T
%?????totalcost=totalcost+para（i1）*p（it）.^2+para（i2）*p（it）+para（i3）*u（it）;??%煤耗成本
%?????totalcost=totalcost+costH（it）;????????????????????????????????%啟動(dòng)成本
%?????totalcost=totalcost+costJ（it）;????????????????????????????????%關(guān)停成本
%?????end
%?end
%%
for?t=1:T
st=st+[sum（p（:t））==PL（1t）];%負(fù)荷平衡約束;
end
%%
for?t=1:T
????for?i=1:gennum
??st=st+[u（it）*limit（i2）<=p（it）<=u（it）*limit（i1）];%機(jī)組出力上下限約束
????end
end
%%?機(jī)組爬坡約束
%按下式進(jìn)行推導(dǎo)編程
%?%啟動(dòng)最大升速率
%?Su=（Pmax+Pmin）/2;
%?%停機(jī)最大降速率
%?Sd=（Pmax+Pmin）/2;
%Ru=Rud;Rd=Rud;
%?%上爬坡約束
%?for?t=2:T
%?st=st+[p（:t）-p（:t-1）<=u（:t-1）.*（Ru-Su）+Su];
%?end
%?%下爬坡約束
%?for?t=2:T
%st=st+[p（:t-1）-p（:t）<=u（:t）.*（Rd-Sd）+Sd];
%?end
%展開表達(dá)式：
for?t=2:T
????for?i=1:gennum
????%?st=st+[-Rud（i1）*u（it）+（u（it）-u（it-1））*limit（i2）-limit（i1）*（1-u（it））<=p（it）-p（it-1）];
????%?st=st+[p（it）-p（it-1）<=Rud（i1）*u（it-1）+（u（it）-u（it-1））*limit（i2）+limit（i1）*（1-u（it））];
????%由于原式可能關(guān)機(jī)以后就無(wú)法再開動(dòng)了，改用下式
????st=st+[p（it-1）-p（it）<=Rud（i1）*u（it）+（1-u（it））*（limit（i2）+l

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????53760??2017-11-13?19:03??電力系統(tǒng)機(jī)組組合優(yōu)化\excel2017.xls

?????文件???????7584??2017-11-24?21:46??電力系統(tǒng)機(jī)組組合優(yōu)化\jizuzuheyouhua.m

?????文件?????184428??2017-11-26?13:35??電力系統(tǒng)機(jī)組組合優(yōu)化\基本要求.docx

?????文件??????34092??2017-11-24?21:53??電力系統(tǒng)機(jī)組組合優(yōu)化\熱備用0.05機(jī)組各時(shí)段最優(yōu)出力圖表.vsdx

?????文件??????69120??2017-11-24?22:16??電力系統(tǒng)機(jī)組組合優(yōu)化\熱備用0.05狀態(tài)下的機(jī)組組合問(wèn)題求解結(jié)果.xls

?????文件??????36685??2017-11-24?22:23??電力系統(tǒng)機(jī)組組合優(yōu)化\熱備用0.2下的機(jī)組最優(yōu)出力.vsdx

?????文件??????69632??2017-11-24?22:21??電力系統(tǒng)機(jī)組組合優(yōu)化\熱備用0.2狀態(tài)下的機(jī)組組合問(wèn)題求解結(jié)果.xls

?????目錄??????????0??2019-04-19?19:20??電力系統(tǒng)機(jī)組組合優(yōu)化

-----------?---------??----------?-----??----

???????????????455301????????????????????8