Matlab中向量场的绘制

% quiver(x,y,u,v)
% x，y是包含坐标位置的矩阵，而u和v则是包含偏导数的矩阵。
% 例如绘制f(x,y)=y-3x-2x^2-3xy-3y^2的方法：
% 先用gradient函数计算偏导数，然后再用quiver绘制向量场


f=@(x,y) -y-3*x-2*x.^2-3*x.*y-y.^2;
[x,y]=meshgrid(-20:0.5:20,-10:0.5:10);
z=f(x,y);
[fx,fy]=gradient(z,0.5);
figure
c=contour(x,y,z);
clabel(c);
hold on;
quiver(x,y,fx,fy);
hold off;

figure,mesh(x,y,z)
hold on;
quiver(x,y,fx,fy);
hold off;