détection de pic

salut
je cherche à détecter le pic d'une courbe.
en fait j'ai pensé à détecter tous les pics de la courbe par une fonction Extrema, ensuite trier le vecteur de sortie.

Code:

---

1 2 3 4 5 6

[xmax,ymax]=extrema(f); % f est la courbe X=sort(xmax,'descend'); X1=X(1); Y=f(X1); figure;plot(f);hold on; plot(X1,Y);

---

mais voyant ici le pic n'est pas bien détecté.(voir figure jointe)
Merci pour un coup de main :oops:

Pourquoi ne pas simplement utiliser la fonction MAX ?

Code:

---

1 2 3 4 5 6 7 8

x = rand(1,100); y = 1:numel(x);   [xmax,idx] = max(x);   figure plot(x,y,'b-',xmax,y(idx),'r*')

---

Citation:

---

Envoyé par Dut

Pourquoi ne pas simplement utiliser la fonction MAX ?

Code:

---

1 2 3 4 5 6 7 8

x = rand(1,100); y = 1:numel(x);   [xmax,idx] = max(x);   figure plot(x,y,'b-',xmax,y(idx),'r*')

---

---

Non, la fonction max (toute seule) n'est pas du tout appropriée dans le cas où le signal est très bruité. Tout simplement parceque le max de signal n'est pas forcément le centre du pic dans ce cas...
J'avais trouvé le code suivant sur internet qui marchait pas trop mal (même très bien) pour mon application ... on peut régler la tolérance au cas ou le code "saute" des positions.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

function pk = lclmax(x,p2) %LCLMAX [ A.K.Booer 3-Dec-1992 ] % % pk = lclmax(x,p) 1-D peak search down columns % using peak half window width of p. % % Returns PK, a BOOLEAN matrix showing local maxima positions. % p = 2*p2 + 1; % full window width [n,m] = size(x); % useful dimensions z = zeros(p2,m); % pre-allocate result i = toeplitz(p:-1:1,p:n); % shift operator % y = zeros(p,(n-p+1)*m); % temporary matrix y(:) = x(i,:); % index into original data matrix ma = max(y); % find maximum in window % pk = [z ; reshape(ma,n-p+1,m) ; z]; % add missing edge elements pk = pk == x; % find matching elements %%

---

pour ce qui est de l'utilisation, la fonction lclmax renvoie je crois un vecteur de la taille de ton vecteur de données experimentales rempli de 1 et de 0 lorsqu'il y a un 1 c'est que tu as un pic tu n'as donc plus qu'a faire une recherche de 1 pour avoir l'indice.
x= tes données
p2 = paramètre de tolérance à tester

Voilà en espérant que ça puisse t'aider... je suis pas très précis mais ça fait longtemps que je n'ai pas utilisé cette fonction et je ne me rappel plus de tout ... à toi de te replonger dedans

merci cher Dut pour le bon placement :oops::oops: aussi pour ton bout de code c'est simple et ça fonctionne ;)
un merci aussi à azabdur , je vais me "replonger" dans ta proposition:aie:
Merci :oops: