Genetic algorithm only works with integer constraints..
7 views (last 30 days)
I am solving a problem using GA, and while my problem formulation does not require to use integer constraints, I found out that using integer constraints is the only way to solve it - and I do not understand why.. I will explain below.
I have 20 optimization variables.
- 9 of them represent switching times, formulated in a non-dimensional form, and bounded by [0 1],
- one is the final time, bounded by [tf1 tf2]
- other 10 variables represent angles, bounded by [0, 2pi]
If I do not introduce any integer constraints, and only use bounds for the optimization variables, then the GA window looks as below. Nothing changes on the plot, and there is no solution if I stop GA.
However, if I introduce integer constraints for the angles, then GA works very well for my problem:
To do this, I bound angles from 0 to 360 (achieving a 1 degree step), or from 0 to 720 (0.5deg step), and so on..
I learnt that no matter how I formulate this problem (different number of variables, times instead of switching times, etc.), this is always the pattern - a portion of the optimization variables have to be integer constrained. If not for angles, then using integer constraints for the times is also working (such that the time variable can take any value with a step of 1 second).
I would like to understand why this is the case. Ideally, I would like not to use integer constraints. Could it be that I need a better computer?
Alan Weiss on 6 Oct 2022
Yakov, your report indicates that the nonlinearly-constrained problem is being solved in the usual way: having very few iterations, because most of the time is taken up with solving subproblem iterations.
For a similar plot in the documentation, see this link.
When you include integer constraints, the solver uses a different algorithm that shows many more iterations. For a simple demonstration of this effect, see the following simple example:
fun = @(x)log(1 + 100*(x(2) - x(1)^2)^2 + (1 - x(1))^2);
opts = optimoptions('ga',PlotFcn='gaplotbestf');
nvar = 2;
lb = [-2 -2];
ub = -lb;
sol1 = ga(fun,nvar,,,,,lb,ub,@nlcon,,opts);
sol2 = ga(fun,nvar,,,,,lb,ub,@nlcon,1,opts);
function [c,ceq] = nlcon(x)
c = x(1)^2 + x(2)^2 - 1;
ceq = ;
MATLAB mathematical toolbox documentation