基于matlab，写代码，用粒子群算法，求函数最大值，Max f(x1,x2)=21.5+x1*sin(4πx1)+x2*sin(20πx2) -2.9<=x1<=12, 4.2<=x2<=5.7 ，要显示出跌代过程图像、收敛过程图像

以下是基于 MATLAB 的粒子群算法求解函数最大值的代码：

% 定义目标函数
fun = @(x) 21.5 + x(1)*sin(4*pi*x(1)) + x(2)*sin(20*pi*x(2));

% 定义问题的边界
lb = [-2.9, 4.2];
ub = [12, 5.7];

% 定义粒子群算法参数
options = optimoptions(@particleswarm, ...
    'SwarmSize', 100, ...
    'MaxIterations', 100, ...
    'FunctionTolerance', 1e-6, ...
    'Display', 'iter');

% 运行粒子群算法
[x, fval, exitflag, output] = particleswarm(fun, 2, lb, ub, options);

% 显示结果
disp(['最大值为：', num2str(-fval)]);
disp(['取得最大值的点为：(', num2str(x(1)), ', ', num2str(x(2)), ')']);

% 显示迭代过程图像
plot(output.fval);
title('目标函数迭代过程');
xlabel('迭代次数');
ylabel('目标函数值');

% 显示收敛过程图像
plot(output.funccount, output.fval);
title('目标函数收敛过程');
xlabel('函数评价次数');
ylabel('目标函数值');

代码中使用 `particleswarm` 函数来实现粒子群算法，其中 `fun` 定义了目标函数，`lb` 和 `ub` 分别定义了问题的边界，`options` 定义了算法的参数。运行算法后，输出结果包括最大值和取得最大值的点。同时，我们也可以使用 `plot` 函数来绘制迭代过程和收敛过程的图像。

注意，粒子群算法的结果可能受到算法参数和初始粒子的影响，因此可能需要多次运行算法来获得更稳定的结果。