Discussion :

[mat_de_lyon]

Bonjour,

Je suis à la recherche de plusieurs maximum dans un même vecteur. Par exemple prenons le vecteur a :
a = [0 1 2 4 6 9 12 11 8 4 2 2 4 6 7 9 7 5 4];
Je souhaite qu'il me sorte les 2 max locaux 12 et 9 avec leurs indices (ce qui doit être simple quand on a les max).
Mais je ne veux pas qu'il me sorte le 11 qui est certes supérieur à 9 mais qui n'est pas un max local.
Une écriture matricielle serait préférable vu qu'il y a beaucoup de max.

Merci d'avance

[Jerome Briot]

Tu devrais trouver la solution dans cette discussion : Détection de pics

[christophe_halgand]

Ma condition serait :

pour i>2
Si x(i) est supérieur à x(i-1) et x(i+1)

et pour i=1
Si x(i) est supérieur à x(i+1)

pour poursuivre ta recherche, tu dois tout simplement faire tourner une boucle sur toutes les données de ton tableau.

[mat_de_lyon]

Merci mais en fait mon problème est pointé dans la réponse de FR119492 : mon signal est bruité et il lui arrive d'avoir des pics parasites du style :
a = [0 1 2 4 6 9 12 11 8 4 5 2 2 4 6 7 9 7 5 4];
J'ai essayé la fonction smooth qui donne un résultat correct mais elle ne filtre pas les "gros" pics parasites.
Existe-t'il une autre fonction du même genre qui pourrait carrement fitter la courbe et de laquelle on pourrait en tirer le max?
Merci

[FR119492]

Salut !

Pour faire un filtrage, c'est-à-dire séparer le signal du bruit, tu dois avoir une définition claire de ce qui est du signal et ce qui est du bruit. Sans une telle définition rigoureuse, tu obtiendras n'importe quoi.

Jean-Marc Blanc

[mat_de_lyon]

Oui... La est tout le problème de l'expérimental. Je sais uniquement que j'ai des cycles qui se reproduisent plus ou moins régulièrement. En admettant que j'ai 1000 points par cycles (dont 1 max uniquement), si il y a petit blocage/bug/je ne sais quoi, je peux me retrouver avec 1 ou 2 points qui apparaitront à matlab comme un max (supérieur aux i-1 et i+1, même si les points sont lissés). Une idée serait de considérer la valeur de ce max par rapport à la valeur max approximative qui serait attendu (si on attend 10 +/- 0,5 et qu'on trouve un max à 4, on le vire). Je viens d'y penser en l'écrivant. J'essaie ça.

Merci

[Dam2227]

Salut

J'ai fait un programme de calibration qui fittait des données expérimentales avec beaucoup de bruit. Je m'en suis sorti en calculant le centre des pics par la méthode des barycentres, en imposant un seuil.

Si tu n'as pas trop de bruit, il y a une autre méthode qui consiste à ne pas considérer chaque points de tes data, mais un point tout les 10, 20, ... points (c'est à toi de voir ça), en répérant le changement de signe de l'intégrale.

Bon courage.