关于滤波器

使用Design Filer进行设计，可以选择FIR和IIR两种。
FIR导出一个参数，比如BBPPFF，然后filter(BBPPFF, 1, dataRaw)进行滤波。
IIR到处两个参数，比如SOS和G，[b a]=SOS2tf(SOS,G); filter(b, a, dataRaw)进行滤波。

用代码设计

Wn1=[20/250 24/250];%10-12
Wn2=[40/250 44/250];%20-22
nf=4; %阶数
[b1,a1] = butter(nf,Wn1,'bandpass'); %巴特沃斯带通
[b2,a2] = butter(nf,Wn2,'bandpass');
C3 = filtfilt(b1,a1,C3);%滤波零相位。

spectrogram时频图函数

[s,f,t] = spectrogram(x,window,noverlap,f,fs)
x是原始信号。
window是短时傅里叶变换的窗。
比如原始信号为1000，窗为2，堆叠为1，那么结果就是999。
比如原始信号为1000，窗为3，堆叠为1，那么结果就是（1000-3+1）/（3-1）=499。
noverlap是堆叠，必须小于window。
f是
fs是采样频率