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

Contribuez Discussion :

[Défis][1] Vecteur colonne ou vecteur ligne ?


Sujet :

Contribuez

  1. #1
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut [Défis][1] Vecteur colonne ou vecteur ligne ?
    ------------------------------------------------------

    Défi : créer une fonction qui, pour un vecteur donné, renvoie ce même vecteur, dans l'ordre inverse et dont chaque valeur est multipliée par 2.

    Contrainte
    : la fonction doit être opérationnelle pour un vecteur colonne comme pour un vecteur ligne.

    ------------------------------------------------------

    La recherche dans la documentation MATLAB avec le mot-clé flip nous donne :
    >> lookfor flip
    FLIPDIM Flip matrix along specified dimension.
    FLIPLR Flip matrix in left/right direction.
    FLIPUD Flip matrix in up/down direction.
    Il n'y a donc pas, a priori de fonction universelle, et il faut donc tester le type de vecteur en début de fonction, comme 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
    function V=test(V)
     
    if size(V,1)==1 % Vecteur ligne type V=[1 2 3]
     
       V = 2*fliplr(V);
       % ou
       % V = 2*flipdim(V,1);
     
    else % Vecteur colonne type V=[1;2;3]
     
       V = 2*flipud(V);
       % ou
       % V = 2*flipdim(V,2);
     
    end
    Testons :
    >> V=[1 2 3]

    V =

    1 2 3

    >> W=test(V)

    W =

    6 4 2

    >> V=[1;2;3]

    V =

    1
    2
    3

    >> W=test(V)

    W =

    6
    4
    2
    Tout est OK... mais est-il possible de faire plus court ?

    Indice : solution en une ligne (déclaration de la fonction mise à part) :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    function W=test(V)
     
    < Votre solution >
    ------------------------------------------------------

    Vous pouvez mettre vos solutions à la suite de ce message

    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  2. #2
    Membre régulier
    Inscrit en
    Décembre 2006
    Messages
    329
    Détails du profil
    Informations forums :
    Inscription : Décembre 2006
    Messages : 329
    Points : 88
    Points
    88
    Par défaut
    voici ma solution
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
     
    function W= test (V)
    W = 2*V(end:-1:1)
    end

  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 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par hanane78 Voir le message
    voici ma solution
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
     
    function W= test (V)
    W = 2*V(end:-1:1)
    end
    C'est en effet effectivement la solution. L'utilisation de l'indexage linéaire permet de s'affranchir de l'orientation du vecteur d'entrée.

    Ceci est relativement utile lorsque l'on développe un code prenant en argument d'entrée un vecteur.

    Dans ce cas, il y a trois solution :
    • soit on demande clairement à l'utilisateur de fournir un vecteur colonne ou un vecteur ligne.
    • soit on teste le type de vecteur en début de code.
    • soit on utilise l'indexage linéaire qui fonctionne aussi bien pour un vecteur colonne que pour un vecteur ligne.


    La dernière solution est donc la plus flexible

    ------------------------------------------------------

    Question subsidiaire
    :

    Est-il possible d'améliorer encore l'efficacité du code ci-dessus en trouvant un moyen de diviser par deux la quantité de mémoire utilisée par la fonction ?
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  4. #4
    Membre éprouvé

    Homme Profil pro
    Développeur informatique
    Inscrit en
    Décembre 2007
    Messages
    979
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 40
    Localisation : France, Hauts de Seine (Île de France)

    Informations professionnelles :
    Activité : Développeur informatique

    Informations forums :
    Inscription : Décembre 2007
    Messages : 979
    Points : 1 256
    Points
    1 256
    Par défaut
    Je me rend moi

    c'est quoi la solution ?
    AlloSchool, votre école sur internet.

  5. #5
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    J'ai oublié de préciser que cette optimisation de mémoire ne peut avoir lieu que si la fonction écrite précédemment est elle-même appelée dans une autre fonction plus globale.

    Indice : il existe trois façon d'appeler une fonction à l'intérieur d'une autre. C'est la technique la plus récente (MATLAB 7 et +) qui nous intéresse ici.

    Le contexte ici serait, par exemple, d'essayer d'optimiser un code qui renvoie un "Out of memory" pour le faire fonctionner. Donc avec un vecteur V de grande taille.
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  6. #6
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  7. #7
    Membre averti

    Inscrit en
    Août 2007
    Messages
    302
    Détails du profil
    Informations forums :
    Inscription : Août 2007
    Messages : 302
    Points : 352
    Points
    352
    Par défaut
    version iterative:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
     
    function b = test(a)
     
    s = numel(a);
    b = a;
     
    i = 1;
    while i <= s 
      b(i) = 2*a(s-i+1);
      i = i+1;
    end;
    version recc:

    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 a = test(b)
        a = b;
        num_el = numel(b);
     
        toto(1);
     
        function toto(i)
            a(i) = b(num_el-i+1);
     
            if i <= num_el/2
                toto(i+1);
            end        
        end
    end

  8. #8
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    La seconde solution proposée utilise la bonne méthode... mais pas complètement. En effet, passer les éléments d'un vecteur un par un réduira la quantité de mémoire nécessaire, mais ralentira très certainement le code si la taille du vecteur augmente. Il nous faut donc trouver comment utiliser cette méthode avec le vecteur complet.

    Il existe bien trois façons d'appeler une fonction sous MATLAB :
    1. en l'écrivant dans une fichier m séparé
    2. en tant que sous-fonction dans une autre fonction
    3. en l'imbriquant dans une autre fonction (fonctions imbriquée - nested functions)


    C'est la dernière qui nous intéresse ici. L'avantage premier des fonctions imbriquées est de pouvoir aisément partager des variables entre elles sans utiliser de techniques de sauvegarde/chargement temporaire.

    Un autre avantage très intéressant est d'éviter la copie temporaire des variables en mémoire si on utilise les fonctions imbriquées sans argument.

    Dans l'exemple suivant, la fonction de la première question a été renommée en
    doubleflip_nf :

    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
    function test_nf
     
    % Augmentez ou diminuez la valeur suivante en fonction 
    % des capacités de votre configuration
    n=15000000; 
     
    V=rand(n,1);
     
    disp(V([1 end])) % Pour vérification
     
    doubleflip_nf;
     
    disp(V([1 end])) % Pour vérification 
     
        function doubleflip_nf
     
            V = 2*V(end:-1:1);
     
        end
     
    end
    Dans ce cas, V n'est pas dupliqué dans doubleflip_nf. Le code utilise donc deux fois moins de mémoire que celui de la première réponse.

    Plus d'information => http://blogs.mathworks.com/loren/200...ariable-scope/
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  9. #9
    Membre averti

    Inscrit en
    Août 2007
    Messages
    302
    Détails du profil
    Informations forums :
    Inscription : Août 2007
    Messages : 302
    Points : 352
    Points
    352
    Par défaut
    Si on cherche a ne pas garder de copie de l'ancienne valeur n'est-il pas mieux de faire:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
     
     
    a = rand(1,1000);
    a = test(a);
     
    ou la function test est definie comme suit:
     
    function W= test (V)
    W = 2*V(end:-1:1)
    end
    D'autre part j'aimerai juste clarifier que le passage d'un argument a une fonction n'implique pas directement la copie de celui-ci. La variable sera copiee que si necessaire puisque matlab implemente la strategie 'Copy-On-Write'.

    Finalement la methode proposee met (a mon avis) en avant une mauvaise pratique de programmation: en particulier la fonction s'applique a une certaine variable en particulier et est donc difficillement reutilisable.

    Norez encore qu'il y a d'autres facon de creer des fonctions a 'effet de bord'. Un bon example est la programmation orientee objet dans MATLAB.

    Salutations,

    Gregoire

  10. #10
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par paradize3 Voir le message
    Si on cherche a ne pas garder de copie de l'ancienne valeur n'est-il pas mieux de faire:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
     
     
    a = rand(1,1000);
    a = test(a);
     
    ou la function test est definie comme suit:
     
    function W= test (V)
    W = 2*V(end:-1:1)
    end
    La quantité de mémoire utilisée sera différente. Dans ce cas, le vecteur temporaire W ne "sert à rien" mais il utilise autant de mémoire que V. Si la taille de V est significative, cela peut entrainer des problèmes de mémoire, problème récurrent sous MATLAB.

    Citation Envoyé par paradize3 Voir le message
    D'autre part j'aimerai juste clarifier que le passage d'un argument a une fonction n'implique pas directement la copie de celui-ci. La variable sera copiee que si necessaire puisque matlab implemente la strategie 'Copy-On-Write'.
    C'est exact. Le terme "V n'est pas dupliqué" dans mon précédent message est assez ambigu. J'aurais dû dire la "quantité de mémoire utilisée par les variables n'est pas doublée"

    Citation Envoyé par paradize3 Voir le message
    Finalement la methode proposee met (a mon avis) en avant une mauvaise pratique de programmation: en particulier la fonction s'applique a une certaine variable en particulier et est donc difficillement reutilisable.
    La méthode décrite ici n'est pas une règle de programmation. Elle permet de contourner un problème particulier (Out of memory), elle s'applique donc bien à un cas particulier. Dans ce sens, elle est parfaitement "réutilisable". A noter, que ce n'est bien entendu pas la seule méthode disponible : Out of memory. Type HELP MEMORY for your options.

    Citation Envoyé par paradize3 Voir le message
    Norez encore qu'il y a d'autres facon de creer des fonctions a 'effet de bord'. Un bon example est la programmation orientee objet dans MATLAB.
    Je ne me prononcerai pas sur ce point, n'ayant jamais utilisé la POO avec MATLAB. Les utilisateurs intéressés pourront consulter ces quelques liens :N'hésitez pas à continuer cette discussion
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  11. #11
    Membre averti

    Inscrit en
    Août 2007
    Messages
    302
    Détails du profil
    Informations forums :
    Inscription : Août 2007
    Messages : 302
    Points : 352
    Points
    352
    Par défaut
    Citation Envoyé par Dut Voir le message
    La quantité de mémoire utilisée sera différente. Dans ce cas, le vecteur temporaire W ne "sert à rien" mais il utilise autant de mémoire que V. Si la taille de V est significative, cela peut entrainer des problèmes de mémoire, problème récurrent sous MATLAB.
    je pensais que a = test(a) (meme variable d'entree que de sortie) resoudrai ce probleme? (je n'ai pas matlab sous la main pour confirmer)

    sur un autre sujet il me semble que faire v = w(end:-1:1) va allouer un tableau temporaire inutile de n elements: c'est l'intervale 'end:-1:1' . Je n'en suis pas 100% sur, peut-être que quelqu'un pourra confirmer.

    Comme je le comprends, en resumé lorsque qu'on écrit:
    où W est un vecteur de N éléments, la mémoire maximum (temporairement utilisée pour pouvoir executer cette instruction) est de:
    N (pour W) + N (pour la version temporaire de W, necessaire à l'execution de l'instruction) + N pour l'intervale = 3*N. Est-ce correct? Une version itérative permetterai une utilisation optimale de mémoire (N éléments).

    Salutations,

    Grégoire

  12. #12
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    J'ai comparé ces deux méthodes :

    Sous-fonction :
    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
    function test_sf
     
    pause(.5)
    disp('start')
     
    % Augmentez ou diminuez la valeur suivante en fonction 
    % des capacités de votre configuration
    n=15000000; 
     
    V=rand(n,1); 
     
    disp(V(end)) % Pour vérification 
     
    pause(2)
     
    disp('DoubleFlip')
     
    V=doubleflip_sf(V);
     
    disp(V(1)) % Pour vérification 
    pause(2)
     
    function W=doubleflip_sf(V)
     
    W = 2*V(end:-1:1);
    pause(2)
    Fonctions Imbriquées (sans argument) :
    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
    function test_nf
     
    pause(.5)
    disp('start')
     
    % Augmentez ou diminuez la valeur suivante en fonction 
    % des capacités de votre configuration
    n=15000000; 
     
    V=rand(n,1); 
     
    pause(2)
     
    disp(V([1 end])) % Pour vérification 
     
    pause(2)
     
    disp('DoubleFlip')
    doubleflip_nf;
    disp(V([1 end])) % Pour vérification 
    pause(2)
     
        function doubleflip_nf
     
            V = 2*V(end:-1:1);
            pause(2)
     
        end
     
    end

    Avec un vecteur V classe Double de taille 15000000x1 => 15000000*8/1024/1024 = 115 Mo (environ)

    J'ai utilisé l'outil MATLAB Monitoring Tool disponible sur le FEX, comme ceci :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    monitormatlab;test_sf;pause(5);test_nf;
    Voici le graphique obtenu :


    Note : il y a un résidu de mémoire MATLAB constant d'environ 10 Mo, le graphique est donc légèrement décalé vers le haut

    Evolution de la mémoire MATLAB (lancement de monitormatlab, appel de test_sf, pause de 5s , appel de test_nf) :

    1. Appel de test_sf, la mémoire grimpe à 115 Mo à la définition de la variable V.
    2. Appel de doubleflip_sf, la mémoire monte à 230 Mo (V+W)
    3. Sortie de doubleflip_sf, la mémoire chute à 115 Mo (destruction de W)
    4. Sortie de test_sf, la mémoire revient à 0 Mo
    5. Palier de 5s à 0 Mo
    6. Appel de test_nf, la mémoire MATLAB grimpe à 115 Mo à la définition de la variable V.
    7. Appel de doubleflip_nf, la mémoire reste à 115 Mo (pas de vecteur temporaire crée)
    8. Sortie de doubleflip_nf, la mémoire reste à 115 Mo
    9. Sortie de test_nf, la mémoire revient à 0 Mo
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  13. #13
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par paradize3 Voir le message
    sur un autre sujet il me semble que faire v = w(end:-1:1) va allouer un tableau temporaire inutile de n elements: c'est l'intervale 'end:-1:1' . Je n'en suis pas 100% sur, peut-être que quelqu'un pourra confirmer.
    Il me semble que non => Is subscript referencing slower than : ?

    Citation Envoyé par Mike Karr
    "A(1:n)" and "I = 1:n; A(I)" have different execution paths,
    and the second can be slower. The reason is that the 1:n is physically
    materialized in the second case, but not in the first. It is something
    we plan to fix.
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  14. #14
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par paradize3 Voir le message
    Une version itérative permetterai une utilisation optimale de mémoire (N éléments).
    Tout à fait, mais le temps d'exécution du code pourrait s'allonger sous MATLAB.

    Dans ce cas, il serait encore plus judicieux de passer en langage bas niveau C/Fortran à l'aide des fichiers MEX.
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  15. #15
    Membre averti

    Inscrit en
    Août 2007
    Messages
    302
    Détails du profil
    Informations forums :
    Inscription : Août 2007
    Messages : 302
    Points : 352
    Points
    352
    Par défaut
    Je pensais que ML arriverait a se 'rendre compte' que la variable d'entree etait la meme que la variable de sortie et donc de ne pas allouer de memoire inutilement.

    Question subsidiaire, est-ce que si on ecrit la fonction en nommant les variables de sortie de la meme facon que les variables d'entrees:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    function V = test(V)
    V = V(end:-1:1)
    et qu'on appelle:

    a = rand(1,100);
    a = test(a)

    matlab arrive-t-il a optimiser la memoire allouee?

    Pour la remarque concernant le fichier MEX, il n'est malheureusement pas possible de modifier une variable d'entree dans un fichier mex, il ne serait donc pas possible de faire des algorithmes 'en place'.

    Salutations,

    Gregoire

  16. #16
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par paradize3 Voir le message
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    function V = test(V)
    V = V(end:-1:1)
    Effectivement, cela résout le problème de mémoire.
    En utilisant la fonction ci-dessus comme sous-fonction dans le test précédent, on obtient ceci :



    Pour conclure, je précise que ce microDéfi repose sur une contribution du FEX :

    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  17. #17
    Membre éprouvé
    Avatar de ol9245
    Homme Profil pro
    Chercheur
    Inscrit en
    Avril 2007
    Messages
    985
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 62
    Localisation : France, Hérault (Languedoc Roussillon)

    Informations professionnelles :
    Activité : Chercheur

    Informations forums :
    Inscription : Avril 2007
    Messages : 985
    Points : 1 158
    Points
    1 158
    Billets dans le blog
    1
    Par défaut alternative avec flipud
    flipud permet d'économiser le calcul du vecteur d'index 1:end

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    >> 2*flipud(1:3)
     
    ans =
     
         2     4     6
     
    >> 2*flipud((1:3)')
     
    ans =
     
         6
         4
         2
    Ceci dit, je ne sais pas si Matlab effectue ce flip "sur place" ou "par copie" et, d'une manière générale, si la solution est vraiment meilleure. On va donc benchmarker (on ne recule devant rien !)

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    function test_flipud
    A = 1:1e6 ;
    tic
    A = flipud(A) ;
    toc
     
    tic
    A = A(end:-1:1) ;
    toc
    end
    ce qui donne :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    Elapsed time is 0.011297 seconds.
    Elapsed time is 0.015921 seconds.
    Gain de temps de 30% : y a pas photo !
    "La vraie grandeur se mesure par la liberté que vous donnez aux autres, et non par votre capacité à les contraindre de faire ce que vous voulez." Larry Wall, concepteur de Perl.

  18. #18
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    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 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Attention de bien comparer la même chose...

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    function test_flipud
    clc
    
    A = (1:1e6).' ;
    tic
    A = flipud(A);
    toc
    
    tic
    A = A(end:-1:1);
    toc
    
    end
    donne :
    Elapsed time is 0.011872 seconds.
    Elapsed time is 0.003905 seconds.
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    function test_fliplr
    clc
    
    A = 1:1e6 ;
    
    tic
    A = fliplr(A) ;
    toc
    
    tic
    A = A(end:-1:1) ;
    toc
    
    end
    donne :
    Elapsed time is 0.013362 seconds.
    Elapsed time is 0.004012 seconds.
    Il suffit d'éditer le code de FLIPUD et FLIPLR pour comprendre
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  19. #19
    Membre éprouvé
    Avatar de ol9245
    Homme Profil pro
    Chercheur
    Inscrit en
    Avril 2007
    Messages
    985
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 62
    Localisation : France, Hérault (Languedoc Roussillon)

    Informations professionnelles :
    Activité : Chercheur

    Informations forums :
    Inscription : Avril 2007
    Messages : 985
    Points : 1 158
    Points
    1 158
    Billets dans le blog
    1
    Par défaut
    Bravo ! ! ! ! Là : Bravo !!
    Je m'incline devant tant de maitrise
    "La vraie grandeur se mesure par la liberté que vous donnez aux autres, et non par votre capacité à les contraindre de faire ce que vous voulez." Larry Wall, concepteur de Perl.

  20. #20
    Membre habitué
    Homme Profil pro
    étudiant ingénieur traitement images médicales
    Inscrit en
    Juin 2012
    Messages
    74
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 35
    Localisation : France, Nord (Nord Pas de Calais)

    Informations professionnelles :
    Activité : étudiant ingénieur traitement images médicales
    Secteur : Santé

    Informations forums :
    Inscription : Juin 2012
    Messages : 74
    Points : 131
    Points
    131
    Par défaut
    Jolie solution, et belle maitrise des rouages de Matlab.

    D'autres défis à proposer? (je viens de m'inscrire sur ce forum, et découvre ces défis solutionnés).

Discussions similaires

  1. vecteur colonne vers ligne
    Par magictol53 dans le forum MATLAB
    Réponses: 2
    Dernier message: 28/10/2013, 13h21
  2. addition de vecteurs colonne et ligne
    Par Dam2227 dans le forum MATLAB
    Réponses: 5
    Dernier message: 18/05/2012, 01h38
  3. [Débutant] transformer des vecteurs colonne en un seul vecteur ligne
    Par sarabouch dans le forum MATLAB
    Réponses: 6
    Dernier message: 14/10/2011, 14h35
  4. transformer un vecteur ligne en vecteur colonne
    Par samiov dans le forum MATLAB
    Réponses: 4
    Dernier message: 26/03/2010, 10h42
  5. transposée vecteur ligne -> vecteur colonne
    Par Pierre845 dans le forum Macros et VBA Excel
    Réponses: 2
    Dernier message: 19/03/2008, 19h29

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