Déterminer les largeurs de pics

Ouvrir le live script

Créez un signal constitué d’une somme de courbes de Gauss. Définissez l’emplacement, la hauteur et la largeur de chaque courbe.

x = linspace(0,1,1000);

Pos = [1 2 3 5 7 8]'/10;
Hgt = [7 6 3 2 2 3]';
Wdt = [3 8 4 3 4 6]'/100;

y = sum(Hgt.*(exp(-((x-Pos)./Wdt).^2)),1);

Mesurez la largeur des pics en utilisant la moitié de la proéminence et la moitié de la hauteur comme référence.

tiledlayout("flow")
nexttile
findpeaks(y,x,Annotate="extents")
title("Half-Prominence Peak Widths")
nexttile
findpeaks(y,x,Annotate="extents",WidthReference="halfheight")
title("Half-Height Peak Widths")

Figure contains 2 axes objects. Axes object 1 with title Half-Prominence Peak Widths contains 4 objects of type line. One or more of the lines displays its values using only markers These objects represent signal, peak, prominence, width (half-prominence). Axes object 2 with title Half-Height Peak Widths contains 6 objects of type line. One or more of the lines displays its values using only markers These objects represent signal, peak, height, width (half-height), border.

Sélectionnez les pics les plus élevés séparés par un écart d’au moins 0,5 unité sur l’axe des x. Mesurez la largeur des pics en utilisant la moitié de la proéminence et la moitié de la hauteur comme référence.

figure
tiledlayout("flow")
nexttile
findpeaks(y,x,MinPeakDistance=0.5,Annotate="extents")
title("Half-Prominence Peak Widths")
nexttile
findpeaks(y,x,MinPeakDistance=0.5,Annotate="extents", ...
    WidthReference="halfheight")
title("Half-Height Peak Widths")

Figure contains 2 axes objects. Axes object 1 with title Half-Prominence Peak Widths contains 4 objects of type line. One or more of the lines displays its values using only markers These objects represent signal, peak, prominence, width (half-prominence). Axes object 2 with title Half-Height Peak Widths contains 6 objects of type line. One or more of the lines displays its values using only markers These objects represent signal, peak, height, width (half-height), border.

Seuls les premier et dernier pics remplissent la condition d’écart minimal. Les largeurs affichées sur le tracé correspondent donc à ces deux pics. Comme l’étendue de chaque pic reste inchangée, la largeur de pic conserve la même valeur indépendamment des conditions spécifiées et du fait que le pic soit sélectionné ou non.