function?[leastCostPt]=ecoDisp（coeffPminPmaxUPload）

%?本函數(shù)用于求解某個時段內(nèi)各投入機(jī)組的最優(yōu)發(fā)電功率和最小損耗=====================>>經(jīng)濟(jì)調(diào)度

%?leastCost是該時段內(nèi)所有機(jī)組的最小發(fā)電成本========>>經(jīng)濟(jì)調(diào)度
%?Pt為該時段內(nèi)所有機(jī)組的最優(yōu)發(fā)電功率

%?coeff為所有機(jī)組的成本系數(shù)矩陣，設(shè)為二次函數(shù)形式[??二次項系數(shù)???一次項系數(shù)????常數(shù)項系數(shù)????]
%?在此假定發(fā)電機(jī)成本函數(shù)為二次函數(shù)，coeff=[c2c1c0]為發(fā)電機(jī)成本函數(shù)的系數(shù)
%?PminPmax為所有機(jī)組的發(fā)電功率下限和上限向量（列向量）
%?U?為由動態(tài)規(guī)劃求解得到的所有機(jī)組在該時段內(nèi)的運(yùn)行狀態(tài)（列向量）
%?Pload為該時段內(nèi)的負(fù)荷（標(biāo)量）

n=size（coeff1）;%?n?為機(jī)組數(shù)包括未投運(yùn)的機(jī)組和已投運(yùn)機(jī)組
c2=zeros（n1）;%?二次項系數(shù)
c1=zeros（n1）;%?一次項系數(shù)
c0=zeros（n1）;%?常數(shù)項系數(shù)
c2=coeff（:1）;%?二次項系數(shù)
c1=coeff（:2）;%?一次項系數(shù)
c0=coeff（:3）;%?常數(shù)項系數(shù)

leastCost=0;
Pt=zeros（n1）;%各機(jī)組在該時段內(nèi)的發(fā)電功率

lambda=10;?%?設(shè)定拉格朗日乘子的初值（可以改變）
epsilon=0.1;?%?允許的誤差
lostPower=0.1;%?負(fù)荷平衡約束條件的偏差
sumPower=0;?%?各機(jī)組的發(fā)電功率之和

maxPower=0;%?投運(yùn)機(jī)組的最大發(fā)電功率之和
for?i=1:n
???maxPower?=?maxPower?+?U（i）*Pmax（i）;?%?計算投運(yùn)機(jī)組的最大發(fā)電功率
end??

if（?Pload<=maxPower）??????????????????%?如果系統(tǒng)負(fù)荷大于投運(yùn)機(jī)組的最大發(fā)電功率，則不必進(jìn)行計算
???whil