%繪制向量場圖
%例一
clear all;clc;
[X,Y] = meshgrid(-2:.2:2,-3:.2:3);
Z = X.*exp(-X.^2 - Y.^2);
[DX,DY] = gradient(Z); %Dx為水平方向上的梯度,第一列元素為原矩陣第二列與第一列元素之差,
%第二列元素為原矩陣第三列與第一列元素之差除以2,以此類推
[DDX,DDY]=gradient(Z,0.1,0.25);%這個就是將默認除以的2,x軸換為0.1,y軸換為0.25.
%[FX,FY]=gradient(F,HX,HY) HX,HY參數表示各方向相鄰兩點的距離
contour(X,Y,Z,7); %7個等級的等值線圖
hold on,colormap hsv
quiver(X,Y,DX,DY,0.8); %繪制梯度(向量)場
%1.2代表所畫出的向量長度的縮放系數
hold off

%例二
n=-2:.15:2;
[X,Y,Z]=peaks(n);
contour(X,Y,Z,10);colormap autumn
[U,V]=gradient(Z);hold on
quiver(X,Y,U,V);hold off
legend('等值線','向量場');

%例三:三維向量場圖繪制
%本例是繪制三維參數曲線的很好的例子
vz = 10;
% Velocity
a = -32;
% Acceleration
t = 0:.1:pi/2;
x = 2*cos(t);y = 3*sin(t);
z = vz*t + 1/2*a*t.^2;
plot3(x,y,z,
'r'
);hold
on
%繪制曲線圖像
u = gradient(x);v = gradient(y);w = gradient(z);
quiver3(x,y,z,u,v,w,0.2,
'b'
);
%繪制切向量
