粒子群入门（二）----------绘制二次函数，并通过粒子群求解y=-(x(1)^2+x(2)^2+x(1)+x(2)+3)；，遇到得一些问题以及解决方案

一、多元函数

1、绘制一元二次函数

1.1代码

clear;clc


%%绘制三维图
X1=-3:0.1:3;
X2=-3:0.1:3;
[X1,X2]=meshgrid(X1,X2);
y=X1.^2+X2.^2+X1+X2+3;
figure(1)
mesh(X1,X2,y)
xlabel('X1');ylabel('X2');zlabel('y');
axis vis3d
hold on;
%% 预设参数
 n=60;
 narvs=2;
 c1=2;
 c2=2;
 w=0.9;
 vmax=[-1.2 1.2];
 x_ub=[3 3];
 x_lb=[-3 -3];
 K=60;
 
 %% 初始化位置和速度
 x=zeros(n,narvs);
 for i=1:narvs
    x(:,i)=x_lb(i)+(x_ub(i)-x_lb(i))*rand(n,1); 
 end
 v=-vmax+2*vmax.*rand(n,narvs);
 
 %% 初始化适应度
 fit=zeros(n,1);
 for i=1:n
    fit(i)=fun3(x(i,:)); 
 end
 pbest=x;
 ind=find(fit==max(fit),1);
 gbest=x(ind,:);
 

 %% 画图
% h=scatter3(x(:,1),x(:,2),fit,'r');
 
%% 迭代
 fitnessbest = ones(K,1);
 for d=1:K
     for i=1:n
         v(i,:)=w*v(i,:) + c1*rand(1)*(pbest(i,:)-x(i,:)) + c2*rand(1)*(gbest-x(i,:));
         for j=