如何用颤动图绘制我的微分方程？

Question

我想求解我的微分方程并绘制速度矢量，但是我遇到了一些麻烦。 我尝试了这个：

syms y(x);
ode = (1+exp(x))*y*diff(y,x)-2*exp(x) == 0;
ySol = dsolve(ode)

[X,Y] = meshgrid(-2:.2:2);
Z = 2*exp(X)/((1+exp(X)).*Y);
[DX,DY] = gradient(Z,.2,.2);

figure
contour(X,Y,Z)
hold on
quiver(X,Y,DX,DY)
hold off

我得到这个错误：

Warning: Matrix is singular to working precision. 
Warning: Contour not rendered for non-finite ZData 

我可能看不到这很简单，但是我只是开始使用Matlab，而我却找不到合适的方法来完成任务。 请帮我...

编辑

正如bconrad所建议的那样，我像这样更改了Z函数：

Z = 2*exp(X)/((1+exp(X)).*Y);

并且以前的错误已修复。 但是，绘制速度矢量的主要目标尚未实现，因为我得到了如下图：

Answer 1

目前尚无法检查，但我认为您需要在该行中逐个元素进行划分。 您在部门上缺少点，请尝试

Z = 2*exp(X)./((1+exp(X)).*Y);

我曾经仔细观察过我的车站。 Pablo提到的零分将inf强制为Z ，因此在缩放向量时（可理解的）颤抖变得混乱了，只是没有显示出来。 尝试以下操作（移除ode部分）：

[X,Y] = meshgrid(-2 : .2 : 2);
Z = 2 * exp(X) ./ ((1 + exp(X)) .* Y);
Z(isinf(Z)) = nan; % To avoid 0-division problems
[DX, DY] = gradient(Z, .2, .2);

figure
contour(X, Y, Z, 30, 'k')
hold on
quiver(X, Y, DX, DY, 6)
hold off

我在这里做了三件事：

1. 添加线Z(isinf(Z)) = nan; 迫使无穷大的quiver实质上忽略了它
2. 在轮廓函数中添加了参数30, 'k'以显示30条线，并使它们变黑（更明显）
3. 在quiver函数中添加了参数6 。 这将覆盖向量的自动长度缩放。

您将需要使用contour和quiver功能中的参数，以使图形随心所欲。

PS：在文件交换中有一个方便的arrow功能，我发现在创建矢量场图时可以更好地控制它。 请参阅https://www.mathworks.com/matlabcentral/fileexchange/278-arrow-评分确实合理。