用R画等高线图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> 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
如何理解f和fa呢，这里用outer函数输出的f是什么啊？
希望大神能指点指点

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏5 回帖

关键词：等高线 contour Levels Level Tour

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后在电脑上执行：

1. panel.3d.contour <- function(x, y, z, rot.mat, distance, nlevels = 20, zlim.scaled, ...) {
   
2.     add.line <- trellis.par.get("add.line")
   
3.     panel.3dwire(x, y, z, rot.mat, distance, zlim.scaled = zlim.scaled, ...)
   
4.     clines <- contourLines(x, y, matrix(z, nrow = length(x), byrow = TRUE), nlevels = nlevels)
   
5.     for (ll in clines) {
   
6.         m <- ltransform3dto3d(rbind(ll$x, ll$y, zlim.scaled[2]), rot.mat, distance)
   
7.         panel.lines(m[1,], m[2,], col = add.line$col, lty = add.line$lty, lwd = add.line$lwd)
   
8.     }
   
9. }

复制代码

library(lattice)
wireframe(f, xlab = "x", ylab = "y", panel.3d.wireframe = "panel.3d.contour", drape = TRUE)

又学了一招

厉害啊 论坛高手真多 会R的也好多啊

cheetahfly 发表于 2018-9-10 13:22
或者可以将下面的代码copy后在电脑上执行：
library(lattice)
wireframe(f, xlab = "x", ylab = "y", pan ...

哈哈哈哈battle起来了

很好的学术范围啊