为什么在过渡带较窄的数字滤波器的输出中出现振铃?


12

我正在为音频的频谱修改类型效果做一些“极限”运算。我正在使用砖墙式滤波器,以及非常窄带的通过和拒绝滤波器(vst插件),我想知道对于线性相位/最小相位滤波器的前置/后置“振铃”,我是否可以做些什么。不幸的是,我必须使用陡峭的eq斜率。我准备使用最小相位,因为它避免了预振铃。

具体来说,我想知道:

  1. 在最小相位滤波器中,输入之后刚好导致脉冲响应振荡的原因是什么?

  2. 这些振荡是否会导致通过陡峭的斜率滤波将通带之前和之后的“振铃”声传给您?

  3. 振荡,因此振铃频率是否总是相同的频率,或者振铃频率是否在某种程度上取决于输入信号?

非常感谢您的专业知识。我期待任何回应。大叔


我已经合并了您的帐户,因此您应该可以在此处的答案下发表评论。我还附加了您的最新帖子作为编辑内容。请删除其中实际上是对一些答复者的评论的部分,并相应地进行评论(现在您可以这样做)
Lorem Ipsum 2012年

Answers:


8

根据OP修订的问题和其他评论进行编辑。

我不同意@JasonR的说法,即滤波器振铃是由于Gibbs现象引起的

如与Jason的答案相关的Wikipedia文章中所述,吉布斯现象是对周期但不连续的信号(例如方波或锯齿波)的傅里叶级数的截断和(前 项)的渐近行为的观察。Wikipedia文章举例说明了方波,它显示出随着越来越多的项(变大),截短的傅立叶和变得越来越接近方波。有迹象表明,发生绕开关瞬间,其中方波的转变从振荡,或反之亦然,但这些越来越小如ñ ñ ñ →交通nnn变大。正如杰森正确指出的那样,振荡幅度变小,频率增加,(观察到的)持续时间也变小。总的来说,截断的傅立叶和看起来像在的极限处收敛到方波。n

吉布斯现象是,观测即使在极限去,n傅立叶级数之和不收敛到高值或者在开关瞬间,其中方波急剧变化值的低值。(会在所有其他时刻发生收敛 )。这 与滤波本身无关,除了可以将截断的傅立叶和视为具有方波输入的理想砖墙低通滤波器的输出之外。如果过滤器截止频率是前 ñ ñn谐波通过不变,高次谐波被阻止,输出是前 项的截短傅立叶和。但是在极限情况下(即发生吉布斯现象的时候),没有滤波器:所有谐波都直接传递到输出而没有任何变化。因此,我不同意滤波器振铃是由于吉布斯现象引起的。n

那么为什么会发生振铃? 所有(非平凡的)滤波器会响起,无论它们是否为砖墙,无论输入信号的形状如何,以及输入是连续的还是陡峭的过渡。原因是,如果输入在被停止的频带中(无论是全部还是大部分)具有能量,则该能量将有效地内部存储在滤波器中,并随着时间的流逝随着带内能量缓慢释放。在大多数情况下,此版本不会引起太多关注,因为它被存在的带内信号的响应淹没了。但是,如果带内信号相对突然改变(或停止),则仍必须释放以前存储的能量,这是在带内信号消失后观察到的振铃。用DSP来讲,FIR滤波器缓冲器即使在信号结束后仍会继续排空,因此即使在信号结束后输出仍会继续。由于锐截止过滤器具有较长的缓冲区(如果需要,则有许多双二阶节),因此清空需要很长时间,并且比更易于使用的过滤器清空要快得多。


1
谢谢迪利普。我看不到您不同意Jason R的帖子。据我了解,你们都说过吉布斯现象确实会引起振铃。您是否同意吉布斯现象相对于其他振铃原因的相对重要性/贡献是振铃的原因?
Dale Newton

请参阅我的修订答案。
Dilip Sarwate 2012年

1
Dilip是正确的,因为您可以从吉布斯现象以外的其他来源的滤波器中获得振铃。如果截断周期性不连续函数的频谱,则结果信号在时域中的振荡就是吉布斯现象。大多数低通滤波器会产生类似的效果,因为它们的脉冲响应通常类似于Sinc函数,其形式与Gibbs现象振荡相同。当您将信号与看起来像正弦的脉冲响应进行卷积时,由于脉冲响应的旁瓣,您经常会看到振铃。
杰森R

@Dilip:我不确定我会说所有非平凡的过滤器都会响。是的,所有非平凡滤波器的脉冲响应都将比脉冲更宽,因此随着时间的流逝,输入信号中的能量会有所散布。但是,这并不意味着您会看到输出振荡。以高斯低通滤波器为例。
杰森R

1
@Dilip:对于这个问题的吉布斯部分,我认为我一点都不反对。我同意您对《任择议定书》问题的回答更为准确。而且,我只听过振荡中使用的“振铃”一词,类似于吉布斯效应。我可能将非振荡效应称为“过冲”,“下冲”或“滤波器瞬变”。不过,这只是术语。我不认为一种解释在客观上是正确的。
杰森R

8

您的观察结果就是吉布斯现象的一个例子。当您使用具有非常陡峭的过渡带的滤波器时,您会观察到输入信号中任何陡峭的过渡(例如,脉冲波形的边界)附近的滤波器输出(或“振铃”)中的振荡。振荡的表观“频率”取决于滤波器的带宽。为您增加滤波器的截止频率,振荡将更加本地化的时间(即“频率较高的”),但峰值超调确实没有变化。上面链接的维基百科文章在一半左右就有很好的解释。


8
  1. 正如杰森(Jason)指出的那样,这里有一个基本的“不确定性原理”:频率非常窄的所有事物的时间都是很宽的,反之亦然。
  2. 如果使用最少的过滤器,则不应有振铃,只有振铃。预振铃仅在线性相位滤波器中发生。预振铃比后振铃更容易听见,因此在这里最小的滤波器往往是更好的选择。它在测量上可能看起来很糟糕,但除非极端,否则由于人耳听觉系统的某些掩盖特性,振铃后的声音不太明显
  3. 它们的振铃通常恰好在滤波器的转折频率处。即2 kHz低通滤波器将产生2 kHz振铃,因此频率是滤波器的函数,而不是内容的函数。内容虽然会有所不同。如果内容少或没有2 kHz,则振铃不会太大。

3

具有陡峭过渡和平坦通带的带通滤波器接近矩形。

一个FT域中的矩形是另一域中的Sinc函数。对于在时域中的矩形窗口在频域中创建频谱“泄漏”而言,这是正确的。或对于频域中的矩形窗口在时域中创建螺旋包。矩形(带宽)越窄,Sinc越宽。(并且Sinc函数在两侧都“响”)。对于一个域中给定的宽度,获得比另一个域中的Sinc更窄的能量范围的唯一方法是使用看起来比矩形更接近高斯的东西,例如没有陡峭的边缘。

现在考虑在一个域中移动该矩形(例如,更改带通滤波器的通带频率)。一个DFT域中的圆形移位是另一域中的线性相位旋转。具有复数共轭的总和可得到真实的响应,而两个相反且快速旋转的复数指数螺旋数据包将成为振铃时域响应。振铃的速度将与带通中心频率相关,而振铃的长度将与带宽的窄度和过渡陡度相关。如果在信封消亡之前,螺旋线旋转了超过半圈,则会出现振铃。使信封在一个域中更快消失的方法是在另一个域中使用更广泛的舍入函数。

第2部分:

如果您使用Remez或Parks-McClellen工具设计滤波器,则最终将得到等波纹响应。一个FT域中的正弦曲线是另一个FT域中的脉冲。因此,频域中的等波纹将是时域中的冲动或“滴答”。该“滴答”将通过频域中纹波的“频率”从脉冲响应的中心移开。Remez设计的滤波器越平坦,纹波越快,脉冲响应中的“滴答声”就越多。那是预响的一部分。使用较不积极的滤波器设计方法可以避免这种情况。


谢谢Hotpaw。由于目前我对DSP的理解非常基础,因此我仅在非常确定的水平上理解您的解释。我不知道用于制作数字滤波器的工具和套件,而我还不了解螺旋包和复杂的共轭物,因为我还不具备数学技能。
Dale Newton

我想问一下:您说的是振铃的速度与滤波器的中心频率有关。这是否意味着您不同意希尔马斯的说法,即振铃发生在转折频率上?还是这两个频率都响?
Dale Newton

矩形的角也可以根据矩形的中心和宽度来定义。
hotpaw2
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.