一种适用于多元离散函数的粒子群算法及matlab实现

%% I. 清空环境
clc
clear

%% II. 绘制目标函数曲线
figure
[x,y] = meshgrid(1:1:100,1:1:100);%x，y 定义域
global z; %全局变量 存放z

for i=1:100
for j = 1:100
z(i,j) = x(i,j)+y(i,j);%z=x+y
end
end

mesh(x,y,z)
hold on

%% III. 参数初始化
c1 = 1.49445;
c2 = 1.49445;

maxgen = 100; % 进化次数
sizepop = 100; %种群规模
%设置最大最小速度，最大最小种群位置
Vmax = 1;
Vmin = -1;
pop1max = 1;
pop1min = -1;
pop2max = 1;
pop2min = -1;

%% IV. 产生初始粒子和速度

for i = 1:sizepop
% 产生一个种群
pop(i,:) = [50rand+51,50rand+51]; %初始种群
V(i,:) = [5rand+51,5rand+51]; %初始化速度
% 计算适应度
fitness(i) = fun(pop(i,:)); %染色体的适应度
end

%% V. 个体极值和群体极值
[bestfitness bestindex] = max(fitness);
zbest = pop(bestindex,:); %全局最佳
gbest =