3-23.m

来自「《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计」· M 代码 · 共 17 行

% 生成二元函数
v = -2:0.2:2;
[x,y] = meshgrid(v);
z = x .* exp(-x.^2 - y.^2);
% 绘制曲面，如图3-4所示
figure(1)
mesh(x,y,z);
[px,py] = gradient(z,.2,.2);
% 求偏导数
figure(2)
contour(v,v,z)
% 绘制等高线，如图3-5所示
hold on
quiver(v,v,px,py)
% 绘制矢量场图，小箭头表示梯度
hold off