Sort of peak analysis

1 次查看（过去 30 天）
Diego R
Diego R 2019-9-6
评论： darova 2019-9-8
I'd like to analyze this graph (is a vector) to have an answer like this:
for 1 to number of peaks:
PEAK_1=[start of peak_1, finish of peak_1, maximum, area under peak_1]
PEAK_2=[start of peak_2, finish of peak_2, maximum, area under peak_2]
...
...
PEAK_N=[start of peak_N, finish of peak_N, maximum, area under peak_N]
I am planning to do a "for" loop checking each value and its predecesor to determine when does a peak start and finish and a while to take the actions. My question is, is there any easier way to determine when does each peak star and finish?
Also, if anybody has any imaginative way better than mine it will be welcome.
  5 个评论
显示 3更早的评论
Image Analyst
Image Analyst 2019-9-8
Diego, are you not seeing all the answers below, in the official "Answers" section rather than up here in the Comments section? If you need more help, then attach your signal in a .mat file with the paper clip icon.

请先登录，再进行评论。

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

采纳的回答

darova
darova 2019-9-7
How i found start and end of each peak
% generate some data
x = linspace(-1,17);
y = sin(x);
y = y.*(y>0);
thresh = 0.1;
ind = find(abs(y) < thresh); % find 'zeros' points
i = find( diff(ind)>1 ); % find breaks
i = [i; i+1]; % add 'end' points for each peak
i = reshape(i,2,[])';
ind = ind(i); % two columns: start - end points
plot(x,y,'.-b')
hold on
plot(x(ind),y(ind),'or')
plot([min(x) max(x)],[1 1]*thresh)
Use finpeaks() to get max value for each peak
Use trapz() to calculate area

更多回答（2 个）

Image Analyst
Image Analyst 2019-9-7
If you have the Image Processing Toolbox, you could use regionprops
props = regionprops(signal < 0, 'PixelIdxList');
It will give you, for each peak, the starting and stopping point. You could also ask for MeanIntensity and Area and multiply those together to get the area under the curve. Let me know if you can't figure it out, and attach your signal in a .mat file.
Fabio Freschi
Fabio Freschi 2019-9-7
If the signal is not too noisy you can work with ‘diff’ and look for the change of sign

类别

在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by