Discussion :

[abel413]

Bonjour

Maintenant dans mon code complet il se trouve que j'appelle la fonction sort à un endroit donné :

Code :

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

xmax_Imp1_sort=sort(xmax_Imp1);

Je suis obligé d'utiliser cette dernière pour trier un vecteur colonne, en l'occurence xmax_Imp1 mais elle prend du temps également et tri simplement par ordre croissant les nombres du vecteur colonne xmax_Imp1. Penses-tu que, de même que tu viens d'écrire ici une fonction pour faire le calcul "sophistiqué" que je demandais, il soit possible d'en faire une similaire qui trie simplement les nombres par ordre croissant (pour un vecteur 1 colonne) et qui soit plus rapide que sort ? Éventuellement, sans boucle for.

merci

[Jerome Briot]

sort fait partie des fonctions optimisées sous MATLAB.

Je doute que tu puisses en améliorer le temps d'exécution.

[abel413]

Effectivement, je me doute bien que s'il existe une telle fonction, il va être difficile de faire mieux. Mais le fait est que dans mon code c'est bien cette fonction qui ralenti le travail. Voici un bout du code :

Code :

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

1
2
3
4
5
6
7
8
9
10
%pour valeurs_Imp2
 
valeurs2=medfilt1(valeurs_Imp2,300);
X_Imp2=find(valeurs_Imp2>valeurs2);
valeurs_Imp2(X_Imp2)=valeurs2(X_Imp2);
[ymax_Imp2,xmax_Imp2,ymin_Imp2,xmin_Imp2] = extrema(valeurs_Imp2);
xmax_Imp2_sort=sort(xmax_Imp2);
ix2=0:floor(size(xmax_Imp2_sort,1)/2)-1;
X2_Imp2_liquide_melange_evenements=[xmax_Imp2_sort(1+2*ix2), xmax_Imp2_sort(2*(ix2+1))];
X2_Imp2_liquide_melange_evenements_1colonne=creerbis(X2_Imp2_liquide_melange_evenements);

La fonction extrema est une fonction qui m'est donnée et que je ne peux pas changer quoi qu'il arrive. La fonction sort est employée une fois ainsi que la fonction creerbis (qui est la variante de Dut pour un autre problème que j'avais et qui est relativement rapide). De plus la fonction medfilt1 n'est pas non plus instantanée mais je ne pense pas que je puisse y faire grand chose.

[Jerome Briot]

Tu peux supprimer le find :

Code :

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

X_Imp2 = valeurs_Imp2>valeurs2;

C'est l'indexage logique

[abel413]

Oui merci Dut, je crois que tu me l'avais déjà signalé.

[abel413]

En fait, je viens d'essayer l'indexage dont tu parles :

Code :

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

X_Imp2 = valeurs_Imp2>valeurs2;

Donc si je comprends bien je stocke dans X_Imp2 1 si la valeurs de valeurs_Imp2 est supérieure à celle de valeurs2 et 0 sinon. En gros les deux blocs suivant sont équivalents :

Code :

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

1
2
3
4
5
6
7
X_Imp2=find(valeurs_Imp2>valeurs2);
valeurs_Imp2(X_Imp2)=valeurs2(X_Imp2);
 
%ou
 
X_Imp2=valeurs_Imp2>valeurs2;
valeurs_Imp2(X_Imp2)=valeurs2(X_Imp2);

Et donc au final, j'aurais le même résultat dans valeurs_Imp2 mais la seconde variante est réputée plus rapide, est-ce ce que tu voulais me faire faire ?