is there a way to obtain peak width after peak detection ?

51 views (last 30 days)
I have used the function findpeaks to detect spikes in a signal
however I have initially used the halfprominence as width-reference
is there a way to obtain widths using half-height after such peak detection, having the peak timestamps and their amplitude values ?

Accepted Answer

Star Strider
Star Strider on 28 May 2019
Not a stupid question at all.
You simply need to ask findpeaks to return them as outputs:
[pks,locs,widths,proms] = findpeaks(PeakSig,x);
See the documentation section on Peak Prominences (link) for an illustration and discussion. The ‘widths’ are defined as the full width at half of the maximum amplitude (FWHM).
  2 Comments
Star Strider
Star Strider on 28 May 2019
I am not certain what you are asking.
The ‘height’ (and related ‘width’) and ‘prominence’ values are different, and only loosely related, at least as I read the documentation. The Prominence (link) property seems difficult to define. I would just use the FWHM value to characterise the peaks.
In other applications, I have used both ‘widths’ and ‘proms’ to select particular peaks, defining a desired peak by ‘(width > wv) & (prom > pv)’, where ‘wv’ and ‘pv’ were values determined elsewhere in the code. Perhaps you can combine them that way.

Sign in to comment.

More Answers (2)

Adam Danz
Adam Danz on 28 May 2019
Edited: Adam Danz on 28 May 2019
Use the "halfheight" width reference. There's an example here: https://www.mathworks.com/help/signal/ref/findpeaks.html#buhd6xj
[pks,locs,w,p] = findpeaks(PeakSig,x,'Annotate','extents','WidthReference','halfheight')
  1 Comment
Adam Danz
Adam Danz on 28 May 2019
LO's comment moved here:
"I know that, thanks. What I wanted to know is whether one could get those values without doing the peak detection again"
Why would you need to run it twice? The widths are returned in the 3rd output.
[pks,locs,w,p] = findpeaks(PeakSig,x,'Annotate','extents','WidthReference','halfheight')

Sign in to comment.


LO
LO on 28 May 2019
I know that, thanks.
what I wanted to know is whether one could get those values without doing the peak detection again
but I guess there is no way. stupid question most likely

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!