ga_f614c.m

其中提到遺傳學的程式碼與應用提供給次淚相向的研究者參考下載

function PI=GA_f614c(x)
global MIN_offset MUL_factor Emin MDL r_deg Init_P Init_V
global V11 V12 V21 V22 W11 W12 b11 b21 c11 FT Mode K a1 a2 Se Sde Su
% Fitness function evaluation of Example 6.1-4

% PenChen Chou, 7-29-2001
Init_P=10;
Init_V=0;
K=10;
a1=2;
a2=3;

V11=GA_round4dp(x(1));
V12=GA_round4dp(x(2));
V21=GA_round4dp(x(3));
V22=GA_round4dp(x(4));
b11=GA_round4dp(x(5));
b21=GA_round4dp(x(6));
W11=GA_round4dp(x(7));
W12=GA_round4dp(x(8));
c11=GA_round4dp(x(9));
Se =GA_round4dp(x(10));
Sde=GA_round4dp(x(11));
Su =GA_round4dp(x(12));

% Get Y
[t, X, Y]=sim(MDL,[0:0.05:FT]);
error=abs(Y(:,1)); 
I1=find(t>5 & t<8);   err1=error(I1);
I2=find(t>9 & t<12); err2=error(I2);
ERR=100*(sum(err1)+sum(err2));
PI=ERR;
if MIN_offset>0
   PI=MIN_offset-PI;
   if PI<0, PI=10; end;
end;
if Emin>=ERR
   Emin=ERR; fprintf('\nERR=%f',ERR);
   %save ex614_data V11 V12 V21 V22 W11 W12 b11 b21 c11 ERR
   %pause(2)
else
   fprintf('>'); 
end;