Genetic Algorithms - Quick Guide
Genetic Algorithms - Quick Guide
Login
___
Home
Coding
HTML Ground CSS Javascript SQL Python Java C C++
Jobs
Whiteboard
Genetic Algorithms - Quick Guide
Tools
Introduction to Optimization
Optimization is the process of making something better. In any process, we have a set of
inputs and a set of outputs as shown in the following figure.
Optimization refers to finding the values of inputs in such a way that we get the “best”
output values. The definition of “best” varies from problem to problem, but in mathematical
terms, it refers to maximizing or minimizing one or more objective functions, by varying the
input parameters.
The set of all possible solutions or values which the inputs can take make up the search
space. In this search space, lies a point or a set of points which gives the optimal solution.
The aim of optimization is to find that point or set of points in the search space.
ADVERTISEMENT
ADVERTISEMENT
GAs were developed by John Holland and his students and colleagues at the University of
Michigan, most notably David E. Goldberg and has since been tried on various optimization
problems with a high degree of success.
In GAs, we have a pool or a population of possible solutions to the given problem. These
solutions then undergo recombination and mutation (like in natural genetics), producing
new children, and the process is repeated over various generations. Each individual (or
candidate solution) is assigned a fitness value (based on its objective function value) and
the fitter individuals are given a higher chance to mate and yield more “fitter” individuals.
This is in line with the Darwinian Theory of “Survival of the Fittest”.
In this way we keep “evolving” better individuals or solutions over generations, till we reach
a stopping criterion.
Genetic Algorithms are sufficiently randomized in nature, but they perform much better
than random local search (in which we just try various random solutions, keeping track of
the best so far), as they exploit historical information as well.
Advantages of GAs
GAs have various advantages which have made them immensely popular. These include −
Does not require any derivative information (which may not be available for many
real-world problems).
Is faster and more efficient as compared to the traditional methods.
Useful when the search space is very large and there are a large number of
parameters involved.
Limitations of GAs
Like any technique, GAs also suffer from a few limitations. These include −
GAs are not suited for all problems, especially problems which are simple and for
which derivative information is available.
Fitness value is calculated repeatedly which might be computationally expensive
for some problems.
Being stochastic, there are no guarantees on the optimality or the quality of the
solution.
If not implemented properly, the GA may not converge to the optimal solution.
ADVERTISEMENT
ADVERTISEMENT
GA – Motivation
Genetic Algorithms have the ability to deliver a “good-enough” solution “fast-enough”. This
makes genetic algorithms attractive for use in solving optimization problems. The reasons
why GAs are needed are as follows −
In computer science, there is a large set of problems, which are NP-Hard. What this
essentially means is that, even the most powerful computing systems take a very long time
(even years!) to solve that problem. In such a scenario, GAs prove to be an efficient tool to
provide usable near-optimal solutions in a short amount of time.
Traditional calculus based methods work by starting at a random point and by moving in
the direction of the gradient, till we reach the top of the hill. This technique is efficient and
works very well for single-peaked objective functions like the cost function in linear
regression. But, in most real-world situations, we have a very complex problem called as
landscapes, which are made of many peaks and many valleys, which causes such methods
to fail, as they suffer from an inherent tendency of getting stuck at the local optima as
shown in the following figure.
Getting a Good Solution Fast
Some difficult problems like the Travelling Salesperson Problem (TSP), have real-world
applications like path finding and VLSI Design. Now imagine that you are using your GPS
Navigation system, and it takes a few minutes (or even a few hours) to compute the
“optimal” path from the source to destination. Delay in such real world applications is not
acceptable and therefore a “good-enough” solution, which is delivered “fast” is what is
required.
ADVERTISEMENT
ADVERTISEMENT
Basic Terminology
Before beginning a discussion on Genetic Algorithms, it is essential to be familiar with
some basic terminology which will be used throughout this tutorial.
Fitness Function − A fitness function simply defined is a function which takes the
solution as input and produces the suitability of the solution as the output. In some
cases, the fitness function and the objective function may be the same, while in
others it might be different based on the problem.
Genetic Operators − These alter the genetic composition of the offspring. These
include crossover, mutation, selection, etc.
ADVERTISEMENT
ADVERTISEMENT
Basic Structure
The basic structure of a GA is as follows −
We start with an initial population (which may be generated at random or seeded by other
heuristics), select parents from this population for mating. Apply crossover and mutation
operators on the parents to generate new off-springs. And finally these off-springs replace
the existing individuals in the population and the process repeats. In this way genetic
algorithms actually try to mimic the human evolution to some extent.
Each of the following steps are covered as a separate chapter later in this tutorial.
A generalized pseudo-code for a GA is explained in the following program −
GA()
initialize population
find fitness of population
Genotype Representation
One of the most important decisions to make while implementing a genetic algorithm is
deciding the representation that we will use to represent our solutions. It has been
observed that improper representation can lead to poor performance of the GA.
In this section, we present some of the most commonly used representations for genetic
algorithms. However, representation is highly problem specific and the reader might find
that another representation or a mix of the representations mentioned here might suit
his/her problem better.
ADVERTISEMENT
ADVERTISEMENT
Binary Representation
This is one of the simplest and most widely used representation in GAs. In this type of
representation the genotype consists of bit strings.
For some problems when the solution space consists of Boolean decision variables – yes
or no, the binary representation is natural. Take for example the 0/1 Knapsack Problem. If
there are n items, we can represent a solution by a binary string of n elements, where the
xth element tells whether the item x is picked (1) or not (0).
For other problems, specifically those dealing with numbers, we can represent the numbers
with their binary representation. The problem with this kind of encoding is that different
bits have different significance and therefore mutation and crossover operators can have
undesired consequences. This can be resolved to some extent by using Gray Coding, as a
change in one bit does not have a massive effect on the solution.
ADVERTISEMENT
ADVERTISEMENT
Integer Representation
For discrete valued genes, we cannot always limit the solution space to binary ‘yes’ or ‘no’.
For example, if we want to encode the four distances – North, South, East and West, we
can encode them as {0,1,2,3}. In such cases, integer representation is desirable.
Permutation Representation
In many problems, the solution is represented by an order of elements. In such cases
permutation representation is the most suited.
A classic example of this representation is the travelling salesman problem (TSP). In this
the salesman has to take a tour of all the cities, visiting each city exactly once and come
back to the starting city. The total distance of the tour has to be minimized. The solution to
this TSP is naturally an ordering or permutation of all the cities and therefore using a
permutation representation makes sense for this problem.
The population is usually defined as a two dimensional array of – size population, size x,
chromosome size.
ADVERTISEMENT
ADVERTISEMENT
Population Initialization
There are two primary methods to initialize a population in a GA. They are −
It has been observed that the entire population should not be initialized using a heuristic,
as it can result in the population having similar solutions and very little diversity. It has
been experimentally observed that the random solutions are the ones to drive the
population to optimality. Therefore, with heuristic initialization, we just seed the population
with a couple of good solutions, filling up the rest with random solutions rather than filling
the entire population with heuristic based solutions.
It has also been observed that heuristic initialization in some cases, only effects the initial
fitness of the population, but in the end, it is the diversity of the solutions which lead to
optimality.
ADVERTISEMENT
ADVERTISEMENT
Population Models
There are two population models widely in use −
Steady State
In steady state GA, we generate one or two off-springs in each iteration and they replace
one or two individuals from the population. A steady state GA is also known as
Incremental GA.
Generational
In a generational model, we generate ‘n’ off-springs, where n is the population size, and the
entire population is replaced by the new one at the end of the iteration.
In most cases the fitness function and the objective function are the same as the objective
is to either maximize or minimize the given objective function. However, for more complex
problems with multiple objectives and constraints, an Algorithm Designer might choose to
have a different fitness function.
In some cases, calculating the fitness function directly might not be possible due to the
inherent complexities of the problem at hand. In such cases, we do fitness approximation
to suit our needs.
The following image shows the fitness calculation for a solution of the 0/1 Knapsack. It is
a simple fitness function which just sums the profit values of the items being picked
(which have a 1), scanning the elements from left to right till the knapsack is full.
However, care should be taken to prevent one extremely fit solution from taking over the
entire population in a few generations, as this leads to the solutions being close to one
another in the solution space thereby leading to a loss of diversity. Maintaining good
diversity in the population is extremely crucial for the success of a GA. This taking up of
the entire population by one extremely fit solution is known as premature convergence and
is an undesirable condition in a GA.
ADVERTISEMENT
ADVERTISEMENT
Consider a circular wheel. The wheel is divided into n pies, where n is the number of
individuals in the population. Each individual gets a portion of the circle which is
proportional to its fitness value.
In a roulette wheel selection, the circular wheel is divided as described before. A fixed point
is chosen on the wheel circumference as shown and the wheel is rotated. The region of the
wheel which comes in front of the fixed point is chosen as the parent. For the second
parent, the same process is repeated.
It is clear that a fitter individual has a greater pie on the wheel and therefore a greater
chance of landing in front of the fixed point when the wheel is rotated. Therefore, the
probability of choosing an individual depends directly on its fitness.
Stochastic Universal Sampling is quite similar to Roulette wheel selection, however instead
of having just one fixed point, we have multiple fixed points as shown in the following
image. Therefore, all the parents are chosen in just one spin of the wheel. Also, such a
setup encourages the highly fit individuals to be chosen at least once.
It is to be noted that fitness proportionate selection methods don’t work for cases where
the fitness can take a negative value.
ADVERTISEMENT
ADVERTISEMENT
Tournament Selection
In K-Way tournament selection, we select K individuals from the population at random and
select the best out of these to become a parent. The same process is repeated for
selecting the next parent. Tournament Selection is also extremely popular in literature as it
can even work with negative fitness values.
Rank Selection
Rank Selection also works with negative fitness values and is mostly used when the
individuals in the population have very close fitness values (this happens usually at the end
of the run). This leads to each individual having an almost equal share of the pie (like in
case of fitness proportionate selection) as shown in the following image and hence each
individual no matter how fit relative to each other has an approximately same probability of
getting selected as a parent. This in turn leads to a loss in the selection pressure towards
fitter individuals, making the GA to make poor parent selections in such situations.
In this, we remove the concept of a fitness value while selecting a parent. However, every
individual in the population is ranked according to their fitness. The selection of the
parents depends on the rank of each individual and not the fitness. The higher ranked
individuals are preferred more than the lower ranked ones.
A 8.1 1
B 8.0 4
C 8.05 2
D 7.95 6
E 8.02 3
F 7.99 5
ADVERTISEMENT
ADVERTISEMENT
Random Selection
In this strategy we randomly select parents from the existing population. There is no
selection pressure towards fitter individuals and therefore this strategy is usually avoided.
Crossover Operators
In this section we will discuss some of the most popularly used crossover operators. It is
to be noted that these crossover operators are very generic and the GA Designer might
choose to implement a problem-specific crossover operator as well.
In this one-point crossover, a random crossover point is selected and the tails of its two
parents are swapped to get new off-springs.
Uniform Crossover
In a uniform crossover, we don’t divide the chromosome into segments, rather we treat
each gene separately. In this, we essentially flip a coin for each chromosome to decide
whether or not it’ll be included in the off-spring. We can also bias the coin to one parent, to
have more genetic material in the child from that parent.
This is commonly used for integer representations and works by taking the weighted
average of the two parents by using the following formulae −
Obviously, if α = 0.5, then both the children will be identical as shown in the following
image.
OX1 is used for permutation based crossovers with the intention of transmitting
information about relative ordering to the off-springs. It works as follows −
Create two random crossover points in the parent and copy the segment between
them from the first parent to the first offspring.
Now, starting from the second crossover point in the second parent, copy the
remaining unused numbers from the second parent to the first child, wrapping
around the list.
Repeat for the second child with the parent’s role reversed.
There exist a lot of other crossovers like Partially Mapped Crossover (PMX), Order based
crossover (OX2), Shuffle Crossover, Ring Crossover, etc.
ADVERTISEMENT
Introduction to Mutation
In simple terms, mutation may be defined as a small random tweak in the chromosome, to
get a new solution. It is used to maintain and introduce diversity in the genetic population
and is usually applied with a low probability – pm. If the probability is very high, the GA gets
reduced to a random search.
Mutation is the part of the GA which is related to the “exploration” of the search space. It
has been observed that mutation is essential to the convergence of the GA while crossover
is not.
Mutation Operators
In this section, we describe some of the most commonly used mutation operators. Like the
crossover operators, this is not an exhaustive list and the GA designer might find a
combination of these approaches or a problem-specific mutation operator more useful.
In this bit flip mutation, we select one or more random bits and flip them. This is used for
binary encoded GAs.
Random Resetting
Random Resetting is an extension of the bit flip for the integer representation. In this, a
random value from the set of permissible values is assigned to a randomly chosen gene.
Swap Mutation
In swap mutation, we select two positions on the chromosome at random, and interchange
the values. This is common in permutation based encodings.
Scramble Mutation
Scramble mutation is also popular with permutation representations. In this, from the
entire chromosome, a subset of genes is chosen and their values are scrambled or
shuffled randomly.
Inversion Mutation
In inversion mutation, we select a subset of genes like in scramble mutation, but instead of
shuffling the subset, we merely invert the entire string in the subset.
Genetic Algorithms - Survivor Selection
The Survivor Selection Policy determines which individuals are to be kicked out and which
are to be kept in the next generation. It is crucial as it should ensure that the fitter
individuals are not kicked out of the population, while at the same time diversity should be
maintained in the population.
Some GAs employ Elitism. In simple terms, it means the current fittest member of the
population is always propagated to the next generation. Therefore, under no circumstance
can the fittest member of the current population be replaced.
The easiest policy is to kick random members out of the population, but such an approach
frequently has convergence issues, therefore the following strategies are widely used.
ADVERTISEMENT
ADVERTISEMENT
For instance, in the following example, the age is the number of generations for which the
individual has been in the population. The oldest members of the population i.e. P4 and P7
are kicked out of the population and the ages of the rest of the members are incremented
by one.
Fitness Based Selection
In this fitness based selection, the children tend to replace the least fit individuals in the
population. The selection of the least fit individuals may be done using a variation of any of
the selection policies described before – tournament selection, fitness proportionate
selection, etc.
For example, in the following image, the children replace the least fit individuals P1 and
P10 of the population. It is to be noted that since P1 and P9 have the same fitness value,
the decision to remove which individual from the population is arbitrary.
When the objective function value has reached a certain pre-defined value.
For example, in a genetic algorithm we keep a counter which keeps track of the
generations for which there has been no improvement in the population. Initially, we set
this counter to zero. Each time we don’t generate off-springs which are better than the
individuals in the population, we increment the counter.
However, if the fitness any of the off-springs is better, then we reset the counter to zero.
The algorithm terminates when the counter reaches a predetermined value.
Like other parameters of a GA, the termination condition is also highly problem specific
and the GA designer should try out various options to see what suits his particular problem
the best.
However, other models of lifetime adaptation – Lamarckian Model and Baldwinian Model
also do exist. It is to be noted that whichever model is the best, is open for debate and the
results obtained by researchers show that the choice of lifetime adaptation is highly
problem specific.
Often, we hybridize a GA with local search – like in Memetic Algorithms. In such cases, one
might choose do go with either Lamarckian or Baldwinian Model to decide what to do with
individuals generated after the local search.
ADVERTISEMENT
ADVERTISEMENT
Lamarckian Model
The Lamarckian Model essentially says that the traits which an individual acquires in
his/her lifetime can be passed on to its offspring. It is named after French biologist Jean-
Baptiste Lamarck.
Even though, natural biology has completely disregarded Lamarckism as we all know that
only the information in the genotype can be transmitted. However, from a computation
view point, it has been shown that adopting the Lamarckian model gives good results for
some of the problems.
In the Lamarckian model, a local search operator examines the neighborhood (acquiring
new traits), and if a better chromosome is found, it becomes the offspring.
Baldwinian Model
The Baldwinian model is an intermediate idea named after James Mark Baldwin (1896). In
the Baldwin model, the chromosomes can encode a tendency of learning beneficial
behaviors. This means, that unlike the Lamarckian model, we don’t transmit the acquired
traits to the next generation, and neither do we completely ignore the acquired traits like in
the Darwinian Model.
The Baldwin Model is in the middle of these two extremes, wherein the tendency of an
individual to acquire certain traits is encoded rather than the traits themselves.
In this Baldwinian Model, a local search operator examines the neighborhood (acquiring
new traits), and if a better chromosome is found, it only assigns the improved fitness to the
chromosome and does not modify the chromosome itself. The change in fitness signifies
the chromosomes capability to “acquire the trait”, even though it is not passed directly to
the future generations.
Effective Implementation
GAs are very general in nature, and just applying them to any optimization problem
wouldn’t give good results. In this section, we describe a few points which would help and
assist a GA designer or GA implementer in their work.
ADVERTISEMENT
ADVERTISEMENT
Switching to rank selection and tournament selection which have more selection
pressure than fitness proportionate selection for individuals with similar fitness.
ADVERTISEMENT
ADVERTISEMENT
Randomization Helps!
It has been experimentally observed that the best solutions are driven by randomized
chromosomes as they impart diversity to the population. The GA implementer should be
careful to keep sufficient amount of randomization and diversity in the population for the
best results.
It may be sometimes useful to hybridize the GA with local search. The following image
shows the various places in which local search can be introduced in a GA.
ADVERTISEMENT
ADVERTISEMENT
Using repair functions which take an infeasible solution and modify it so that the
violated constraints get satisfied.
Not allowing infeasible solutions to enter into the population at all.
ADVERTISEMENT
ADVERTISEMENT
Schema Theorem
Researchers have been trying to figure out the mathematics behind the working of genetic
algorithms, and Holland’s Schema Theorem is a step in that direction. Over the year’s
various improvements and suggestions have been done to the Schema Theorem to make it
more general.
In this section, we don’t delve into the mathematics of the Schema Theorem, rather we try
to develop a basic understanding of what the Schema Theorem is. The basic terminology
to know are as follows −
The schema theorem states that this schema with above average fitness, short defining
length and lower order is more likely to survive crossover and mutation.
Building Blocks are low order, low defining length schemata with the above given average
fitness. The building block hypothesis says that such building blocks serve as a foundation
for the GAs success and adaptation in GAs as it progresses by successively identifying and
recombining such “building blocks”.
Wolpert and Macready in 1997 published a paper titled "No Free Lunch Theorems for
Optimization." It essentially states that if we average over the space of all possible
problems, then all non-revisiting black box algorithms will exhibit the same performance.
It means that the more we understand a problem, our GA becomes more problem specific
and gives better performance, but it makes up for that by performing poorly for other
problems.
ADVERTISEMENT
ADVERTISEMENT
It should be kept in mind that the standard issue like crossover, mutation, Lamarckian or
Darwinian, etc. are also present in the GBML systems.
In this section, we list some of the areas in which Genetic Algorithms are frequently used.
These are −
Neural Networks − GAs are also used to train neural networks, particularly
recurrent neural networks.
Parallelization − GAs also have very good parallel capabilities, and prove to be very
effective means in solving certain problems, and also provide a good area for
research.
Image Processing − GAs are used for various digital image processing (DIP) tasks
as well like dense pixel matching.
Vehicle routing problems − With multiple soft time windows, multiple depots and a
heterogeneous fleet.
Scheduling applications − GAs are used to solve various scheduling problems as
well, particularly the time tabling problem.
Parametric Design of Aircraft − GAs have been used to design aircrafts by varying
the parameters and evolving better solutions.
DNA Analysis − GAs have been used to determine the structure of DNA using
spectrometric data about the sample.
Traveling salesman problem and its applications − GAs have been used to solve
the TSP, which is a well-known combinatorial problem using novel crossover and
packing strategies.
Get Started
Advertisements
Tutorials Point is a leading Ed Tech company
striving to provide the best learning material on
technical and non-technical subjects.
Affiliates eBooks
Contact Us
Tutorials Point India Private Limited, Incor9
Building, Kavuri Hills, Madhapur, Hyderabad,
Telangana - 500081, INDIA