我喜欢设计自己的过滤器,那么试试这个:
Fs = 44100;%采样频率(赫兹)
Fn = f / 2;%奈奎斯特频率(赫兹)
Wp = (5 [5] + 2.2 e + 3) / Fn;%阻带频率(正常)
Ws = (-10 [10] + 2.2 e + 3) / Fn;%通带频率(正常)
Rp = 1;%通带波纹(dB)
Rs = 50;%阻带波动(dB)
[n, Ws] = cheb2ord (Wp、Ws Rp, Rs);%过滤器订单
(z,磷、钾)= cheby2 (Ws n, Rs,“停止”);%滤波器设计,Sepcify Bandstop
(sosn, gn) = zp2sos (z,磷、钾);%转换为稳定Second-Order-Section
图(3)
freqz (Fs sosn 2 ^ 16日)%滤波德图
filtered_signal = filtfilt (gn, sosn original_signal);%过滤信号