用频率采样法设计用频率采样法设计FIR滤波器滤波器
窗函数法和频率采样法是两种较为典型的FIR数字滤波器设计方法。目前,相关的《数字信号处理》教科书对窗
函数法设计FIR滤波器进行了较为详细的论述,但对用频率采样法设计FIR滤波器这部分内容讲解得不够细致,
让初学的学生感到难以理解。针对用频率采样法设计FIR滤波器的相关问题进行了较为深入的探讨,并结合实例
借助Matlab软件进行了仿真和验证。仿真结果表明,选择合适的过渡采样点和滤波器长度,可以有效地控制阻
带衰减、过渡带宽及计算复杂度。
有限长脉冲响应(FIR)
1 设计原理及滤波器性能分析设计原理及滤波器性能分析
频率采样法是从频域出发,对给定的理想滤波器的频响 进行N点等间隔采样,即 ,然后
以此Hd(k)作为实际FIR滤波器的频率特性采样值H(k),即令:
由DFT定义,可以用这N个频域的采样值H(k)来惟一确定FIR的单位脉冲响应h(n),即:
下面对设计出的滤波器 ,这个 将逼近理想滤波器的频响 。 的内插公式为:
从式(5)可以看到,在各频率采样点上,设计的滤波器,实际的频率响应严格地与理想滤波器的频率响应数值相等,即
。但是在采样点之间的频率响应是由各采样点的加权内插函数叠加而
形成的,因而有一定的逼近误差。该误差大小取决于理想频率响应的形状,理想频响特性变化越平缓,内插值越接近理想值,
逼近误差越小;反之,如果采样点之间的理想频响特性变化越陡,则内插值与理想值之间的误差越大,因而在理想滤波器不连
续点的两边,就会产生尖峰,而在通带和阻带就会产生波纹。用频率采样法设计的实际滤波器频率响应如图1所示。由图1可
知,实际滤波器的阻带衰减取决于内插函数第一旁瓣幅度值的大小,其大小决定了所设计的滤波器的阻带性能。
2 线性相位条件线性相位条件
FIR滤波器的最大优点是严格的线性相位特性。下面讨论为实现线性相位,在频域内采样得到的Hd(k)应满足什么条件。FIR
滤波器具有线性相位的条件是h(n)是实序列,且满足h(n)=±h(N-1-n),即h(n)关于 对称,其中N为滤波器的长
度。以第一类线性相位条件h(n)=h(N-1-n)(偶对称)为基础来推导频域采样Hd(k)满足的条件。
理想滤波器的频率响应可以表示为:
为实现第一类线性相位条件,相位函数θ(ω)和幅度函数Hg(ω)应分别满足: