Discussion :

[fatenov]

Salut,

Je cherche à écrire un bout de code qui me detecte les pics d'un séquencement donné, j'ai pensé à fixer un seuil et regarder les montés et les descentes par ce seuil, mais cette méthode ne permettra de détecter les pics en dessous du seuil (même si leur valeur est faible, ils restent toujours des pics) et aussi si le sommet d'un pic, lui même, est divisé en deux pics successif, il ne sera compté que pour un seul.

Enfin, cette méthode ne donne pas de bon résultats, et par suite je cherche de l'aide pour détecter efficacement les pics.

Merci pour toute proposition et à bientôt.

[WooD*]

tout dépend de comment tu définis tes "pics", mais si tu donne comme valeur seuil la valeur moyenne de ton signal à la fonction que tu imagine tu devrais trouver ton bonheur nan? peut-être qu'un exemple concret(une image ) pourrait aider à mieux comprendre ce que tu recherche

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
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);
 
figure
plot(x,y,'b-',x(idx),y(idx),'r*')

[fatenov]

J'ai pensé dès le début à mettre cette image jointe, mais dsl pour sa qualité du graphe

Voila, les pic detecté sont: 1, 2 et 3.
Mais, en oubliant le seuil, je veux aussi detecter les pics 6, 7 et 8.
De même, le pic 3 contient lui même deux pics qui doivent être distinqués.

J'espère que c'est plus clair maintenant, merci d'essayer.

[fatenov]

Merci bien Dut, ça me donne exactement ce que je cherchai.
Mais j'aimerai que tu m'expliques un peu ces quelques lignes, le principe (mathématique) de cette detection. Car ça me semble intéressant.

[FR119492]

Salut !

Il faudrait savoir ce qui est un pic et ce qui n'en est pas un. Les deux bouts de programmes suivants représentent la même fonction avec deux résolutions différentes. Le nombre de pics est-il le même?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
n=3;
x=1e6;
m=100;
d=1e-7*x;
r=x*ones(1,n);
p=poly(r);
xmin=x-m*d;
xmax=x+m*d;
xtab=linspace(xmin,xmax,2*m+1);
ztab=xtab-x;
ytab=polyval(p,xtab);
plot(ztab,ytab)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
n=3;
x=1e6;
m=200;
d=.5e-7*x;
r=x*ones(1,n);
p=poly(r);
xmin=x-m*d;
xmax=x+m*d;
xtab=linspace(xmin,xmax,2*m+1);
ztab=xtab-x;
ytab=polyval(p,xtab);
plot(ztab,ytab)

Un maximum, ce que tu appelles un pic, correspond à un changement de signe de la dérivée. Or le calcul numérique de la dérivée d'un signal bruité n'a souvent aucun sens.

Jean-Marc Blanc

[Jerome Briot]

Un pic est un maximum donc la valeur qui le précède et celle qui le suit sont inférieures.

Il suffit donc de faire la différence des valeurs dans un sens puis dans l'autre et de prendre les valeurs pour lesquels le signe de cette différence vaut +1 dans les deux cas.

On ajoute le -inf pour identifier directement le bon indice. Car DIFF retourne un vecteur de taille n-1

Par exemple, pour un vecteur [1 2 4 3], le maximum vaut 4 est se trouve à la 3ème place :

Code :

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

1
2
3
4
5
6
7
8
y=[1 2 4 3]
t=diff([-inf y]);
sa=sign(t)
 
t=diff([-inf y(end:-1:1)]);
 
sb=sign(t);
sb=sb(end:-1:1)

y =

1 2 4 3


sa =

1 1 1 -1


sb =

-1 -1 1 1


idx =

0 0 1 0

[fatenov]

Merci Dut, ton idée est bien claire maintenant.

[Jerome Briot]

Une autre méthode (quasi similaire) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
y=rand(1,30);
 
figure
plot(y)
 
s=sign(diff([-inf y]));
 
idxmax=strfind(s,[1 -1]);
idxmin=strfind(s,[-1 1]);
 
hold on 
 
plot(idxmax,y(idxmax),'ro',idxmin,y(idxmin),'go')

[fatenov]

Salut Dut,


Oui ça donne le même résultat et c'est bien le même principe.

Dans le 1er code,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
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);
 
figure
plot(x,y,'b-',x(idx),y(idx),'r*')

On peut détecer les creus par cette ligne de code:

Code :

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

idx_creu=(sa==-1 & sb==-1);

Ce code détecte bien les maximas au début et à la fin, mais pas les creux , je pense faire :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
y=rand(1,30);
 
figure
plot(y)
 
s=sign(diff([+inf y +inf]));
 
idxmax=strfind(s,[1 -1]);
idxmin=strfind(s,[-1 1]);
 
hold on 
 
plot(idxmax,y(idxmax),'ro',idxmin,y(idxmin),'go')

Il détecte bien les creux au début et fin.
Mais est ce qu'on ne pourra pas avoir quelque chose qui détecte les deux à la fois?

Merci et à bientôt.

[Jerome Briot]

Mais est ce qu'on ne pourra pas avoir quelque chose qui detecte les deux à la fois?

Comme ceci ?

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,30);
 
figure
plot(y)
 
s=sign(diff([+inf y +inf]));
 
idx=(s(2:end).*s(1:end-1))==-1;
 
hold on 
 
plot(find(idx),y(idx),'ro')

[fatenov]

Il ne détecte pas les pics (maximas) aux bornes.

[Jerome Briot]

Je ne suis pas sûr de bien comprendre ce que tu cherches, mais de toute façon, la première et la dernière valeur sont des valeurs particulières pour lesquelles il faudra faire un test à part avec respectivement, la valeur suivante et la valeur précédente, non ?

[fatenov]

Oui c'est bien le cas, d'ailleur, on peut les négliger si on s'intéresse à des signaus de taille très grande, en éliminant les bornes.
Mais j'ai voulu savoir s'il y a une méthode pour tout détecter qui nous permet d'être comme ça:

Merci et à bientôt.

[zaiani]

bonjour,
pourquoi ne pas utilisé la fonction findpeaks

Code :

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

[peak locs]= findpeaks(y)

dont peak donne les maximums des pics et locs les indices des pics
tu peux essayé ça pour un vecteur et ça marche.

[Jerome Briot]

findpeaks

C'est bien entendu la fonction à utiliser en priorité... quand on a la chance de posséder la Signal Processing Toolbox

[acx01b]

salut

la référence dans le domaine ce sont les articles de Karin Dressler pour le mirex

celui ci par exemple : http://www.music-ir.org/evaluation/m...y/dressler.pdf

[chester2518]

thanks for your help DUT !!! it is an amazing simple IDEa