ga_rule.m

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

% GA_rule.m file
% Find scaling factors Se, Sde and Su and rules for a FLC.

% PenChen Chou, 5-9-2002

clear all
ans=1;
who
fprintf('>>>Start this program [GA_rule.m] now. Please wait!\n');
pause(1)
global MIN_offset MUL_factor
%**********************************************************/
%       User can modify the following in blocks
%**********************************************************/
global MDL Se Sde Su FT IN_main Emin f0408A f0408 f0408C ruleList
MIN_offset=125e9; MUL_factor=1; IN_main=0; FT=0.4;
MDL='MD04081B'; Emin=MIN_offset; eval(MDL)
f0408A=readfis('f0408A');
obj_fcn = 'GA_frule';	% Objective function
range = [ones(1,9)   0.01 0.001 20 0 0 0   
         5*ones(1,9)  10   10   220 50 50 50];	% Range of the input variables
IC=[]; elite=1; 
gen_no=80; popuSize=45;
bit_n=100; xover_rate=0.8; mutate_rate=0.1;
%************************************************************
% 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 the final plot and data
%Se=0.5;Sde=0.1;Su=100;para=[Se Sde Su];
%Se=0.312;Sde=0.0343;Su=38.24;para=[Se Sde Su];
IN_main=1; load garule % Get the final x and ruleList
GA_frule(x);