24小时热门版块排行榜

返回列表

zongshengl

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 5218.9
散金: 144
红花: 1
帖子: 151
在线: 103.2小时
虫号: 1285718
注册: 2011-05-04
性别: GG
专业: 运筹与管理

[求助] 遗传算法整数规划中变量取整的设置

本人在用遗传算法求解非线性整数规划问题，关于变量整数的约束前面看到一个有关的帖子，我发现中间有错误（红色部分），但是不知道怎么修改，请指点。
源程序如下：
function [x,fval] = gainteger_demo
% Fitness function and numver of variables
fitnessFcn = @(x) norm(x);
numberOfVariables = 15;
% If decision variables are bounded provide a bound e.g, LB and UB.
LB = -5*ones(1,numberOfVariables);
UB = 5*ones(1,numberOfVariables);
Bound = [LB;UB]; % If unbounded then Bound = []
% Create an options structure to be passed to GA
% Three options namely 'CreationFcn', 'MutationFcn', and
% 'PopInitRange' are required part of the problem.
options = gaoptimset('CreationFcn',@int_pop,'MutationFcn',@int_mutation, ...
      'PopInitRange',Bound,'Display','iter','StallGenL',40,'Generations',150, ...
      'PopulationSize',60,'PlotFcns',{@gaplotbestf,@gaplotbestindiv});
[x,fval] = ga(fitnessFcn,numberOfVariables,options);
%---------------------------------------------------
% Mutation function to generate childrens satisfying the range and integer
% constraints on decision variables.
function mutationChildren = int_mutation(parents,options,GenomeLength, ...
      FitnessFcn,state,thisScore,thisPopulation)
shrink = .01;
scale = 1;
scale = scale - shrink * scale * state.Generation/options.Generations;
range = options.PopInitRange;
lower = range(1,

;
upper = range(2,

;
scale = scale * (upper - lower);
mutationPop =    length(parents);
% The use of ROUND function will make sure that childrens are integers.
mutationChildren =    repmat(lower,mutationPop,1) +    ...
      round(repmat(scale,mutationPop,1) .* rand(mutationPop,GenomeLength));
% End of mutation function
%---------------------------------------------------
function Population = int_pop(GenomeLength,FitnessFcn,options)
totalpopulation = sum(options.PopulationSize);
range = options.PopInitRange;
lower= range(1,

;
span = range(2,

- lower;
% The use of ROUND function will make sure that individuals are integers.
Population = repmat(lower,totalpopulation,1) + ...
round(repmat(span,totalpopulation,1) .* rand(totalpopulation,GenomeLength));
% End of creation function

回复此楼