1 clear;clc;close all 2 [X,Y] = meshgrid(-2:.2:2); % 產生網格數據X和Y 3 Z = X.*exp(-X.^2 - Y.^2); % 計算網格點處曲面上的Z值 4 [DX,DY] = gradient(Z,0.2,0.2); % 計算曲面上各點處的梯度 5 figure; 6 contour(X,Y,Z) ; % 繪制等高線 7 title('等高線') 8 9 figure; 10 quiver(X,Y,DX,DY) ; % 繪制梯度場 11 title('梯度場') 12 13 figure; 14 contour(X,Y,Z) ; % 繪制等高線 15 hold on 16 quiver(X,Y,DX,DY) ; % 繪制梯度場 17 h = get(gca,'Children'); % 獲取當前axes對象的所有子對象的句柄 18 set(h, 'Color','k'); % 設置當前axes對象的所有子對象的顏色為黑色