带通滤波算法 python

1、打开matlab，在命令行窗口输入fdatool，回车即可。

2、以FIR窗函数带通滤波器为例。

• 滤波器类型中有lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的FIR滤波器,本文选择带通和FIR,函数设置为窗函数。其他窗函数可以自己选择。
• design method（设计方法）选项，包括iir滤波器的butterworth（巴特沃思）法、chebyshev typei（切比雪夫i型）法、chebyshev type ii(切比雪夫i型）法、elliptic（椭圆滤波器）法和fir 滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法。
• filter order(滤波器阶数）选项，定义滤波器的阶数，包括 specify order（指定阶数）和minimum
  order（最小阶数)。在specify order中填入所要设计的滤波器的阶数（n阶滤波器，specify
  order=n-1)，如果选择minimum order则matlab根据所选择的滤波器类型自动使用最小阶数。例如：50阶应该设置为49
  
  3、frenquency specifications选项设置。
• frenquency specifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filtertype选项和 design method选项决定，例如bandpass（带通）滤波器需要定义fstop1（下阻带截止频率）、fpass1(通带下限截止频率）、fpass2（通带上限截止频率）、fstop2(上阻带截止频率），而lowpass（低通）滤波器只需要定义fstop1、fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。
• magnitude specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义wstop1（频率fstop1处的幅值衰减）、wpass（通带范围内的幅值衰减）、wstop2（频率fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。
  
  4、生成脚本代码函数。
  设计完滤波器后，一定要点击下方的Design filter按钮，然后点击file里面生成。
  
  有三种类型，可生成.m脚本文件，可直接调用，具体需要哪一种自己测试。