基于遗传算法的微电网负荷优化，基于遗传算法的微电网负荷均衡优化（代码完整，数据齐全）

function [F,x] = fun(E,umax,umin,Size,BPoutput) CodeL=10; nx = size(umax,2) ;%变量个数 x = zeros(Size,nx);%存储变量 n0 = 0.98;%系数1 n1 = 0.95;%系数2 F = zeros(Size,1);%适应度值 for s=1:1:Size %解码，二进制转换为十进制 m=E(s,:); for jj = 1:nx y1 = 0; m1=m((jj-1)*CodeL+1:1:jj*CodeL); for ii=1:1:CodeL y1=y1+m1(ii)*2^(ii-1); end x1=(umax(1)-umin(1))*y1/1023+umin(1); x(s,jj) = x1; end end pdt = zeros(Size,nx);%蓄热罐状态 pct = zeros(Size,nx);%制热设备 % pt = zeros(Size,nx); for s = 1:Size%求解ST for jj = 1:nx if jj >=8&&jj<20&&BPoutput(jj)>x(s,jj-1)%约束条件1 x(s,jj) = 0; elseif jj >=8&&jj<20&&BPoutput(jj)<=x(s,jj-1)%约束条件2 x(s,jj) = x(s,jj-1)-BPoutput(jj); elseif jj ==1&&x(s,24)<max(x(s,1:jj))-5%约束条件3 1 x(s,jj) = x(1,24)+5*n1; elseif jj<8&&jj>1&&x(s,24)<max(x(s,:))-5%约束条件3 2-7 x(s,jj) = x(s,24)+5*n1; elseif jj<8&&jj>1&&x(s,24)>max(x(s,1:jj))-5%约束条件4 2-7 x(s,jj) =max(x(s,1:jj)); elseif jj>20 &&x(s,24)>max(x(s,1:jj))-5%约束条件4 21-24 x(s,jj) = max(x(s,1:jj)); elseif jj>20&&x(s,24)<max(x(s,1:jj))-5 %约束条件3 21-24 x(s,jj) = x(s,jj-1)+5*n1; elseif jj ==1&&x(s,24)>=max(x(s,1:jj))-5 %约束条件4 1 x(s,jj) = max(x(s,:)); end end end pt = x; for s=1:1:Size for jj = 1:nx %%pdt if jj >=8&&jj<20&&BPoutput(jj)<=pt(s,jj-1) %约束条件1 pdt(s,jj) = 0; elseif jj >=8&&jj<20&&BPoutput(jj)>=pt(s,jj-1)&&pt(s,jj-1)>=0 %约束条件2 pdt(s,jj) = (BPoutput(jj)-pt(s,jj-1))/n0; else %约束条件3 pdt(s,jj) = BPoutput(jj)/n0; end %%pct if jj>1%大于1的情况 if jj<8&&pt(s,jj-1)<max(pt(s,:))%约束条件1 2-7 pct(s,jj) = min(15,(max(pt(s,:))-pt(s,jj-1))/n1); elseif jj<=24&&jj>=20&&pt(s,jj-1)<max(pt(s,:))%约束条件1 21-24 pct(s,jj) = min(15,(max(pt(s,:))-pt(s,jj-1))/n1); else %约束条件2 pct(s,jj) = 0; end elseif jj ==1; if pt(s,24)<max(pt(s,:))%约束条件1 pct(s,jj) = (max(pt(s,:))-pt(s,24))/n1; else %约束条件2 pct(s,jj) = 0; end end end end %%计算目标函数 for s=1:1:Size for jj = 1:nx if pdt(s,jj)*n0+pct(s,jj)*n1<BPoutput(jj)%不满足热负荷平衡约束 F(s) = F(s)+2000; break elseif pdt(s,jj)>25%不满足直接制热器电功率上下限约束 F(s) = F(s)+2000; break elseif pct(s,jj)>20%不满足蓄热装置容量约束 F(s) = F(s)+2000; break else%满足所有的约束条件 F(s) = F(s)+pdt(s,jj)*0.4883+(pdt(s,jj)+pct(s,jj))*0.1; end end end F = 1./F; end