資源簡介
風電光伏蓄電池最優組合運行代碼

代碼片段和文件信息
function?F=ff(kx)
global?CW?CS?CB?CR?PEENS?T?PL?PW?Ppv?EBmin?EB?EB0?PB?fyz?h
e=14500;%WG單價
f=1100;%PV單價
l=580;%B單價
r0=0.0022;
mm=10;
uw=160;
upv=180;
ub=30;
CW(k)=0;
CW(k)=CW(k)+x(k1)*(e*r0*power(1+r0mm)/(power(1+r0mm)-1)+uw);
CS(k)=0;
CS(k)=CS(k)+2*x(k2)*(f*r0*power(1+r0mm)/(power(1+r0mm)-1)+upv);
CB(k)=0;
CB(k)=CB(k)+2*x(k3)*(l*r0*power(1+r0mm)/(power(1+r0mm)-1)+ub);
b=25;
ee=0.4883;%單位電價
kcoe=b*ee;%賠償系數
%計算PEENS;
t=1;
PEENS0=PL(t)-PW(kt)-Ppv(kt)-(EB0-EBmin);
if?PEENS0<0
????PEENS0=0;
end
PEENS(k)=PEENS0;
for?t=2:T
????PEE(kt)=PL(t)-PW(kt)-Ppv(kt)-(EB(kt-1)-EBmin);
????if?PEE(kt)<0
????????PEE(kt)=0;
????end
????PEENS(k)=PEENS(k)+PEE(kt);
end
CR(k)=kcoe*PEENS(k);
%能量平衡約束罰函數
h(k)=abs(EB(k1)-EB0+PB(k1));
for?t=2:T
????h(k)=h(k)+abs(EB(kt)-EB(kt-1)+PB(kt));
end
F=CW(k)+CS(k)+CB(k)+CR(k)+fyz*h(k);
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????924??2010-05-24?12:59??程序\ff.m
?????文件???????9831??2010-05-24?13:27??程序\zhuchengxu.m
?????目錄??????????0??2010-06-23?16:38??程序
-----------?---------??----------?-----??----
????????????????10755????????????????????3
- 上一篇:labview 自動售貨機
- 下一篇:自考計算機網絡名詞解釋大集合
評論
共有 條評論