Производные данные о рассеянии в Matlab

Question

У меня есть данные, собранные из сети GPS. Мои данные состоят из двух массивов станций (lat, long) и еще двух массивов, заполненных вертикальной и горизонтальной скоростью каждой станции. Мой сценарий для интерполирования является:

clear all; clc; format compact 
load('lat_long_Ve_Vn.mat'); 
x = 34.5:0.1:42; 
y = 19:0.1:28.5; 
[Xq,Yq] = meshgrid(x,y); 
Ve_i = griddata(lat,long,Ve,Xq,Yq); 
Vn_i = griddata(lat,long,Vn,Xq,Yq); 

Я получаю интерполированные данные для каждого узла на моей сетке с двумя векторами, Ve_i и Ve_n Я хочу, чтобы вычислить следующие производные, но я понятия не имею о том, как это сделать ,

Следует отметить, что Vx это мой Ve_i и Vy это мой Vn_i, и я не имею математическую формулу, чтобы я мог рассчитать производные MuPAD. Любая идея о том, как это сделать?

Answer 1

Если вы используете Градиент:

[Vxx Vxy] = gradient(Vx); 

[Vyx Vyy] = gradient(Vy);