contour是把第三维的f值画到x-y平面，形成等高线图，在x-y平面上看就是相同的f值会在一个圈内，形成像年轮一样的结构。这里是指把x,y值代人cos(y)/(1+x^2))公式得到对应的f矩阵值，相当于外积。可理解为x,y在第三维f对应的所有值。fa 就好理解了，t是translocation矩阵转置，行和列兑换，兑换后维度相同采能做矩阵加减法。你要补点线性代数，知道矩阵内积和外积是什么东东！

matlab 代码

X = -pi:0.1:pi;
Y = X;
[X,Y] = meshgrid(X,Y);
f = cos(Y)./(1+X.^2);
meshc(X,Y,f);
fa = (Z-Z\')/2; meshc(X,Y,fa);

试试R代码 x <- seq(-pi, pi, length = 50)
y <- x
f <- outer(x, y, function(x, y)cos(y)/(1+x^2))
contour(x, y, f)
contour(x, y, f, nlevels = 45)
fa <- (f-t(f))/2
p1 <- plot_ly(x=x, y = y, z = f) %>% add_surface()
p2 <-plot_ly(x=x, y = y, z = fa) %>% add_surface()
还是Matlab简单

用plotly画3D图，google  了解一下
install.packages(plotly)
library(plotly)

“还是Matlab简单？”
不存在的

y <- x <- seq(-pi, pi, length = 50)
f <- outer(x, y, function(x, y) cos(y) / (1 + x ^ 2))
lattice::wireframe(f, xlab = "x", ylab = "y", drape = TRUE)

或者可以将下面的代码copy后在电脑上执行：
library(lattice)
wireframe(f, xlab = "x", ylab = "y", panel.3d.wireframe = "panel.3d.contour", drape = TRUE)