在 R 中編碼數值導數

Question

我對 R 完全陌生，我正在努力編寫代碼來查找向量場的數值導數。 我有兩個矩陣 U 和 V，例如，

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

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

這些矩陣（顯然不是實際值）表示 2D 風矢量場的分量。 我想對 2 個向量分量du/dy和dv/dx的數值導數進行編碼。 我不知道如何在 R 中執行此操作。 請幫忙。 如果已經回答了這個問題，請提前道歉。

Answer 1

您正在尋找的是diff() function。 您可以使用 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

希望這有幫助。