1

我在这里的 .mat 文件中有两个变量:https: //www.yousendit.com/download/UW13UGhVQXA4NVVQWWNUQw

testz 是累积距离的向量(以米为单位,单调且有规律地增加)

testSDT 是使用距离向量和速度向量生成的积分(累积)声波传播时间(以毫秒为单位)的向量(有一个创建间隔传播时间的中间步骤)

由于速度是一个连续可变的函数,因此得到的间隔行程时间以及积分行程时间都是非整数且大小可变

我想要的是以固定的时间间隔(例如 1 ms、2 ms、...、n ms)重新采样距离向量

困难在于最大行程时间 994.6659 小于 2 个向量中的样本数,因此使用 interp1 并不简单。IE:

X=testSDT -> 1680 个样本

Y=testz -> 1680 个样本

XI=[1:1:994] -> 994 个样本

这是我想出的代码。这是一个有效的代码,我认为这还不错。

%% Initial chores
M=fix(max(testSDT));
L=(1:1:M);     

%% Create indices
% this loops finds the samples in the integrated travel time vector
% that are closest to integer milliseconds and their sample number
for i=1:M
  [cl(i) ind(i)] = min(abs(testSDT-L(i)));
  nearest(i) = testSDT(ind(i)); 
end

%% Remove duplicates
% this is necessary to remove duplicates in the index vector (happens in this test). 
% For example: 2.5 ms would be the closest to both 2 ms and 2 ms
[clsst,ia,ic] = unique(nearest);
idx=(ind(ia));

%% Interpolation
% this uses the index vectors to resample the depth vectors at
% integer times
newz=interp1(clsst,testz(idx),[1:1:length(idx)],'cubic')';

据我所知,这段代码存在一个问题:我依靠向量 idx 作为我的 XI 进行插值。向量 idx 比向量 ind 短 1 个样本(删除了一个副本)。

因此,我的新时间将停止一毫秒。这是一个非常小的问题,不太可能出现重复,但我想知道是否有人可以想到解决方法,或者完全解决问题的不同方法。

谢谢

4

1 回答 1

1

如果我对你的理解正确,你想推断出那个额外的点。您可以通过多种方式做到这一点,一种是将额外的点添加到 interp1 行。如果您有一些希望跟踪数据的功能,您可以通过将其拟合到数据然后获得该额外点或使用类似的工具来使用它fnxtr

但是由于您使用该线路的方式,我无法理解您想要什么。您使用的第三个参数 [1:1:length(idx)],只是系列 [1 2 3 ...],通常在插值时,使用一些x_i兴趣点向量,尽管我怀疑您的兴趣点恰好是整数系列1:length(idx),您想要的只是[1:length(idx) xi]xi那个额外的点 x 轴值在哪里。

编辑:

L而不是循环只是从and产生矩阵形式testSDT,然后矩阵运算在执行以下操作时会更快一些min(abs(...

  MM=ones(numel(testSDT),1)*L;
  TT=testSDT*ones(1,numel(L));
  [cl ind]=(min(abs(TT-MM)));
  nearest=testSDT(ind);
于 2013-01-01T21:22:13.017 回答