sympy设置隐式图范围（实体截断）

Question

我无法设置范围以使圆圈不被切断。我希望有人能帮助我。

编辑后的代码：

我想在 g 线处镜像点“V”。

from sympy import symbols
from sympy.plotting import plot
from sympy.plotting import plot_implicit
from sympy.geometry import Circle
from sympy import solve
import numpy as np   

x, y = symbols(['x', 'y'])
g = 4*x-3*y+1
V = np.array([5,2])
n = np.array([4, -3])
abs_n = np.linalg.norm(n)
hnf_g = abs(g/abs_n)
mirrored_V = V - 2*hnf_g.evalf(subs={x:V[0], y:V[1]}) * n/abs_n
y_g = solve(g,y)

p = plot(y_g[0], show=False, xlim=(-10,10), ylim=(-10,10))
p.extend(plot_implicit(Circle(V,1).equation(), var_start_end_x=(x,-10,10), var_start_end_y=(y,-10,10)))
p.extend(plot_implicit(Circle(mirrored_V,1).equation(), var_start_end_x=(x,-10,10), var_start_end_y=(y,-10,10)))
p.show()

新结果：

旧：

这是整个图的结果：

如果我只看一下隐含的情节：

Answer 1

从

plot_implicit

的文档中，可以通过设置范围来完成此操作，如第三个示例所示，复制如下。

p2 = plot_implicit(Eq(x**2 + y**2, 3), (x, -3, 3), (y, -3, 3))

但这需要一个

Eq

对象，不幸的是，它不能从

Circle

构造用于这些目的 - 可能是由于传递

bool

而不是

bool

元组的错误的结果 - 或者我我缺少一些东西。如果是前者的话

Eq(Circle(Point(0,0),1).equation())

是不够的，你必须显式地构造圆方程，即

Eq((x - 5)**2 + (y - 2)**2 - 1, 1)

---

完整示例

from sympy import symbols
from sympy.plotting import plot
from sympy.plotting import plot_implicit
from sympy.geometry import Circle
from sympy import solve
import numpy as np   

x, y = symbols(['x', 'y'])
g = 4*x-3*y+1
V = np.array([5,2])
n = np.array([4, -3])
abs_n = np.linalg.norm(n)
hnf_g = abs(g/abs_n)
mirrored_V = V - 2*hnf_g.evalf(subs={x:V[0], y:V[1]}) * n/abs_n
y_g = solve(g,y)

p = plot(y_g[0], show=False, xlim=(-10,10), ylim=(-10,10))
p.extend(plot_implicit(Eq((x - V[0])**2 + (y - V[1])**2, 1), 
                       (x,-10,10), (y,-10,10), show=False))
p.extend(plot_implicit(Eq((x - mirrored_V[0])**2 + (y - mirrored_V[1])**2, 1), 
                       (x,-10,10), (y,-10,10), show=False))
p.show()