getting the area under the peaks

Question

Rashid Hussein 2017-11-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/367783-getting-the-area-under-the-peaks

评论： Star Strider 2017-11-19

I have a graph that contain many peaks , but I want to get the sum( area under peak) of counts (y-axis) in every peaks as shown in the image that I uploaded , thanks in advanced

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2017-11-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/367783-getting-the-area-under-the-peaks#answer_291838

在 MATLAB Online 中打开

I would first identify the low points (‘valleys’) by inverting your signal and then using the Signal Processing Toolbox findpeaks function to locate them. Then integrate between those points to identify the peaks.

Example —

t = linspace(0, 10*pi, 250);                        % Create Data
s = sin(t) + 0.2*cos(5*t) + 2;                      % Create Data
figure(1)
plot(t, -s)
grid
[Vlys, Idx] = findpeaks(-s, 'MinPeakHeight',-1.1);
Area = cumtrapz(t,s);
idxvct = [1 Idx length(s)];
for k1 = 1:length(idxvct)-1
    PkAreas(k1) = Area(idxvct(k1+1)) - Area(idxvct(k1));
end

The ‘PkAreas’ vector will have the areas of the peaks.

You will have to adapt this to your data. Be sure to plot the negative of your data first so you can see and set the appropriate arguments for findpeaks.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Star Strider 2017-11-19

My pleasure.

My code creates the cumulative integral, identifies the valleys, then uses those indices to separate and calculate the areas for each peak in the ‘PkAreas’ vector.

My code worked correctly with the data I created to test it, because I do not have your data to test it with. You have included my code, although you commented it out.

请先登录，再进行评论。