pso_common.m

PSO算法

function fout=pso_common[functionname,D,varargin]
%%
%out是函数输出值；
%functionname是待求的函数名称，是字符变量；
%Dimen是待求函数的维数；
%varargin是永久变量。在本程序中它包括 mxvel、VarRange、minormax、p、psoseedval
if nargin<2
    error('输入变量不够。');
end

if nargin==2 %只确定了functionname及D，还要定义其他默认参数
    mv       =4;
    VarRange =[ones(D,1)*(-100),ones(D,1)*100;];
    minormax =0;
    Param    =[];
elseif nargin==3 %确定了functionname、D及mv，还要定义其他默认参数
    mv       =varargin{1};
    VarRange =[ones(D,1)*(-100),ones(D,1)*100;];
    minormax =0;
    Param    =[];
elseif nargin==4 %确定了functionname、D、mv及VarRange，还要定义其他默认参数
    mv       =varargin{1};
    VarRange =varargin{2};
    minormax =0;
    Param    =[];
elseif nargin==5 %还要确定p[]
    mv       =varargin{1};
    VarRange =varargin{2};
    minormax =varargin{3};
    Param    =[];
elseif nargin==6
    mv       =varargin{1};
    VarRange =varargin{2};
    minormax =varargin{3};
    Param    =varargin{4};
elseif nargin==7
    mv       =varargin{1};
    VarRange =varargin{2};
    minormax =varargin{3};
    Param    =varargin{4};
    seedval  =varargin{5};
else
    error('输入变量错误!')
end
%%
%下面设置Param的默认值及在具体应用中它的值。
Pdefault     =[24,2000,2,2,0.9,0.4,1500,1e-25,250,nan,0];
Plength      =length(Param);
Param        =[Param,Pdefault(Plength+1:end)];

ps =Param(1);
me      =Param(2);
acc1    =Param(3);
acc2    =Param(4);
iwt1    =Param(5);
iwt2    =Param(6);
iwte    =Param(7);
ergrd   =Param(8);
ergrdep =Param(9);
errgoal =Param(10);
seedpos =Param(11);
% 上面设置好了默认的Param参数的默认值及在具体应用中它的值，并且给它的每个元素限定了含义。

%%
% 给出速度范围
if length(mv)==1
    velmaskmin =-mv*ones(ps,D);
    velmaskmax =mv*ones(ps,D);
elseif length(mv)==D
    velmaskmin =repmat(forcerow(-mv),ps,1);
    velmaskmax =repmat(forcerow(mv),ps,1);
else 
    error('Max velocity must be either a scalar or the same length as problem dimension D');
end
%%
%给出微粒的位置范围
posmaskmin  =repmat(VarRange(1:D,1)',ps,1);
posmaskmax  =repmat(VarRange(1:D,2)',ps,1);
%%
message = sprintf('PSO: %%g/%g iterations, GBest = %%20.20g.\n',me);
%%
%初始化开始
pos(1:ps,1:D) = normmat(rand([ps,D]),VR',1);

























    return