粒子群算法实验报告

Windows

Matlab

PSO

ParticleSwarmOptimization

(Evolu2tionaryAlgorithm-EA

(PSO)

(evolutionarycomputation)1995Eberhart

Kennedy

PSO

(crossover)

PSO

(mutation)

.m

matlab

clearall;

cic;

formatlong;

%——

c1=1.4902;

c2=1.4901;

w=0.7281;

MaxDT=1000;

D=5;

N=40;

eps=10A(-6);

%——

fori=1:N

forj=1:D

x(ij)=randn;%

end

end

%——

% 1

% 2

%

%

%

,(1

v(ij)=randn;%

PiPg fori=1:N

p(i)=function(x(i,:));y(i,:)=x(i,:);

end

Solution=

0.576475699699576

10cos(2八rj+10) ( 1

Rastrigrin

pg=x(1,:); %Pg matlab

fori=2:N

iffunction(x(i,:))<function(pg)pg=x(i,：)；

end

%——

fort=1:MaxDT

fori=1:N

v(i,:)=w*v(i,:)+crrand*(y(i,:)-x(i5:))+c2*rand*(pg-x(i,:));x(i5:)=x(i5:)+v(i,:);

iffunction(x(i,:))<p(i)

p(i)=function(x(i,:));

y(i,：)=x(i,:);

end

ifp(i)<function(pg)

pg=y&：);

end

end

Pbest(t)=function(pg);end

% disp(disp(' ')

Solution=pg'

plot(Solution)

disp(' ')

Result=function(pg)disp(

Rastrigins.m Rastrigins

function[out]=Rastrigin(x)

x=-5.12:0.01:5.12;

cos_in=cos(2*pi*x);

out=sum((x.A2-10*cos_in+10),2);

Rastrigrin f

V=0tH /(jt*)=0,EZo1维Rastrigrin 4

(2)

Griewank

Griewank

Griewank.m

Dx=length(in(1,:));

tlenx=length(in(:,1));

ifisempty(D)|D~=Dx|tlen~=tlenx

D=Dx;%dimensionofprob

tlen=tlenx;%howmanyseparatestatesd=repmat([1:D],tlen,1);sqrtd=sqrt(d);

end

dat1=sum([(in-100).a2],2)./4000;

dat2=prod((cos((in-100)./sqrtd)),2);

out=datl-dat2+1;

functionGriewankmatlab

J/=0,1

4,7

4.7Griewankp

(4.10)|

2 a Giwmik

(b)伞局:

Solution=

1.0e+002*

0.089393835024922

1.355075075547163

0.945667645113542

1.188118773920475

0.929927307049068

Result=

2.497904795831135

41 」 I I 」 j I 」

1 16 2 2& 3 35 4 45 5

⑶ 3 Foxhole Foxhole.m

function[out]=Foxhole(in)%x=in(,1);%y=in(,2);term_sum=0;x=in(:,1);y=in(:,1);

a{1}=[..•

32-1601632

32-1601632

32-1601632

32-16

32

-16

16

16

32

32

];

啊..

32-32--32-

32 3

-16-16-16-16-16；

00000;...

1616161616;...

3232323232;...

]；

term_sum=O;

forj=1:numel((a{1}))

ax=a{1}(j);

ay=a{2}(j);

term_sum=(x-sx).a6+(y-ay).A6;term_sum=term_sum+1.0/(j+term_sum);endout=.002+term_sum;

•••

Gricwank

Z

tmk

Z

4.7

八二0,

IUIU-J

Solution=

31.999503320926419

6.742047876319