优化算法——粒子群算法(PSO)

时间：2022-12-11 07:34:26

一、粒子群算法的概述

    粒子群算法(PSO)属于群智能算法的一种，是通过模拟鸟群捕食行为设计的。假设区域里就只有一块食物（即通常优化问题中所讲的最优解），鸟群的任务是找到这个食物源。鸟群在整个搜寻的过程中，通过相互传递各自的信息，让其他的鸟知道自己的位置，通过这样的协作，来判断自己找到的是不是最优解，同时也将最优解的信息传递给整个鸟群，最终，整个鸟群都能聚集在食物源周围，即我们所说的找到了最优解，即问题收敛。

二、粒子群算法的流程

    粒子群算法通过设计一种无质量的粒子来模拟鸟群中的鸟，粒子仅具有两个属性：速度优化算法——粒子群算法(PSO)和位置优化算法——粒子群算法(PSO)，速度代表移动的快慢，位置代表移动的方向。每个粒子在搜索空间中单独的搜寻最优解，并将其记为当前个体极值优化算法——粒子群算法(PSO)，并将个体极值与整个粒子群里的其他粒子共享，找到最优的那个个体极值作为整个粒子群的当前全局最优解优化算法——粒子群算法(PSO)，粒子群中的所有粒子根据自己找到的当前个体极值优化算法——粒子群算法(PSO)和整个粒子群共享的当前全局最优解优化算法——粒子群算法(PSO)来调整自己的速度和位置。粒子群算法的思想相对比较简单，主要分为：1、初始化粒子群；2、评价粒子，即计算适应值；3、寻找个体极值优化算法——粒子群算法(PSO)；4、寻找全局最优解优化算法——粒子群算法(PSO)；5、修改粒子的速度和位置。下面是程序的流程图：
优化算法——粒子群算法(PSO)
 （PSO流程）
下面我们具体解释下流程图里面的每一个步骤：

1、初始化

   首先，我们需要设置最大的速度区间，防止超出最大的区间。位置信息即为整个搜索空间，我们在速度区间和搜索空间上随机初始化速度和位置。设置群体规模优化算法——粒子群算法(PSO)。

2、个体极值与全局最优解

   个体极值为每个粒子找到的历史上最优的位置信息，并从这些个体历史最优解中找到一个全局最优解，并与历史最优解比较，选出最佳的作为当前的历史最优解。

3、更新速度和位置的公式

   更新公式为： 优化算法——粒子群算法(PSO)
 优化算法——粒子群算法(PSO)
 其中，优化算法——粒子群算法(PSO)称为惯性因子，优化算法——粒子群算法(PSO)和优化算法——粒子群算法(PSO)称为加速常数，一般取优化算法——粒子群算法(PSO)。优化算法——粒子群算法(PSO)表示区间优化算法——粒子群算法(PSO)上的随机数。优化算法——粒子群算法(PSO)表示第优化算法——粒子群算法(PSO)个变量的个体极值的第优化算法——粒子群算法(PSO)维。优化算法——粒子群算法(PSO)表示全局最优解的第优化算法——粒子群算法(PSO)维。

4、终止条件

有两种终止条件可以选择，一是最大代数：优化算法——粒子群算法(PSO)；二是相邻两代之间的偏差在一个指定的范围内即停止。我们在实验中选择第一种。

