Discussion :

[Laulau29]

Bonjour,

Je dois traiter une matrice dont les éléments ont la forme suivante :
1 1 1 1 2 2 3 3 3 3 4 4 2 2 2 2 2 2
je voudrais pouvoir obtenir une matrice contenant le numéro (1, 2 3 et 4) avec sa première position et sa dernière.
Pour mon exemple, je devrais obtenir :
1 1 4
2 5 6
3 7 10
4 11 12
2 13 18.
J'ai une solution en effectuant une boucle for, mais j'aurai voulu savoir si une solution plus rapide existe.

Merci

[phryte]

Bonjour,

Une première approche :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
clear
x=[1 1 1 1 2 2 3 3 3 3 4 4 2 2 2 2 2 2];
[C,IA,IC]=unique(x)
D=diff(x)
F=find(D~=0)
R=[C' IA F']

[Laulau29]

J'avais essayé avec la fonction unique mais je perds les données de fin de la matrice (retour à 2)

[phryte]

Bonjour,

Un peu tiré par les cheveux.
Sans doute pas robuste.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
clear
x=[1 1 1 1 2 2 3 3 3 3 4 4 2 2 2 2 2 2];
xx=[x max(x)+1];
[C,IA,~]=unique(xx);
IA(end)=IA(end)-C(end)-1;
D=diff([[xx]]);
F=find(D~=0);
R=[x(F)' IA F']

[Jerome Briot]

Comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
M = [1 1 1 1 2 2 3 3 3 3 4 4 2 2 2 2 2 2]
 
idxA = [1 diff(M)]~=0;
 
idxB = find(idxA);
 
idxC = [idxB(2:end)-1 numel(M)];
 
[M(idxA(:)).' idxB(:) idxC(:)]

Ce qui donne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
ans =
 
     1     1     4
     2     5     6
     3     7    10
     4    11    12
     2    13    18

[Laulau29]

Comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
M = [1 1 1 1 2 2 3 3 3 3 4 4 2 2 2 2 2 2]
 
idxA = [1 diff(M)]~=0;
 
idxB = find(idxA);
 
idxC = [idxB(2:end)-1 numel(M)];
 
[M(idxA(:)).' idxB(:) idxC(:)]

Ce qui donne :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
ans =
 
     1     1     4
     2     5     6
     3     7    10
     4    11    12
     2    13    18

Merci pour cette solution.
Mon problème est résolu