關(guān)于濾波器
使用Design Filer進(jìn)行設(shè)計(jì),可以選擇FIR和IIR兩種优妙。
FIR導(dǎo)出一個(gè)參數(shù)助琐,比如BBPPFF,然后filter(BBPPFF, 1, dataRaw)進(jìn)行濾波舅巷。
IIR到處兩個(gè)參數(shù)羔味,比如SOS和G,[b a]=SOS2tf(SOS,G); filter(b, a, dataRaw)進(jìn)行濾波钠右。
用代碼設(shè)計(jì)
Wn1=[20/250 24/250];%10-12
Wn2=[40/250 44/250];%20-22
nf=4; %階數(shù)
[b1,a1] = butter(nf,Wn1,'bandpass'); %巴特沃斯帶通
[b2,a2] = butter(nf,Wn2,'bandpass');
C3 = filtfilt(b1,a1,C3);%濾波零相位赋元。
spectrogram時(shí)頻圖函數(shù)
[s,f,t] = spectrogram(x,window,noverlap,f,fs)
x是原始信號(hào)。
window是短時(shí)傅里葉變換的窗飒房。
比如原始信號(hào)為1000搁凸,窗為2,堆疊為1狠毯,那么結(jié)果就是999护糖。
比如原始信號(hào)為1000,窗為3嚼松,堆疊為1嫡良,那么結(jié)果就是(1000-3+1)/(3-1)=499。
noverlap是堆疊献酗,必須小于window寝受。
f是
fs是采樣頻率