基于MATLAB实现的PSO和POWELL相结合的快速搜索算法，有很强的实用性+使用说明文档.rar

function [pa,mi]=PSO(R,F) %粒子群算法 %初始化 D=3; %维数 ps=20; %种群规模 VRmin=ones(D,1)*-20; %最小速度 VRmax=ones(D,1)*20; %最大速度 VR=[VRmin,VRmax]; %minmax = 1; pos=40*rand(ps,D)-20; %随机产生初始位置 vel=8*rand(ps,D)-4; %产生随机速度 % ps=15; %种群规模 % VRmin=ones(D,1)*-10; %最小速度 % VRmax=ones(D,1)*10; %最大速度 % VR=[VRmin,VRmax]; % %minmax = 1; % pos=20*rand(ps,D)-10; %随机产生初始位置 % vel=4*rand(ps,D)-2; %产生随机速度 %一些参数 maxinterations=20; %最大迭代次数 iw=1; %固定权重 iw1 = 0.9; % 最大惯性权重 iw2 = 0.4; iwe = 15; ac1=2; ac2=2; flagg=1; % ergrd=1e-5; % ergrdep=5; %？？？？ % mv=4;% ergrd=1e-3; ergrdep=5; %？？？？ mv=4;% %初始个体极值 pbest=pos; %求初始全局极值 for i=1:ps p=pos(i,:); %第i个粒子位置 out(i)=pv_mi(p(1),p(2),p(3),R,F); %求函数值 end pbestval=out; %每个粒子当前函数值 [gbestval,idx]=max(pbestval); %全局最优函数值 gbest=pbest(idx,:); %全局极值 tr(1)=gbestval; %保存当前全局最优函数值 % start PSO iterative procedures cnt=0; % counter used for updating display according to df in the options cnt2=0; % counter used for the stopping subroutine based on error convergence %开始迭代 for i=1:maxinterations for j=1:ps if flagg==1 % randomization control, one random set for each particle at each epoch rannum1=rand(1); rannum2=rand(1); end p=pos(j,:); %第i个粒子位置 out(j)=pv_mi(p(1),p(2),p(3),R,F); %求函数值 e(j)=out(j); %更新pbest if pbestval(j)<=e(j) pbestval(j)=e(j); pbest(j,:)=pos(j,:); end %更新gbest [iterbestval,idx1]=max(pbestval); if gbestval<=iterbestval gbestval=iterbestval; gbest=pbest(idx1,:); end tr(i+1)=gbestval; te=i; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%更新速度，位置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if i<=iwe iwt(i)=((iw2-iw1)/(iwe-1))*(i-1)+iw1; %惯性权重 else iwt(i)=iw2; end % iwt(i)=1; %%%%%%%%%%速度%%%%%%%%%%%%%%%%%%%% if flagg==2 % 粒子的每个参数的随机数不同 for dimcnt=1:D rannum1=rand(1); rannum2=rand(1); vel(j,dimcnt)=iwt(i)*vel(j,dimcnt)... +ac1*rannum1*(pbest(j,dimcnt)-pos(j,dimcnt))... +ac2*rannum2*(gbest(1,dimcnt)-pos(j,dimcnt)); end else % 粒子的每个参数的随机数相同 vel(j,:)=iwt(i)*vel(j,:)... +ac1*rannum1*(pbest(j,:)-pos(j,:))... +ac2*rannum2*(gbest(1,:)-pos(j,:)); end % 固定权重 % vel(j,:)=iw*vel(j,:)... % +ac1*rannum1*(pbest(j,:)-pos(j,:))... % +ac2*rannum2*(gbest(1,:)-pos(j,:)); %%%%%%%%%%位置%%%%%%%%%%%%%%%%%%%%%%% pos(j,:)=pos(j,:)+vel(j,:); %%%%%%%%%速度和位置范围％％％％％％％％％ for dimcnt=1:D if vel(j,dimcnt)>mv vel(j,dimcnt)=mv; end if vel(j,dimcnt)<-mv vel(j,dimcnt)=-mv; end if pos(j,dimcnt)>=VR(dimcnt,2) pos(j,dimcnt)=VR(dimcnt,2); end if pos(j,dimcnt)<=VR(dimcnt,1) pos(j,dimcnt)=VR(dimcnt,1); end end end %结束一次粒子循环 temp=gbest'; fprintf('%f,%f,%f,%f\n',temp(1),temp(2),temp(3),gbestval); Y(i)=gbestval; X(i)=i; % 收敛条件 tmp1=abs(tr(i)-gbestval); if tmp1>ergrd cnt2=0; elseif tmp1<=ergrd cnt2=cnt2+1; if cnt2>=ergrdep break end end end %迭代结束 fprintf('total interations:%d\n',i); %OUT=[gbest';gbestval]; pa=gbest; mi=gbestval; %plot(X,Y);