在序列中寻找零岛

[poster]

假设您的序列很长。查找序列全为零的时间间隔（或更准确地说，序列下降到接近零的值abs(X) =的间隔达到某个指定值（例如3 ），并在所有这些间隔组合中返回值的索引：

indices = [3 4 5 6 14 15 16]; 最后一部分与上一个问题有关：

MATLAB：从开始/结束索引列表创建向量化数组

这是我到目前为止的内容：

sig = [1 1 0 0 0 0 1 1 1 1 1 0 1 0 0 0 1 1 1 1 1 1 1 1 0 0 1 1 1 0]; len = length(sig); thresh = 3; %# align the signal with itself successively shifted by one %# v will thus contain 1 in the starting locations of the zero interval v = true(1,len-thresh+1); for i=1:thresh v = v & ( sig(i:len-thresh+i) == 0 ); end %# extend the 1's till the end of the intervals for i=1:thresh-1 v(find(v)+1) = true; end %# get the final indices v = find(v); 我正在寻求矢量化/优化代码，但是我对其他解决方案持开放态度。我必须强调空间和时间效率非常重要，因为我正在处理大量的长生物信号。



回答：

从给定的vector sig开始，我将按照以下这些步骤以向量化的方式解决您的问题：

• 首先，对向量进行阈值处理以获得零和一的向量tsig （信号绝对值下降到接近零的零，其他地方为零）：
  
  tsig = (abs(sig) >= eps); %# Using eps as the threshold
• 接下来，使用函数DIFF和FIND查找每个零字符串的开始索引，结束索引和持续时间：
  
  dsig = diff([1 tsig 1]); startIndex = find(dsig < 0); endIndex = find(dsig > 0)-1; duration = endIndex-startIndex+1;
• 然后，找到持续时间大于或等于某个值（例如您的示例中为3）的零字符串：
  
  stringIndex = (duration >= 3); startIndex = startIndex(stringIndex); endIndex = endIndex(stringIndex);
• 最后，使用我对所链接问题的答案中的方法来生成最终的索引集：
  
  indices = zeros(1,max(endIndex)+1); indices(startIndex) = 1; indices(endIndex+1) = indices(endIndex+1)-1; indices = find(cumsum(indices));

更多&回答...