基于粒子群算法的线性规划问题求解matlab程序

基于粒子群算法的线性规划问题求解matlab程序

1 基本粒子群算法流程

粒子群算法基于“种群”和“进化”的概念，通过个体间的协作与竞争，实现复杂空间最优解的搜索，其流程如下：

（1）初始化粒子群，包括群体规模 NN，每个粒子的位置 xixi 和速度 vivi。

（2）计算每个粒子的适应度值 fit[i]fit[i]。

（3）对每个粒子，用它的适应度值 fit[i]fit[i] 和个体极值 pbest(i)pbest(i) 比较。如果 fit[i]>pbest(i)fit[i]>pbest(i)，则用 fit[i]fit[i]替换掉 pbest(i)pbest(i)。

（4）对每个粒子，用它的适应度值 fit[i]fit[i] 和全局极值 gbestgbest 比较。如果 fit[i]>gbestfit[i]>gbest 则用 fit[i]fit[i]替换掉 gbestgbest。

（5）迭代更新粒子的速度 vivi 和位置 xixi。

（6）进行边界条件处理。

（7）判断算法终止条件是否满足：若是，则结束算法并输出优化结果；否则返回步骤（2）。

基本粒子群算法的算法流程图如下图所示：

2 线性规划算例

3 粒子群算法求解结果

1）迭代曲线

2）求解

4 matlab程序

1）主函数

% 粒子群算法版 
clc;
clear;
close all;

%% 算法参数

nVar=3;              %变量个数
VarMin=ones(1,3)*0; % 变量下限
VarMax=ones(1,3)*15; % 变量上限
MaxIt=60;      % 最大迭代次数
nPop=100;        % 种群规模

%% 计算
[ bestPosition, fitValue ] = PSOFUN( @fun_objective,nVar,VarMin,VarMax,MaxIt,nPop );

%% 输出
display(['粒子群算法获得的最优变量为 : ', num2str(bestPosition)]);
display([' 获得的最优目标函数: ', num2str(fitValue)]);
。。。。。。。略
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22