function?[leastCostPt]=ecoDisp（coeffPminPmaxUPload）

%?本函數用于求解某個時段內各投入機組的最優發電功率和最小損耗=====================>>經濟調度

%?leastCost是該時段內所有機組的最小發電成本========>>經濟調度
%?Pt為該時段內所有機組的最優發電功率

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

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

leastCost=0;
Pt=zeros（n1）;%各機組在該時段內的發電功率

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

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

if（?Pload<=maxPower）??????????????????%?如果系統負荷大于投運機組的最大發電功率，則不必進行計算
???whil