使用`Plots.jl`在2D空间中绘制2d函数作为表面

Question

使用Plots.jl绘图时出现以下问题。我喜欢绘制rosenbrock函数

rosenbrock(x) = (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2

作为表面，需要输入2d Tuple{Float64,Float64}。

我能想到的是以下内容：

using Plots
gr()

rosenbrock(x) = (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2

ts = linspace(-1.0, 1.0, 100)
x = ts
y = map(rosenbrock, [(x, z) for (x,z) in zip(ts,ts)])
z = map(rosenbrock, [(x, y) for (x,y) in zip(ts,ts)])
# plot(x, x, z)
plot(x, y, z, st = [:surface, :contourf])

产生这个情节：

我认为我搞砸了一些方面，但我看不出我的错误。

我是否必须嵌套y和x的映射计算才能获得结果？

Answer 1

在我找到Rosenbrock函数的快速调查之后，如果我错了，请纠正我，但你需要指定y - 你不应该将它嵌套在z或类似的东西中。 其他人尝试了here显示的相同内容，但使用了Plots

解决方案如 Patrick Kofod Mogensen

所做

using Plots

function rosenbrock(x::Vector)
  return (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2
end

default(size=(600,600), fc=:heat)
x, y = -1.5:0.1:1.5, -1.5:0.1:1.5
z = Surface((x,y)->rosenbrock([x,y]), x, y)
surface(x,y,z, linealpha = 0.3)

这导致

旁注

我很高兴我搜索了这个，因为我一直在寻找除PyPlot之外的Julia 3D绘图仪（因为为我的程序的用户设置可能有点麻烦）这甚至看起来更好，图像可以旋转。