信号与系统分析——宗伟 8.ppt

已知模拟滤波器的传输函数Ha(s)为 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 解 首先将Ha(s)写成部分分式： 再通过 转换到z平面上，得到： 下面分析模拟频率Ω和数字频率ω之间的关系。 令s=jΩ,z=e jω，并代入(6.4.3)式中，有 例6.4.1试分别用脉冲响应不变法和双线性不变法将图6.4.4所示的RC低通滤波器转换成数字滤波器。 解 首先按照图6.4.4写出该滤波器的传输函数Ha(s)为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 下面我们总结利用模拟滤波器设计IIR数字低通滤波器的步骤。 (1)确定数字低通滤波器的技术指标：通带截止频率ωp、通带衰减αp、阻带截止频率ωs、阻带衰减αs。 (2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。 (3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。 (4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。 通过以上分析可知，对 加矩形窗处理后， 和原理想低通 差别有以下两点： (1)在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度，近似等于 主瓣宽度，即4π/N。 (2)通带内增加了波动，最大的峰值在ωc-2π/N处。阻带内产生了余振，最大的负峰在ωc+2π/N处。 在主瓣附近， 可近似为 下面介绍用窗函数设计FIR滤波器的步骤。 (1)根据技术要求确定待求滤波器的单位取样响应 。如果给出待求滤波器的频响为 ，那么单位取样响应用下式求出： 例如，理想低通滤波器如(7.2.1)式所示，求出单位取样响应hd(n)如(7.2.2)式，重写如下： (2)根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。设待求滤波器的过渡带用Δω表示，它近似等于窗函数主瓣宽度。 (3) 计算滤波器的单位取样响应 (4)验算技术指标是否满足要求。设计出的滤波器频率响应用下式计算： 例7.3.1 利用频率采样法设计线性相位低通滤波器，要求截止频率ωc=π/2rad，采样点数N=33，选用h(n)=h(N-1-n)情况。 解 用理想低通作为逼近滤波器。 8.4 FIR数字滤波器的设计 1．FIR数字滤波器的性质 FIR DF 系统单位响应为： FIR数字滤波器的系统函数： 对于模拟理想低通滤波器，其频率特性具有线性性质；若FIR滤波器的 满足对称条件，则其具有严格线性的相频特性，即相位特性是频率 的线性函数。 FIR数字滤波器在一定的条件下，具有严格线性的相频特性。 系数 表示任意频率信号通过该系统产生相同的时间延迟。 FIR数字滤波器具有线性相位的充分必要条件是 单位响应 的序列必须以 为对称中心，时间延迟 为长度 的一半，即 FIR数字滤波器的系统函数： 判断某数字滤波器频率响应是否具有线性相位的方法： 1）直接求出频率响应判断； 2）利用FIR DF 具有线性相位的充要条件判断。 FIR 数字滤波器还有另外一个重要的性质，即由于系统函数只有零极点，系统必然稳定且是因果系统。 2．窗函数法 设希望设计的滤波器传输函数为 , 是与其对应的单位脉冲响应，因此 （1）由 时延 截断 逼近 设计过程为： （2）截断 窗函数的长度为 ，即单位响应的长度，则 （3）时延 矩形窗对理想低通 幅度特性的影响 根据频率采样定理，hM(k)与hd(k)应满足如下关系： 3．频域取样法 窗函数法设计 滤波器是从时域出发，频域取样法时从频域出发，直接对理想滤波器的频域特性 进行等间隔采样，然后求出相应的数字滤波器的系统函数 即为给定的系统频率特性的设计指标样本，在一个周期 内对理想滤波器的频率特性作N点均匀抽样，得 , 对 进行 可得到相应得系统单位响应 实际滤波器得系统函数和频率特性为 所设计得滤波器得频率特性为 式中， 式中， 为内插函数； 为权，且 由此，在采样点 上， 与理想值 严格相等，在非采样点上， 为N个以相应插值函数为权得采样值之和。 内插函数 具有线性相位特性， 也具有线性相位。 为使 具有线性相位，