✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。
🍎个人主页：Matlab科研工作室
🍊个人信条：格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
为了解决现有冷热电联供型综合能源系统大多只单一考虑系统机组投资成本或系统环境污染，影响系统整体优化运行的问题，以系统经济性和环保性为目标，对冷热电联供系统进行研究分析。构建含燃气轮机、燃气锅炉、电制冷机等机组 的冷热电联供系统优化模型并建立约束条件。结果表明粒子群算法能够同时兼顾系统的经济性和环保性，使系统运行更加优化，为之后的能源供给系统的规划提供前期依据。
⛄ 部分代码
function result=fitness(x,k)global PV;global WT;global P_load;global R_load;global L_load;global G_price;global C_grid_sell;%卖电电价gas_price=0.175; %气价nGT_e=0.29; %燃气轮机转化效率nGT_h=0.61;nGB_h=0.85; %热锅炉的热效率nGE_e=0.35; %天燃气内燃机转化效率nGE_h=0.7;COP_EC=5; %电制冷机的制冷效率COP_AC=1.2;%吸收制冷机的制冷效率n_hs=0.9; %余热回收系统的效率n_hr=0.85; %换热器的效率for j=1:120 if j<25 PGT(j)=x(j); elseif j>24&&j<49 PGB(j-24)=x(j); elseif j>48&&j<73 PGE(j-48)=x(j); elseif j>72&&j<97 PEC(j-72)=x(j); elseif j>96 Grid(j-96)=x(j); endendPGTe=PGT*nGT_e; %燃气轮机转化PGTh=PGT*nGT_h;PGBh=PGB*nGB_h; %热锅炉转化PGEe=PGE*nGE_e; %天燃气内燃机转化PGEh=PGE*nGE_h;PE= PV'+WT'+ Grid+PGTe+PGEe-PEC; %电网总出力P_hs=(PGEh+PGTh)*n_hs+PGBh; %余热回收系统的热量P_hr=P_hs*n_hr; %过换热器提高的热量 总热量Q_AC=P_hs*COP_AC; %吸收制冷机的冷量Q_EC=PEC*COP_EC; %电制冷机的冷量Q_all=Q_AC+Q_EC; %总冷量sum_delt_P=0;sum_delt_R=0;sum_delt_L=0;delt_P=(P_load-PE);delt_R=(R_load-P_hr);delt_L=(L_load-Q_all);sum_delt_P=sum(delt_P);sum_delt_R=sum(delt_R);sum_delt_L=sum(delt_L);BT1=delt_P;BT2=delt_R;BT3=delt_L;SOC=0; dsoc=0;for i=1:24 SOC=SOC+BT1(i); if SOC>1000 dsoc=dsoc+(SOC-1000); end if SOC<-1000 dsoc=dsoc+abs((SOC+1000)); end endr1=0;r2=0;r3=0;if(sum_delt_P<2000) r1=1.0;elseif(sum_delt_P>2000&&sum_delt_P<=3000) r1=5;elseif(sum_delt_P>3000&&sum_delt_P<=4000) r1=10;elseif(sum_delt_P>4000&&sum_delt_P<=6000) r1=100;else r1=1000;endif(sum_delt_R<1000) r2=1.0;elseif(sum_delt_R>1000&&sum_delt_R<=2000) r2=5;elseif(sum_delt_R>2000&&sum_delt_R<=3000) r2=10;elseif(sum_delt_R>3000&&sum_delt_R<=4000) r2=20;else r2=1000;endif(sum_delt_L<3000) r3=1.0;elseif(sum_delt_L>3000&&sum_delt_L<=4000) r3=5;elseif(sum_delt_L>4000&&sum_delt_L<=8000) r3=10;elseif(sum_delt_L>8000&&sum_delt_L<=10000) r3=50;else r3=2000;endC_all=0;for i=1:24 C_PEC(i)= G_price(i)*PEC(i); %电制冷机买电成本 if Grid(i)>0 C_grid(i)= G_price(i)*Grid(i); %电网购电成本 else C_grid(i)= C_grid_sell*Grid(i);%电网售电收益 endendC_all=gas_price*(sum( PGT)+sum( PGB)+sum( PGE))+sum( C_grid)++sum( C_PEC)+0.8*sum(abs(BT1)+abs(BT2)+abs(BT3));% result=C_all+r1*sum_delt_P+r2*sum_delt_R+r3*sum_delt_L;result=C_all+10*dsoc;
⛄ 运行结果
