matlab 双峰,如何在MATLAB中区分双峰和单峰阵列？

我发现

PEAKDET的功能是相当可靠和快速的,虽然它是循环的.它不需要对噪声数据进行预平滑,而是发现差异大于参数delta的局部最大和最小极值.

由于PEAKDET从左到右运行,有时会错过正确的站点上的高峰.为了避免它,我更喜欢运行它两次：

%# some data

n = 100;

x = linspace(0,3*pi,n);

y = sin(x) + rand(1,n)/5;

%# run peakdet twice left-to-right and right-to-left

delta = 0.5;

[ymaxtab, ymintab] = peakdet(y, delta, x);

[ymaxtab2, ymintab2] = peakdet(y(end:-1:1), delta, x(end:-1:1));

ymaxtab = unique([ymaxtab; ymaxtab2],'rows');

ymintab = unique([ymintab; ymintab2],'rows');

%# plot the curve and show extreme points based on number of peaks

plot(x,y)

hold on

if size(ymaxtab,1) == 2 && size(ymintab,1) == 1 %# if double peak

plot(ymintab(:,1),ymintab(:,2),'r.','markersize',30)

elseif size(ymaxtab,1) == 1 && size(ymintab,1) == 0 %# if single peak

plot(ymaxtab(:,1),ymaxtab(:,2),'r.','markersize',30)

else %# if more (or less)

plot(ymintab(:,1),ymintab(:,2),'r.','markersize',30)

plot(ymaxtab(:,1),ymaxtab(:,2),'r.','markersize',30)

end

hold off