Global Optimization Toolbox 

This example shows how to perform a multiobjective optimization using multiobjective genetic algorithm function gamultiobj in Global Optimization Toolbox.
On this page… 

Simple Multiobjective Optimization Problem 
Simple Multiobjective Optimization Problem
gamultiobj can be used to solve multiobjective optimization problem in several variables. Here we want to minimize two objectives, each having one decision variable.
min F(x) = [objective1(x); objective2(x)] x
where, objective1(x) = (x+2)^2  10, and objective2(x) = (x2)^2 + 20
% Plot two objective functions on the same axis x = 10:0.5:10; f1 = (x+2).^2  10; f2 = (x2).^2 + 20; plot(x,f1); hold on; plot(x,f2,'r'); grid on; title('Plot of objectives ''(x+2)^2  10'' and ''(x2)^2 + 20''');
The two objectives have their minima at x = 2 and x = +2 respectively. However, in a multiobjective problem, x = 2, x = 2, and any solution in the range 2 <= x <= 2 is equally optimal. There is no single solution to this multiobjective problem. The goal of the multiobjective genetic algorithm is to find a set of solutions in that range (ideally with a good spread). The set of solutions is also known as a Pareto front. All solutions on the Pareto front are optimal.
We create a MATLABfile named simple_multiobjective.m:
function y = simple_multiobjective(x) y(1) = (x+2)^2  10; y(2) = (x2)^2 + 20;
The Genetic Algorithm solver assumes the fitness function will take one input x, where x is a row vector with as many elements as the number of variables in the problem. The fitness function computes the value of each objective function and returns these values in a single vector output y.
To use the gamultiobj function, we need to provide at least two input arguments, a fitness function, and the number of variables in the problem. The first two output arguments returned by gamultiobj are X, the points on Pareto front, and FVAL, the objective function values at the values X. A third output argument, exitFlag, tells you the reason why gamultiobj stopped. A fourth argument, OUTPUT, contains information about the performance of the solver. gamultiobj can also return a fifth argument, POPULATION, that contains the population when gamultiobj terminated and a sixth argument, SCORE, that contains the function values of all objectives for POPULATION when gamultiobj terminated.
FitnessFunction = @simple_multiobjective; numberOfVariables = 1; [x,fval] = gamultiobj(FitnessFunction,numberOfVariables);
Optimization terminated: average change in the spread of Pareto solutions less than options.TolFun.
The X returned by the solver is a matrix in which each row is the point on the Pareto front for the objective functions. The FVAL is a matrix in which each row contains the value of the objective functions evaluated at the corresponding point in X.
size(x) size(fval)
ans = 18 1 ans = 18 2
Constrained Multiobjective Optimization Problem
gamultiobj can handle optimization problems with linear inequality, equality, and simple bound constraints. Here we want to add bound constraints on simple multiobjective problem solved previously.
min F(x) = [objective1(x); objective2(x)] x
subject to 1.5 <= x <= 0 (bound constraints)
where, objective1(x) = (x+2)^2  10, and objective2(x) = (x2)^2 + 20
gamultiobj accepts linear inequality constraints in the form A*x <= b and linear equality constraints in the form Aeq*x = beq and bound constraints in the form lb < x < ub. We pass A and Aeq as matrices and b, beq, lb, and ub as vectors. Since we have no linear constraints in this example, we pass [] for those inputs.
A = []; b = []; Aeq = []; beq = []; lb = 1.5; ub = 0; x = gamultiobj(FitnessFunction,numberOfVariables,A,b,Aeq,beq,lb,ub);
Optimization terminated: maximum number of generations exceeded.
All solutions in X (each row) will satisfy all linear and bound constraints within the tolerance specified in options.TolCon. However, if you use your own crossover or mutation function, ensure that the new individuals are feasible with respect to linear and simple bound constraints.
gamultiobj can accept one or more plot functions through the options argument. This feature is useful for visualizing the performance of the solver at run time. Plot functions can be selected using gaoptimset. The help for gaoptimset contains a list of plot functions to choose from.
Here we use gaoptimset to create an options structure to select two plot functions. The first plot function is GAPLOTPARETO, which plots the Pareto front (limited to any three objectives) at every generation. The second plot function is GAPLOTSCOREDIVERSITY, which plots the score diversity for each objective. The options structure is passed as the last argument to the solver.
options = gaoptimset('PlotFcns',{@gaplotpareto,@gaplotscorediversity});
gamultiobj(FitnessFunction,numberOfVariables,[],[],[],[],lb,ub,options);
Optimization terminated: maximum number of generations exceeded.
Vectorizing Your Fitness Function
Consider the previous fitness functions again:
objective1(x) = (x+2)^2  10, and objective2(x) = (x2)^2 + 20
By default, the gamultiobj solver only passes in one point at a time to the fitness function. However, if the fitness function is vectorized to accept a set of points and returns a set of function values you can speed up your solution.
For example, if the solver needs to evaluate five points in one call to this fitness function, then it will call the function with a matrix of size 5by1, i.e., 5 rows and 1 column (recall that 1 is the number of variables).
Create a MATLABfile called vectorized_multiobjective.m:
function scores = vectorized_multiobjective(pop) popSize = size(pop,1); % Population size numObj = 2; % Number of objectives % initialize scores scores = zeros(popSize, numObj); % Compute first objective scores(:,1) = (pop + 2).^2  10; % Compute second obective scores(:,2) = (pop  2).^2 + 20;
This vectorized version of the fitness function takes a matrix 'pop' with an arbitrary number of points, the rows of 'pop', and returns a matrix of size populationSizebynumberOfObjectives.
We need to specify that the fitness function is vectorized using the options structure created using gaoptimset. The options structure is passed in as the ninth argument.
FitnessFunction = @(x) vectorized_multiobjective(x); options = gaoptimset('Vectorized','on'); gamultiobj(FitnessFunction,numberOfVariables,[],[],[],[],lb,ub,options);
Optimization terminated: average change in the spread of Pareto solutions less than options.TolFun.