Another Particle Swarm Toolbox: rastriginsfcn(x) - File Exchange

Code covered by the BSD License

Highlights from
Another Particle Swarm Toolbox

...
... Find the minimum of a function using Particle Swarm Optimization.
... Particle swarm optimization for binary genomes.
...
... This is like "soft" boundaries, except that some kind of penalty value
ackleysfcn(x) Ackley's Function.
binarytestfcn(x)
dejongsfcn(x)
dropwavefcn(x)
evolutioncomplete(options,sta... Plays a notification when genetic algorithm finishes. Requires a
griewangksfcn(x) Template for writing custom test/demonstration functions for psodemo.
langermannsfcn(x)
nonlinearconstrdemo(x) Nonlinear constraints demo with Rosenbrock's function. Select this
psoboundsabsorb(state,Aineq,b...
psocalculatepenalties(state) Calculates the penalty to apply for out-of-bounds particles when using
psodemo(DemoMode) Runs the PSO on a few demonstration functions, which should be located
psoiterate(options,state,flag... Updates swarm positions and velocities. Called to iterate the swarm from
psooptimset(varargin) Creates an options structure for pso.
psoplotbestf(options,state,fl... Plots the best, mean, and worst scores of particle swarm.
psoplotscorediversity(options... Plots a histogram containing the best and mean scores of particle swarm.
psoplotswarm(options,state,fl... Plots the positions of particle swarm.
psoplotswarmsurf(options,stat... Shows the evolution of a 2-variable population over a known surface. This
rastriginsfcn(x) Takes row inputs. If not row, will attempt to correct it.
rosenbrocksfcn(x) Takes row inputs. If input is not row, will attempt to correct it.
schwefelsfcn(x) Schwefel's function.
templatefcn(x) Template for writing custom test/demonstration functions for psodemo.
testfcn1(x)
View all files

from Another Particle Swarm Toolbox by Sam
Implementation of a PSO algorithm with the same syntax as the Genetic Algorithm Toolbox.

rastriginsfcn(x)

function y = rastriginsfcn(x)
% Takes row inputs. If not row, will attempt to correct it.
% Syntax:
% y = rastriginsfcn(x)
% options = rastringinsfcn('init')

if strcmp(x,'init')
    y.Aineq = [] ; y.bineq = [] ;
    y.Aeq = [] ; y.beq = [] ;
    y.LB = [-4,-4] ; y.UB = [4,4] ;
    y.nonlcon = [] ;
    y.options.PopInitRange = [-4, -4; 4, 4] ;
    y.options.KnownMin = [0,0] ;
    y.options.PopulationSize = 100 ;
    y.options.Vectorized = 'on' ;
    y.options.ConstrBoundary = 'absorb' ;
%     y.options.InitialPopulation = rand(100,2) ;
%     y.options.InitialVelocities = repmat([-1,0],100,1) ;
else
%     x = reshape(x,1,[]) ;
%     y = 10*size(x,2) + x*x' - 10*sum(cos(2*pi*x),2) ;
    
    y = 10*size(x,2) + sum(x.^2,2) - 10.*sum(cos(2*pi*x),2) ;
end

Highlights from Another Particle Swarm Toolbox

Highlights from
Another Particle Swarm Toolbox