VectorGA: A Vectorized Implementation of a Genetic Algorithm in Matlab

VectorGA is a vectorized implementation of a genetic algorithm in the Matlab programming language.

Matlab is optimized for performing operations on arrays. Loops, especially nested loops, tend to run slowly in Matlab. It is possible to significantly improve the performance of Matlab programs by converting loops into array operations. This process is called vectorization. Matlab provides a rich set of functions and many expressive indexing schemes that make it possible to vectorize code. Such code not only runs faster, it is also shorter, and simpler to understand and change (provided that you know a little about Matlab of course).

Genetic Algorithms that are implemented in C/C++ or Java typically have multiple nested loops. Therefore direct ports of such implementations to Matlab will run very slowly. Many of the nested loops found in a typical GA implementation have been eliminated from VectorGA. The resulting code is short, fast and simple. It is a neat coincidence when the constructs of a programming language match a programming task so well that a program can be written this succintly.

VectorGA is proof that Matlab is a useful language for the rapid prototyping of Genetic Algorithms. This, in addition to Matlab’s extensive data visualization capabilities, make Matlab an extremely useful platform for the experimental analysis of GAs.

VectorGA has been created and tested under Matlab 7 (R14). Out of the box it evolves a population against the one-max fitness function. The royal-roads fitness function has also been included but is not currently being called. If you find vectorGA useful or find any bugs please let me know.

Enjoy!

Original vectorGA release site

The latst version of VectorGA is SpeedyGA

VectorGA: A Vectorized Implementation of a Genetic Algorithm in Matlab

6 thoughts on “VectorGA: A Vectorized Implementation of a Genetic Algorithm in Matlab

  1. Masaya says:

    It is very useful.
    I found a small bug.Please check it.

    line 43
    %%%%%%%%%%%%%%%%%%%%%%%%
    if size(nargin==1)
    rand(‘state’,varargin{:});
    else
    rand(‘state’,sum(100*clock));
    end
    %%%%%%%%%%%%%%%%%%%%%%%%

                    ↓

    %%%%%%%%%%%%%%%%%%%%%%%%
    if size(varargin)==1
    rand(‘state’,varargin{:});
    else
    rand(‘state’,sum(100*clock));
    end
    %%%%%%%%%%%%%%%%%%%%%%%%

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s