Discussion :

[cailloup]

Bonjour, ne sachant pas trop quoi taper pour chercher la réponse à ma question je poste un nouveau sujet

Pour éviter une boucle for, je me demande s'il est possible de remplir des cellules par la multiplication d'une matrice par des nombres:

Code :

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

1
2
3
4
b=cell(4,1);
a=[1 2 3];
b{1:4}=a*{1:4}
??? Undefined function or method 'mtimes' for input arguments of type 'cell'.

Cette syntaxe ne fonctionne pas par exemple
Je souhaiterai arriver à ça:

Code :

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

1
2
3
4
5
6
7
8
b{1}=a*1;b{2}=a*2;b{3}=a*3;b{4}=a*4
 
b = 
 
    [1x3 double]
    [1x3 double]
    [1x3 double]
    [1x3 double]

Sauriez-vous comment faire si c'est possible?
Merci!

NB: a est en réalité une matrice 3D

[le fab]

quelque chose comme:

Code :

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

b =cellfun(@(x) a*x,mat2cell((1:4).',ones(1,4)),'UniformOutput',false)

[cailloup]

Excellent, merci!

Résolu

[cailloup]

Après avoir implémenté l'astuce dans mon calcul, je suis passé d'un temps de calcul de 8s avec la boucle for à 7.2s avec cellfun, je m'attendais à bien mieux

Si je double la longueur de ma boucle for je passe à 16s et 14.4s pour cellfun

Théoriquement avec Matlab et le calcul matriciel on pourrait avoir un temps de calcul qui ne dépend que très peu de la longueur des vecteurs non ?

Si quelqu'un a une autre proposition qui pourrait accélérer ce type de calcul, je suis preneur!

[Invité]

Bonjour,

Le mécanisme interne de la fonction cellfun est très proche de la boucle : il applique la fonction dont on passe le handle à chaque cellule. Donc c'est normal que les temps soient proches.

Théoriquement avec Matlab et le calcul matriciel on pourrait avoir un temps de calcul qui ne dépend que très peu de la longueur des vecteurs non ?

Le calcul matriciel va être en général plus rapide, mais va consommer plus de mémoire. Par contre, le temps de calcul sera toujours fortement dépendant de la taille des données.

Il sera difficile je pense de faire mieux... Voici une petite amélioration en utilisant arrayfun :

Code :

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

1
2
b = 1:4;
arrayfun(@(x) a*x,b,'UniformOutput',false)

[cailloup]

même résultat niveau temps de calcul
mais l'écriture est un peu moins lourde

Je resterai la dessus alors
Merci pour l'aide