%function?[leastCostPt?U]=UC（?coeff?Pmin?Pmax?Ploadlambda）

%?leastCost為機(jī)組的最小發(fā)電成本，是某一時間段內(nèi)所有機(jī)組的發(fā)電成本=========>>經(jīng)濟(jì)調(diào)度
%?Pt為機(jī)組在各時段的發(fā)電功率，行對應(yīng)一個機(jī)組，列對應(yīng)某一時間段。最終需要的是Pt而不是P因?yàn)镻t才是經(jīng)濟(jì)調(diào)度的結(jié)果，滿足負(fù)荷平衡約束
%?U為機(jī)組在各時段的運(yùn)行狀態(tài)，行對應(yīng)一個機(jī)組，列對應(yīng)某一時間段。若其值為1，則表示機(jī)組投運(yùn)；若值為0，則表示機(jī)組停機(jī)

%?coeff為各機(jī)組的成本系數(shù)矩陣，設(shè)為二次函數(shù)形式[??二次項(xiàng)系數(shù)???一次項(xiàng)系數(shù)????常數(shù)項(xiàng)系數(shù)????]
%?PminPmax為各機(jī)組的發(fā)電功率下限和上限向量（列向量）
%?Pload為各時段的負(fù)荷向量（行向量）
%?lambda為拉格朗日系數(shù)初值

coeff=[0.00128??6.48??459;
???????0.00194??7.85??310;
???????0.00482??7.97??78];?%?每一行對應(yīng)一個機(jī)組
Pmin=[150?;100?;50];
Pmax=[600?;400?;200];
Pload=[850];
n=size（coeff1）;%?n?為機(jī)組數(shù)
T=size（Pload2）;%?T?為需要計算的總時間段數(shù)
lambda=zeros（1T）;
??????
%????????num=nargin;
%????????if（?num==6）
%??????????????lambda=zeros（1T）;
%????????end
???????
???????
???????q=0;?????????????????%?對偶問題
???????minCost=zeros（n1）;??%對應(yīng)每個機(jī)組在所有時間段內(nèi)的發(fā)電成本=========>>動態(tài)規(guī)劃
???????P=zeros（nT）;????????%對應(yīng)每個機(jī)組=============>>動態(tài)規(guī)劃
???????U=zeros（nT）;????????%?每個機(jī)組在所有時間段內(nèi)的運(yùn)行狀態(tài)，若其值為1，則表示機(jī)組投運(yùn)；若值為0，則表示機(jī)組停機(jī)
???????lostPower=zeros（1T）;%?各時段內(nèi)的系統(tǒng)發(fā)電功率缺額=========>>動態(tài)規(guī)劃結(jié)束后
???????
???????J=0;?????????????????%?原問題
???????leastCost=zeros（1T）;%對應(yīng)每個時間段內(nèi)所有機(jī)組的發(fā)電成本========>>經(jīng)濟(jì)調(diào)度
???????Pt=zeros（nT）;???????%對應(yīng)每個時間段==========>>經(jīng)濟(jì)調(diào)度
???????lackCost=10000;??????%?當(dāng)機(jī)組容量不能滿足經(jīng)濟(jì)調(diào)度時，強(qiáng)制設(shè)定此時的總發(fā)電成本為10000
???????flag=0;
???????
???????iter=1;??????????????%?計算迭代次數(shù)
???????epsilon=0.05;????????%?設(shè)定容許誤差
???????gap=0.05;????????????%?設(shè)定對偶間隙初值
????????????????????
?????while（?abs（gap）>=epsilon?）?
???
??????????for?i=1:n???%?n為機(jī)組數(shù)
????????????????[?minCost（i）P（i:）U（i:）?]?=?dp（?coeff（i:）Pmin（i）Pmax（i）Tlambda）;%?利用動態(tài)規(guī)劃求解每個機(jī)組在所有時段的運(yùn)行狀態(tài)和最優(yōu)發(fā)電功率
??????????end
??????????%=========================================================================>>
??????????%?本部分用于求對偶問題最優(yōu)值q
??????????sumPower=zeros（1T）;%?系統(tǒng)各時段的總發(fā)電功率
??????????for?t=1:T
??????????????for?i=1:n
??????????????????sumPower（t）=sumPower（t）?+?U（it）*P（it）;%?求各時段內(nèi)的發(fā)電功率
???????????????end
???????????end????
??????????lostPower=Pload-sumPower;%?求得各時段的系統(tǒng)發(fā)電功率缺額
??????????
??????????q=0;
??????????for?i=1:n
???????????????q?=?q?+?minCost（i）;???%?疊加各機(jī)組的最小發(fā)電成本
??????????end
??????????for?t=1:T
???????????????q?=?q?+?lambda（t）*lostPower（t）;?%?疊加拉格朗日乘子項(xiàng)???
??????????end
??????????%=========================================================================>>
??????????%本部分用于求原問題最優(yōu)值J
??????????????for?t=1:T
????????????????????for?i=1:n
?????????????????????????if?U（it）==1???%??判斷某時段內(nèi)是否有機(jī)組投運(yùn)
?????????????????????????????flag=1;
?????????????????????????end????
????????????????????end??
????????????????????if?flag==1???%?如果所有機(jī)組停運(yùn)則不必求解
?????????????????????????[leastCost（t）Pt（:t）]?=?ecoDisp