三、实验

    我们选择的测试函数是：Griewank。其基本形式如下： 优化算法——粒子群算法(PSO)
 图像为： 优化算法——粒子群算法(PSO)
 （Griewank函数图像） 在实验中我们选择的维数是20；MATLAB程序代码如下： 主程序： [plain] view plain copy
  1. c1=2;%学习因子  
  2. c2=2;%学习因子  
  3. Dimension=20;  
  4. Size=30;  
  5. Tmax=500;  
  6. Velocity_max=1200;%粒子最大速度  
  7.   
  8. F_n=2;%测试函数名  
  9.   
  10. Fun_Ub=600;%函数上下界  
  11. Fun_Lb=-600;  
  12. Position=zeros(Dimension,Size);%粒子位置  
  13. Velocity=zeros(Dimension,Size);%粒子速度  
  14. Vmax(1:Dimension)=Velocity_max;%粒子速度上下界  
  15. Vmin(1:Dimension)=-Velocity_max;  
  16. Xmax(1:Dimension)=Fun_Ub;%粒子位置上下界，即函数自变量的上下界  
  17. Xmin(1:Dimension)=Fun_Lb;  
  18. [Position,Velocity]=Initial_position_velocity(Dimension,Size,Xmax,Xmin,Vmax,Vmin);  
  19.   
  20. Pbest_position=Position;%粒子的历史最优位置，初始值为粒子的起始位置，存储每个粒子的历史最优位置  
  21. Gbest_position=zeros(Dimension,1);%全局最优的那个粒子所在位置，初始值认为是第1个粒子  
  22.   
  23. for j=1:Size  
  24.     Pos=Position(:,j);%取第j列，即第j个粒子的位置  
  25.     fz(j)=Fitness_Function(Pos,F_n,Dimension);%计算第j个粒子的适应值  
  26. end  
  27. [Gbest_Fitness,I]=min(fz);%求出所有适应值中最小的那个适应值，并获得该粒子的位置  
  28. Gbest_position=Position(:,I);%取最小适应值的那个粒子的位置，即I列  
  29.   
  30. for itrtn=1:Tmax  
  31. time(itrtn)=itrtn;  
  32.   
  33. Weight=1;  
  34. r1=rand(1);  
  35. r2=rand(1);  
  36. for i=1:Size  
  37.    Velocity(:,i)=Weight*Velocity(:,i)+c1*r1*(Pbest_position(:,i)-Position(:,i))+c2*r2*(Gbest_position-Position(:,i));  
  38. end  
  39. %限制速度边界  
  40. for i=1:Size  
  41.     for row=1:Dimension  
  42.         if Velocity(row,i)>Vmax(row)  
  43.             Veloctity(row,i)=Vmax(row);  
  44.         elseif Velocity(row,i)<Vmin(row)  
  45.             Veloctity(row,i)=Vmin(row);  
  46.         else  
  47.         end  
  48.     end  
  49. end  
  50.   
  51. Position=Position+Velocity;  
  52.   
  53. %限制位置边界  
  54. for i=1:Size  
  55.     for row=1:Dimension  
  56.         if Position(row,i)>Xmax(row)  
  57.             Position(row,i)=Xmax(row);  
  58.         elseif Position(row,i)<Xmin(row)  
  59.             Position(row,i)=Xmin(row);  
  60.         else  
  61.         end  
  62.     end  
  63. end  
  64.   
  65.   for j=1:Size  
  66.      P_position=Position(:,j)';%取一个粒子的位置  
  67.      fitness_p(j)=Fitness_Function(P_position,F_n,Dimension);  
  68.      if fitness_p(j)< fz(j) %粒子的适应值比运动之前的适应值要好，更新原来的适应值  
  69.          Pbest_position(:,j)=Position(:,j);  
  70.          fz(j)=fitness_p(j);  
  71.      end  
  72.      if fitness_p(j)<Gbest_Fitness  
  73.          Gbest_Fitness=fitness_p(j);  
  74.      end  
  75.   end  
  76.   [Gbest_Fitness_new,I]=min(fz);%更新后的所有粒子的适应值，取最小的那个，并求出其编号  
  77.    Best_fitness(itrtn)=Gbest_Fitness_new; %记录每一代的最好适应值  
  78.    Gbest_position=Pbest_position(:,I);%最好适应值对应的个体所在位置  
  79. end  
  80. plot(time,Best_fitness);  
  81. xlabel('迭代的次数');ylabel('适应度值P_g');  

初始化： [plain] view plain copy
  1. function [Position,Velocity] = Initial_position_velocity(Dimension,Size,Xmax,Xmin,Vmax,Vmin)  
  2.   for i=1:Dimension  
  3.       Position(i,:)=Xmin(i)+(Xmax(i)-Xmin(i))*rand(1,Size); % 产生合理范围内的随机位置，rand(1,Size)用于产生一行Size个随机数  
  4.       Velocity(i,:)=Vmin(i)+(Vmax(i)-Vmin(i))*rand(1,Size);  
  5.   end  
  6. end  

适应值计算： [plain] view plain copy
  1. function Fitness=Fitness_Function(Pos,F_n,Dimension)  
  2.  switch F_n  
  3.     case 1  
  4.         Func_Sphere=Pos(:)'*Pos(:);  
  5.         Fitness=Func_Sphere;  
  6.     case 2  
  7.         res1=Pos(:)'*Pos(:)/4000;  
  8.         res2=1;  
  9.         for row=1:Dimension  
  10.             res2=res2*cos(Pos(row)/sqrt(row));  
  11.         end  
  12.         Func_Griewank=res1-res2+1;  
  13.         Fitness=Func_Griewank;  
  14. end  

最终的收敛曲线： 优化算法——粒子群算法(PSO)
 （收敛曲线）
标签：

相关文章