ga_f362.m

书籍代码:遗传演算法原理与应用_活用MATLAB(Source Code)

function PI=GA_f362(x)
global MIN_offset Se Sde dSu MDL Emin Final_time
% Fitness function evaluation of Example 3.6-2

% PenChen Chou, 7-21-2001

Se =GA_round4dp(x(1));
Sde=GA_round4dp(x(2));
dSu =GA_round4dp(x(3));
[T, X, Y]=sim(MDL,[0 Final_time]);
in=Y(:,1); y=Y(:,2); STP=sum(Y(:,3)); 
err=Y(:,4); N=length(err); ERR=err(N);
if STP>0, ERR=Emin+1000; end;
PI=ERR;
if MIN_offset>0
   PI=MIN_offset-PI;
   if PI<0, PI=10; end;
end;
if Emin>ERR
   Emin=ERR; 
   save xxx Se Sde dSu ERR
   fprintf('\nSe=%f;Sde=%f;dSu=%f;ERR=%f;\n',Se,Sde,dSu,ERR);
else
   fprintf('>'); 
end;