700字范文 > 【MATLAB教程案例11】基于PSO粒子群优化算法的函数极值计算matlab仿真及其他应用

【MATLAB教程案例11】基于PSO粒子群优化算法的函数极值计算matlab仿真及其他应用

时间：2020-09-02 15:54:34

FPGA教程目录

MATLAB教程目录

---------------------------------------------------------------------------------------------------------------

1.软件版本

2.PSO粒子群优化的理论知识

3.完整核心代码学习

3.1通过PSO计算函数最小值

3.2通过PSO计算函数最大值

4.操作步骤与仿真结论

4.1通过PSO计算函数最小值

4.2通过PSO计算函数最大值

5.其他应用分析

1.软件版本

MATLABa或者高级版本

2.PSO粒子群优化的理论知识

PSO是粒子群优化算法（Particle Swarm Optimization）的英文缩写，是一种基于种群的随机优化技术，由Eberhart和Kennedy于1995年提出。粒子群算法模仿昆虫、兽群、鸟群和鱼群等的群集行为，这些群体按照一种合作的方式寻找食物，群体中的每个成员通过学习它自身的经验和其他成员的经验来不断改变其搜索模式。

在PSO算法中，用粒子来描述鸟类个体，每个粒子在搜索空间中搜索，飞行过程为搜索过程，飞行速度可以根据粒子历史最优位置和种群历史最优位置动态调整，粒子当前位置就是对应优化问题的一个候选解，称为个体极值，粒子群中的最优的个体极值称为全局最优解。不断地迭代，粒子更新速度和我位置，最终满足条件的全局最优解即为所求目标。PSO的基本步骤如下：

（1）初始化：迭代次数，目标函数的自变量个数，粒子的最大速度，位置信息，我们在速度区间和搜索空间上随机初始化速度和位置，粒子群规模，每个粒子随机初始化一个飞翔速度。此外还有粒子位置信息和速度信息的最小值最大值约束范围。

（2）适应度函数的定义。

适应度函数为粒子群算法的求解目标函数，比如我们要计算某个函数f(x)的最小值，那么可以定义适应度函数为f(x)，x为粒子的位置，即搜索最优的粒子位置x，使得适应度函数f(x)达到最小值。

（3）粒子更新。

3.完整核心代码学习

3.1通过PSO计算函数最小值

首先，定义一个适应度函数程序：func_fitness.m

function [minf,x,y] = func_fitness(x)minf = (x-5)^2;x=x;

这个函数表达式为：

然后定义粒子群优化算法程序：

参数初始化

%%c1 = 2;c2 = 2;%x的范围min1 = -100;max1 = 100;%速度的范围vmin =-1;vmax = 1;%粒子数目Pop = 100;BsJ = 1;%迭代次数tmax = 100;x= zeros(1,Pop);va= zeros(size(x)); x_best = zeros(1,Pop);

粒子群更新

%速度1设置va(1,i) = va(1,i) + c1*rand(1)*(x_best(1,i)-x(1,i)) + c2*rand(1)*(Tx_best-x(1,i));%更新x(1,i) = x(1,i) + va(1,i);

最优值输出

[BsJ,x(1,i)] = func_fitness(x(1,i)); if BsJ<BsJi(i)BsJi(i) = BsJ;x_best(1,i) = x(1,i);endif BsJi(i)<minJiminJi= BsJi(i);Tx_best = x(1,i);end

完整程序如下：

clc;clear;close all;warning off;rng('default');%%c1 = 2;c2 = 2;%x的范围min1 = -100;max1 = 100;%速度的范围vmin =-1;vmax = 1;%粒子数目Pop = 100;BsJ = 1;%迭代次数tmax = 100;x= zeros(1,Pop);va= zeros(size(x)); x_best = zeros(1,Pop);%粒子初始化for i=1:Popx(1,i) = rand(1)*(max1-min1)+min1;x_best(1,i) = rand(1)*(max1-min1)+min1;[BsJ,x(1,i)] = func_fitness(x(1,i)); BsJi(i) = BsJ;va(1,i) =(vmax-vmin)*rand(1)+vmin;end[minJi,index]= min(BsJi);Tx_best = 10;for t=1:tmaxttime(t) = t;for i=1:Poprng(i);%N%速度1设置va(1,i) = va(1,i) + c1*rand(1)*(x_best(1,i)-x(1,i)) + c2*rand(1)*(Tx_best-x(1,i));%更新x(1,i) = x(1,i) + va(1,i);%变量1的限制if x(1,i) >= max1x(1,i) = max1;endif x(1,i) <= min1x(1,i) = min1;end [BsJ,x(1,i)] = func_fitness(x(1,i)); if BsJ<BsJi(i)BsJi(i) = BsJ;x_best(1,i) = x(1,i);endif BsJi(i)<minJiminJi= BsJi(i);Tx_best = x(1,i);endendJibest(t) = minJi;endTx_bestfigure;plot(Jibest,'b','linewidth',1);xlabel('迭代次数');ylabel('J');grid on

3.2通过PSO计算函数最大值

首先，定义一个适应度函数程序：func_fitness.m

function [minf,x] = func_fitness(x)minf = 1/(1+abs(x));x=x;

这个函数表达式为：

完整程序如下：

clc;clear;close all;warning off;rng('default');%%c1 = 2;c2 = 2;%x的范围min1 = -1000;max1 = 1000;%速度的范围vmin =-1;vmax = 1;%粒子数目Pop = 100;BsJ = 1;%迭代次数tmax = 200;x= zeros(1,Pop);va= zeros(size(x)); x_best = zeros(1,Pop);%粒子初始化for i=1:Popx(1,i) = rand(1)*(max1-min1)+min1;x_best(1,i) = rand(1)*(max1-min1)+min1;[BsJ,x(1,i)] = func_fitness(x(1,i)); BsJi(i) = BsJ;va(1,i) =(vmax-vmin)*rand(1)+vmin;end[maxJi,index]= max(BsJi);Tx_best = 0;for t=1:tmaxttime(t) = t;for i=1:Poprng(i);%N%速度1设置va(1,i) = va(1,i) + c1*rand(1)*(x_best(1,i)-x(1,i)) + c2*rand(1)*(Tx_best-x(1,i));%更新x(1,i) = x(1,i) + va(1,i);%变量1的限制if x(1,i) >= max1x(1,i) = max1;endif x(1,i) <= min1x(1,i) = min1;end [BsJ,x(1,i)] = func_fitness(x(1,i)); if BsJ>BsJi(i)BsJi(i) = BsJ;x_best(1,i) = x(1,i);endif BsJi(i)>maxJimaxJi= BsJi(i);Tx_best = x(1,i);endendJibest(t) = maxJi;endTx_bestfigure;plot(Jibest,'b','linewidth',1);xlabel('迭代次数');ylabel('J');grid on