IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

Parcours d'une matrice


Sujet :

MATLAB

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre expérimenté Avatar de LordPeterPan2
    Inscrit en
    Février 2007
    Messages
    247
    Détails du profil
    Informations forums :
    Inscription : Février 2007
    Messages : 247
    Par défaut Parcours d'une matrice
    Bonjour,

    Je voudrais trouver une manière assez rapide pour parcourir une matrice n x m de la façon suivant (dans l'exemple n = 6 et m=4) :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
     
          |   1   2   6    7 |
          |   3   5   8   14 |
          |   4   9  13   15 |
          |  10  12  16   21 |
          |  11  17  20   22 |
          |  18  19  23   24 |
    Ma première idée est de créer la matrice ci dessous pour ensuite me servir de l'indexage des matrices, mais là encore ce n'est pas évident (en ne faisant pas de boucle for).

    Pour info je cherche à refaire l'algorithme de codage JPEG.

  2. #2
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 317
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 317

  3. #3
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 317
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 317
    Par défaut
    Sinon en fouillant dans mes archives personnelles

    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
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    function [r,c]=zigzag(n,m)
    % ZIGZAG...(Complete H1 line)
    %   ZIGZAG...(Complete file description)
    %
     
    %   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
    %   Contact: dutmatlab@yahoo.fr
    %   Revision: 1.0 (15-Feb-2007)
    %   Comments:
    %
     
    % Check the number of input arguments
    error(nargchk(1,2, nargin))
     
    if nargin==1
        m=n;
    end
     
     
    idx=reshape(1:n*m,n,m);
    idx=idx(:,end:-1:1);
     
    r=zeros(1,numel(idx));
     
    flag=0;
    x=1;
    j=1;
     
    if n>m
     
        idx=rot90(idx,3);
        idx=idx(:,end:-1:1);
        [n,m]=size(idx);
    end
     
    if n==m
        dim=min([1:m-1 m m-1:-1:1],n);
    else
        dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
    end
     
     
    for k=m-1:-1:1-m
     
    	if isempty(diag(idx,k))
     
            break
    	end
     
    	if ~flag
            r(x:x+dim(j)-1)=diag(idx,k);
    	else
            r(x+dim(j)-1:-1:x)=diag(idx,k);
    	end
     
    	flag=~flag;
    	x=x+dim(j);
        j=j+1;
     
    end
     
    if n>m
        r=r(end:-1:1);
    end
     
    if nargout==2
        [r,c]=ind2sub(size(idx),r);
    end
    A utiliser comme ceci :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    M = [1   2   6    7
         3   5   8   14
         4   9  13   15
        10  12  16   21
        11  17  20   22
        18  19  23   24];
     
    M = M.';
     
    idx = zigzag(size(M,1),size(M,2));
     
    M(idx)
    Le code n'est pas vraiment achevé mais c'est l'idée...

  4. #4
    Membre expérimenté Avatar de LordPeterPan2
    Inscrit en
    Février 2007
    Messages
    247
    Détails du profil
    Informations forums :
    Inscription : Février 2007
    Messages : 247
    Par défaut
    Merki bicoup,

    Je préfère (à priori) la seconde solution qui ne se limite pas aux matrices carrées

  5. #5
    Membre expérimenté Avatar de LordPeterPan2
    Inscrit en
    Février 2007
    Messages
    247
    Détails du profil
    Informations forums :
    Inscription : Février 2007
    Messages : 247
    Par défaut


    En faite çà fait un zigzag dans l'autre sens :

    Avec ta macro la matrice (pour une 3 x 5) est lu comme suit :

    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
     
    >> A=reshape(1:15,3,5)
     
    A =
     
         1     4     7    10    13
         2     5     8    11    14
         3     6     9    12    15
     
    >> zigzag(3,5)
     
    ans =
     
      Columns 1 through 14
     
         1     2     4     7     5     3     6     8    10    13    11     9    12    14
     
      Column 15
     
        15
    Alors que je voudrai :

    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
     
    >> A=reshape(1:15,3,5)
     
    A =
     
         1     4     7    10    13
         2     5     8    11    14
         3     6     9    12    15
     
    >> zigzag(3,5)
     
    ans =
     
      Columns 1 through 14
     
         1     4     2     3     5     7    10     8     6     9    11    13    14    12
     
      Column 15
     
        15

  6. #6
    Membre expérimenté Avatar de LordPeterPan2
    Inscrit en
    Février 2007
    Messages
    247
    Détails du profil
    Informations forums :
    Inscription : Février 2007
    Messages : 247
    Par défaut
    Merci Dut, en modifiant ton code (y avait pas de copyright ) j'ai réussi à faire ce que je voulais, je met le code aux cas où çà intéresse quelqu'un d'autre :

    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
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127
    128
    129
    130
    131
    132
    133
    134
    135
    136
    137
    138
    139
    140
    141
    142
    143
    144
    145
    146
    147
    function [vect_elt]=zigzag(A);
    
    [n,m] = size(A) ;
    vect_elt = A(zigzag_aux(n,m));
    
    function [r,c]=zigzag_aux(n,m)
    % ZIGZAG...(Complete H1 line)
    %   ZIGZAG...(Complete file description)
    %
    
    %   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
    %   Contact: dutmatlab@yahoo.fr
    %   Revision: 1.0 (15-Feb-2007)
    %   Comments:
    %
    
    % Check the number of input arguments
    error(nargchk(1,2, nargin))
    
    if nargin==1
        m=n;
    end
    
    idx=reshape(1:n*m,n,m);
    idx=idx(:,end:-1:1);
    
    r=zeros(1,numel(idx));
    
    flag=0;
    x=1;
    j=1;
    flip = 0 ;
    
    if n>m
        idx=rot90(idx,3);
        idx=idx(:,end:-1:1);
        [n,m]=size(idx);
        flip = 1;
    end
    
    if n==m
        dim=min([1:m-1 m m-1:-1:1],n);
    else
        dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
    end
    
    
    for k=m-1:-1:1-m
    
    	if isempty(diag(idx,k))
            
            break
    	end
    
        % 	if ~flag
    	if flag
            r(x:x+dim(j)-1)=diag(idx,k);
    	else
            r(x+dim(j)-1:-1:x)=diag(idx,k);
    	end
        
    	flag=~flag;
    	x=x+dim(j);
        j=j+1;
        
    end
    
    if flip
        r=r(end:-1:1);
    end
    
    if nargout==2
        [r,c]=ind2sub(size(idx),r);
    endfunction [vect_elt]=zigzag(A);
    
    [n,m] = size(A) ;
    vect_elt = A(zigzag_aux(n,m));
    
    function [r,c]=zigzag_aux(n,m)
    % ZIGZAG...(Complete H1 line)
    %   ZIGZAG...(Complete file description)
    %
    
    %   Author: Jérôme Briot, Matlab 6.1.0.450 (R12.1)
    %   Contact: dutmatlab@yahoo.fr
    %   Revision: 1.0 (15-Feb-2007)
    %   Comments:
    %
    
    % Check the number of input arguments
    error(nargchk(1,2, nargin))
    
    if nargin==1
        m=n;
    end
    
    idx=reshape(1:n*m,n,m);
    idx=idx(:,end:-1:1);
    
    r=zeros(1,numel(idx));
    
    flag=0;
    x=1;
    j=1;
    flip = 0 ;
    
    if n>m
        idx=rot90(idx,3);
        idx=idx(:,end:-1:1);
        [n,m]=size(idx);
        flip = 1;
    end
    
    if n==m
        dim=min([1:m-1 m m-1:-1:1],n);
    else
        dim=min([1:n repmat(n,1,m-n-1) n:-1:1],n);
    end
    
    
    for k=m-1:-1:1-m
    
    	if isempty(diag(idx,k))
            
            break
    	end
    
     	if flag
            r(x:x+dim(j)-1)=diag(idx,k);
    	else
            r(x+dim(j)-1:-1:x)=diag(idx,k);
    	end
        
    	flag=~flag;
    	x=x+dim(j);
        j=j+1;
        
    end
    
    if flip
        r=r(end:-1:1);
    end
    
    if nargout==2
        [r,c]=ind2sub(size(idx),r);
    end


+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. Parcours d'une matrice triangulaire superieure
    Par proc02 dans le forum Langage
    Réponses: 9
    Dernier message: 09/02/2014, 13h05
  2. concernant parcours d'une matrice
    Par zied_m dans le forum VB.NET
    Réponses: 2
    Dernier message: 25/01/2012, 14h14
  3. parcours d'une matrice
    Par doudi20 dans le forum Général Java
    Réponses: 5
    Dernier message: 17/03/2009, 15h51
  4. Parcours d'une matrice/tableau à deux dimensions
    Par yal001 dans le forum Mathématiques
    Réponses: 5
    Dernier message: 20/11/2008, 14h59
  5. Parcours d'une matrice(Gantt)
    Par asmiti dans le forum Algorithmes et structures de données
    Réponses: 19
    Dernier message: 17/01/2008, 10h40

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo