Genetic Algorithm options 'PopInitRange' not working

2013 6 11
2 个回答
回答已采纳
5 次查看（30 天）

0 个投票

I am using GAOPTIMSET to establish an upper and lower bound for my initial population in a GA optimisation.
I have abbreviated my inputs as below:
InitPopRange = [0.0076 0.0030 0.0011 0.0004 0.0068 0.0004; 0.0303 0.0068 0.0149 0.0049 0.0118 0.0303];
Options = gaoptimset('PopInitRange' , InitPopRange);
The optimisation kept crashing in my fitness function just after the very first individual had been created. During debugging, I noticed that my first individual had the following values:
x = [0.2508 0.0000 0.0000 0.2132 0.0970 1.0000];
Why is it that every single value of x falls outside InitPopRange? This is what has caused my fitness function to crash.
Obviously the work-around solution for me is to create my own initial population using a random number generator and apply my own limits. I'm just curious to know why the GA has seemingly ignored my option for initial population range. Any ideas?

5 个评论
显示 3更早的评论 隐藏 3更早的评论

Walter Roberson
Walter Roberson 2013-6-11
Which MATLAB version are you using?
Matt J
Matt J 2013-6-11
编辑：Matt J 2013-6-11
Did you actually pass "Options" to GA? We need to see more code, including the actual function call to GA.
Craig Cowled
Craig Cowled 2013-6-11
编辑：Craig Cowled 2013-6-11
Version = Matlab R2013a.
MattJ, good question. I used debug mode to check that the Options were correctly passed into the GA and it was.
My code is really very long, so I'll just provide the relevant sections ...
function [x, Frequ, CombMAC, PrepParameters, objFunVal] = GAtoFEMupdate(para)
% A lot of code
InitPopRange = [0.0076 0.0030 0.0011 0.0004 0.0068 0.0004; ...
0.0303 0.0068 0.0149 0.0049 0.0118 0.0303];
GA_Options = gaoptimset('PopInitRange', InitPopRange); %omit other opt
[x, Frequ, CombMAC, PrepParameters, objFunVal] = ...
NestedGA(A, b, LB, UB, IntCon, GA_Options);
end % GAtoFEMupdate()
function [x, Frequ, CombMAC, PrepParameters, objFunVal] = NestedGA(A, b, LB, UB, IntCon, GA_Options)
% some code
[x, objFunVal] = ga(@St7Fitness, nvars, A, b, [], [], LB, UB, [], ...
IntCon, GA_Options);
function objFunVal = St7Fitness(x)
% All my fitness func code
end % St7Fitness()
end % NestedGA
Matt J
Matt J 2013-6-11
I noticed that my first individual had the following values:
What method did you use to observe the individuals? How do you know this individual was the "first"?
Even if you do succeed in getting the Initial Population that you want, what is to prevent the population from evolving at later iterations into a region which causes your fitness function to crash?
Because this is a Genetic Algorithm, the fitness function doesn't have to be smooth/continuous, so I think one normally defines the fitness function to return a safe value everywhere.
Craig Cowled
Craig Cowled 2013-6-11
MattJ, I was in debug mode when I observed the value of the first individual in my population. I placed a break within my fitness function so that I could catch the first value of 'x'.
I take your point about constraining the fitness function so that it only evaluates safe values. I have allowed values of x to vary between 0 and 1. In my case, some of the values can be zero and the fitness function will still evaluate, however, the fitness function becomes unstable when too many values of x are very close to zero (I don't know what combinations though). I don't want to move my lower bounds up because it is quite possible that my solution has some values close to zero.
In any case, I have a work-around. I guess I was just curious to know why PopInitRange didn't seem to work.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

Alan Weiss
Alan Weiss 2013-6-11

2 个投票

I believe that the problem is that you have linear constraints. By default, ga uses the gacreationlinearfeasible function to create individuals in that case, and gacreationlinearfeasible does not pay attention to PopInitRange when there are linear constraints. To see how it works, look at this example.
You can choose to use the gacreationuniform creation function, which does exactly what you think it should when you pass in an initial range (it creates individuals uniformly between the bounds given in the initial range). To use this function,
options = gaoptimset('CreationFcn',@gacreationuniform,...);
Alan Weiss
MATLAB mathematical toolbox documentation

更多回答（1 个）

Craig Cowled
Craig Cowled 2013-6-12

0 个投票

Alan, thank you for such a good explanation. Cleared it up for me nicely.

类别

在 帮助中心 和 File Exchange 中查找有关 Genetic Algorithm 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by