一、多元函数
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=