#r
Вопрос:
Я совершенно новичок в R и изо всех сил пытаюсь написать код для поиска числовых производных векторных полей. У меня есть две матрицы U и V, например,
U lt;- matrix(runif(9), nrow = 3, ncol = 3, byrow = T)
V lt;- matrix(runif(9), nrow = 3, ncol = 3, byrow = T)
Эти матрицы (очевидно, не фактические значения) представляют компоненты векторного поля 2D ветра. Я хотел бы закодировать числовые производные 2-х векторных компонент du/dy
и dv/dx
. Я понятия не имею, как это сделать в R. Пожалуйста, помогите. Заранее извините, если на этот вопрос уже был дан ответ.
Ответ №1:
То, что вы ищете, — это diff()
функция. Вы можете эффективно применить его к измерению матрицы с помощью приложения
U lt;- 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
Надеюсь, это помогло.