零一矩阵中零范围的指数

时间:2011-03-26 21:41:31

标签: arrays algorithm matlab

我正在使用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 

等等。

当然,我可以通过编写一个简单的循环来实现。我想知道是否有更优雅(矢量化)的方法来解决这个问题。我正在考虑像前缀一样的东西,但是,现在没有运气。

谢谢,

阿尼尔。

3 个答案:

答案 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)。