Another Particle Swarm Toolbox: psoplotswarmsurf(options,state,flag) - File Exchange

Code covered by the BSD License

Highlights from
Another Particle Swarm Toolbox

... This is like "soft" boundaries, except that some kind of penalty value
...
... Particle swarm optimization for binary genomes.
...
... Find the minimum of a function using Particle Swarm Optimization.
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.
psoboundsabsorb(state,Aineq,b...
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.
psoplotscores(options,state,f... Plots 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.

psoplotswarmsurf(options,state,flag)

function state = psoplotswarmsurf(options,state,flag)
% Shows the evolution of a 2-variable population over a known surface. This
% plotting function is useful for observing the performance of the swarm
% over a given 2-D search space. It is used when running demonstration
% cases of pso.

notinf = isfinite(state.Score') ;
% notinf = ~state.OutOfBounds ;
if isfield(state,'Penalty')
    particleHeights = state.Score(notinf)' ; %- state.Penalty(notinf)' ;
else
    particleHeights = state.Score(notinf)' ;
end

if strcmpi(flag(1:4),'init') % Initialize
    delete(findobj(gca,'-regexp','Tag','*Locations'))
    
    line(state.Population(notinf,1),state.Population(notinf,2),...
        particleHeights,...
        'Color',0.75*ones(1,3),...
        'Marker','.',...
        'LineStyle','none',...
        'Tag','Initial Locations') ;
    
    % Set reasonable axes limits
    % ---------------------------------------------------------------------
    xlim([options.PopInitRange(1,1) options.PopInitRange(2,1)])
    ylim([options.PopInitRange(1,2) options.PopInitRange(2,2)])
    view(3)
    % ---------------------------------------------------------------------
    
    if ~strcmpi(options.DemoMode,'off')
        set(state.hfigure,'Units','normalized',...
            'Position',[0.2, 0.3, 0.6, 0.5])
        overlaysurface(state.fitnessfcn,options) ;
    end % if ~isempty
    
    line(state.Population(notinf,1),...
        state.Population(notinf,2),...
        state.Score(notinf)',...
        'LineStyle','none',...
        'Marker','.',...
        'Color','blue',...
        'Tag','Swarm Locations') ;
    set(gca,'Tag','Swarm Plot 3D',...
        'NextPlot','add',...
        'Projection','perspective')
    rotate3d(gca,'on')
elseif strcmpi(flag(1:4),'iter') % Iterate
    set(findobj(gca,'Tag','Swarm Locations','Type','line'),...
        'XData',state.Population(notinf,1),...
        'YData',state.Population(notinf,2),...
        'ZData',particleHeights)
end

Highlights from Another Particle Swarm Toolbox

Highlights from
Another Particle Swarm Toolbox