火电成本包括运行成本、碳捕集成本,其中火电运行成本由运行维护成本和发电煤耗成本构成,发电煤耗与其出力关系: 2 F = aP +bP + c,F为机组煤耗量(kg/h),P 为机组出力/MW;运行维护成本按照0.5 倍煤耗成本考虑,碳捕集成本取决于碳排放量及碳捕集单价,火电机组相关参数如附表1 所示,电煤价格为700 元/t。风电成本仅考虑运维成本,相关参数如附表2。
储能成本由投资成本、运维成本构成,相关参数如附表3。注:在计算每天成本时,需将投资成本平摊至每天,即平均每天投资成本=总投资成本/运行年限/365 天。
弃风损失按0.3 元/kWh 计算,失负荷损失按8 元/kWh 计算。
理清了问题结构还是好做的
已经写出论文和代码,文件见文末:
function [xm,fv]=AsyLnCPSO(fitness,N,c1max,c1min,c2max,c2min,w,M,D)
%待优化的目标函数:fitness
%粒子数目:N
%学习因子1的最大值:c1max
%学习因子1的最小值:c1min
%学习因子2的最大值:c2max
%学习因子2的最小值:c2min
%惯性权重:w
format long;
for i=1:N
for j=1:D
x(i,j)=randn; %随机初始化位置
v(i,j)=randn; %随机初始化速度
end
end
for i=1:N
p(i)=fitness(x(i,:));
y(i,:)=x(i,:);
end
pg=x(N,:); %pg为全局最优
for i=1:(N-1)
if fitness(x(i,:))
function [xv,fv]=GMGA(fitness,a,b,NP,NG,Pc,Pm,alpha,Pbm,eps)
%待优化目标函数:fitness
%自变量下界:a
%自变量下界:b
%种群个体数:NP
%最大进化代数:NG
%杂交概率:Pc
%变异概率:Pm
%密集因子:alpha
%大变异概率:Pbm
%自变量离散精度:eps
%目标函数取最小值时的自变量值:xm
%目标函数的最小值:fv
L=ceil(log2((b-a)/eps+1)); %根据离散精度,确定二进制编码需要的码长
x=zeros(NP,L);
for i=1:NP
x(i,:)=Initial(L); %种群初始化
fx(i)=fitness(Dec(a,b,x(i,:),L)); %个体适应值
end
for k=1:NG
sumfx=sum(fx); %所有个体适应值之和
Px=fx/sumfx; %所有个体适应值的平均值
[fmax,xmax]=max(fx); %群体适应值最大值
if kfv
fv=fitx; %取个体中最好值作为最终结果
xv=Dec(a,b,x(i,:),L);
end
end
function result=Initial(length) %初始化函数
for i=1:length
r=rand();
result(i)=round(r);
end
全部代码和论文https://mianbaoduo.com/o/bread/YpqcmJhs



