在Julia中绘制代数曲线

Question

我想在Julia中可视化一些代数曲线

我有多项式：

f1=(x^4+y^4-1)(x^2+y^2-2)+x^5y

f2 = x^2+2xy^2-2y^2-1/2

我想绘制V(f1)和V(f2)，这样我就可以看到它们的公共交叉点。我已经尝试在Gadfly.jl中使用轮廓图，但它似乎只允许我一次绘制一条曲线。有没有一种方法可以在Gadfly.jl中绘制两条曲线，或者在另一个Julia软件包中绘制？

这是我到目前为止所拥有的。

​

Answer 1

Gadfly使用了一个方便的组合项：Layers

https://gadflyjl.org/stable/man/compositing/#Layers

它们可以作为plot_name.layers通过绘图自由访问，并且可以手动附加(例如，使用append!(p.layers, new_layer) )。个人最喜欢的做法是在调用plot()之前构建两个层，并在plot()函数中实现任何必要的图形标签：

using Gadfly

pol_one = layer(z=(x,y) -> (x^4 + y^4 - 1) * (x^2+y^2-2) + x^5 * y,
               xmin=[-2], xmax=[2], ymin=[-2], ymax=[2],
               Geom.contour(levels=[0;]))

pol_two = layer(z=(x,y) -> x^2 + 2x*y^2 - 2y^2 - 1/2,
               xmin=[-2], xmax=[2], ymin=[-2], ymax=[2],
               Geom.contour(levels=[0;]))

plot(p_layer, q_layer, Guide.xlabel("x"), Guide.ylabel("y"))

这将产生以下图：

​