使用 for 循环在 MATLAB 中对错误事件进行排序？

Question

我在 MATLAB 中的一行脚本存在一些问题，该脚本应该对错误检测到的事件进行分类。所以基本上，我有这两个变量spikemax 和spikethresh，它们包含我的数据跟踪的索引。 spikemax 在我的数据跟踪中包含正确的时间点索引，但 spikethresh 仍然有太多值，我需要对其中的一些进行排序。 规则应该是，spikethresh 的每个值总是在下一个spikethresh 之前跟一个spikemax。我添加了一个图像，其中显示了红色和绿色的点，我的两个变量最终在数据跟踪中指示。

我的想法是运行一个 for 循环，其条件是从 spikethresh 中删除值以更新变量。在代码中，我添加了一些示例值，这些值简单地显示了什么不起作用。 这里的输出是[15029;15055;15066;15339]，但我想得到

spikethresh = [15029;15066;15339]

如何修复我的代码以产生所需的输出？

这是我的代码：

% Example variables:
spikethresh = [15029; 15038;15055; 15060; 15066; 15339];
spikemax= [15034; 15072; 15344];

spikethresh_corr = spikethresh ;

for i = 1:length(spikemax)-1
    if (spikethresh(i+1)-spikethresh(i)) < (spikemax(i)-spikethresh(i))
        spikethresh_corr(i) = [] ; 
    end
    spikethresh = spikethresh_corr ;
end 


% Results in spikethresh = [15029;15055;15066],
% but I am trying to get [15029;15066;15339]

Answer 1

您的代码中有两个问题：

1. (spikethresh(i+1)-spikethresh(i)) < (spikemax(i)-spikethresh(i)) 并不能真正测试您的状况。你要测试的是spikethresh(i+1) < spikemax(i)。也就是说，如果spikethresh 中的下一个值小于spikemax 中的当前值，则我们不需要spikethresh 中的当前值。

2. 您的测试正好是 2 次（length(spikemax)-1 是 2）。但是您可能需要删除更多元素。对于每个i，您必须重复应用您的测试，直到它不再满足（在当前spikemax(i) 之前没有更多值要删除）。

这就是您的代码的样子：

spikethresh = [15029; 15038; 15055; 15060; 15066; 15339];
spikemax = [15034; 15072; 15344];

for i = 1:length(spikemax)-1
   while spikethresh(i+1) < spikemax(i)
      spikethresh(i) = [];
   end
end

这导致spikethresh = [15029; 15066; 15339]。