Assegna il conteggio sequenziale per le corse numeriche

Question

Desidero assegnare un valore numerico cumulativo per le esecuzioni sequenziali in un vettore binario. Quello che ho è

x = [0 0 0 1 1 0 1 1 1 0 1 0 0 0 0 0 0], 

e quello che vorrei è

y = [1 2 3 1 2 1 1 2 3 1 1 1 2 3 4 5 6]. 

La soluzione con somma/cumSum/unico/trovare gamma di funzioni me allude. Qualsiasi aiuto sarebbe molto apprezzato.

Answer 1

Ecco un modo:

a = arrayfun(@(x)(1:x), diff(find([1,diff(x),1])), 'uni', 0); 
[a{:}] 

L'idea è quella di generare un elenco dei 'tirature', vale a dire [3,2,1,3,1,1,6] nel tuo caso, poi basta concatenare un gruppo di vettori che contano ad ogni valore in tale elenco , ovvero cat(2, 1:3, 1:2, 1:1, 1:3.... Io uso arrayfun come collegamento per riapplicare l'operatore : e quindi utilizzare l'elenco separato da virgole che {:} restituisce come collegamento per la concatenazione.

Answer 2

(Non è un one-liner, ahimè ...):

F = find(diff(x))+1; 
y = ones(size(x)); 
y(F) = y(F)-diff([1,F]); 
y = cumsum(y); 

primo luogo, trovare tutte le posizioni in x dove c'è un cambiamento; quindi costruisci un vettore di 1 dove devi sottrarre la lunghezza di ogni segmento continuo. Alla fine, prendi il cumsum di esso.

Answer 3

Creare una matrice sparse tale che ogni corsa è su una colonna diversa, e poi fare la somma cumulativa:

t = sparse(1:numel(x), cumsum([1 diff(x)~=0]), 1); 
y = nonzeros(cumsum(t).*t).'; 

---

Utilizzare accumarray con una funzione personalizzata per generare ogni aumento modello in una cella diversa , e poi concatenare tutte le cellule:

y = accumarray(cumsum([1 diff(x)~=0]).', 1, [], @(x) {1:numel(x).'}); 
y = [y{:}];