风电光伏蓄电池最优组合运行代码.rar

风电光伏蓄电池最优组合运行代码

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