La traduction de cette page n'est pas à jour. Cliquez ici pour voir la dernière version en anglais.

Déterminer les largeurs de pics

Créez un signal constitué d’une somme de courbes en cloche. 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 = [4 4 2 2 2 3];
Wdt = [3 8 4 3 4 6]/100;

for n = 1:length(Pos)
    Gauss(n,:) =  Hgt(n)*exp(-((x - Pos(n))/Wdt(n)).^2);
end

PeakSig = sum(Gauss);

Tracez chacune des courbes ainsi que leur somme.

plot(x,Gauss,'--',x,PeakSig)
grid

Figure contains an axes object. The axes object contains 7 objects of type line.

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

findpeaks(PeakSig,x,'Annotate','extents')

Figure contains an axes object. The axes object 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).

Mesurez à nouveau les largeurs, cette fois en utilisant la moitié de la hauteur comme référence.

findpeaks(PeakSig,x,'Annotate','extents','WidthReference','halfheight')
title('Signal Peak Widths')

Figure contains an axes object. The axes object with title Signal 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.