今天给大侠带来FIR数字滤波器设计,由于篇幅较长,分三篇。今天带来第二篇,FIR数字滤波器设计基础,包括FIR数字滤波器的特点、线性相位条件以及基本结构。话不多说,上货。
数字滤波器的输入输出均为数字信号,信号通过数字滤波器后,可以改变频率成分的相对比例或滤除某些频率成分。数字滤波器可以分为IIR数字滤波器和FIR数字滤波器。
本篇只介绍FIR数字滤波器的设计,可以根据所给定的频率特性直接设计FIR数字滤波器。FIR数字滤波器在保证幅度特性满足要求的同时,能够做到严格的线性特性。
本篇采用了窗函数法、频率采样法以及基于firls函数和remez函数的最优化方法设计FIR滤波器。对FIR滤波器进行了详细的理论分析,并且对应于每种方法都给出了设计实例。通过编写MATLAB语言程序,运行程序,得到幅频和相频特性图。
对于窗函数和firls函数设计的滤波器,还通过建立Simulink系统模块进行仿真,观察滤波器滤波情况。
FIR滤波器在保证幅度特性的同时,很容易做到严格的线性相位特性。
在数字滤波器中,FIR滤波器的最主要特点是没有反馈回路,故不存在不稳定的问题。
同时,在幅度特性可以任意设置的同时,保证了精确的线性相位。
稳定和线性相位是FIR滤波器的突出优点。
另外还有以下特点:
设计方式是线性的;
硬件容易实现;
滤波器过渡过程具有有限区间;
相对IIR滤波器而言,阶次较高,其延迟也要比同样性能的IIR滤波器大得多。
设滤波器单位脉冲响应的长度为N,系统函数为关系式(2-1),如下:
由此式可见,H(z)是图片的(N-1)次多项式,它在Z平面上有(N-1)个零点,原点z=0是(N-1)阶重极点,位于r =1的单位圆内,系统永远稳定。稳定性和线性相位特性是FIR滤波器的突出优点。
FIR滤波器的设计任务是选择有线长度的h(n),使传输函数
满足要求。
线性相位条件为关系式(2-2),如下:
对于长度为N的h(n),传输函数为关系式(2-3),如下:
上式中,图片称为幅度特性,图片称为相位特性。线性相位是指相位函数图片满足如下特性:
或
是起始相位,
为常数,一般称第一种情况为第一类线性相位,称第二种情况为第二类线性相位。
满足第一类线性相位的充要条件是:h(n)为实序列,并且对(N-1)/2偶对称,
即:
满足第二类线性相位的充要条件是:h(n)为实序列,并且对(N-1)/2奇对称。
即:
。
FIR滤波器的基本结构有以下几种:直接型、级联型、线性相位型、频率采样型。
设FIR滤波器的单位冲击响应h(n)为一个长度为N的序列,则滤波器系统函数为关系式(2-4),如下所示:
表示这一系统输入输出关系的差分方程为关系式(2-5),如下所示:
直接由差分方程可得出对应的网络结构如图2-1所示:
图2-1 FIR滤波器的直接型结构
直接型结构的优点:简单直观,乘法运算量较少。
缺点:调整零点较难。
当需要控制滤波器的传输零点时,可将H(z)分解为实系数二阶因子的乘积形式,见关系式(2-6),如下所示:
上式中,
为
的 z变换,
为实数。级联型结构如图2-2所示:
图2-2 FIR滤波器的级联型结构
该结构的优点:调整零点比直接型方便。
缺点:图片中的系数比直接型多,因而需要的乘法器多。当图片的阶次高时,也不易分解。
FIR滤波器的线性相位结构有偶对称和奇对称,不论图片为偶对称还是奇对称都有:
当N为偶数时,系统函数为关系式(2-7),如下所示:
当N为奇数时,系统函数为关系式(2-8),如下所示:
对这两种情况,都可以用FIR直接型实现,其信号流图如图2-3所示。
(a)N为偶数
(b)N为奇数
图2-3 线性相位型结构
这种结构在本质上是直接型,但乘法次数比直接型省了一半。
频率采样型结构是一种用系数将滤波器参数化的一种实现结构。一个有限长序列可以由相同长度频域采样值惟一确定。
系统函数在单位圆上作N等分取样就是单位取样相应h(n)的离散傅里叶变换。H(k)与系统函数之间的关系可用内插公式表示,为关系式(2-9),如下所示:
上式中
,
这样,图片是由梳状滤波器和N个一阶网络的并联结构进行级联而成的,其网络结构(信号流图)如图2-3所示。是一个梳状网络,其零点为:
刚好和极点一样,等间隔地分布在单位圆上。理论上,极点和零点相互抵消,保证了网络的稳定性。
图2-5 FIR滤波器的频率采样结构
频率采样结构的优点:
1)在频率采样点
,
,只要调整H(k)就可以有效地调整频响特性。
2)只要h(n)长度N相同,对于任何频响,其梳状滤波器部分和N个一阶网络部分完全相同,只是各支路增益H(k)不同。相同部分便于标准化、模块化。
缺点:
1)寄存器长度都是有限的,零、级点可能不能正好抵消,造成系统不稳定。
2)当N很大时,其结构很复杂,需要的乘法器和延时单元很多。
第二篇就到这里,下一篇带来第三篇,FIR数字滤波器设计,包括窗函数法设计FIR滤波器、频率采样法设计FIR滤波器以及基于firls函数和remez函数的最优化方法设计FIR滤波器。
完
后续会持续更新,带来Vivado、 ISE、Quartus II 、candence等安装相关设计教程,学习资源、项目资源、好文推荐等,希望大侠持续关注。
江湖偌大,继续闯荡,愿大侠一切安好,有缘再见!
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。