我想计算连续整数序列的所有长度并将它们作为向量返回。例如,考虑向量:
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]
我怎样才能做到这一点?
我想计算连续整数序列的所有长度并将它们作为向量返回。例如,考虑向量:
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]
我怎样才能做到这一点?
这应该可以解决问题:
y = diff(find(diff([nan ; x(:) ; nan]) ~= 1))
内部diff
查找不是+1 的步骤(序列中断),find
确定相应的位置(索引),外部diff
计算序列长度作为序列中断位置之间的差异。nan
s 用于确保找到向量开头的序列和结尾的序列,方法是通过诱导不同于diff
1 的值。
A. Donda 回答的一个小变种:
diff
检测大于 1 的差异。这会在每次运行结束1
时给出一个值。cumsum
)为每次运行分配不同的数字标签。累加是向后进行的,因为1
第 1 步中的值是在每次运行结束时,而不是在开始时。histc
。代码:
y = [diff(x)>1 1]; %// step 1
y = cumsum(fliplr(y)); %// step 2
y = fliplr(histc(y, 1:y(end))); %// step 3
这等于 的运行长度。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);
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)
这是经过测试的并且有效。