xianxingguihua.rar_优化 工程问题

[x,y]=meshgrid(-2:0.1:2,-1:0.1:3); z=100*(y-x.^2).^2+(1-x).^2; figure(1); mesh(x,y,z) figure(2); contour(x,y,z,20) drawnow hold on plot(-1.2,2,'o'); text(-1.2,2,'start point') plot(1,1,'o') text(1,1,'solution')