Discussion :

[cgoury]

Bonjour,

j'ai cherché dans le forum mais n'ai rien trouvé à ce sujet.

J'ai 6 matrices A, B, C, D, E et F à 2 dimensions et je cherche à générer une matrice G qui soit la moyenne terme à terme des 6 premières matrices.
J'utilise une double boucle for qui créé pour chaque élément (i,j) une sous matrice formée des 6 élément A(i,j), B(i,j), C(i,j), D(i,j), E(i,j) et F(i,j), en fait la moyenne qui est stockée dans G(i,j).
De plus, comme j'ai des éléments non numériques Nan dans chacune des 6 matrices je suis quasi contrainte d'utiliser la fonction nanmean.

Y'a-t-il un moyen plus élégant et surtout plus rapide qu'une boucle for pour faire ce calcul? J'ai un temps de calcul proche de 3minutes...


D'avance merci beaucoup pour votre aide
Celine

[Jerome Briot]

À tester :

Code :

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

1
2
G = zeros(size(A,1), size(A,2));
G(:) = nanmean([A(:) B(:) C(:) D(:) E(:) F(:)],2);

[cgoury]

Génial! fonctionne parfaitement!
Je ne connaissais pas cette syntaxe, merci!

[Jerome Briot]

Un peu de lecture complémentaire : MATLAB : Qu'est-ce que l'indexage linéaire ? (s'applique aussi à Scilab et à GNU Octave)

[cgoury]

merci pour le lien super intéressant, je vais pouvoir optimiser pleins de choses...

Une autre question: pour quelle raison quand je remplace la fonction nanmean précédemment utilisée, par nanmedian qui a exactement la même syntaxe, mon programme tourne indéfiniment ?

[Jerome Briot]

pour quelle raison quand je remplace la fonction nanmean précédemment utilisée, par nanmedian qui a exactement la même syntaxe, mon programme tourne indéfiniment ?

Aucune idée. Peut être un bug ?