ga_gen_p43.asv

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

% GA_gen_p43.m file
% Fin min of a function of 2 varaibles
% See P.43 of Gen & Cheng's Genetic Algorithms & Engineering Design
% Wiley Interscience. 1997.

% PenChen Chou, 2003-3-12

clear all
ans=1;
who
fprintf('>>>Start this program [GA_gen_p43.m] now. Please wait!\n');
pause(1)
global MIN_offset MUL_factor LOCUS x_data y_data
MIN_offset=250; LOCUS=1; x_data=[]; y_data=[];
MUL_factor=1;
%**********************************************************/
%       User can modify the following in blocks
%**********************************************************/
obj_fcn = 'GA_gen_fp43';	% Objective function
range = [-5 5
    -5 5];	% Range of the input variables
IC=[]; elite=1; 
gen_no=500; popuSize=50; bit_n=40;
xover_rate=0.8; mutate_rate=0.05;

x=-5:0.1:5; y=x; [X,Y]=meshgrid(x,y);
c1=20;c2=0.2;c3=2*pi;efactor=2.71282;
Z=;
figure(2);contour(X,Y,Z)
%************************************************************
% call GA_genetic
tic
[popu, popu_real, fcn_value, upper, average,...
     lower, BEST_popu, popuSize, gen_no, para, best_pi,...
     bit_n, xover_rate, mutate_rate]=...
GA_genetic(obj_fcn, range, IC, elite, gen_no, popuSize,...
     bit_n, xover_rate, mutate_rate);
t=toc/60;
fprintf('==>  Computation time is (%.2f) minutes.\n\n',t);            
            
% Plot locus of optimization process
LAST=length(x_data);
figure(2); hold on
plot(x_data(1),y_data(1),'x');
plot(x_data,y_data);
plot(x_data(LAST),y_data(LAST),'o');
%plot(0,0,'*');   % The best position
hold off
title('*=BEST, x=Starting point, o=Ending point')