2013DSP实验讲义1.pptVIP

窗函数法设计FIRDF的具体步骤 FIRDF的设计方法： 如果希望得到的滤波器的理想频率响应为 ，那么FIR滤波器的设计就在于寻找一个传递函数 去逼近 ，逼近方法有三种： 窗函数设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等波纹逼近） 窗函数法设计FIRDF的思路 窗函数设计法是从单位脉冲响应序列着手，使h(n)逼近理想的单位脉冲响应序列hd(n)。 理想频响 通常具有片段常数特性，因此hd(n)是无限长非因果序列，不能作为FIRDF的单位脉冲响应。 窗函数设计法就是截取hd(n)为有限长的一段因果序列，并用合适的窗函数进行加权后作为FIRDF的单位脉冲响应h(n)。 截取长度和加权窗函数的类型直接影响逼近精度。 窗函数法设计FIRDF的具体步骤 （1）构造希望逼近的频率响应函数，以低通线性FIRDF为例，选择截止角频率为ωc，延时为 的线性相位理想低通滤波器，即 （2）求出hd(n)，即对上式进行IFT得到 窗函数法设计FIRDF的具体步骤 （3）加窗得到FIRDF的单位脉冲响应h(n) 式中w(n)为窗函数，长度为N。如果要设计严格线性相位FIRDF，则要求h(n)关于(N-1)/2点偶对称，而hd(n)是关于 点偶对称的，所以延时： ，同时w(n)也要关于(N-1)/2点偶对称。 矩形窗设计FIRDF的时频域示意图 矩形窗的卷积过程（P95的图4.5来说明） 加矩形窗引起的幅度逼近误差 在 处形成过渡带，宽为 ， 近似等于WRg(ω)的主瓣宽度。在 处幅度衰减6dB（幅度下降1/2）。 在通带和阻带均有波纹，由WRg(ω)的旁瓣引起，旁瓣幅度相对主瓣越大，波纹幅度就越大，旁瓣越多，波纹起伏密度就越大。 截取长度N的变化会影响过渡带的宽度变化，而不会影响通带和阻带波纹的最大起伏幅度。 截取长度N的改变不能改变主瓣与旁瓣的比例关系，只能改变WRg (ω)的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大起伏幅度永远为8.95%，这种现象称为吉布斯（Gibbs）效应。 举例：矩形窗设计的6dB截止?c=π/2 FIRDF的幅度响应 0 0.25 0.5 0.75 1 -40 -30 -21 -10 0 N=15 N=31 通过窗函数形式改善波纹幅度 波纹幅度的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。 改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求： ①窗谱主瓣宽度要窄，以获得较陡的过渡带； ②相对于主瓣幅度，旁瓣要尽可能小，使能量集中在主瓣中，这样就 可以减小肩峰（通带波纹）和余振（阻带波纹），提高通带平稳性和阻带衰减。 但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 常用窗函数 矩形窗 汉宁窗（升余弦窗） 三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为 。 常用窗函数 汉明窗（改进的升余弦窗） 对汉宁窗的改进，在主瓣宽度相同的情况下，旁瓣进一步减小，可使99.96%的能量集中在窗谱的主瓣内。 常用窗函数 布莱克曼窗（三阶升余弦窗） 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，为 。 几种常用窗函数时域图形 四种常用窗函数的频谱N=51，A=20lg|W(ω)/W(0)| 四种常用窗函数在同一设计指标下的滤波器的频率特性 四种常用窗函数的基本参数 实验步骤 第一步 明确并理解设计指标： 使用窗函数法设计数字低通滤波器，数字频率域指标为： 通带截止频率：0.25π，通带最大衰减：0.2db， 阻带截止频率：0.325π，阻带最小衰减：40db 第二步 理解设计指标： 数字频率域指标到模拟频率域指标的转化（通过采样率相互联系），设采样率为8KHz（对应2π），则模拟频率域指标为： 通带截止频率：1KHz，通带最大衰减：0.2db， 阻带截止频率：1.3KHz，阻带最小衰减：40db 1KHz频点的衰减不大于0.2dB(对数值0 ~ - 0.2dB) 1.3KHz频点的衰减不小于40dB(对数值小于-40dB)。 实验步骤 第三步 确定设计参数： 根据阻带最小衰减40db , 选择汉明窗 根据通带0.2db截止频率为0.25π，取6dB截止频率 ，可取初始值 可