我有一些时间/频率数据,我尝试使用Matlab的interp2
函数进行插值。数据[F,T,data]
是从Matlab的不同Matlab例程中获得的,如果您感兴趣,可以使用频谱图。
[~,F,T,data] = spectrogram(...)
data = 10*log10(data);
我可以使用surf
绘制数据。我相信数据很好。然而,插入数据似乎是一个问题。即使使用interp2(F,T,data,F,T)
(实际上没有插值)也会出现以下错误。
这里出了什么问题?
我有我在这里使用的数据:https://www.dropbox.com/s/zr5zpfhp6qyarzw/test.mat
interp2(F,T,data,f,t)
Error using griddedInterpolant
The grid vectors do not define a grid of points that match the given values.
Error in interp2>makegriddedinterp (line 228)
F = griddedInterpolant(varargin{:});
Error in interp2 (line 128)
F = makegriddedinterp({X, Y}, V, method,extrap);
>> size(F),size(T),size(data),size(f),size(t)
ans =
129 1
ans =
1 52
ans =
129 52
ans =
200 1
ans =
1 121
答案 0 :(得分:1)
问题是你应该交换F
和T
:
interp2(T,F,data,t,f);
第一个参数对应于矩阵的列,第二个参数对应的行记录为here:
如果X和Y是网格向量,则V必须是包含的矩阵 长度(Y)行和长度(X)列。
作为替代,您可以采用data
的转置:
interp2(F,T,data',f,t);
推理(奇怪的)争论顺序
interp2(X,Y,V,Xq,Yq)
被解释为函数的插值,由矩阵 V表示,即样本值。问题是函数/矩阵的参数/索引是以相反的顺序提供的:
matrix(row, column)
与
function(x,y)
x
(第一个参数)通常表示水平轴,因此与column
({1}}和y
的{{1}}(第二个参数)参数和同义词相对应。