digtlal signal processing —— using matlab - read.ppt

Digtlal Signal Processing—— Using MATLAB 第七章 FIR滤波器设计 数字频率w的概念 定义： 其中：Ω=2πf为模拟角频率 T:抽样时间间隔，fs:抽样频率 所以数字滤波器设计必须给出抽样频率 数字频率的2π等价于模拟抽样频率Ωs=2πfs 按照Nyquist抽样定理，基带信号的频率特性只能限于|w|ws/2=π的范围 数字滤波器幅度响应(1) 数字滤波器幅度响应(2) 7.1 概论 滤波器设计：给定技术要求设计系统 设计步骤： 确定技术要求：由具体应用条件决定 提供一个逼近要求的滤波器的表述 根据表述实现滤波器 下面讨论时我们均假设技术要求已知 7.1.1 技术要求的给定 幅度要求： 绝对指标要求：对幅度响应|H(ejw)|给出要求 相对指标要求：以分贝dB形式给出 相位要求：线性相位 一、绝对指标要求(1) 绝对指标(2) 频带 [0,wp] 称为 通带passband, δ1 是在理想通带响应上可以接受的容度（或波纹） 频带[ws,pi] 称为 阻带stopband, δ2 是相应的阻带容度（或波纹） 频带[wp, ws] 称为 过渡带transition band, 在这个频带内幅度响应不作要求 二、 相对指标要求(1) 相对指标(2) Rp:以dB计的通带波纹 As:以dB计的阻带衰减 两种指标之间的关系： 三、为什么只讨论低通滤波器(LPF) 上述指标都是针对低通滤波器的 其他类型的频率选择性滤波器(如高通或带通)也能给出类似要求 滤波器设计最重要的参数是频带容限和频带边缘频率 四、技术指标举例 设计一个低通滤波器，它具有一个通带 [0,wp] ，通带内频带容限为δ1(或Rp，单位 dB)，一个阻带[ws,pi]，阻带内容度为δ2(或As,单位dB) 最后求得结果是得出滤波器的系统函数H(z)或差分方程 五、FIR滤波器的优点 相位响应可以真正线性 系统绝对稳定，设计相对容易 高效实现 可用DFT实现 实际应用时，我们感兴趣的是线性相位的FIR滤波器 六、线性相位响应的优点 设计问题中仅有实数运算 时延固定，没有时延失真 对长为M的滤波器，运算次数只有M/2量级 7.2 线性相位FIR滤波器性质 包括脉冲和频率响应的形状，系统函数零点的位置 设h(n)是长为M的脉冲响应，0≤n≤M-1，则 线性相位的脉冲响应形状(1) 因为频率响应函数具有线性相位 线性相位的脉冲响应形状(1) 线性相位的脉冲响应形状(2) 第二类线性相位满足条件 线性相位的脉冲响应形状(2) 对应频率响应特性H(ejw) 将M为奇和偶数结合对称和反对称的情况， 得到四种类型的线性FIR滤波器 对应每种类型其频率响应特性都有独特性质，令 例：设脉冲响应为h(n)={1，1，1，1}，求出并画出频率响应 解：频率响应函数为 I类线性相位：对称脉冲响应，M为奇数 这种情况下，beta=0，alpha=(M-1)/2是整数 h(n)=h(M-1-n)， 0≤n≤M-1 II类线性相位:对称脉冲响应，M为偶数 这种情况下，beta=0，alpha=(M-1)/2不是整数 h(n)=h(M-1-n)， 0≤n≤M-1 III类线性相位：反对称脉冲响应，M为奇数 这种情况下，beta=pi/2，alpha=(M-1)/2是整数 h(n)=-h(M-1-n)， 0≤n≤M-1 IV类线性相位：反对称脉冲响应，M为偶数 这种情况和II类似，有 MATLAB实现 Hr_type1:求I类线性相位的Hr(w) 调用格式：[Hr,w,a,L]=Hr_type1(h) Hr_type2:求II类线性相位的Hr(w) 调用格式：[Hr,w,b,L]=Hr_type2(h) Hr_type3:求III类线性相位的Hr(w) 调用格式：[Hr,w,c,L]=Hr_type3(h) Hr_type4:求IV类线性相位的Hr(w) 调用格式：[Hr,w,d,L]=Hr_type4(h) 小结 了解了线性相位FIR滤波器的各种特性，便可根据实际需要选择合适的FIR滤波器，同时设计时要遵循有关约束条件。 如：第3、4种情况，对于任何频率都有固定的π/2相移，一般微分器及90°相移器采用这两种情况，而选频性滤波器则用第1、2种情况。 (1)设计线性相位的低通Digtal Filter 从幅度特性考虑，只能选择第1种或第2种 第一种： 第二种 (2)设计线性相位的高通DF 从幅度特性看，可用第一种或第四种 第一种 第四种 (3)设计线性相位的带阻DF 从幅度特性考虑，只能选