MATLAB绘制等高线和梯度场

clear;clc;close all
[X,Y] = meshgrid(-2:.2:2); % 产生网格数据X和Y
Z = X.*exp(-X.^2 - Y.^2); % 计算网格点处曲面上的Z值
[DX,DY] = gradient(Z,0.2,0.2); % 计算曲面上各点处的梯度
figure;
contour(X,Y,Z) ; % 绘制等高线
title('等高线')
 
figure;
quiver(X,Y,DX,DY) ; % 绘制梯度场
title('梯度场')
 
figure;
contour(X,Y,Z) ; % 绘制等高线
hold on
quiver(X,Y,DX,DY) ; % 绘制梯度场
h = get(gca,'Children'); % 获取当前axes对象的所有子对象的句柄
set(h, 'Color','k'); % 设置当前axes对象的所有子对象的颜色为黑色