Discussion :

[pachalcs]

Bonjour,

J’ai un petit problème, en fait je travaille sur un programme écrit en matlab.
Et je voudrai optimiser le temps de calcul qu’effectuent ces lignes de codes suivantes

Code :

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

1
2
3
4
while t<=length(x)
   tab=[tab (x(t)+y([1:length(y)]))'];
   t=t+1;
end

Concrètement j’ai deux tableaux d’entiers :

Par exemple x = [1 2 6] et y = [2 1 5]
Et le résultat qui sera dans Tab est [3 2 6 4 3 7 8 7 11].
[3 2 6]= 1 + [2 1 5] et [4 3 7]= 2 + [2 1 5] et [8 7 11]= 6 + [2 1 5]
C'est-à-dire qu’on prend chaque valeur de x, on le somme avec toutes les valeurs de y et le résultat est concaténé au fur et à mesure dans Tab.

De plus le traitement de données que j’effectue s’effectue sur quelques dizaine de milliers de points. Par conséquent je voudrai que quelqu’un m’aide pour une vectorisation possible.

Je vous remercie d’avance,

[Merel]

Salut, tu peux déjà utiliser la préallocation de mémoire...

preallocation

[Jean Dumoncel]

Salut,

Code :

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

tab = repmat(y,length(x),1)+repmat(x(:),1,size(y,2))

après il ne reste plus qu'à mettre le résultat dans un vecteur ligne.

[pachalcs]

Merci Magellan, cela m'a permis de reduire mon temps de calcul de plus de 10 sur de grandes séries de données.

Bonne journée

[pachalcs]

Bonjour

Au fait j'ai ce code ci-dessus et je voudrai le remplacer par une autre tournure en utilisant une méthode de vectorisation.

Code :

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

1
2
3
4
for i=1:length(Tab_u)
   index_valeurs=find(abs(tab-Tab_u(i))<eps);
   tabpr_U(i)=sum(tab_pr(index_valeurs));
end

En effet je veux récupérer les index des occurences de chaque élément de Tab_u dans tab.
Par la suite ces index sont utilisés pour construire un autre vecteur (tabpr_U).

Merci d'avance,

[Jean Dumoncel]

et avec l'exemple que je t'ai donné, ça ne te donne pas des idées?

propose quelque chose (même si c'est faux), on t'aidera à le corriger.

[pachalcs]

Bonjour,

En effet je suis toujours bloqué sur le même probleme depuis vendredi et jusqu'à présent je suis ren réflexion. Si je n'y arrive pas, je vous ferai savoir par contre vous auriez pas une petite piste sans pour autant me donner la solution complete si vous l'avez.

Merci