目录
主要内容
该模型构建了含冷、热、电负荷的冷热电联供型综合能源系统优化调度模型，设备包括燃气轮机、电制冷机、锅炉、风光出力以及购售电交易，以用户购电购热冷量的成本、CCHP收益以及碳排放成本为目标，采用权值多目标方式将多目标转化为单目标求解，同时将约束采用罚函数形式进行表达，实现系统的经济运行，程序采用matlab编写，注释清晰，方便学习。
部分代码
PGnom1=1200;PGnom2=2800;%电功率基准值为1200和2800
Copac=0.80;Copec=3;
%-------------%
for t=1:24
ygg(1,t)=(Aa+Bb*(X((t-1)*D+1)/PGnom1)+Cc*(X((t-1)*D+1)/PGnom1)^2+Dd*(X((t-1)*D+1)/PGnom1)^3)/100; %X((t-1)*D+1)在对应的该时刻的该燃气发电机的电功率
ygg(2,t)=(Aa+Bb*(X((t-1)*D+2)/PGnom2)+Cc*(X((t-1)*D+2)/PGnom2)^2+Dd*(X((t-1)*D+2)/PGnom2)^3)/100;
%ygg对应的是第i台燃气发电机的发电效率
end
%-------------%
for t=1:24
%R是燃气机发电模型的λ
FG(t)=R*(X((t-1)*D+1)/ygg(1,t)+X((t-1)*D+2)/ygg(2,t)); %燃气发电机在单时段的燃气耗量
%R疑似是锅炉模型的λ
FBoi(t)=R*(X((t-1)*D+3)+X((t-1)*D+4))/yBoi; %锅炉在单时段消耗的燃料量
FE(t)=FG(t)+FBoi(t); %t时刻的燃气量
if X((t-1)*D+5)>=0
Rr(t)=rElc(t);
PG(t)=X((t-1)*D+5); %向主电网交互的电量，此时是买电
else
Rr(t)=RElc(t);
PG(t)=0; %向主电网交互的电量，此时是卖电
end
CE(t)=Rr(t)*X((t-1)*D+5); %单时段与主电网间能量交互成本
HJ(t)=ue*PG(t)+uf*FE(t); %二氧化碳的排放量
%---单时段功率平衡约束---%
h(t)=abs(X((t-1)*D+1)+X((t-1)*D+2)+Pwind(t)+Ppv(t)-X((t-1)*D+6)+X((t-1)*D+5)-PD(t)-X((t-1)*D+8));%电负荷平衡约束
end
%---用的燃气量不同，得到的价格也不同---%
if sum(FE)<250
rLrg=5.257;
elseif 250<=sum(FE)<=4167
rLrg=3.25;
else
rLrg=2.814;
end
%-----------------------------------%
f1=rLrg*sum(FE)+sum(CE); %燃气成本和向电网购电或是售电后的花费
f2=W*sum(HJ); %碳排放消费费用结果一览
