How to get the longest consecutive values in a column vector and the position at which it starts

Question

Meghana Balasubramanian 2019-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/477407-how-to-get-the-longest-consecutive-values-in-a-column-vector-and-the-position-at-which-it-starts

回答： Benjamin Hezrony 2024-10-6

采纳的回答： Bruno Luong

Hello,

Suppose i have a single column vector A'=[0 0 0 1 0 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 1]

I only want the longest consecutive values of 1's and display only that.

I'd really appreciate any help!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bruno Luong 2019-8-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/477407-how-to-get-the-longest-consecutive-values-in-a-column-vector-and-the-position-at-which-it-starts#answer_388866

编辑：Bruno Luong 2019-8-24

在 MATLAB Online 中打开

A=[0 0 0 1 0 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 1]'
i=reshape(find(diff([0;A;0])~=0),2,[]);
[lgtmax,jmax]=max(diff(i));
istart=i(1,jmax);
lgtmax % length of the longest sequence of 1s
istart % where it starts

3 个评论
显示 1更早的评论隐藏 1更早的评论

Bruno Luong 2019-8-24

编辑：Bruno Luong 2019-8-24

在 MATLAB Online 中打开

The terminologies are mine it doesn't matter for MATLAB syntax.

Actually my code is compact but difficult to understand, even for people who are initiated.

If you want to understand, you can split into multiple small and basic commands to analyze it and read the doc of corresponding command and experiment with your own examples.

[0;A;0]

Pad 0s at the header and trailer of A.

The command

diff([0;A;0])~=0

returns a logical array with TRUE whene there is a transitions (0->1) or (1->0).

And the TRUE positions must come in pairs: alternate in order of 0->1, 1->0, 0->, 1->0, etc ... since I take care to pad the arrays with 0s in both ends.

So by reshaping

i=reshape(... ,2,[]);

I'll get two-row array, the first row contains index positions of 0->1, and the second row of 1->0. You'll need to know about MATLAB major-column storage scheme for array to fully understand such trick.

Therefore the difference of the two columns is the the length of the consecutive 1s in A with this command:

diff(i);

By doing

[lgtmax,jmax]=max(diff(i));

I search for the longest sequence of 1s, and the number of this sequence is in retreived in jmax.

The last command

istart=i(1,jmax)

just maps back the starting position of this sequence in A (since I pad 0 in the head, it's actually the position of the first 1).

Laurent Haret 2022-5-16