關(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是采樣頻率