draw_fun.m

采用标准群粒子优化计算方法

flag=8;
switch flag
    case 1
        x=-10:0.2:10;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=X.^2+Y.^2;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Spherical函数');
    case 2
        %x=-10:0.2:10;
        x=-1:0.02:1;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=100*(X.^2-Y).^2+(1-X).^2;
        surf(X,Y,Z);
        colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Rosenbrock函数');
    case 3
        %x=-10:0.2:10;
        x=-3.5:0.05:3.5
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=(sin(sqrt(X.^2+Y.^2)).^2-0.5)./((0.001*(X.^2+Y.^2)+1).^2)-0.5;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Shaffer’s f6函数');
    case 4
        %x=-5:0.1:5;
        x=-1:0.02:1;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=X.^2-10*cos(2*pi*X)+10+Y.^2-10*cos(2*pi*Y)+10;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Rastrigrin函数');
    case 5
        %x=-1000:20:1000;
        x=-10:0.2:10;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=(1/4000)*(X.^2+Y.^2)-cos(X./sqrt(1)).*cos(Y./sqrt(2))+1;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Griewank函数');
    case 6
       %x=-50:1:50;
        x=-5:0.1:5;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=-20*exp(-0.2*sqrt((X.^2+Y.^2)./2))-exp((cos(2*pi*X)+cos(2*pi*Y))./2)+20+exp(1);
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Ackley函数');
    case 7
        x=-10:0.2:10;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=X.^2+2*Y.^2
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Axis Parallel hyper-ellipsoid函数');
    case 8
        %x=-500:10:500;
        x=-600:4:-300;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=X.*sin(sqrt(abs(X)))+Y.*sin(sqrt(abs(Y)))+418.9829*2;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Schwefel函数');
    otherwise
        %x=-1000:20:1000;
        x=-600:4:-300;
        y=x;
        [X,Y]=meshgrid(x,y);
        Z=X.*sin(sqrt(abs(X)))+Y.*sin(sqrt(abs(Y)))+418.9829*2;
        surf(X,Y,Z);
       % colormap();
        xlabel('x1');
        ylabel('x2');
        zlabel('f(x1,x2)');
        %title('Schwefel函数');
end