测试函数shubert(十)自适应变异粒子群求解方法

测试函数shubert（十）

shubert函数属于周期性多峰函数，如图1所示拥有多个全局最优值，如图二所示在一个周期内只有一个全局最优值，局部最优解较多，适合测试算法的收敛性能，粒子群算法是一种收敛速度较快的算法，运算速度快，但是粒子群算法容易陷入局部最优，有些时候会导致收敛慢，或者不收敛，本文用标准粒子群算法进行改进，加入自适应变异，用自适应变异粒子群AMPSO求解，如有疑问，欢迎大家留言交流！
函数图像如下


MATLAB编程shubert代码如下：

clc
clear
close all
x = -2:0.1:2;
y  = -2:0.1:2;
x = -10:0.1:10;
y  = -10:0.1:10;
[x,y] = meshgrid(x,y);
[m,n] = size(x);
z = zeros(m,n);
for ii = 1:m
    for jj = 1:n
        xx = [x(ii,jj) y(ii,jj)];
        z(ii,jj) = shubertfun(xx);
    end
end
figure
surf(x,y,z)
xlabel('x1')
ylabe