数字滤波器设计浙江大学光电科学与工程学院.ppt

73 第六讲 数字滤波器设计 一个窗函数的性能主要取决于它的两个参数，即主瓣宽度和相对旁瓣级。 相对旁瓣级是最大旁瓣与主瓣以dB为单位的幅度差值。 最大通带偏移和最小阻带值之间的距离近似等于窗的主瓣宽度，过渡带宽小于主瓣宽度。 为了保证从通带快速过渡到阻带，窗函数应该有一个非常小的主瓣宽度，另一方面，为了减小通带和阻带波纹，旁瓣下的面积也要求非常小。遗憾的是，这两个要求是相互矛盾的。 矩形窗、Hann窗、Hamming窗的波纹与滤波器长度和截止频率无关，为一常数。 过渡带宽近似为： c为一常数 Rectangular Window Hamming Window ht(n) Ht(ejw) Hanning Window Blackman Window ht(n) Ht(ejw) 3.32?/M 54.5dB 42.7dB 8?/(2M+1) Hamming 5.56?/M 75.3dB 58.1dB 12?/(2M+1) Blackman 3.11?/M 43.9dB 31.5dB 8?/(2M+1) Hann 0.92?/M 20.9dB 13.3dB 4?/(2M+1) Rectangular Transition bandwidth Minimum stopband attenuation Relative sidelobe level Main lobe width Type of windows 固定窗函数的特性 固定窗函数设计的滤波器的波纹是固定的 一些可变窗函数提供了额外的参数以控制波纹 多尔夫-切比雪夫窗 凯泽窗——使用最广泛 Tk(x)为k阶切比雪夫多项式 I(x)为修正0阶贝塞尔函数 7.7.5 可调窗函数 另一种消除Gibbs现象的方法是修正数字滤波器的频率响应的指标，从而使该滤波器在通带和阻带之间出现一个过渡带，使各频带之间平滑过渡 0 -? ? ?p ?s -?p -?s 7.7.6 具有平滑过渡带的FIR滤波器的冲激响应 h[n]的N点DFT可由其DTFT进行等间隔的N点抽样得到? 数字滤波器的冲激响应序列h[n]可以利用它的频率样本上的IDFT来计算——频率抽样法 7.8 频率抽样法 由于用来使计算机设计的滤波器与要求的频率响应之间的误差最小化的迭代最优化技术的提出，出现了很多基于该最优化技术的滤波器设计方法 需设计的滤波器传输函数 期望的频率响应 某种准则下的逼近 7.9 数字滤波器的计算机辅助设计 在商业上有很多可供选择的软件包，它们使得数字滤波器的计算机设计相当简单——FDA Tools@Matlab 设计一个低通IIR数字滤波器G(z)，性能指标如下: 7.3 低通IIR滤波器设计 将数字截止频率预畸变为模拟截止频率 k：过渡比 k1：分辨率 巴特沃兹逼近的阶数估计： 频率解归一化： 三阶归一化低通Butterworth传输函数: 第一种方法分为以下几个步骤： 7.4 高通、带通和带阻IIR数字滤波器设计 Step 1：用式 预畸变所求数字滤波器GD(z)的数字频率指标，从而得到一个等价的模拟滤波器HD(s) 的频率指标。 Step 2：选取一种合适的频率变换，将HD(s)的频率指标转换成原型低通滤波器HLP(s)的频率指标。（附录B） Step3:设计模拟低通滤波器HLP(s).（附录A） Step4: 用步骤2中频率变换的逆变换将传输函数HLP(s)转换为HD(s) 。（附录B） Step5:对传输函数HD(s)进行双线性变换，从而得到所求的数字IIR传输函数GD(z) 第二种方法分为以下几个步骤： 步骤1：用式 预畸变所求数字滤波器GD(z)的数字频率指标，从而得到一个等价的模拟滤波器HD(s)的频率指标。 步骤2：选取一种合适的频率变换，将HD(s)的频率指标转换成原型低通滤波器HLP(s)的频率指标。 步骤3：设计模拟低通滤波器HLP(s)。 步骤4：对传输函数HLP(s)进行双线性变换，将其转换为IIR数字滤波器的传输函数GLP(z)。 步骤5：选取一种合适的谱变换将GLP(z)转换成所求的数字传输函数GD(z)。 GD(z)的数字频率指标 与HD(s)相同类型的模拟滤波器频率指标 原型低通滤波器HLP(s) 的频率指标 Designing of HLP(s) HD(s) GD(z) GLP(z) GD(z) 谱变换可以用来将给定的低通数字IIR传输函数GL(z)转换成另一个低通、高通、带通或带阻滤波器的数字传输函数GD(z) . 7.5 IIR