Discussion :

[jena]

Le rôle de la fonction "findpeaks" est de trouver des crêtes avec l'interpolation quadratique facultative [ K, v]=(x, m, w)

Entrées: X est le signal d'entrée
M est mode:
' q ' exécute l'interpolation quadratique
' v ' trouve des vallées au lieu des crêtes
' W ' est la tolérance de largeur; on éliminera une crête s' il y a une crête plus élevée dans des échantillons de +-w

Sorties: K sont les locations maximaux dans X

Code :

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

1
2
[k,v]=findpeaks(x,'q',60);
figure(1);plot((fe*(1:nfft-1)/nfft),x), hold on , plot(k,v,'r*'),hold off

Dans mon cas; le valeur k n’est pas sur le pic de la figure 1. car on a plot (abscisse,x) alors que dans la fonction findpeaks on a que x donc il n’existe pas de correspondance entre les valeurs des crêtes détectés par la fonction findpeaks et les pics donnés par la figura 1.
Ma question existe t’il une solution.
Merci

[Jerome Briot]

As-tu déjà lu cette discussion Détection de pics ?

[jena]

Salut; j'ai bien créer les croix sur les pics dans ma figure à l'aide de votre fonction "détection pic".
Mais, je veux voir sur ma figure des croisx que sur les pics qui ont une amplitude >0.6

Code :

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

1
2
3
4
5
 
if y(idx)>0.6
figure
plot(x,y,'b-',x(idx),y(idx),'r*')
end

ce bout de code avec cette condition n' a pas donné de solution:
J'espère que je suis clair
Merci A+

[Jerome Briot]

voir sur ma figure des croisx que sur les pics qui ont une amplitude >0.6

Il suffit d'ajouter une condition supplémentaire :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
y = rand(1,15);
x = 1:15;

sa = sign(diff([-inf y]));
sb = sign(diff([-inf y(end:-1:1)]));
sb = sb(end:-1:1);

idx = (sa==1 & sb==1) & y>0.6;

figure
plot(x,y,'b-',x(idx),y(idx),'r*')

[jena]

Salut;
Maintenant je veux que le croix ne s'affiche dans la figure que sur le premier pic.
donc qu'est ce que je doit ajouter comme condition ou bien dans le plot
Merci

[Jerome Briot]

Il suffit de ne prendre que le premier élément contenu dans le vecteur idx.

Code :

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

idx = find((sa==1 & sb==1) & y>0.6,1,'first');