global M;

creatematrix(2);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画ackley图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ackley x from[-5 5]

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% temp1=x.^2+y.^2;

% temp2=cos(2*pi*x)+cos(2*pi*y);

% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的ackley图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%Rotated ackley x from[-5 5

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=x.^2+y.^2;

% temp2=cos(2*pi*x)+cos(2*pi*y);

% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画cigar图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%cigar

x=-5:0.01:5;

[x,y]=meshgrid(x);

z=x.^2+(10^4)*y.^2;

axis([-5,5,-5,5]);

meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的cigar图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rotated cigar x from [-5 5]

% x=-7:0.01:7;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=x.^2+(10^4)*y.^2;

% axis([-7,7,-7,7]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画ellipse图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%ellipse x from [-5 -5]

% x=-5:0.01:5;

% [x,y]=meshgrid(x);

% z=x.^2+10.^4*(y.^2);

% axis([-5,5,-5,5]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画griewangk图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%griewangk

% x=-600:1:600;

% [x,y]=meshgrid(x);

% temp1=(x.^2+y.^2)/4000;

% temp2=cos(x)*cos(y/sqrt(2));

% z=temp1-temp2+1;

% axis([-600,600,-600,600]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的griewangk图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated griewangk

% x=-600:1:600;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=(x.^2+y.^2)/4000;

% temp2=cos(x)*cos(y/sqrt(2));

% z=temp1-temp2+1;

% axis([-600,600,-600,600]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画nocon_rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%nocon_rastrigin

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% temp1=0;

% temp2=0;

% z=0;

% if abs(x)<1/2

% temp1=x;

% else

% temp1=round(2.*x)/2;

% end

% if abs(y)<1/2

% temp2=y;

% else

% temp2=round(2.*y)/2;

% end

% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的nocon_rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated nocon_rastrigin

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=0;

% temp2=0;

% z=0;

% if abs(x)<1/2

% temp1=x;

% else

% temp1=round(2.*x)/2;

% end

% if abs(y)<1/2

% temp2=y;

% else

% temp2=round(2.*y)/2;

% end

% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rastrigin

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% temp1=x.^2-10*cos(2*pi.*x);

% temp2=y.^2-10*cos(2*pi.*y);

% z=temp1+temp2+20;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的rastrigin图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated rastrigin

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% temp1=x.^2-10*cos(2*pi.*x);

% temp2=y.^2-10*cos(2*pi.*y);

% z=temp1+temp2+20;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画rosenbrock图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%rosenbrock

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=100*(y-x.^2).^2+(x-1).^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画schwefel图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%schwefel

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的schwefel图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated schwefel

% x=-15:0.05:15;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;

% axis([-15,15,-15,15]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画sphere图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%sphere

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=x.^2+y.^2;

% axis([-10,10,-10,10]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画sumcan图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%sumcan

% x=-0.16:0.0005:0.16;

% [x,y]=meshgrid(x);

% z=1./(10.^(-5)+abs(x)+abs(x+y));

% z=log10(z);

% axis([-0.16,0.16,-0.16,0.16]);

% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画tablet图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%tablet

% x=-10:0.05:10;

% [x,y]=meshgrid(x);

% z=10.^4*x.^2+y.^2;

% axis([-10,10,-10,10]);

% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画Weierstrass图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Weierstrass

% x=-1:0.005:1;

% [x,y]=meshgrid(x);

% z=0;

% for k=0:20

% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));

% end

% temp=0;

% for k=0:20

% temp=temp+0.5^k*cos(2*pi*3^k*0.5);

% end

% z=z-2*temp;

% axis([-1,1,-1,1]);

% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%画旋转的Weierstrass图。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Rotated Weierstrass

% x=-1.5:0.005:1.5;

% [x,y]=meshgrid(x);

% for i=1:size(x,1)

% for j=1:size(y,1)

% p=[x(i,j),y(i,j)]';

% x(i,j)=M(1,:)*p;

% y(i,j)=M(2,:)*p;

% end

% end

% z=0;

% for k=0:20

% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));

% end

% temp=0;

% for k=0:20

% temp=temp+0.5^k*cos(2*pi*3^k*0.5);

% end

% z=z-2*temp;

% axis([-1.5,1.5,-1.5,1.5]);

% meshc(x,y,z);

%
xlabel('x');

ylabel('y');

zlabel('f(x,y)');

colormap jet;

% clear;

%saveas(gcf,'temp.jpg')