MATLAB－PSO工具箱

2017-01-16 by:CAE仿真在線來源:互聯網

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

%%%%%%%%%%% 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<jhzhu99@gmail.com>
% 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<RAND>
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畫出的圖形

MATLAB-PSO工具箱