基于多目标粒子群算法冷热电联供综合能源系统运行优化（matlab代码）

主要内容   

该模型构建了含冷、热、电负荷的冷热电联供型综合能源系统优化调度模型，设备包括燃气轮机、电制冷机、锅炉、风光出力以及购售电交易，以用户购电购热冷量的成本、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);  %碳排放消费费用

  结果一览   

下载链接