新建文件夹 (2).rar_matlab_optimization_ourselvesiyk_slowzk2_粒子群算法

clear all clc x0=-1:0.02:1; y0=-1:0.02:1; [X,Y]=meshgrid(x0,y0); Z=fun_zhl_pso2(X,Y); surf(X,Y,Z) view(35,30) % hold on % [T1,T2]=find(Z==max(Z(:))) % % % xmax=x0(T1); % % ymax=y0'(T2); % % zmax=fun_zhl_pso2(xmax,ymax); % % plot3(xmax,ymax,zmax,'r*') % plot3(X(T1'),Y(T2'),fun_zhl_pso2(X(T1'),Y(T2')),'r*') % % hold off % % %%　粒子群初始化 xMax=1; xMin=-1; vMax=0.2; dimSize=2; popSize=100; kmax=120; pBest=zeros(popSize,dimSize); gBest=zeros(popSize,dimSize); x=2*rand(popSize,dimSize)-1; v=2*rand(popSize,dimSize)-1; F=zeros(popSize,1); F1=zeros(popSize,1); mov1 = avifile('wave1.avi')%产生空的avi文件 for k=1:kmax for popIndex = 1:popSize F(popIndex ,1)=fun_zhl_pso2(x(popIndex ,1),x(popIndex ,2)); end % %%粒子速度更新 w =0.9-k*(0.9-0.4)/kmax ; %最大加权因子- (最大加权因子- 最小加权因子)*当前世代数/总世代数; subtract1 = pBest - x; subtract2 = gBest - x; v =w.*v+ 2.*rand().*subtract1 + 2.*rand().*subtract2 ; v(find(v>vMax))=vMax; v(find(v<-vMax))=-vMax; x=x+v; x(find(x>xMax))=xMax; x(find(x<xMin))=xMin; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 产生动画 figure(k+1) hold on surf(X,Y,Z) plot3(x(:,1),x(:,2),fun_zhl_pso2(x(:,1),x(:,2)),'*') grid on view(35,30) hold off FP= getframe(gcf);%获取图形句柄 mov1 = addframe(mov1,FP);%向avi文件中添加图形帧 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for popIndex = 1:popSize F1(popIndex )=fun_zhl_pso2(x(popIndex ,1),x(popIndex ,2)); end k1=find(F1>F); %找到F1>F的行序号 pBest(k1,1:dimSize)=x(k1,1:dimSize);%更新 pBest k2=min(find(F1==max(F1))); %找到F1为最大时的行序号之一 if max(F1)>max(F) %更新 gBest for popIndex = 1:popSize gBest(popIndex ,1:dimSize)=x(k2,1:dimSize); end end end mov1 = close(mov1);%关闭avi文件