如何在Julia上绘制矢量场?

Gus*_*uri 4 plot julia

我是Julia的新手,我尝试制作以下文章的代码如何在Julia中绘制矢量场?,但是没有用,所以,我想知道是否可以使用“ Plots”软件包进行绘图以及如何进行?这对我的研究非常重要。

附:有人给我下面的代码,但是,实际上,我不知道为什么不起作用:

using Plots
gr(size=(600,400))

function example()
  X = linspace(-2, 2, 100)
  Y = linspace(-2, 2, 100)
  f(x, y) = x^3 - 3x + y^2
  contour(X, Y, f)

  x = linspace(-2, 2, 11)
  y = linspace(-2, 2, 11)
  df(x, y) = [3x^2 - 3; 2y] / 25
  quiver!(x, y', quiver=df, c=:blue)

  xlims!(-2, 2)
  ylims!(-2, 2)
  png("example")
end

example()
Run Code Online (Sandbox Code Playgroud)

crs*_*nbr 6

正如评论中已经提到的那样,您应该提供错误消息,否则人们将不得不猜测您的代码出了什么问题。

但是,在您的情况下,我想我可以猜到:)

在Julia 1.0上,以下作品:

using Plots
gr(size=(600,400))

function example()
  X = range(-2, stop=2, length=100)
  Y = range(-2, stop=2, length=100)
  f(x, y) = x^3 - 3x + y^2
  contour(X, Y, f)

  x = range(-2, stop=2, length=11)
  y = range(-2, stop=2, length=11)
  df(x, y) = [3x^2 - 3; 2y] / 25
  quiver!(x, y', quiver=df, c=:blue)

  xlims!(-2, 2)
  ylims!(-2, 2)
  png("example")
end

example()
Run Code Online (Sandbox Code Playgroud)

并给出以下输出

在此处输入图片说明

请注意,我只是改变了所有出现的linspacerange(-2, stop=2, length=X),因为linspace功能已被否决的朱莉娅0.7。