必威体育Betway必威体育官网
当前位置:首页 > IT技术

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

时间:2019-08-14 14:41:03来源:IT技术作者:seo实验室小编阅读:68次「手机版」
 

PSO

一、粒子群算法的概述

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

二、粒子群算法的流程

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

(PSO流程)

下面我们具体解释下流程图里面的每一个步骤:

1、初始化

 首先,我们需要设置最大的速度区间,防止超出最大的区间。位置信息即为整个搜索空间,我们在速度区间和搜索空间上随机初始化速度和位置。设置群体规模

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

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

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

 更新公式为:

其中,称为惯性因子,称为加速常数,一般取表示区间上的随机数。表示第个变量的个体极值的第维。表示全局最优解的第维。

4、终止条件

有两种终止条件可以选择,一是最大代数:;二是相邻两代之间的偏差在一个指定的范围内即停止。我们在实验中选择第一种。

三、实验

我们选择的测试函数是:Griewank。其基本形式如下:

图像为:

(Griewank函数图像)

在实验中我们选择的维数是20;Matlab程序代码如下:

主程序:


 
  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');

  82.  

初始化:


 
  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

适应值计算:


 
  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)

1. 概念   粒子群优化算法(PSO:Particle swarm optimization) 是一种进化计算技术(evolutionary computation)。 源于对鸟群捕食

爱普生epson lq-630k打印机没有反应,发送成功,就是不

epson lq-630k打印机系统:windows7旗舰版,刚买的台式机。故障现象:驱动安装成功,测试没有文档挂起,打印机没有反应,XP系统可以打印。故

Linux下安装EPSON L310打印机驱动

Linux下安装EPSON L310打印机驱动 在这里http://download.ebz.epson.net/dsc/search/01/search/searchModule搜索l310查到如下结

rhapsody

实现了一些基本的容器,如list vector rbtree Ttree deque hashtable heap avltree 一些常用的算法quick sort ,heap sort,  bina

爱普生Epson WFM5693墨仓式打印机怎么安装驱动?

爱普生Epson WFM5693墨仓式打印机驱动想要安装驱动,该怎么安装呢?下面我们就来看看详细的教程。一、扫描功能驱动安装1、首先到下载

分享到:

栏目导航

推荐阅读

热门阅读