请你修改以下代码，并给出修改后的完整代码：clc
clear
tic %tic用来保存当前时间，而后使用toc来记录
程序完成时间
%% 参数设置
%---燃气轮机2台、余热回收锅炉1台、吸收式制冷机1台---%
PGmax=[1200,2800]; %燃气发电机最大电功率
PGmin=[300,300]; %燃气发电机最小电功率
PGnom=[1200,2800]; %燃气发电机额定电功率
%这里的Aa,Bb,Cc,Dd分别对应四台不同的燃气发电机的发电效率
Aa=8.93
5; Bb=33.1
57; Cc=-27.081; Dd=17.989;
QHrsnom=1000; %余热回收
系统的额定功率
QAcnom=400; %吸收式制冷机的额定功率
yHrs=0.60; %热回收效率
Copac=0.80; %吸收制冷能效比
%---燃气锅炉1台、电制冷机1台---%
Qbolnom=1000; %燃气锅炉额定功率
QEcnom=300; %电制冷机的额定功率
Copec=3; %电制冷能效比
%----上网购电参数-------%
Pgridmax=600;
Pgridmin=400;
% ----------------蓄电池组参数-------------------%
SBatmax=180; %蓄电池保持稳定运行时的最大容量
SBatmin=40; %蓄电池保持稳定运行时的最小容量
SBatnom=200; %蓄电池在协调周期内的额定容量
RBatcha=0.20;%蓄电池组l在t时段的最大充电率
RBatdis=0.40;%蓄电池组l在t时段的最大放电率
dBat=0.04; %蓄电池组的自放电率
yBatcha=0.9
5; yBatdis=0.9
5; %充放电效率
% 风力发电机的功率
Pwind=[693,718,883,698,888,683,723,69
5,678,721,
568,
583,614,693,678,703,692,693,808,994,869,813,873,808];
% 光伏发电机的功率
Ppv=[0,0,0,0,0,0,16
50,24
50,32
50,33
50,3400,37
50,34
50,32
50,3200,2400,2100,1300,0,0,0,0,0,0];
% 负荷功率
PD=[2800,2700,3000,3800,4600,4600,
5200,
5400,
5800,6300,7400,8700,9700,10000,10100,10300,9000,7000,6700,
5900,4
500,3000,2700,2800]; %为一天中不同时段的
系统需要的预测电功率
QDtherm=[1300,1400,1360,1700,1600,2000,3100,34
50,3900,4400,4300,4
500,
5400,
5000,
5100,49
50,4600,4710,4400,4000,36
50,3
570,2900,2420]; %
系统在一天中需要的热负荷
QDcool=[23
50,2800,28
50,3400,3800,4000,4700,43
50,
5100,7
500,7600,84
50,8700,80
50,7700,74
50,7200,66
50,7300,
5700,
5100,44
50,3600,3080]; %
系统在一天中需要的冷负荷
%% 粒子群算法参数设置
w=0.7; %惯性权重
c1=2; %个体
c2=2; %群体
N=100; %粒子个数
SBat=zeros
(N,2
5); %蓄电池
QAc=zeros
(N,24
); %吸收式制冷机
QEc=zeros
(N,24
); %电制冷机
QHrs=zeros
(N,24
); %余热回收锅炉
yg=zeros
(N,48
); %燃气轮机
k=1;
D=8;
Tmax=1000; %最大迭代次数
%% 主循环开始
for i=1:N
SBat
(i,19
)=63.
5209380080000; %初始储能量，单位为kW.h,
end
for i=1:N %分时段初始化满足约束的初始值
for t=1:24
%燃气轮机出力初始化
for j=1:2
X
(i,
(t-1
)*D
+j
)=unifrnd
(PGmin
(j
),PGmax
(j
));
V
(i,
(t-1
)*D
+j
)=unifrnd
(PGmin
(j
)-PGmax
(j
),PGmax
(j
)-PGmin
(j
));
end
%余热回收锅炉提供热出力的初始化
for j=1:2
yg
(i,
(t-1
)*2
+j
)=
(Aa
+Bb*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))+Cc*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^2
+Dd*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^3
)/100;
end
QHrs(i,t
)=
(X
(i,
(t-1
)*D
+1
)*
(1-yg
(i,
(t-1
)*2
+1
))+X
(i,
(t-1
)*D
+2
)*
(1-yg
(i,
(t-1
)*2
+2
)))*yHrs; %热回收
系统在单时段提供的总热功率
if
QHrs(i,t
)>1000
QHrs(i,t
)=1000;
end
if
QHrs(i,t
)-QDtherm
(t
)>=0
X
(i,
(t-1
)*D
+3
)=0;
V
(i,
(t-1
)*D
+3
)=0;
else
X
(i,
(t-1
)*D
+3
)=unifrnd
(abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
))),QDtherm
(t
));
V
(i,
(t-1
)*D
+3
)=unifrnd
(abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
)))-QDtherm
(t
),QDtherm
(t
)-abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
))));
end
%燃气锅炉提供的热功率用于冷负荷部分量的初始化
X
(i,
(t-1
)*D
+4
)=unifrnd
(0,min
((Qbolnom-X
(i,
(t-1
)*D
+3
)),QDcool
(t
)/Copac
));
V
(i,
(t-1
)*D
+4
)=unifrnd
(-min
((Qbolnom-X
(i,
(t-1
)*D
+3
)),QDcool
(t
)/Copac
),min
((Qbolnom-X
(i,
(t-1
)*D
+3
)),QDcool
(t
)/Copac
));
%热回收
系统用于提供冷功率初始化
for j=1:2
yg
(i,
(t-1
)*2
+j
)=
(Aa
+Bb*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))+Cc*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^2
+Dd*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^3
)/100;
end
QHrs(i,t
)=
(X
(i,
(t-1
)*D
+1
)*
(1-yg
(i,
(t-1
)*2
+1
))+X
(i,
(t-1
)*D
+2
)*
(1-yg
(i,
(t-1
)*2
+2
)))*yHrs; %热回收
系统在单时段提供的总热功率
if
QHrs(i,t
)>
QHrsnom
QHrs(i,t
)=
QHrsnom;
end
X
(i,
(t-1
)*D
+7
)=min
(max
(0,
QHrs(i,t
)-QDtherm
(t
)),QDcool
(t
)/Copac
);
V
(i,
(t-1
)*D
+7
)=0;
%电制冷机制冷所消耗的电功率初始化
QAc
(i,t
)=X
(i,
(t-1
)*D
+7
)*Copac; %吸收制冷机在单时段提供的冷功率
X
(i,
(t-1
)*D
+8
)=max
(0,min
((QDcool
(t
)-QAc
(i,t
))/Copec,2*QEcnom/Copec
));
V
(i,
(t-1
)*D
+8
)=0;
X
(i,
(t-1
)*D
+4
)=abs
(min
(0,
(QAc
(i,t
)+ X
(i,
(t-1
)*D
+8
)*Copec-QDcool
(t
))));
V
(i,
(t-1
)*D
+4
)=0;
%微网在单个时段内与电网交互的电量初始化
X
(i,
(t-1
)*D
+5)=unifrnd
(-Pgridmin,Pgridmax
);
V
(i,
(t-1
)*D
+5)=unifrnd
(-Pgridmin-Pgridmax,Pgridmax
+Pgridmin
);
%蓄电池组的充放电功率初始化
X
(i,
(t-1
)*D
+6
)=unifrnd
(max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
)),min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
)));
V
(i,
(t-1
)*D
+6
)=unifrnd
(max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
))-min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
)),min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
))-max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
)));
%--------------------%
if X
(i,
(t-1
)*D
+6
)>=0
SBat
(i,t
+1
)= SBat
(i,t
)*
(1-dBat
)+X
(i,
(t-1
)*D
+6
)*yBatcha; %蓄电池组在单时段时的剩余容量
else
SBat
(i,t
+1
)= SBat
(i,t
)*
(1-dBat
)+X
(i,
(t-1
)*D
+6
)/yBatdis; %蓄电池组在单时段时的剩余容量
end
if SBat
(i,t
+1
)>SBatmax
SBat
(i,t
+1
)=SBatmax;
end
if SBat
(i,t
+1
)<SBatmin
SBat
(i,t
+1
)=SBatmin;
end
end
end
WEI=24*8;
for i=1:N
p
(i
)=fitness
(X
(i,:
),WEI
);
y
(i,:
)=X
(i,:
);
end
pg = X
(N,:
); %Pg为全局最优
for i=1:
(N-1
)
if fitness
(X
(i,:
),WEI
)<fitness
(pg,WEI
)
pg=X
(i,:
);
end
end
for K=1:Tmax
if ~mod
(K,10
) %K与10作除法后的余数
fprintf
('当前迭代次数:%d\n',K
);
end
for i=1:N
V
(i,1:WEI
)=w*V
(i,1:WEI
)+c1*rand*
(y
(i,1:WEI
)-X
(i,1:WEI
))+c2*rand*
(pg-X
(i,1:WEI
)); %粒子速度更新
for t=1:24
%燃气轮机出力速度越限处理
for j=1:2
if V
(i,
(t-1
)*D
+j
)>0.
5*
(PGmax
(j
)-PGmin
(j
))
V
(i,
(t-1
)*D
+j
)=0.
5*
(PGmax
(j
)-PGmin
(j
));
end
if V
(i,
(t-1
)*D
+j
)<0.
5*
(PGmin
(j
)-PGmax
(j
))
V
(i,
(t-1
)*D
+j
)=0.
5*
(PGmin
(j
)-PGmax
(j
));
end
end
%燃气锅炉提供的热功率用于热负荷部分量的速度越限处理
for j=1:2
yg
(i,
(t-1
)*2
+j
)=
(Aa
+Bb*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))+Cc*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^2
+Dd*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^3
)/100;
end
QHrs(i,t
)=
(X
(i,
(t-1
)*D
+1
)*
(1-yg
(i,
(t-1
)*2
+1
))+X
(i,
(t-1
)*D
+2
)*
(1-yg
(i,
(t-1
)*2
+2
)))*yHrs; %热回收
系统在单时段提供的总热功率
if
QHrs(i,t
)>
QHrsnom
QHrs(i,t
)=
QHrsnom;
end
if
QHrs(i,t
)-QDtherm
(t
)>=0
if V
(i,
(t-1
)*D
+3
)>0
V
(i,
(t-1
)*D
+3
)=0;
end
if V
(i,
(t-1
)*D
+3
)<0
V
(i,
(t-1
)*D
+3
)=0;
end
else
if V
(i,
(t-1
)*D
+3
)>0.
5*
(QDtherm
(t
)-abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
))))
V
(i,
(t-1
)*D
+3
)=0.
5*
(QDtherm
(t
)-abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
))));
end
if V
(i,
(t-1
)*D
+3
)<0.
5*
(abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
)))-QDtherm
(t
))
V
(i,
(t-1
)*D
+3
)=0.
5*
(abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
)))-QDtherm
(t
));
end
end
%锅炉提供的热功率用于冷负荷部分量的速度越限处理
if V
(i,
(t-1
)*D
+4
)>0
V
(i,
(t-1
)*D
+4
)=0;
end
if V
(i,
(t-1
)*D
+4
)<0
V
(i,
(t-1
)*D
+4
)=0;
end
%微网在单个时段与电网交互的电量的速度越限处理
if V
(i,
(t-1
)*D
+5)>0.
5*
(Pgridmax
+Pgridmin
)
V
(i,
(t-1
)*D
+5)=0.
5*
(Pgridmax
+Pgridmin
);
end
if V
(i,
(t-1
)*D
+5)<0.
5*
(-Pgridmin-Pgridmax
)
V
(i,
(t-1
)*D
+5)=0.
5*
(-Pgridmin-Pgridmax
);
end
%蓄电池组的充放电功率的速度越限处理
if V
(i,
(t-1
)*D
+6
)>0.
5*
(min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
))-max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
)))
V
(i,
(t-1
)*D
+6
)=0.
5*
(min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
))-max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
)));
end
if V
(i,
(t-1
)*D
+6
)<0.
5*
(max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
))-min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
)))
V
(i,
(t-1
)*D
+6
)=0.
5*
(max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
))-min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
)));
end
%热回收
系统用于提供冷功率的热功率的速度越限处理
if V
(i,
(t-1
)*D
+7
)>0
V
(i,
(t-1
)*D
+7
)=0;
end
if V
(i,
(t-1
)*D
+7
)<0
V
(i,
(t-1
)*D
+7
)=0;
end
%电制冷机制冷所消耗的电功率的速度越限处理
QAc
(i,t
)=X
(i,
(t-1
)*D
+7
)*Copac; %吸收制冷机在单时段提供的冷功率
if V
(i,
(t-1
)*D
+8
)>0
V
(i,
(t-1
)*D
+8
)=0;
end
if V
(i,
(t-1
)*D
+8
)<0
V
(i,
(t-1
)*D
+8
)=0;
end
end
X
(i,1:WEI
)=X
(i,1:WEI
)+V
(i,1:WEI
); %粒子位置更新
for t=1:24
%燃气发电机出力位置越限处理
for j=1:2
if X
(i,
(t-1
)*D
+j
)>PGmax
(j
)
X
(i,
(t-1
)*D
+j
)=PGmax
(j
);
V
(i,
(t-1
)*D
+j
)=-rand*V
(i,
(t-1
)*D
+j
);
end
if X
(i,
(t-1
)*D
+j
)<=PGmin
(j
)
X
(i,
(t-1
)*D
+j
)=0;
V
(i,
(t-1
)*D
+j
)=0;
end
end
%锅炉提供的热功率用于热负荷部分量的位置越限处理
for j=1:2
yg
(i,
(t-1
)*2
+j
)=
(Aa
+Bb*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))+Cc*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^2
+Dd*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^3
)/100;
end
QHrs(i,t
)=
(X
(i,
(t-1
)*D
+1
)*
(1-yg
(i,
(t-1
)*2
+1
))+X
(i,
(t-1
)*D
+2
)*
(1-yg
(i,
(t-1
)*2
+2
)))*yHrs; %热回收
系统在单时段提供的总热功率
if
QHrs(i,t
)>
QHrsnom
QHrs(i,t
)=
QHrsnom;
end
if
QHrs(i,t
)-QDtherm
(t
)>=0
if X
(i,
(t-1
)*D
+3
)>0
X
(i,
(t-1
)*D
+3
)=0;
V
(i,
(t-1
)*D
+3
)=-rand*V
(i,
(t-1
)*D
+3
);
end
if X
(i,
(t-1
)*D
+3
)<0
X
(i,
(t-1
)*D
+3
)=0;
V
(i,
(t-1
)*D
+3
)=-rand*V
(i,
(t-1
)*D
+3
);
end
else
if X
(i,
(t-1
)*D
+3
)>QDtherm
(t
)
X
(i,
(t-1
)*D
+3
)=QDtherm
(t
);
V
(i,
(t-1
)*D
+3
)=-rand*V
(i,
(t-1
)*D
+3
);
end
if X
(i,
(t-1
)*D
+3
)<abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
)))
X
(i,
(t-1
)*D
+3
)=abs
(min
(0,
QHrs(i,t
)-QDtherm
(t
)));
V
(i,
(t-1
)*D
+3
)=-rand*V
(i,
(t-1
)*D
+3
);
end
end
%锅炉提供的热功率用于冷负荷部分量的位置越限处理
if X
(i,
(t-1
)*D
+4
)>min
((Qbolnom-X
(i,
(t-1
)*D
+3
)),QDcool
(t
)/Copac
)
X
(i,
(t-1
)*D
+4
)=min
((Qbolnom-X
(i,
(t-1
)*D
+3
)),QDcool
(t
)/Copac
);
V
(i,
(t-1
)*D
+4
)=-rand*V
(i,
(t-1
)*D
+4
);
end
if X
(i,
(t-1
)*D
+4
)<0
X
(i,
(t-1
)*D
+4
)=0;
V
(i,
(t-1
)*D
+4
)=-rand*V
(i,
(t-1
)*D
+4
);
end
%微网在单个时段内与电网交互的电量的位置越限处理
if X
(i,
(t-1
)*D
+5)>Pgridmax
X
(i,
(t-1
)*D
+5)=Pgridmax;
V
(i,
(t-1
)*D
+5)=-rand*V
(i,
(t-1
)*D
+5);
end
if X
(i,
(t-1
)*D
+5)<
(-Pgridmin
)
X
(i,
(t-1
)*D
+5)=
(-Pgridmin
);
V
(i,
(t-1
)*D
+5)=-rand*V
(i,
(t-1
)*D
+5);
end
%蓄电池组的充电功率的位置越限处理
if X
(i,
(t-1
)*D
+6
)>min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
))
X
(i,
(t-1
)*D
+6
)=min
(SBatnom*RBatcha,SBatmax-SBat
(i,t
));
V
(i,
(t-1
)*D
+6
)=-rand*V
(i,
(t-1
)*D
+6
);
end
if X
(i,
(t-1
)*D
+6
)<max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
))
X
(i,
(t-1
)*D
+6
)=max
(-SBatnom*RBatdis,SBatmin-SBat
(i,t
));
V
(i,
(t-1
)*D
+6
)=-rand*V
(i,
(t-1
)*D
+6
);
end
%热回收
系统用于提供冷功率的热功率的位置越限处理
for j=1:2
yg
(i,
(t-1
)*2
+j
)=
(Aa
+Bb*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))+Cc*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^2
+Dd*
(X
(i,
(t-1
)*D
+j
)/PGnom
(j
))^3
)/100;
end
QHrs(i,t
)=
(X
(i,
(t-1
)*D
+1
)*
(1-yg
(i,
(t-1
)*2
+1
))+X
(i,
(t-1
)*D
+2
)*
(1-yg
(i,
(t-1
)*2
+2
)))*yHrs; %热回收
系统在单时段提供的总热功率
if
QHrs(i,t
)>
QHrsnom
QHrs(i,t
)=
QHrsnom;
end
if X
(i,
(t-1
)*D
+7
)>min
(max
(0,
QHrs(i,t
)-QDtherm
(t
)),QDcool
(t
)/Copac
)
X
(i,
(t-1
)*D
+7
)=min
(max
(0,
QHrs(i,t
)-QDtherm
(t
)),QDcool
(t
)/Copac
);
end
if X
(i,
(t-1
)*D
+7
)<min
(max
(0,
QHrs(i,t
)-QDtherm
(t
)),QDcool
(t
)/Copac
)
X
(i,
(t-1
)*D
+7
)=min
(max
(0,
QHrs(i,t
)-QDtherm
(t
)),QDcool
(t
)/Copac
);
end
%电制冷机制冷所消耗的电功率的位置越限处理
QAc
(i,t
)=X
(i,
(t-1
)*D
+7
)*Copac; %吸收制冷机在单时段提供的冷功率
if X
(i,
(t-1
)*D
+8
)>max
(0,min
((QDcool
(t
)-QAc
(i,t
))/Copec,2*QEcnom/Copec
))
X
(i,
(t-1
)*D
+8
)=max
(0,min
((QDcool
(t
)-QAc
(i,t
))/Copec,2*QEcnom/Copec
));
V
(i,
(t-1
)*D
+8
)=-rand*V
(i,
(t-1
)*D
+8
);
end
if X
(i,
(t-1
)*D
+8
)<max
(0,min
((QDcool
(t
)-QAc
(i,t
))/Copec,2*QEcnom/Copec
))
X
(i,
(t-1
)*D
+8
)=max
(0,min
((QDcool
(t
)-QAc
(i,t
))/Copec,2*QEcnom/Copec
));
V
(i,
(t-1
)*D
+8
)=-rand*V
(i,
(t-1
)*D
+8
);
end
if X
(i,
(t-1
)*D
+4
)>abs
(min
(0,
(QAc
(i,t
)+ X
(i,
(t-1
)*D
+8
)*Copec-QDcool
(t
))))
X
(i,
(t-1
)*D
+4
)=abs
(min
(0,
(QAc
(i,t
)+ X
(i,
(t-1
)*D
+8
)*Copec-QDcool
(t
))));
end
if X
(i,
(t-1
)*D
+4
)<abs
(min
(0,
(QAc
(i,t
)+ X
(i,
(t-1
)*D
+8
)*Copec-QDcool
(t
))))
X
(i,
(t-1
)*D
+4
)=abs
(min
(0,
(QAc
(i,t
)+ X
(i,
(t-1
)*D
+8
)*Copec-QDcool
(t
))));
end
if X
(i,
(t-1
)*D
+6
)>=0
SBat
(i,t
+1
)= SBat
(i,t
)*
(1-dBat
)+X
(i,
(t-1
)*D
+6
)*yBatcha; %蓄电池组在单时段时的剩余容量
else
SBat
(i,t
+1
)= SBat
(i,t
)*
(1-dBat
)+X
(i,
(t-1
)*D
+6
)/yBatdis; %蓄电池组在单时段时的剩余容量
end
if SBat
(i,t
+1
)>SBatmax
SBat
(i,t
+1
)=SBatmax;
end
if SBat
(i,t
+1
)<SBatmin
SBat
(i,t
+1
)=SBatmin;
end
end
if fitness
(X
(i,:
),WEI
)<p
(i
)
p
(i
)=fitness
(X
(i,:
),WEI
);
y
(i,:
)=X
(i,:
);
end
if p
(i
)<fitness
(pg,WEI
)
pg=y
(i,:
);
end
end
Gbest
(K
)=fitness
(pg,WEI
);
end
P_mt = pg
(1,1:24
); %燃气轮机
P_mt1 = pg
(1,
(24
+1
):24*2
);
P_mt2 = pg
(1,
(24*2
+1
):24*3
);%燃气锅炉
P_mt3 = pg
(1,
(24*3
+1
):24*4
);
P_mt4 = pg
(1,
(24*4
+1
):24*
5);
P_buy = pg
(1,
(24*
5+1
):24*6
); %上网交互
P_discharge = pg
(1,
(24*6
+1
):24*7
); %电储能
C_ec = pg
(1,
(24*7
+1
):24*8
); %电制冷机
扫一扫
981
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
