0

我想计算连续整数序列的所有长度并将它们作为向量返回。例如,考虑向量: x = [1 2 3 4 6 8 9 10 12 13];

长度将是:

length([1 2 3 4]) = 4;
length([6]) = 1;
length([8 9 10]) = 3;
length([12 13]) = 2;  

所以,我想要生成的结果是:

y = [4 1 3 2] 

我怎样才能做到这一点?

4

4 回答 4

3

这应该可以解决问题:

y = diff(find(diff([nan ; x(:) ; nan]) ~= 1))

内部diff查找不是+1 的步骤(序列中断),find确定相应的位置(索引),外部diff计算序列长度作为序列中断位置之间的差异。nans 用于确保找到向量开头的序列和结尾的序列,方法是通过诱导不同于diff1 的值。

于 2015-01-21T15:54:24.113 回答
2

A. Donda 回答的一个小变种:

  1. 用于diff检测大于 1 的差异。这会在每次运行结束1时给出一个值。
  2. 向后累积(使用cumsum)为每次运行分配不同的数字标签。累加是向后进行的,因为1第 1 步中的值是在每次运行结束时,而不是在开始时。
  3. 用 计算运行长度histc

代码:

y = [diff(x)>1 1];               %// step 1
y = cumsum(fliplr(y));           %// step 2
y = fliplr(histc(y, 1:y(end)));  %// step 3
于 2015-01-21T17:40:15.870 回答
1

这等于 的运行长度x - (1:length(x))因此你可以使用:

runLengths = @(x) diff([0, reshape(find(x(1:end-1)~=x(2:end)),1,[]), numel(x)]);
sequenceCounts = @(x) runLengths(x(:)-(1:numel(x)).');
result = sequenceCounts(x);
于 2015-01-21T15:51:42.887 回答
0
x = [1 2 3 4 6 8 9 10 12 13];

consecs = []; % empty vector to store answers
consec = 1; % initialize
for I=2:length(x)
    if x(I)==x(I-1)+1 % this one is consecutive with previous
        consec = consec+1;
    else % this one starts a new set of consecutives
        consecs(end+1) =  consec;
        consec = 1;
    end
end
consecs(end+1)=consec; % remember to include the last consecutives
display(consecs)

这是经过测试的并且有效。

于 2015-01-21T16:06:52.233 回答