数字滤波器的原理
1. 基本的系统分析
下面是一个基本的延时器,作用就是将输入的信号进行一拍的延时。所有的信号都能通过该系统,除了有一拍的延时外,并不会对信号有其他影响。
根据系统的转移函数,可以得到对应的旋转矢量,可以看到,旋转矢量其实就是一个在单位圆上逆时针旋转的单位旋转矢量。可以把系统的转移函数从幅度和相位两个部分组合起来看。可以看到,系统的幅度响应为1,相位响应呈线性。
2. 基本系统的分析方法
下面是一个基本的系统结构,其差分方程如下。系统的作用就是将输入的信号与其延时一拍的信号进行累加作为输出,若给系统一个单位冲击响应 [1,0,0],得到的输出将是 [1 ,1, 0];
对差分方程两边进行Z变换很容易就得出系统的转移函数。
2.1 对信号进行幅频响应和相频响应分析
对该系统的转移函数稍微进行处理,即可得到幅频响应和相频响应:
幅频响应如下:
相频响应如下:
其中ω是一个归一化频率,其是信号的频率相较于采样频率的归一化。
幅频响应的特性曲线如下:
可以把该傅里叶变换看成一个在单位圆上顺时针旋转的一个旋转矢量,旋转矢量在 [0, π) 内模长由2 逐渐减小至0,角度从0减小到-π/2;在 [π, 2π] 内模长由0逐渐增加至2,角度从π/2减小至0;
2.2 使用零点极点图粗略描绘幅频响应
将Z变换表示成分数形式们可以得到一个极点Z=0,一个零点Z=π;
这时候幅频响应就有如下表达式: 幅频响应就是|1+z| 和|z|这两个向量的模长比。
在单位圆上画出极点和零点的表示图,通过这个确定零点极点地位置,同时还能得到模长地信息,因此也就可以初略地画出幅频响应。
3. 信号经过系统
现在若让一个信号经过上面地系统,可以发现系统会对原始信号产生影响。
原始的两个信号是两个不同频率的单音信号的叠加,将该信号通过上面提出的最简单的一个系统,来观察其变化。
从频谱上看,该信号经过该系统过后,低频的幅度相较于原始的信号放大了,高频的信号相较于原来的信号被抑制了。这也就说明了这个最简单的滤波器系统具有一定的滤波能力。
4.滤波器的基本概念
滤波器对通过系统信号会对它的幅度和相位进行改变。滤波器主要有以下4类,分别对应的让不同频率的信号通过。
4.1 设计一个低通滤波器
低通滤波器就是能够让低频的信号通过,对到幅频特性上,就是在低频部分为1,在高频部分为0。根据对称性,很容易就能够得到滤波器在时域上表达,可以看出是一个Sa函数。
上面的是一个理想的滤波器,可以看到在时域上需要无穷多个点才能够实现一个理想的低通滤波器。在实际中这是不可能实现的。
因为信号的能力主要集中在主瓣,因此可以保留低通滤波器频谱的其主瓣即可。采用一个矩形脉冲序列和Sa函数序列进行相乘,就能得到一个实际中的滤波器的表达。在时域进行相乘,在频域就进行卷积,因此可以看到实际的低通滤波器其不再是平坦的,而是有一定的波动。
上面这个给原始信号乘以一个矩形脉冲序列的操作也叫做加窗。在实际中,可以选择不同的窗函数,根据窗函数的不同,能够得到不同的滤波器。
4.2 滤波器的性能指标
一个实际中的滤波器的幅频响应如下:可以把其分为通道,阻带和过渡带,在通带中,信号具有一定的波纹,该波纹对应着在通带内的最大的增益。通带的最大频率对应着一个通带的边缘,该处有一个具体增益值。
在阻带部分,信号不能超过阻带的波纹的增益。通带与阻带之间为过渡带。
在滤波器设计中,通常有一个-3dB带宽,该带宽的含义是,信号的功率谱密度下降为原来的1/2,也就是信号的能量在该频率下只有原来的一半,因此该频率通常也被叫做截至频谱。
参考: