我正在使用Matlab进行我的一个项目。从某种程度上说,我实际上已经陷入了困境。尝试在谷歌搜索,但是,没有太大的成功。
我有一个0和1的数组。类似的东西:
A = [0,0,0,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,0,0,0,0];
我想提取一系列指标:[x_1, x_2, x_3, x_4, x_5, ..]
这样x_1
是第一个零范围的起始索引。 x_2
是第一个零范围结束的索引。
x_3
是第二个零范围的起始索引。 x_4
是第二个零范围结束的索引。
对于上面的例子:
x_1 = 1, x_2 = 3
x_3 = 9, x_4 = 10
等等。
当然,我可以通过编写一个简单的循环来实现。我想知道是否有更优雅(矢量化)的方法来解决这个问题。我正在考虑像前缀一样的东西,但是,现在没有运气。
谢谢,
阿尼尔。
答案 0 :(得分:3)
diff功能对于这类东西很有用而且很快。
temp = diff(A);
Starts = find([A(1) == 0, temp==-1]);
Ends = find([temp == 1,A(end)==0])
编辑:修复了gnovice捕获的结束计算中的错误。
答案 1 :(得分:2)
零之前没有其他零:A==0 & [true A(1:(end-1))~=0]
Zeros后面没有其他零:A==0 & [A(2:end)~=0 true]
使用以下每个加find
来获取零运行的开始和结束。然后,如果您真的希望它们在您描述的单个向量中,则将它们交错。
答案 2 :(得分:1)
如果你想在上面描述的单个向量中得到结果(即x = [x_1 x_2 x_3 x_4 x_5 ...]
),那么你可以使用函数DIFF执行二阶差分并找到大于0:
x = find(diff([1 A 1],2) > 0);
修改强>
以上内容适用于每个零字符串中至少 2个零的情况。如果您在A
中出现单个零,则可以修改上述零以处理它们:
diffA = diff([1 A 1],2);
[~,x] = find([diffA > 0; diffA == 2]);
在这种情况下,单个零值将在x
中创建重复索引(即如果A
以单个零开头,那么x(1)
和x(2)
都将是1)。