Hi,
I was looking for the exact same thing (I want to integrate the signal covered by the horizontal yellow segments, to calculate the corresponding percentage of the total energy of the signal contained in each peak). And I find a way to do exactly this (skip the bullet points if you just want the result without the investigations) :
- Open the findpeaks function (using debug mode and ctrl+D with the text cursor somewhere on the function name in my code) and search for "annotate" with ctrl+F. It will quickly lead you to lines 187-190 :
hAxes = plotSignalWithPeaks(x,y,iPk);
if strcmp(annotate,'extents')
plotExtents(hAxes,x,y,iPk,bPk,bxPk,byPk,wxPk,refW);
end
which seems to be where it plots the so called "horizontal yellow segments".
- It seems that the parameter wxPk contains what we want. The function returns wPk as 3-rd output, which is pretty similar. To verify, you can place a breakpoint on the line 189 (the one calling plotExtents() ), and lauch the findpeaks function to get paused on this line. Then copy it in your command window, the figure will appear and the values in wxPk will describe exactly the bounds of each yellow segment.
Once this investigation done, it turns out that all you need to do is the following :
_______________________________
Since you aren't allowed to modify directly this function to add wkPk as an output, just copy its code in a new function file that you can rename (my_findpeaks.m for example). There you just have to modify the header :
function [Ypk,Xpk,Wpk,Ppk] = my_findpeaks(Yin,varargin)
in :
function [Ypk,Xpk,Wpk,Ppk,wxPk] = my_findpeaks(Yin,varargin)
This allows you to get the wanted parameter when calling the function in your code :
[pks,loc,w,p,wxPk] = my_findpeaks(...).
And wxPk(n,:) will give you the two boundaries of the n-th peak.
_______________________________
Done !
Regards.
