Attempting to find data points of two different matrices based off of time, but outcome is that I am getting the same value repeatedly

Question

Krispy Scripts 2016-10-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/309329-attempting-to-find-data-points-of-two-different-matrices-based-off-of-time-but-outcome-is-that-i-am

评论： dpb 2016-11-1

I have one matrix with data in the first column and time stamps in the second column (datamatrix.mat). The next matrix contains spiketimes (spiketimematrix.mat). I want to get the data point in the first column of first matrix that is the closest time point corresponding to the spike times in spiketimematrix.mat. For example, the first spiketime is 166.1670, which corresponds to the closet time point of 166.1696 and corresponds with the data point 2.5281. I have tried using the following code:

for k=1:length(spiketimesmatrix)
  dataout(k)=datamatrix(find(spiketimesmatrix(k)>=datamatrix(:,2),1));
end

However, I am getting the first value of =1.8577 repeatedly and I am not sure why?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2016-10-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/309329-attempting-to-find-data-points-of-two-different-matrices-based-off-of-time-but-outcome-is-that-i-am#answer_240877

在 MATLAB Online 中打开

'Cuz you're not restricting the search to only begin with the last found location and by default find will find the first location that satisfies the condition. If the data are ordered in ascending order, then once you have a location you're looking for that is greater, then any other that is that larger or greater will also satisfy...

Two ways to fix--

Use the 'last' optional argument and search for '<=' instead of '>=' which will give you the last value that isn't over, or
Use the search for a minimum difference instead, which actually is what your topic says you're looking for, anyway. Note the search for the first/last over/under a given value isn't necessarily the closest, only the one before the sign change of the difference.

for k=1:length(spiketimesmatrix)
  dataout(k)=datamatrix(find(min(abs(spiketimesmatrix(k)-datamatrix(:,2))),1));
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

dpb 2016-10-26

编辑：dpb 2016-10-26

在 MATLAB Online 中打开

No, but got to make sure you don't try to index using the longer into the shorter...here's a working sample written without find...

>> spike=linspace(0,100,13).';    % a sample list of times
>> data(:,2)=linspace(1,155,6).'; % different length data array
>> data(:,1)=1:length(data);      % just an id to look at
>> spike
spike =
       0
3333
6667
0000
3333
6667
0000
3333
6667
0000
3333
6667
0000
>> data
data =
0000    1.0000
8000   31.8000
6000   62.6000
4000   93.4000
2000  124.2000
0000  155.0000
>> for i=1:length(spike)
     [~,ix]=min(abs(spike(i)-data(:,2)));  % find nearest location
     disp([spike(i) ix data(ix,:)])  % show results
  end
   1     1     1
3333    1.0000    1.0000    1.0000
6667    2.0000    2.0000   31.8000
0000    2.0000    2.0000   31.8000
3333    2.0000    2.0000   31.8000
6667    2.0000    2.0000   31.8000
0000    3.0000    3.0000   62.6000
3333    3.0000    3.0000   62.6000
6667    3.0000    3.0000   62.6000
0000    3.0000    3.0000   62.6000
3333    4.0000    4.0000   93.4000
6667    4.0000    4.0000   93.4000
0000    4.0000    4.0000   93.4000
>>

Hopefully that mimics your situation well enough you can adapt it...note I made the granularity purposefully large to illustrate the min() selection of which point chosen...

Krispy Scripts 2016-11-1