第7章_FIR数字滤波器的设计方法.ppt

第7章FIR数字滤波器的设计方法 7 1线性相位FIR数字滤波器的特点7 2FIR数字滤波器的窗函数设计法7 3频率抽样法设计FIR数字滤波器7 4IIR和FIR数字滤波器的比较 7 1线性相位FIR数字滤波器的特点 7 1 1FIR数字滤波器具有线性相位的条件 对于离散线性时不变系统 系统频率响应可表示为 其中称为滤波器的幅度响应 称为相位响应 幅度响应反映了滤波器对信号通过该滤波器后各频率成分衰减的强弱 而相位响应则反映了各频率成分通过滤波器时间上延时的多少 滤波器具有线性相位是指相位响应与频率呈线性关系 满足 或者 其中 都是常数 当 也标志滤波器具有线性相位 对于FIR滤波器 设其单位冲激响应为h n 长度为N 则 对应的系统函数为 下面讨论有限长单位冲激响应h n 为实序列 并关于 偶对称或者奇对称两种情况的相位特性 1 有限长单位冲激响应h n 为偶对称 令 则有 因此 系统的频率响应 幅度响应 相位响应 2 有限长单位冲激响应h n 为奇对称 幅度响应 相位响应 结论 当FIR滤波器的单位冲激响应h n 为实序列 具有偶对称性或者奇对称 那么该FIR滤波器具有线性相位特性 7 1 2线性相位FIR数字滤波器的幅度特点 1 h n 偶对称 N为奇数 由于N为奇数 中间项为 其余项偶对称 关于偶对称 可实现任意形式滤波器 2 h n 偶对称 N为奇数 FIR滤波器的幅度响应如右图所示 幅度响应具有以下特点 1 当时 对呈奇对称 对呈偶对称 因此 具有h n 偶对称 N为偶数的FIR滤波器不能用于高通滤波器或者带阻滤波器 3 h n 奇对称 N为奇数 幅度响应具有以下特点 1 当时 对呈奇对称 因此 具有h n 奇对称 N为奇数的FIR滤波器只能实现带通滤波器 4 h n 奇对称 N为偶数 幅度响应具有以下特点 1 当时 对呈偶对称 因此 具有h n 奇对称 N为偶数的FIR滤波器不能实现低通 带阻滤波器 对 呈奇对称 7 1 3线性相位FIR数字滤波器零点分布特点 线性相位FIR滤波器的系统函数 线性相位FIR数字滤波器零点分布特点 互为共轭对称 线性相位FIR滤波器零点分布图 7 2FIR滤波器的窗函数设计法 7 2 1设计原理 1 根据滤波器技术要求直接设计一个理想数字滤波器的频率响应 2 求具有上述频率响应的滤波器对应的冲激响应序列理想频率响应的滤波器具有矩形频率特性 因而其时域响应为无限长的 3 将求得的冲激响应序列截短 并使之具有对称性 下面以低通滤波器为例说明设计过程 理想低通数字滤波器的频率响应为 式中表示群延时 其单位冲激响应为 从上式可以看出 理想低通滤波器的单位冲激响应是无限长非因果序列 且具有偶对称性 其中心点在如下图所示 理想单位冲激响应 用一个有限长的窗函数序列 比如说矩形窗 来截取 截取后的实际单位冲激响应h n 如下图所示 它是一个以 为对称点的对称图形 实际单位冲激响应 这时FIR滤波器的频率响应为 从上式中可以看出 FIR滤波器的频率响应实际上是的傅里叶级数的有限项之和 所以窗函数法也称为傅里叶级数法 显然 N越大 和的差别就越小 但所需的计算量就越大 7 2 2窗函数设计法的截断效应 根据卷积定理可知 时域相乘则频域相卷 实际FIR滤波器的频率特性应为 式中是窗函数w n 的傅里叶变换 对于理想低通滤波器 为了保证截取的h n 中心对称 需保证其相移满足 实际滤波器的频率响应为 说明 1 加窗后的滤波器仍然具有线性相位 2 加窗后滤波器的幅度特性等于理想低通滤波器的幅度特性与矩形窗的幅度特性的卷积 即 理想滤波器幅度特性 窗函数幅度特性 卷积结果 卷积情况 卷积情况 卷积情况 加窗处理的影响 1 理想频率特性的边沿形成过渡 宽度近似为窗函数频率响应的主瓣宽度 2 在截止频率的两边 出现最大肩峰 幅度取决于窗函数幅度特性旁瓣的相对幅度 波动的次数取决于旁瓣的个数 3 改变截取长度N 仅能改变过渡带的宽度 不能改变的肩峰相对值 称为吉布斯现象 7 2 3常用窗函数及其选择原则 1 主瓣宽度窄 以获得较陡的过渡带 窗函数选择要求 2 最大旁瓣相对主瓣值尽可能小 改善通带平稳度和增大阻带衰减 常用窗函数 1 矩形窗 窗函数 频率响应 2 三角窗 窗函数 频率响应 3 汉宁窗 升余弦窗 窗函数 频率响应 时 频率响应的幅度函数近似为 4 海明窗 改进升余弦窗 窗函数 频率响应 时 频率响应的幅度函数近似为 为了便于比较和设计时参考 下表归纳常用窗函数的基本参数 7 2 4窗函数设计步骤 1 首先由给定滤波器的理想频率响应 利用傅里叶拟变换 求理想滤波器的单位冲激响应即 2 根据过渡带以及阻带衰减的要求 选择窗函数的形式 并估计窗口的大小 3 计算滤波器的单位冲激响应 4 检验设计出的数字滤波器是否满足技术指标 例7 1设计一个线性相位FIR数字低通滤波器 通带截至频率 阻带起始频率 阻带衰减不小于50dB 抽样频率 解 1 由模拟滤波器指标求数字滤波器指标 通带截至频率 阻带起始频率 2 确定理想线性相位数字低通滤波器的频率响应 其对应的数字频率 利用序列傅里叶变换拟变换 可求得理想低通滤波器的单位冲激响应 3 由阻带衰减来确定窗形状 由过渡带确定N 由于 查表可知 选择海明窗其阻带最小衰减53dB满足要求 所要求的过渡带宽 由于海明窗过渡带宽满足 所以 4 由海明窗表达式确定FIR滤波器的h n 海明窗函数序列为 5 由h n 求 检验各项指标是否满足要求 如不满足可改变N 或改变窗函数形状来重新计算 7 3频率抽样法设计FIR滤波器 7 3 1设计原理 2 根据傅立叶变换的特性 N点频响采样通过傅立叶反变换成N点时域序列 即为FIR滤波器时域冲激响应 1 对理想滤波器频响进行N点采样 设序列h n 长度为M H z Z h n 在单位圆上等间隔抽样N点 得到 只要抽样点数N与序列长度M满足 则M个频域抽样值H k 能够不失真地代表该有限长序列h n 由此可见 只要知道FIR数字滤波器的单位冲激响应h n 的长度N 频率响应函数在上的N点等间隔抽样就可以确定滤波器的单位冲激响应h n 或者系统函数H z 其滤波器的频率响应为 因此 频率抽样法就是根据频域抽样理论 将待设计的滤波器的频率响应 在到之间进行等间隔抽样N点 得到 即 然后再利用 求得H k 和h n 以及 7 3 2频率抽样法的优化设计 在频率采样点上 实际滤波器频率响应与理想滤波器频率响应的数值严格相等 在频率采样点之间 实际滤波器频率响应由各采用点的加权插值形成 与理想滤波器频率响应存在误差 理想滤波器频率响应变化越平缓 插值结果与实际频响曲线越接近 增加采样点可以减小插值误差 优化措施 在实际的设计过程中 通常在理想频率响应不连续点的边缘增加一些过渡的抽样点 一般增加1 3个 7 3 3线性相位与抽样的约束 1 频率抽样的两种方法 第一种 起始点在 第二种 起始点在 2 线性相位的约束 1 第一类线性相位FIR滤波器 h n h N 1 n N为奇数 2 第二类线性相位FIR滤波器 h n h N 1 n N为偶数 幅度值 相位值 幅度值 相位值 3 第三类线性相位FIR滤波器 h n h N 1 n N为奇数 幅度值 相位值 4 第四类线性相位FIR滤波器 h n h N 1 n N为偶数 幅度值 相位值 7 3 4频率抽样法的设计步骤 1 根据理想滤波器性能指标 计算在通带和阻带中的抽样点数 确定滤波器冲激响应的对称性 2 根据单位冲激响应的对称性 计算抽样点的幅度 和相位 3 由求 4 检验是否满足设计要求 若不满足 修改抽样点参数 重复 1 2 3 步计算直至满足要求 例7 2利用频率抽样法设计线性相位低通滤波器 单位冲激响应满足h n h N 1 n 幅度相应为 要求截止频率 抽样点数N 33 讨论过渡点数以及过渡点的取值大小对滤波器的影响 解 1 首先确定截止频率所处的位置 因为抽样点数N 33 抽样间隔 所以 即截止频率位于第8个抽样与第9个抽样点之间 取 2 确定各抽样点的幅度和相位大小由于属于第一类线性相位FIR滤波器设计 所以得 3 利用IDFT求FIR滤波器的单位冲激响应 由抽样的幅度和相位值求理想滤波器频率响应的抽样点 然后求所设滤波器的单位冲激响应 4 最后计算所设FIR滤波器的频率响应 在第2步中可以设置过渡点 例如在本题中 可把第9点和第24点作为不为零的过渡点 7 4IIR和FIR滤波器的比较 1 IIR利用模拟滤波器设计结果 简单 方便 有简单的解析式 不容