Précédent   Forum des professionnels en informatique > Environnements de développement > MATLAB > Contribuez
Contribuez Proposez vos tutoriels, FAQ, sources, astuces MATLAB
Partagez cette discussion sur d'autres réseaux sociaux : Viadeo Twitter Google Facebook Digg Delicious MySpace Yahoo
Réponse Actualité déjà publiée
 
Outils de la discussion
Publicité
Vieux 02/01/2008, 15h25   #1 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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 :
Citation:
>> 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 :
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 :
Citation:
>> 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 :
function W=test(V)

< Votre solution >
------------------------------------------------------

Vous pouvez mettre vos solutions à la suite de ce message


Dernière modification par Dut ; 15/01/2008 à 11h08.
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 15/01/2008, 13h09   #2 (permalink)
Membre Confirmé
 
Date d'inscription: décembre 2006
Localisation: salé(Maroc)
Messages: 252
Par défaut

voici ma solution
Code :
function W= test (V)
W = 2*V(end:-1:1)
end
hanane78 est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 15/01/2008, 14h01   #3 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
Par défaut

Citation:
Envoyé par hanane78 Voir le message
voici ma solution
Code :
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 ?
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)

Dernière modification par Dut ; 15/01/2008 à 16h39.
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 17/01/2008, 10h27   #4 (permalink)
Membre émérite
 
Avatar de mr_samurai
 
Date d'inscription: décembre 2007
Localisation: ile-de-france
Âge: 26
Messages: 942
Par défaut

Je me rend moi

c'est quoi la solution ?
mr_samurai est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 17/01/2008, 11h42   #5 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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.
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 18/01/2008, 09h10   #6 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
Par défaut

Un dernier indice => http://blogs.mathworks.com/loren/200...ariable-scope/
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 21/01/2008, 04h28   #7 (permalink)
Membre Confirmé
 
Date d'inscription: août 2007
Messages: 293
Par défaut

version iterative:

Code :
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 :
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

Dernière modification par paradize3 ; 21/01/2008 à 04h47.
paradize3 est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 21/01/2008, 08h48   #8 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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 :
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/
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 00h45   #9 (permalink)
Membre Confirmé
 
Date d'inscription: août 2007
Messages: 293
Par défaut

Si on cherche a ne pas garder de copie de l'ancienne valeur n'est-il pas mieux de faire:

Code :
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
paradize3 est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 08h57   #10 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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 :
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
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 11h01   #11 (permalink)
Membre Confirmé
 
Date d'inscription: août 2007
Messages: 293
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:
Code :
 
W = W(end:-1:1)
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
paradize3 est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 11h40   #12 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
Par défaut

J'ai comparé ces deux méthodes :

Sous-fonction :
Code :
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 :
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 :
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
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)

Dernière modification par Dut ; 01/02/2008 à 09h08.
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 11h45   #13 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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.
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 23/01/2008, 12h02   #14 (permalink)
Dut
Responsable C
 
Avatar de Dut
 
Date d'inscription: novembre 2006
Localisation: Toulouse, France
Âge: 33
Messages: 10 016
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.
__________________
Mes articles et critiques de livre MATLAB
.
Langages : MATLAB - C | Outils : MATLAB R2009a - Code::Blocks (svn 5966) - GCC 4.4.1 | OS : Kubuntu 9.10 - Windows XP SP2
.

• 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)
Hjartað Hamast (Bamm Bamm Bamm) (traduction française : Le cœur bat) (enfin je crois... je parle pas islandais !)
Dut est déconnecté   Envoyer un message privé Réponse avec citation
Vieux 01/02/2008, 01h30   #15 (permalink)
Membre Confirmé
 
Date d'inscription: août 2007
Messages: 293
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 :
function V = test(V)
V = V(end:-1:1)
et qu'on appelle:

Citation:
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

Dernière modification par Dut ; 01/02/2008 à 08h55.
paradize3 est déconnecté   Envoyer un message privé Réponse avec citation
NEWS MATLABFAQ MATLABTutoriels MATLABLivres MATLABSources MATLAB

Réponse Actualité déjà publiée

Précédent   Forum des professionnels en informatique > Environnements de développement > MATLAB > Contribuez



Outils de la discussion

Règles de messages
Vous ne pouvez pas créer de nouvelles discussions
Vous ne pouvez pas envoyer des réponses
Vous ne pouvez pas envoyer des pièces jointes
Vous ne pouvez pas modifier vos messages

Les balises BB sont activées : oui
Les smileys sont activés : oui
La balise [IMG] est activée : oui
Le code HTML peut être employé : non
Trackbacks are non
Pingbacks are non
Refbacks are non



Fuseau horaire GMT +1. Il est actuellement 18h19.


Vos questions techniques : forum d'entraide MATLAB - Publiez vos articles, tutoriels et cours
et rejoignez-nous dans l'équipe de rédaction du club d'entraide des développeurs francophones
Nous contacter - Hébergement - Participez - Copyright © 2000-2010 www.developpez.com - Legal informations.