（MATLAB）一维信号峰值检测findpeaks

by HPC_ZY

一维信号波峰波谷检测，峰值、峰宽、到达时间计算。MATLAB自带的findpeaks函数能完成这一任务，在这里分享其基本用法与特殊用法。

基本介绍

关于波峰，峰值，峰宽

findpeaks函数用法

基本用法

1. 无参数，简单输出

x = 0:0.1:2*pi;
y0 = sin(x);

% 这里的locs对应y0索引
[pks,locs] = findpeaks(y0);
plot(y0, '.-'),hold on
plot(locs, pks, 'ro'),hold off

2. 指定x坐标

x = 0:0.1:2*pi;
y0 = sin(x);

% 这里的locs对应实际x坐标
[pks, locs] = findpeaks(y0, x);
plot(x, y0, '.-'),hold on
plot(locs, pks, 'ro'),hold off

3. 绘制结果
   如果不设置输出，就会自动绘制结果

x = 0:0.1:2*pi;
y0 = sin(x);

findpeaks(y0, x);

4. 完整输出+完整绘制

x = 0:0.1:2*pi;
y0 = sin(2*x);

% w为半(凸/峰)宽，p为凸高
[pks, locs, w, p] = findpeaks(y0, x);
% 加入如下参数即可绘制所有结果
findpeaks(y0, x, 'Annotate','extents');

特殊用法

为了简便，后续例子不再输出结果，直接绘制。

1. 最小峰值高度

x = 0:0.1:2*pi;
y0 = sin(2*x)+x/2;

MPH = 2;
% 仅搜索高度MPH以上的峰
findpeaks(y0, 'MinPeakHeight', MPH, 'Annotate', 'extents'); 

2. 指定最小峰凸高度

x = 0:0.1:2*pi;
y0 = sin(3*x)+((x-pi)/2).^2;

MPP = 1;
% 仅搜索高度MPH以上的峰
findpeaks(y0, 'MinPeakProminence', MPP, 'Annotate', 'extents'); 

3. 指定相邻最小差
   指的是波峰点需要比左右相邻点高出一定值

x = 0:0.1:2*pi;
y0 = sin(x).^20;

TH = 0.05;
% 仅搜索从两侧峰值垂直下降超过MPP的峰
findpeaks(y0, 'Threshold', TH, 'Annotate', 'extents'); 

可以看出第一个峰的左邻点与峰值差太小

4. 指定最小半宽

x = 0:0.1:2*pi;
y0 = sin(x.^2);

MINW = 5; 
% 搜索半宽大于MINW的峰
findpeaks(y0, 'MinPeakWidth', MINW, 'Annotate', 'extents'); 

5. 指定最大半宽

x = 0:0.1:2*pi;
y0 = sin(x.^2);

MAXW= 5; 
% 搜索半宽小于MAXW的峰
findpeaks(y0, 'MaxPeakWidth', MAXW, 'Annotate', 'extents'); 

其他

1. 后续还有，可能更新。
2. 求峰谷只需要对信号取负，计算结果再取负就可以了。