python - 绘制非多项式多元隐式方程

Question

我目前正在尝试用 matplotlib 绘制以下方程：

-4xy + x² - y² + 4(x³y - xy³) = 0

我不想解决这个可怕的方程，以便通过简单地分配x = np.linspace(-5,5)然后将其放入我的方程来绘制它。

我想知道是否可以直接绘制解决方案matplotlib.pyplot？我能够使用 Geogebra 快速快速地绘制函数，所以我想这也应该可以使用 python。但是，我还没有发现任何其他问题。有没有人有任何想法？

标签： pythonmatplotlibmathplotgraph

Answer 1

一种方法是使用 Sympy 的plot_implicit函数。

代码：

from sympy import plot_implicit, symbols, Eq

x, y = symbols('x y')

plot_implicit(Eq(-4*x*y+x**2-y**2+4*(x**3-x*y**3), 0),
              adaptive=False,
              points=1000)

输出：

在引擎盖下，这使用网格来决定是否在图上的每个点上绘制函数。这也是 Geogebra 将在幕后做的事情。要在 vanilla matplotlib 中应用这种方法，我们可以使用等高线图：

代码：

import numpy as np
import matplotlib.pyplot as plt

# Plot axes in middle
fig, ax = plt.subplots()
ax.spines['left'].set_position('center')
ax.spines['bottom'].set_position('center')
ax.spines['right'].set_color('none')
ax.spines['top'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.yaxis.set_ticks_position('left')

# Set up mesh grid
x = y = np.linspace(-1.0, 1.0, 100)
X, Y = np.meshgrid(x,y)

# Plot contour
ax.contour(X, Y, -4*X*Y+X**2-Y**2+4*(X**3-x*Y**3), [0])
plt.show()

输出：