xfcn6a.m

来自「其中提到遺傳學的程式碼與應用提供給次淚相向的研究者參考下載」· M 代码 · 共 26 行

function PI = xfcn6a(x)

global KAV KNA KVF Mass KBEMF Emin COL DZO
KAV=x(1); KNA=x(2); KVF=x(3); Mass=x(4); KBEMF=x(5);
COL=x(6); DZO=x(7);
KBEMF=0; %Mass=1.5;
MDL='hong2000';
[T X Y]=sim(MDL,[0 1.5]);
E=Y(:,1); STP=sum(Y(:,2));
if STP>=1
   Err=10000;
   fprintf('>');
else   
   Err=E(length(E));
   if Emin>Err
      fprintf('S');
      Emin=Err;
      save xxx
      fprintf('\nKAV=%.2f; KNA=%.2f; KVF=%.6f;\n',KAV,KNA,KVF);
      fprintf('Mass=%.6f; KBE=%.2f; \n',Mass,KBEMF);
      fprintf('COL=%.6f; DZO=%.2f; Err=%10.6f;\n',COL,DZO,Err);
   end;
end;   
PI=10000-Err;