有限脉冲响应数字滤波器设.......pptVIP

第七章 有限脉冲响应数字滤波器的设计 三 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数综合表示如下： 表明 如果Z=Zi是H(z)的零点，则Z=Zi-1也是H(z)的零点。 h(n)为实序列，零点必定共轭成对。则Zi*和(Zi-1)*也是H(z)的零点；即H(z)的零点必定互为倒数的共轭对。 三 线性相位FIR滤波器零点分布特点 零点分布的几种情况 (1)一般情况：零点是互为倒数的两组共轭对，确定了一个零点，就可确定其它三个。 (2)Zi在|Z|=1且不在实轴上，零点是一组共轭对，因为共轭对的倒数是它们本身。 (3)Zi在实轴上且不在|Z|=1上，零点是互为倒数两个实数； (4)Zi在实轴上且在|Z|=1上，则零点只有一个，位于Z=1或Z=?1。 四 线性相位FIR滤波器网络结构 1、当N为偶数 四 线性相位FIR滤波器网络结构 四 线性相位FIR滤波器网络结构 2、N为奇数时 四 线性相位FIR滤波器网络结构 例:某FIR DF的系统函数H(z)＝1＋2z-1+4z-2+2z-3+z–4，试求h(n)的表达式，判断是否具有线性相位，画出该滤波器流图的直接结构和线性相位结构形式。 7.2 用窗函数法设计FIR滤波器 一、设计过程 7.2 用窗函数法设计FIR滤波器 二、设计思想 1、Hd(ejω) Hd(ejω)为希望设计的滤波器传输函数，hd(n)是与其对应的单位脉冲响应 ： 问题 如果Hd(ejω)是逐段恒定的，即边界频率处有断点，则hd(n)是无限时宽，且是非因果序列，这样的系统是不能实现的。 7.2 用窗函数法设计FIR滤波器 例：理想低通滤波器的传输函数Hd(ejω)为 7.2 用窗函数法设计FIR滤波器 2、加窗截断 加窗后的单位响应为： h(n)= hd(n)?RN(n) 用一个有限长的序列h(n)去代替一个无限长的序列hd(n) 问题： 产生误差； 时域中是截断处理； 频域中滤波器的通带和阻带中有波动。 H(ω)形成的过程 H(w)与理想低通Hd(w)差别 H(w)在 w=wC附近形成过渡带，过渡带宽度B=4?/N，近似于矩形序列幅度谱RN(w)的主瓣宽度； 通带和阻带内有波动，最大的峰值在 w=wC?2?/N 处，最大的负峰值在 w=wC+2?/N 处。幅度谱RN(w)波动越快(N加大)，通带、阻带内波动越快，其旁瓣的大小直接影响H(w)波动的大小。 减小吉布斯效应措施 增加N值 减小过渡带宽度，由于主瓣与旁瓣幅度也增加，且主瓣和旁瓣的相对值不变， H(w)的波动幅度没有改变。 带内最大肩峰比H(0)高8.95%，阻带最大负峰比零值超过8.95% 。使阻带最小的衰减只有21dB。 谱间干扰未减小，波动更明显，因此加大N并不是减少吉布斯效应的有效方法； 改善窗函数的形状 对于窗函数的形状主要考虑以下2点因素： 尽量减小主瓣宽度，以获得较窄的过渡带； 尽量使窗函数的最大副瓣相对于主瓣要小，使设计出来的滤波器幅度特性中肩峰和余振较小，阻带衰减较大。 几种常见的窗函数 矩形窗 (Rectangle Window) wR(n)=RN(n) 其频率响应为 三角窗 (Bartlett Window)，巴特利特窗 汉宁(Hanning)窗——升余弦窗 哈明窗 (Hamming Window)—改进的余弦窗 布莱克曼窗 (Blackman Window):二阶升余弦窗 凯塞-贝塞尔窗 (Kaiser-Basel Window) 凯塞-贝塞尔窗 (Kaiser-Basel Window) 六种窗函数的基本参数 窗函数选择方法 根据过渡带与阻带衰减的要求，选择满足条件的窗函数形式，并估计窗口长度N。原则是保证阻带衰减的前提下，尽量选主瓣窄的窗函数。 3、计算滤波器的单位脉冲响应h(n) h(n)= hd(n)?w(n)。其中 w(n)是上面选择好的窗函数，hd(n)与w(n)都应满足线性相位要求。 4．验证技术指标是否满足要求 已设计出的滤波器的频率响应 。验算H(ejw)是否满足设计要求，若H(ejw)不满足要求，重复上面2，3，4过程。 窗函数法优缺点： 1、优点 从时域出发的设计方法，设计简单，方便，实用。 2、缺点 要求用计算机实现，边界频率不易控制。 例：用矩形窗设计法设计一个FIR线性相位低通滤波器，已知ωc =0.5π，N=21。 ?＝(N-1)/2，画出h(n)和20lg|H(ω)/H(0)|曲线，再计算正、负肩峰的位置和过渡带宽带度。 加矩形窗：h(n)=hd(n).RN(n)