【導(dǎo)讀】算法通過(guò)粒子群中粒子間的合作在復(fù)雜的搜索空間中找到全局的最優(yōu)解。注，已經(jīng)成為發(fā)展最快的優(yōu)化算法之一。論文介紹了粒子濾波算法的原理、特。圍繞粒子濾波算法的原理特點(diǎn)和參數(shù)的設(shè)置進(jìn)行了各方面的闡述和論證，參數(shù)，論證評(píng)估了各種改進(jìn)粒子濾波算法的方法的性能和可行性。濾波算法的研究和應(yīng)用提出了一些建議和展望。

　　

【正文】 39。 fv = fitness(pg)。 r=(dot(xm,xm))^。 程序 ： function [xm,fv,r] = LnCPSO(fitness,N,cmax,cmin,w,M,D) format long。 %179。245。202。188。187。175。214。214。200。186。181。196。184。246。204。229。 for i=1:N for j=1:D x(i,j)=randn。 %203。230。187。179。245。202。188。187。175。206。187。214。195。 v(i,j)=randn。 %203。230。187。179。245。202。188。187。175。203。217。182。200。 end end 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 37 %207。200。188。198。203。227。184。247。184。246。193。163。211。181。196。202。202。211。166。182。200。163。172。178。162。179。245。202。188。187。175。Pi186。205。Pg for i=1:N p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end pg = x(N,:)。 %Pg206。170。200。171。190。214。238。211。197。 for i=1:(N1) if fitness(x(i,:))fitness(pg) pg=x(i,:)。 end end %189。248。200。235。214。247。210。170。209。 187。163。172。176。180。213。213。185。171。202。189。210。192。180。206。181。252。180。 for t=1:M c = cmax (cmax cmin)*t/M。 for i=1:N 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 38 v(i,:)=w*v(i,:)+c*rand*(y(i,:)x(i,:))+c*rand*(pgx(i,:))。 x(i,:)=x(i,:)+v(i,:)。 if fitness(x(i,:))p(i) p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end if p(i)fitness(pg) pg=y(i,:)。 end end Pbest(t)=fitness(pg)。 end xm = pg39。 fv = fitness(pg)。 r=(dot(xm,xm))^。 程序 8 異步變化學(xué)習(xí)因子的粒子濾波算法： function [xm,fv,r] = AsyLnCPSO(fitness,N,c1max,c1min,c2max,c2min,w,M,D) 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 39 format long。 %179。245。202。188。187。175。214。214。200。186。181。196。184。246。204。229。 for i=1:N for j=1:D x(i,j)=randn。 %203。230。187。179。245。202。188。187。175。206。187。214。195。 v(i,j)=randn。 %203。230。187。179。245。202。188。187。175。203。217。182。200。 end end %207。200。188。198。203。227。184。247。184。246。193。163。211。181。196。202。202。211。166。182。200。163。172。178。162。179。245。202。188。187。175。Pi186。205。Pg for i=1:N p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end pg = x(N,:)。 %Pg206。170。200。171。190。214。238。211。197。 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 40 for i=1:(N1) if fitness(x(i,:))fitness(pg) pg=x(i,:)。 end end %189。248。200。235。214。247。210。170。209。 187。163。172。176。180。213。213。185。171。202。189。210。192。180。206。181。252。180。 for t=1:M c1 = c1max (c1max c1min)*t/M。 c2 = c2max (c2max c2min)*t/M。 for i=1:N v(i,:)=w*v(i,:)+c1*rand*(y(i,:)x(i,:))+c2*rand*(pgx(i,:))。 x(i,:)=x(i,:)+v(i,:)。 if fitness(x(i,:))p(i) p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 41 end if p(i)fitness(pg) pg=y(i,:)。 end end Pbest(t)=fitness(pg)。 end xm = pg39。 fv = fitness(pg)。 r=(dot(xm,xm))^。 程序 9：二階粒子濾波算法 function [xm,fv,r] = SecPSO(fitness,N,w,c1,c2,M,D) format long。 %179。245。202。188。187。175。214。214。200。186。181。196。184。246。204。229。 for i=1:N for j=1:D x(i,j)=randn。 %203。230。187。179。245。202。188。187。175。206。187。214。195。 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 42 xl(i,j)=randn。 v(i,j)=randn。 %203。230。187。179。245。202。188。187。175。203。217。182。200。 end end %207。200。188。198。203。227。184。247。184。246。193。163。211。181。196。202。202。211。166。182。200。163。172。178。162。179。245。202。188。187。175。Pi186。205。Pg for i=1:N p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end pg = x(N,:)。 %Pg206。170。200。171。190。214。238。211。197。 for i=1:(N1) if fitness(x(i,:))fitness(pg) pg=x(i,:)。 end 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 43 end %189。248。200。235。214。247。210。170。209。 187。163。172。176。180。213。213。185。171。202。189。210。192。180。206。181。252。180。 for t=1:M for i=1:N v(i,:)=w*v(i,:)+c1*rand*(y(i,:)2*x(i,:)+xl(i,:))+c2*rand*(pg2*x(i,:)+ xl(i,:))。 xl(i,:) = x(i,:)。 x(i,:)=x(i,:)+v(i,:)。 if fitness(x(i,:))p(i) p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end if p(i)fitness(pg) pg=y(i,:)。 end 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 44 end Pbest(t)=fitness(pg)。 end xm = pg39。 fv = fitness(pg)。 r=(dot(xm,xm))^。 程序 10：二階振蕩粒子濾波算法 function [xm,fv,r] = SecVibratPSO(fitness,N,w,c1,c2,M,D) format long。 %179。245。202。188。187。175。214。214。200。186。181。196。184。246。204。229。 for i=1:N for j=1:D x(i,j)=randn。 %203。230。187。179。245。202。188。187。175。206。187。214。195。 xl(i,j)=randn。 v(i,j)=randn。 %203。230。187。179。245。202。188。187。175。203。217。182。200。 end 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 45 end %207。200。188。198。203。227。184。247。184。246。193。163。211。181。196。202。202。211。166。182。200。163。172。178。162。179。245。202。188。187。175。Pi186。205。Pg for i=1:N p(i)=fitness(x(i,:))。 y(i,:)=x(i,:)。 end pg = x(N,:)。 %Pg206。170。200。171。190。214。238。211。197。 for i=1:(N1) if fitness(x(i,:))fitness(pg) pg=x(i,:)。 end end %189。248。200。235。214。247。210。170。209。 187。163。172。176。180。213。213。185。171。202。189。210。192。180。206。181。252。180。 for t=1:M for i=1:N 安徽財(cái)經(jīng)大學(xué)管理科學(xué)與工程學(xué)院本科畢業(yè)論文 46 phi1 = c1*rand()。 phi2 = c2*rand()。 if t M/2 ks1 = (2*sqrt(phi1)1)*rand()/phi1。 ks2 = (2*sqrt(phi2)1)*rand()/phi2。 else ks1 = (2*sqrt(phi1)1)*(1 + rand())/phi1。 ks2 = (2*sqrt(phi2)1)*(1 + rand())/phi2。 end v(i,:)=w*v(i,:)+phi1*(y(i,:)(1+ks1)*x(i,:)+ks1*xl(i,:))+ ... phi2*(pg(1+ks2)*x(i,:)+ks1*xl(i,:))。 xl(i,:) = x(i,:)。 x(i,:)=x(i,:)+v(i,:)。 if fitness(x(i,:))p(i) p(i)=fitnes