unc_n2_braninplot.m

来自「it contains many classic Test Problems f」· M 代码 · 共 23 行

% teix=0.1;
% [X,Y] = meshgrid(-5:teix:10,0:teix:15);
% Z=(Y-1.275*X.*X/pi/pi+5*X/pi-6).^2+10*(1-0.125/pi)*cos(X)+10;
% mesh(X,Y,Z)
% axis([-5 10 0 15 0 250]);
% title('unc#n2#branin');
% xlabel('x');
% ylabel('y');
% zlabel('z')


%contour
teix=0.01;
[X,Y] = meshgrid(-5:teix:10,0:teix:15);
Z=(Y-1.275*X.*X/pi/pi+5*X/pi-6).^2+10*(1-0.125/pi)*cos(X)+10;
v=[  0 1 2 5 10 50 60];
[c,h]=contour(X,Y,Z,v);
clabel(c,h);
xlabel('x');
ylabel('y');
title('unc#n2#branin');