数字滤波器,带通滤波器设计

巴特沃思（Butterworth）模拟／数字滤波器带通滤波器设计

格式：[b,a] = butter(n,Wn)

[b,a] = butter(n,Wn,'ftype')

[b,a] = butter(n,Wn,'s')

[b,a] = butter(n,Wn,'ftype','s')

[z,p,k] = butter(...)

[A,B,C,D] = butter(...)

说明：butter函数可以设计模拟或数字的低通、高通、带通和带阻Butterworth滤波器。Butterworth滤波器可以使通带内的幅频响应最大程度地平坦，但这也使得它的通带到阻带的过渡过程较慢。在这方面Chebyshev滤波器和椭圆滤波器性能较好。

在设计数字滤波器时，butter函数中的参数Wn与besself函数有很大的区别，它是一个相对量，其定义区间为Wn∈[0,1]，其中1对应于0.5fs，fs为采样频率（单位Hz）；在设计模拟滤波器时，Wn采用真实频率，单位为Hz。

[b,a] = butter(n,Wn)返回截止频率为Wn（单位为弧度/秒）的n阶Butterworth

数字低通滤波器，b、a分别为滤波器传递函数的分子和分母系数向量（降幂排列）。

当Wn为二元向量，即Wn = [W1 W2] (W1

举例：信号采样率fs为10000Hz，带通为50—600Hz

WW=fs/2;

Wn=[50/WW 600/WW]; 通带为50Hz--600Hz

[b,a] = butter(3,Wn);

Original_Signal=filter(b,a,Original_Signal);