%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%% PSO工具箱函数 %%%%%%%%%%%%%%%

pso PSO主程序

psoopt PSO参数项

psooptimset　 PSO参数项设置

psooptimget　 PSO参数项获取

psoplot　　　 PSO图示

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

各函数调用格式

% [x,fval,output]=pso(fun,nvars);

% [x,fval,output]=pso(fun,nvars,options);

% optoins = psodefaultopt;

% options = psodefaultopt('default')

% options = psooptimset(options,defaultopt);

% options = psooptimset(options,value);

% options = psooptimset(options,name,value);

% options = psooptimset(options,name1,value1,name2,value2,...);

% options = psooptimget('default');

% [name,values]=psooptimget(options);

% value = psooptimget(options,name);

% psoplot(output,options)

% h = psoplot(output,options)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

以下给出PSO主函数程序代码

function [x,fval,output]=pso(fun,nvars,options)

% PSO

% [x,fval,output]=pso(fun,nvars);

% [x,fval,output]=pso(fun,nvars,options);

% fun: function handle

% nvars: dimension of problem

% options: pso options setting

%

% call subfunction: bestinitialize, updatepbest, updategbest,

% updateparticle, chkerrgoal, limitpos, limitvel

% See also: psoopt, psoopitmset, psooptimget, readjust, psoplot

% Author: J.H.Zhu

% Date: Nov. 1, 2010.

% Copyright reserved by Author

% Check number of input arguments

error(nargchk(1,3,nargin));

if nargin<3

options

=[]; end

defaultopt = psoopt('default');

% Take defaults for parameters that are not in options

structure

options = psooptimset(options,defaultopt);

rand('state',options.RandState);

%#ok

popsize = options.PopulationSize;

poprange = options.PopRange;

partvel = options.ParticleVelocity;

% initialize population of particles and their velocities

pos = readjust(rand(popsize,nvars),poprange,1);

vel = readjust(rand(popsize,nvars),partvel,1);

% initial partial and global best pos

[pbest,pbestval,gbest,gbestval]=bestinitialize(fun,pos,options);

output.gbest(1,:)=gbest;

output.gbestval(1)=gbestval;

output.counter=0;

counter2=0;

for i = 1:options.MaxEpoch

for j

=1:options.PopulationSize

funval(j) = feval_r(fun,pos(j,:));

% update pbest

[pbest(j,:),pbestval(j)]=...

updatepbest(pbest(j,:),pbestval(j),pos(j,:),funval(j),options);

% update gbest

[gbest,gbestval]=...

updategbest(gbest,gbestval,pbest,pbestval,options);

% update velocities and positions

[pos(j,:),vel(j,:)]=...

updateparticle(i,pos(j,:),vel(j,:),gbest,pbest(j,:),options);

end % end

popsize

% check the condition of end loop

output.gbest(i,:)=gbest;

output.gbestval(i)=gbestval;

output.counter=i;

[ExitFlag,counter2]=chkerrgoal(output,counter2,options);

% figure plots

if

~isempty(options.PlotFcn) &&

options.Display

feval_r(options.PlotFcn,output,options);

end

if

(~ExitFlag)

x = gbest;

fval = gbestval;

break;

end

end %% end maxepoch

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

psoplot画出的图形