Discussion :

[miss_angel]

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 :

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

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

[Jerome Briot]

Pourquoi ne pas simplement utiliser la fonction MAX ?

Code :

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

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*')

[azabdur]

Pourquoi ne pas simplement utiliser la fonction MAX ?

Code :

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

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 :

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

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

[miss_angel]

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