Question

如何用Matlab画出z=tanx/(tanx-tany)的向量场？

我尝试了如下代码：
[X,Y] = meshgrid(-80:2:80,-80:2:80);
Z =tand(X)./(tand(X)-tan(Y));
[DX,DY] = gradient(Z,2,2); %梯度
% contour(X,Y,Z) %等高线
hold on
quiver(X,Y,DX,DY) %向量场
xlabel('x');
ylabel('y');
colormap hsv
hold off
但结果不对，正确的结果应该如图，求帮忙。

Answer 1

主要问题是某些点的梯度是发散的，结果算出来DX,DY里面有一堆Inf。画图的时候，Inf的箭头很长，所以把其它地方的箭头都变成一个点了。解决方法是把Inf都替换成NaN，代码如下(缩进的地方是修改过的代码)：

[X,Y] = meshgrid(-80:2:80,-80:2:80);
Z =tand(X)./(tand(X)-tand(Y));
[DX,DY] = gradient(Z,2,2); %梯度
　　DX(abs(DX)==Inf)=NaN;
　　DY(abs(DY)==Inf)=NaN;
% contour(X,Y,Z) %等高线
hold on
　　quiver(X,Y,DX,DY,10) %向量场
xlabel('x');
ylabel('y');
colormap hsv
hold off