问R中的数值导数编码
EN

Stack Overflow用户

提问于 2021-11-17 15:17:12

回答 1查看 43关注 0票数 0

我完全是R的新手，我正在努力编写一个代码来寻找向量场的数值导数。我有两个矩阵U和V，例如，

U <- matrix(runif(9), nrow = 3, ncol = 3, byrow = T)

V <- matrix(runif(9), nrow = 3, ncol = 3, byrow = T)

这些矩阵(显然不是实际值)表示二维风矢量场的分量。我想编码两个向量分量du/dy和dv/dx的数字导数。我不知道如何在R中做到这一点。请帮助。如果这个问题已经回答了，请提前道歉。

回答已采纳

发布于 2021-11-17 16:53:08

您要查找的是diff()函数。可以使用apply将其有效地应用于矩阵的某个维度

U <- matrix(runif(9), nrow = 3, ncol = 3, byrow = T) #Your wind component

apply(U,2,diff) #change the '1' by '2' to apply diff over the other dim

希望这能帮上忙。

票数 0

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/70007008

复制

相似问题

问R中的数值导数编码EN