x <- seq(1, 5, by = 0.1)
y <- seq(1, 5, by = 0.1)
df <- expand.grid(x = x, y = y)
z <- 310 * exp(-((df$x - 3)^2 + (df$y - 4)^2)/(2*0.5^2))
df$z <- z
这里使用expand.grid函数生成一个x、y的网格,步长为0.1,并给网格中每一个点计算对应的z值。
2.使用geom_tile和geom_contour绘制z图像:
library(ggplot2)
p1 <- ggplot(df, aes(x = x, y = y, fill = z)) + geom_tile() +
scale_fill_gradient(low = "white", high = "red")
p2 <- ggplot(df, aes(x = x, y = y, z = z)) + geom_contour() +
scale_fill_gradient(low = "white", high = "red")
这里使用ggplot2包中的geom_tile和geom_contour函数分别绘制z的图像并设置了颜色刻度范围为从白色到红色。
对z在x,y平面的单调性进行描述:
根据z的表达式与实际图像,我们可以发现z在x,y平面上形成了一个高斯分布的山峰形状。当x和y的值在山峰中心(3,4)处时,z取得最大值310。当x和y的值距离山峰中心越远时,z的值会逐渐变小,即z随着x和y的增大而减小,因此z在x,y平面上是单调递减的。