matlab - 从 Matlab/Octave 中的向量中识别(和删除)序列

Question

我正在尝试从 Matlab(或 Octave)中的数字向量中删除任何长度为 3 或更长的序列。例如，给定向量 dataSet，

dataSet = [1 2 3 7 9 11 13 17 18 19 20 22 24 25 26 28 30 31];

删除所有长度为 3 或更长的序列将产生 prunedDataSet:

prunedDataSet = [7 9 11 13 22 28 30 31 ];

我可以暴力破解一个解决方案，但我怀疑有一种更简洁(也许更有效)的方法可以使用向量/矩阵运算来解决这个问题，但我总是对某些东西是否产生索引或所述索引处的值感到困惑.有什么建议吗？

这是我想出的蛮力方法:

dataSet = [1 2 3 7 9 11 13 17 18 19 20 22 24 25 26 28 30 31];
benign = [];
for i = 1:size(dataSet,2)-2;
    if (dataSet(i) == (dataSet(i+1)-1) && dataSet(i) == dataSet(i+2)-2);
        benign = [benign i ] ;
    end;
end;

remove = [];
for i = 1:size(benign,2);
    remove = [remove benign(i) benign(i)+1 benign(i)+2 ];
end;

remove = unique(remove);

prunedDataSet = setdiff(dataSet, dataSet(remove));

Answer 1

这是一个使用 DIFF 的解决方案和 STRFIND

%# define dataset
dataSet = [1 2 3 7 9 11 13 17 18 19 20 22 24 25 26 28 30 31];

%# take the difference. Whatever is part of a sequence will have difference 1
dds = diff(dataSet);

%# sequences of 3 lead to two consecutive ones. Sequences of 4 are like two sequences of 3
seqIdx = findstr(dds,[1 1]);

%# remove start, start+1, start+2
dataSet(bsxfun(@plus,seqIdx,[0;1;2])) = []
dataSet =

     7     9    11    13    22    28    30    31