转自:https://blog.csdn.net/it_beecoder/article/details/78681517
原波形 为250HZ的正弦波型,findpeaks 寻找峰值函数
pks = findpeaks(data)
[pks,locs] = findpeaks(data) ------pks 对应峰值,locs 对应峰值位数
[...] = findpeaks(data,'minpeakheight',mph)----mph 设定峰值的最小高度
[...] = findpeaks(data,'minpeakdistance',mpd)----mpd 设定两峰值间的最小间隔数
原始波形
波峰查找和描绘
>> [maxv,maxl]=findpeaks(m,'minpeakdistance',100); %maxv峰峰值点 maxl:峰峰值点对应的位置 最小间隔=100
>> hold on
>> plot(t,m) %绘制原波形
>> plot(maxl,maxv,'*','color','R'); %绘制最大值点
通过增加 最小间隔条件 可以过滤掉单周期内的干扰极值 mpb的值 建议取 周期/1.5 具体看自己的波形做调整
波谷查找和描绘
由于findpeaks不带直接查找波谷,因此需要用些方法来查找,网上的diff方法,不适用于单周期内有干扰极值的情况。
因此我还是采用了findpeaks
首先 将原波形数据 “取反” m=原始数据(数组类型)
nm=3000-m ; % 3000是原始数据中的最大值限度
取反后的波形如图
然后再使用findpeaks 找出取反后的波形中的峰峰值的位置即可
程序:
nm=3000-m(t);% 数组整体减
[minv,minl]=findpeaks(nm,'minpeakdistance',100); %maxv峰峰值点 maxl:峰峰值点对应的位置
hold on
plot(t,m)
plot(minl,m(minl),'*','color','G'); %绘制最大值点
实现的波谷波形
最后的波峰波谷波形
联系客服