我从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计算衍生物。关于如何做的任何想法?
答案 0 :(得分:0)
如果您使用Gradient:
[Vxx Vxy] = gradient(Vx);
[Vyx Vyy] = gradient(Vy);