方法1 findpeaks
t=A;
x=B;
y=detrend(x);
plot(t,y);
grid; hold
on
[K,V]=findpeaks(x);
thr=1000;
vdex=find(V>thr);
V1=V(vdex);
K1=K(vdex);
delta=10;if (a(i)<0.5)
a(i)=0;
new(vdex) = 0.5*( y(vdex-delta)+ y(vdex+delta) ) ;
plot(t(K1),V1,'r.','linewidth',5);
title(['信号提取极值,阈值thr=' num2str(thr)]);
方法2 diff
for i=1:3000if (a(i)<0.5)a(i)=0;endend
b=find(diff(sign(diff(a)))<0)+1;% 求解极值点序号plot(t, a,t(b), a(b),'r*') % 画出去噪声后的波形，由于滤波的关系，前一秒钟的信号异常，因此只画1s以后的数据xlabel('t(s)')ylabel('Voltage(v)')axis([1 6 0 2]);title('Final Signal')grid on
方法3 localMaximum
a=mydata(:,2);a=detrend(a);
peaks = localMaximum(a,37000);