绘制 3D 椭圆抛物面的问题

Question

我希望使用 R 绘制 3D 函数。例如，采用 f(x,y) = (????−2????−1)^2 + (3????+????−2)^2 给出的椭圆抛物面。这是我尝试过的：

require(lattice)
x <- seq(-10, 10, by=0.5)
y <- seq(-10, 10, by=0.5)
g <- expand.grid(x = x, y = y)
g$z <- (x-2*y-1)^2 + (3*x-y-2)^2
wireframe(z ~ x * y, g, drape = TRUE,
          aspect = c(1,1), colorkey = TRUE)`

这是输出

但是，这是 f 的“真实”图：

我尝试更改x 和y 的定义，但无济于事。我还尝试了 emdbook 包中的 curve3d() 函数。看起来更糟。

Answer 1

你乘错了x和y。你需要使用g里面的那些：

g$z <- with(g, (x-2*y-1)^2 + (3*x-y-2)^2)
wireframe(z ~ x * y, g, drape = TRUE,
          aspect = c(1,1), colorkey = TRUE)