Discussion :

[titi_du_34]

Bonjour à tous,

Voilà, je souhaiterais obtenir la matrice A_ordre depuis la matrice A comme le montre l'exemple suivant

Code :

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

1
2
3
 
A=[1 1 1 1 1 2 2 2  3 4 4 9 9 9;    -3 1 5 0 -4  8 7 -2  5  9 12  2 7 -1]';
A_ordre=[1 2 3 4 9 1 2 4 9 1 2 9 1 1;   -3 8 5 9 2 1 7 12 7 5 -2 -1 0 -4]';

C'est à dire ranger les valeurs 1 1 1 1 1 2 2 2 3 4 4 9 9 9 de telle manière à ce que les valeurs doublées se répètent seulement lorsque toutes les valeurs sont passées une seule fois, soit 1 2 3 4 9 1 2 4 9 1 2 9 1 1.
Bien entendu, il faut que l'autre partie des coordonnées associées à ces valeurs (-3 1 5 0 -4 8 7 -2 5 9 12 2 7 -1) suivent la nouvelle matrice ordonnée (-3 8 5 9 2 1 7 12 7 5 -2 -1 0 -4).

Auriez vous des idées car là j'avoue que je sèche ?

Encore merci pour votre aide ?

Cdt,

Titi le petit poussin

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
A = [1 1 1 1 1 2 2 2 3 4 4 9 9 9 ; -3 1 5 0 -4 8 7 -2 5 9 12 2 7 -1];
 
flag = false(1,size(A,2));
 
A_ordre = zeros(size(A));
 
k = 1;
flag(k) = true;
A_ordre(:,k) = A(:,k);
 
while ~all(flag==true)
 
    for n = 2:size(A,2)
 
        if flag(n)
            continue
        end
 
        if A(1,n)==A_ordre(1,k) && ~all(A(1,n)==A(1,~flag))
            continue
        end
 
        k = k+1;
 
        A_ordre(:,k) = A(:,n);
 
        flag(n) = true;
 
    end
 
end

[titi_du_34]

Mais c'est impeccable , c'est exactement ça que je cherchais à faire. Je pensais que ça aurait pu être évident à réaliser avec les fonctions basiques de Matlab, (sort, diff, find) mais finalement ça n'est pas si évident et je n'ai pas trouvé de fonctions permettant de réaliser une telle opération..

Donc MERCI, MERCI, et REMERCIE,

Passe de joyeuses fêtes de fin d'année au passage,

Cdt,

Titi le petit poussin

[phryte]

Bonjour,

Un essai avec "unique" :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
A=[1 1 1 1 1 2 2 2 3 4 4 9 9 9;-3 1 5 0 -4 8 7 -2 5 9 12 2 7 -1];
AB=A;
[C1,IA1,IC1] =unique(AB(1,:));
A1=[C1;A(2,IA1)];
AB(1,IA1)=0;
while any(AB(1,:)) ~= 0
    [C,IA,IC]=unique(AB(1,:));
    A1=[[A1(1,:) C(2:end)];[A1(2,:) A(2,IA(2:end))]];
    AB(1,IA)=0;
end
A1

[titi_du_34]

Bonjour,

Tout d'abord bonne année
Merci de vous être penché sur le problème énoncé. Effectivement votre méthode marche tout aussi bien.
Je suis content de voir que vos deux méthodes marchent. La seule chose que je trouve un peu dommage, c'est qu'à ce jour je n'ai pas vu de commande directes faites par Mathworks qui permettent d'avoir ce type résultat.
Résultat des courses on est obligé de passer par des boucles pour résoudre le problème donné, or nous savons que les boucles while ou encore for ne sont malheureusement plus adaptées lorsqu'on a énormément de données à traiter.

[phryte]

Bonjour,
On n'est pas obligé de passer par des boucles. Il y a la récursivité pour les calculs répétitifs.
Cela n'est certainement pas plus rapide mais plus élégant tout en prenant plus de mémoire de stockage intermédiaire.
On peut sans doute encore optimiser ce code !

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
function R=rangement(A,init)
persistent A1 AB
if init ==1
    AB=A;
    [C1,IA1,IC1] =unique(AB(1,:));
    A1=[C1;A(2,IA1)];
    AB(1,IA1)=0;
end
if any(AB(1,:)) ~= 0
    [C,IA,IC]=unique(AB(1,:));
    A1=[[A1(1,:) C(2:end)];[A1(2,:) A(2,IA(2:end))]];
    AB(1,IA)=0;
    R=rangement(AB,0);
end
R=A1;

Appel :

Code :

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

1
2
A=.....;
R=rangelent(A,1)

[petchi]

Un petit "one liner" des familles (uniquement compatible depuis l'introduction de la version 2016b et de l’extension implicite):

Code :

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

1
2
3
4
 
A=[1 1 1 1 1 2 2 2 3 4 4 9 9 9; -3 1 5 0 -4 8 7 -2  5 9 12 2 7 -1]
[~,ind] = sort(sum(triu(A(1,:)==A(1,:)')))
A_ordre = A(:,ind)

Possible d'optimiser la consommation de mémoire en utilisant une boucle pour le calcule de

Code :

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

sum(triu(A(1,:)==A(1,:)'))

mais du coup plus de one liner