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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????53760??2017-11-13?19:03??電力系統機組組合優化\excel2017.xls

?????文件???????7584??2017-11-24?21:46??電力系統機組組合優化\jizuzuheyouhua.m

?????文件?????184428??2017-11-26?13:35??電力系統機組組合優化\基本要求.docx

?????文件??????34092??2017-11-24?21:53??電力系統機組組合優化\熱備用0.05機組各時段最優出力圖表.vsdx

?????文件??????69120??2017-11-24?22:16??電力系統機組組合優化\熱備用0.05狀態下的機組組合問題求解結果.xls

?????文件??????36685??2017-11-24?22:23??電力系統機組組合優化\熱備用0.2下的機組最優出力.vsdx

?????文件??????69632??2017-11-24?22:21??電力系統機組組合優化\熱備用0.2狀態下的機組組合問題求解結果.xls

?????目錄??????????0??2019-04-19?19:20??電力系統機組組合優化

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

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