隐式绘图的同集范围(实体截断)

2024-04-25 12:00:54 发布

您现在位置:Python中文网/ 问答频道 /正文

我没能设置好范围,这样圆就不会被切断。我希望有人能帮助我。你知道吗

编辑的代码:

我想在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()

新结果:

New result

旧的:

这是整个绘图的结果: whole plot

如果我只看一下隐含的情节: only implicit plot with the circle


Tags: fromimportplotvarnpabsplottingarray
1条回答
网友
1楼 · 发布于 2024-04-25 12:00:54

^{}的文档中,可以通过如下复制的第三个示例中所示的设置范围来实现这一点。你知道吗

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()

相关问题 更多 >

    热门问题