finding accurate number of peaks and eleminate unwanted peaks

3 次查看（过去 30 天）
sarmad m
sarmad m 2017-3-20
编辑： Greg Dionne 2017-3-22
HI
I'm trying to find peaks in this data :
the actual data has 10 peaks as shown in the image below (peaks are in rectangle)
. I applied sgf filter got the peaks below
order = 7;
framelen =15;
x=data;
lx = 34;
sgf = sgolayfilt(x,order,framelen);
plot(x);
hold on;
plot(sgf);
sgf=-sgf;
[pks_smoothed,locs_smoothed,widths_smoothed,proms_smoothed] = findpeaks(sgf,'MinPeakProminence',0.05);
pks_smoothed = -pks_smoothed;
plot(locs_smoothed,pks_smoothed,'g*');
text(locs_smoothed+.02,pks_smoothed,num2str((1:numel(pks_smoothed))'))
it detects many unwanted peaks how can I find the required peaks accurately and get rid of unwanted peaks ?

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

回答（1 个）

Greg Dionne
Greg Dionne 2017-3-20
Try:
findpeaks(sgf,'MinPeakProminence',1.2);
  3 个评论
显示 1更早的评论
Adam
Adam 2017-3-21
编辑：Adam 2017-3-21
Well, you can just set the 'NPeaks' property to 10, but whether it picks the 10 you want or not is another matter. You may still have to tweak around with some of the other options for findpeaks.
Greg Dionne
Greg Dionne 2017-3-22
编辑：Greg Dionne 2017-3-22
An alternative is to perhaps take the output of your widths_smoothed and prom_smoothed and come up with some heuristic that matches the most of the peaks you care about then filter out the rest based upon that somehow. Some of your peaks look kind of "flat". Maybe use a combination of medfilt1() and diff() to identify the edges.

请先登录，再进行评论。

类别

在 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