eldga1.m

基于遗传算法的经济调度处理,word文档有简要说明,并且有算例

   function[ F P1 Pl]=eldga(x)
   global data B B0 B00 Pd


 x=abs(x);
n=length(data(:,1));
for i=1:n-1
    if x(i)>1;
        x(i)=1;
    else
    end
    P(i)=data(i+1,4)+x(i)*(data(i+1,5)-data(i+1,4));
end
% Demand (MW)
Pd=150;
B11=B(1,1);
B1n=B(1,2:n);
Bnn=B(2:n,2:n);
A=B11;
BB1=2*B1n*P';
B1=BB1-1+B0(1);
C1=P*Bnn*P';
C=Pd-sum(P)+C1+B00+B0(2:n)*P';
x1=roots([A B1 C]);
% x=.5*(-B1-sqrt(B1^2-4*A*C))/A
 x=abs(min(x1));
 if x>data(1,5)
     x=data(1,5);
 else
 end
   if x<data(1,4)
x=data(1,4);
   else
 end

P1=[x P];
for i=1:n
   F1(i)=data(i,1)* P1(i)^2+data(i,2)*P1(i)+data(i,3);
end
Pl=P1*B*P1'+B0*P1'+B00;
 lam=abs(sum(P1)-Pd-P1*B*P1'-B0*P1'-B00);
 F=sum(F1)+10000*lam;