用MATLAB编制PSO算法程序求解简单的函数极值问题(2)

while k<=MaxNum

for i=1:particlesize

for j=1:narvs

f(i)=fitness(x(i,j));

end

if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置 personalbest_faval(i)=f(i);

personalbest_x(i,:)=x(i,:);

end

end

[globalbest_faval i]=min(personalbest_faval);

globalbest_x=personalbest_x(i,:);

for i=1:particlesize %更新粒子群里每个个体的最新位置

v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))... +c2*rand*(globalbest_x-x(i,:));

for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度

if v(i,j)>vmax;

v(i,j)=vmax;

elseif v(i,j)<-vmax;

v(i,j)=-vmax;

end

end

x(i,:)=x(i,:)+v(i,:);

end

if abs(globalbest_faval)<E0,break,end

k=k+1;

end

Value1=1/globalbest_faval-1; Value1=num2str(Value1);

% strcat指令可以实现字符的组合输出

disp(strcat('the maximum value','=',Value1));

%输出最大值所在的横坐标位置

Value2=globalbest_x; Value2=num2str(Value2);

disp(strcat('the corresponding coordinate','=',Value2)); x=-5:0.01:5;

y=2.1*(1-x+2*x.^2).*exp(-x.^2/2);

plot(x,y,'m-','linewidth',3);

holdon;

plot(globalbest_x,1/globalbest_faval-1,'kp','linewidth',4);

legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on; toc;

运行结果为：