% main5_2
clc
clear
gen=zeros(5,12);
price_e=zeros(2,24);
y=zeros(5,24);
gen = [
1 335 125 1169.88 2.6998 0.00753 10 -1 -1 14 14 500;
2 232 150 779.64 0.9286 0.0156 10 -1 -1 10 10 300;
3 260 50 636.46 3.826 0.0089 10 -1 -1 15 15 500;
4 440 160 669.12 7.9215 0.0000017 10 -1 -1 25 25 603;
5 250 130 590.72 2.4435 0.00906 10 -1 -1 12 12 500;
];
price_f=3;
%%energy price 各时段的电价
%T price_energy($/MWh)
price_p=[ 20 20 20 20 20 24 26 28 30 32 34 36 38 40 42 48 48 48 48 48 48 48 48 20;
20 20 20 20 20 24 26 28 30 32 34 36 38 40 42 48 48 48 48 48 48 48 48 20;
20 20 20 20 20 24 26 28 30 32 34 36 38 40 42 48 48 48 48 48 48 48 48 20;
20 20 20 20 20 24 26 28 30 32 34 36 38 40 42 48 48 48 48 48 48 20 20 20;
20 20 20 20 20 24 26 28 30 32 34 36 38 40 42 48 48 48 48 48 48 20 20 20;
];
Ng=5;
NT=24;
%机组初始状态
statini=gen(:,9);
Lon=zeros(5,24);
for i=1:Ng
if statini(i)<0
statini(i)=0;
else
statini(i)=1;
end
a=price_f.*gen(i,4);
b=price_f.*gen(i,5);
c=price_f.*gen(i,6);
fst=gen(i,12);%启动燃料
Pgmax=gen(i,2);
Pgmin=gen(i,3);
ep=0.001;
lam_P=30*ones(Ng,NT); %拉格朗日乘子
P(i,:)=(price_p(i,:)-lam_P(i,:)-b)./c/2;
P1(i,:)=(price_p(i,:)-b)./c/2;
k=1;
for n=1:NT
%修正系数——机组P的约束
while(k>ep)
step=1/sum(0.5./c);
x=zeros(1,NT);
Stat(i,n)=1;
if P(i,n)>Pgmax
lam_P(i,n)=lam_P(i,n)+step*(P(i,n)-Pgmax);
x(n)=1;
end
if P(i,n)<Pgmin
lam_P(i,n)=lam_P(i,n)+step*(P(i,n)-Pgmin);
x(n)=1;
end
P(i,n)=(price_p(i,n)-lam_P(i,n)-b)./c/2;
P1(i,n)=(price_p(i,n)-b)./c/2;
Fon(i,n)=(price_p(i,n)*P(i,n)-a-b.*P(i,n)-c.*P(i,n).^2)*Stat(i,n)+lam_P(i,n)*P(i,n)*Stat(i,n);%基于机组组合的拉格朗日值
Lon(i,n)=(price_p(i,n)*P1(i,n)-a-b.*P1(i,n)-c.*P1(i,n)^2)*Stat(i,n)+lam_P(i,n)*P1(i,n)*Stat(i,n);%基于经济分配的拉格朗日值
k=abs((Lon(i,n)-Fon(i,n))/Lon(i,n));
if (P(i,n)<Pgmax)|(P(i,n)>Pgmin)|(P(i,n)==Pgmax)|(P(i,n)==Pgmin)
break
end
end
y(i,n)=price_p(i,n)*P(i,n)-a-b.*P(i,n)-c.*P(i,n).^2;%机组的盈利
if y(i,n)>0
Stat(i,n)=1;
else if y(i,n)<0
Stat(i,n)=0;
end
end
end
end
y
Stat
P
Fon
评论0