pso.asv

粒子群优化算法用于天线阵列方向图综合

%nominal PSO
%no neibour best effct, no changable w,c1,c2
%fitness

%initialize
name='De Jong';
if strcmp(name,'De Jong') 
    itmax=5000;
    M=20;
    D=2;
    Vmax=2.048;
    Vmin=-2.048;
    Xmax=2.048;
    Xmin=-2.048;
    Opt=[1;1];
    Optvalue=3905.3;
elseif strcmp(name,'Goldstein')
    itmax=5000;
    M=20;
    D=2;
    Vmax=2;
    Vmin=-2;
    Xmax=2;
    Xmin=-2;
    Opt=[0;-1];
    Optvalue=3;   
elseif strcmp(name,'Branin')
     itmax=1000;
    M=20;
    D=2;
    Vmax=5;
    Vmin=-5;
    Xmax=10;
    Xmin=-5;
    Opt=[-22/7;12.275];
    Optvalue=0.3977272;   
elseif strcmp(name,'Martin')
     itmax=1000;
    M=20;
    D=2;
    Vmax=10;
    Vmin=0;
    Xmax=10;
    Xmin=0;
    Opt=[5;5];
    Optvalue=0;   
elseif strcmp(name,'Rosenbrock')
     itmax=1000;
    M=20;
    D=2;
    Vmax=1.2;
    Vmin=-1.2;
    Xmax=1.2;
    Xmin=-1.2;
    Opt=[1;1];
    Optvalue=0;   
elseif strcmp(name,'Hyper sphere')
    itmax=1000;
    M=20;
    D=6;
    Vmax=5.12;
    Vmin=-5.12;
    Xmax=5.12;
    Xmin=-5.12;
    Opt=[0;0;0;0;0;0];
    Optvalue=0;    
end
X=zeros(itmax,M,D);
FX=zeros(itmax,M);
V=zeros(itmax,M,D);  
Pbest=zeros(itmax,M,D);
Gbest=zeros(itmax,D);
FPbest=zeros(itmax,M);
FGbest=zeros(itmax,1);
for i=1:M
    X(1,i,:)=(Xmax-Xmin)*rand(D,1)+Xmin.*ones(D,1);
    FX(1,i)=Fitness(name,X(1,i,:));
    V(1,i,:)=(Vmax-Vmin)*rand(D,1)+Vmin.*ones(D,1);
    Pbest(1,i,:)=X(1,i,:);
    FPbest(1,i)=FX(1,i);
end
[FGbest(1),id]=max(FPbest(1,:));
Gbest(1,:)=Pbest(1,id,:);
w=0.9;
c1=2;
c2=2;
    
Threshold=0.001;
Xbest=zeros(D,1);
Xbest=X(1,1,:);
FXbest=FX(1,1);
    
trial=1;
AvgTr=zeros(trial,1);
%main body
for j=1:trial
for t=2:itmax
    for i=1:M
       for d=1:D
           V(t,i,d)=w*V(t-1,i,d)+c1*rand()*(Pbest(t-1,i,d)-X(t-1,i,d))+c2*rand()*(Gbest(t-1,d)-X(t-1,i,d));
           V(t,i,d)=min([Vmax,max([Vmin,V(t,i,d)])]);
           X(t,i,d)=min([Xmax,max([Xmin,X(t-1,i,d)+V(t,i,d)])]);
       end
       FX(t,i)=Fitness(name,X(t,i,:));
       if FX(t,i)>FPbest(t-1,i)
           Pbest(t,i,:)=X(t,i,:);
           FPbest(t,i)=FX(t,i);
           if FX(t,i)>FGbest(t-1)
               Gbest(t,:)=X(t,i,:);
               FGbest(t)=FX(t,i);
           end
       end  
    end
    if abs(FGbest(t)-Optvalue)<=Optvalue*Threshold
        Xbest=Gbest(t,:);
        FXbest=FGbest(t);
        break;
    end
end
AvgTr(j)=t;
Xbest=Gbest(t,:);
FXbest=FGbest(t);
end
AvgIter=mean(AvgTr);
%plot result