滤波器的设计方法概要.ppt

2. 三角形（Bartlett）窗 w(n)的傅里叶变换为 （7-88） (7-89) 近似结果在N1 时成立。 此时，主瓣宽度为8π/N, 比矩形窗主瓣宽度增加一倍， 但旁瓣却小很多。 3. 汉宁（Hanning）窗 汉宁窗又称升余弦窗。 (7-90) 利用傅里叶变换特性，可得 (7-91) 当N1 时，N-1≈N, 所以窗函数的幅频函数为 （7-92） 这三部分之和，使旁瓣互相抵消，能量更集中在主瓣，它的最大旁瓣值比主瓣值约低31dB。但是代价是主瓣宽度比矩形窗的主瓣宽度增加一倍， 即为 8π/N。 4. 海明（Hamming）窗 海明窗又称改进的升余弦窗。 把升余弦窗加以改进， 可以得到旁瓣更小的效果， 窗形式为 (7-93) w(n)的频率响应的幅度特性为 (7-94) 与汉宁窗相比，主瓣宽度相同，为 8π/N，但旁瓣又被进一步压低， 结果可将99.963%的能量集中在窗谱的主瓣内，它的最大旁瓣值比主瓣值约低41dB。 5. 布拉克曼（Blackman）窗 布拉克曼窗又称二阶升余弦窗。 为了进一步抑制旁瓣，对升余弦窗函数再加上一个二次谐波的余弦分量， 变成布拉克曼窗，故又称二阶升余弦窗。 (7-95) w(n)的频率响应的幅度特性为 (7-96) 图 7-21 五种常用的窗函数 图 7-22 图 7-21的各种窗函数的傅里叶变换（N=51），A=20 lg|W(ω)/W(θ)| (a) 矩形窗; (b) 巴特利特窗（三角形窗）; (c) 汉宁窗; (d) 海明窗; (e) 布拉克曼窗 图 7-23 理想低通滤波器加窗后的幅度响应（N=51）, A=20lg|H(ω)/H(0)| (a) 矩形窗; (b) 巴特利特窗（三角形窗）; (c) 汉宁窗; (d) 海明窗; (e) 布拉克曼窗 6. 凯塞（Kaiser）窗 这是一种适应性较强的窗，其窗函数的表示式为 0≤n≤N-1 （7-97） 式中，I0(x)是第一类变形零阶贝塞尔函数，β是一个可自由选择的参数。 ? 零阶贝塞尔函数的曲线如图7-24所示。 图 7-24 零阶贝塞尔函数 图7-25 凯塞窗函数 表7-4 凯塞窗的性能 虽然凯塞窗看上去没有初等函数的解析表达式。但是，在设计凯塞窗时，对零阶变形贝塞尔函数可采用无穷级数来表达 （7-98） 这个无穷级数，可用有限项级数去近似，项数多少由要求的精度来确定。 因而采用计算机是很容易求解的。 表7-5 六种窗函数基本参数的比较 把变换式（7-41）和变换式（7-42）结合起来，可得到直接从模拟低通原型变换成数字带阻滤波器的表达式，也就是直接联系s与z之间的变换公式 经推导后得 （7-43） 式中： 根据双线性变换，模拟带阻频率与数字带阻频率之间的关系仍为 定义: （7-44） （7-45） （7-46） 式中，Ω0 为带阻滤波器阻带的几何对称中心角频率; B为带阻滤波器的阻带宽度，它与低通原型中的截止频率Ωc成反比。  设数字带阻的中心频率为ω0，数字带阻滤波器的上、下边带的截止频率分别为ω2和ω1， 则将式（7-42）代入式（7-43）、 式（7-43），可得: （7-47） （7-48） 考虑到模拟带阻到数字带阻是阻带中心频率相对应的映射关系，则有 （7-49） 将式（7-44）、 式（7-45）和式（7-46）代入式（7-40）及式（7-41），并应用一些标准三角恒等式运算后可得: 所以，在设计时，要给定中心频率和带宽或者是中心频率和边带频率， 利用式（7-48）和式（7-49）来确定D1和E1两常数， 然后利用式（7-40）的变换，把模拟低通系统函数一步变成数字带阻系统函数 （7-50） 式中，Ha(s)为模拟低通原型传递函数。  可以看出，数字带阻滤波器的极点数（或阶数）将是模拟低通滤波器极点数的两倍。 下面来讨论模拟低通滤波器与数字带阻滤波器频率之间的关系。令s=jΩ，z=ejω代入式（7-40），经推导后可得 (7-51) 其变换关系曲线如图7-18所示，其映射关系为: Ω=0→ω=0, ω=π Ω=±∞→ω=ω0 也就是说，低通滤波器的通带（Ω=0 附近）映射到带阻滤波器的阻带范围之外（ω=0, π），低通滤波器的阻带（Ω=±∞）映射到带阻滤波器的阻带上（ω=ω0 附近）。 如图7-18 变换关系曲线 表 7-3 利用双线性变换法从截止频率为Ωc