exn527a.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 14 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-2-7a
% 用梯度图形分析空间曲面的极值点
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月

clear
[x,y]=meshgrid([-5:0.5:3],[-3:0.5:5]);
z=x.^3-y.^3+3*x.^2+3*y.^2-9*x;
contour(x,y,z,20), colorbar         % 等高线绘制
hold on,grid on
plot([-3,1,1,-3],[0,0,2,2],'o')     % 画出极点
[px,py] = gradient(z,.2,.2);		% 以0.2步长求z1的梯度的x,y分量
quiver(x,y,px,py)                   % 画出梯度向量图
set(gcf,'color','w')