汉明数字滤波器中的吉布斯现象

信息处理 傅里叶级数 吉布斯现象
2022-02-17 09:30:56

在汉明的“数字滤波器”中,有一个神秘的部分,他描述了吉布斯现象如何被视为两个函数中心之间的位移,因为它们被卷积在一起。这是在第 3 版的第 112 - 113 页上。

在此过程中,他表明截断傅立叶级数与将傅立叶系数与矩形函数相乘相同。然后他继续证明具有 2N+1 个系数为 1 的函数是:

h(θ)=sin(2N+1/2)θsin(θ/2)

我很困惑:我认为矩形的频率响应是一个 sinc 函数,而他在另一页展示了这一点(当他推导出 Lanczos 平滑因子时)。

谁能帮我解决这个问题?

1个回答

我认为主要问题是你超越了自己。你可能记得或在某处读到过傅里叶变换rect函数是一个sinc功能。这是真的; 但是,他在本节中没有提到傅立叶变换!事实上,他所做的并不是傅立叶变换。

他在本节中所做的是将任何周期函数表示为傅里叶级数:

(1) f(θ)=a02+k=0(akcoskθ+bksinkθ)

这里的关键是这个函数不包含每个频率。它包含频率 0(直流)和k2π在哪里k=1,2,3,.... 这实际上非常重要。这个函数总是周期性的,周期为2π. 系数akbk几乎是傅立叶变换的采样值(但不完全是,为什么?家庭作业!:p)。你有很多这些系数。在本书的后面,您将了解到,当您在一个空间(频率或时间)中采样数据时,它必然会使另一个空间(时间或频率)中的计数器部分具有周期性。

在第一种情况下,他进行了 Lanczos 平滑推导,通过在函数中运行一个矩形窗口(与 rect 卷积)来平均函数。他所展示的是,毫不奇怪,系数乘以这个项:

(2) sinkπ/NkπN

当然,这对你来说应该很熟悉,因为它是sinc功能。但是,您缺少的是k是离散的!它实际上是一个采样版本sinc功能。

实际上,他将函数与rect并表明所得傅立叶级数的系数(松散地读作:采样傅立叶变换)是采样的sinc功能。那里并不奇怪。Convolution thm 表示时间上的卷积变成了频率上的乘法。傅里叶变换,你知道的,rectsinc, 所以卷积由rect是乘以sinc在频率空间。

在下一节中,他做了一些不同的事情。他采用傅里叶级数(阅读:采样傅里叶变换)并删除所有较高频率的系数。实际上,他是在进行傅里叶变换,乘以一个rect,然后对其进行采样。为简单起见,他将所有没有被丢弃的傅里叶系数设置为1.

他剩下的是这样的:

(3) h(θ)=sin(N+1/2)θsin(θ/2)

你问,为什么这不是一个sinc功能?你现在能回答吗?

快速的答案是因为在频域中应用的不仅仅是截断,它是截断和采样算子您所知道的是,当您在频域中截断(即乘以矩形)时,时域会与sinc(通过卷积 thm 和傅里叶变换rect),但这是没有采样的。

至于为什么公式看起来像它的样子,有两种看待它的方法。他展示的第一个是你可以将傅里叶级数相加NN,这就是你得到的。

第二个,更深刻,可能会在以后出现,是(*)当你在一个空间(比如频率)中对一个函数进行采样时,另一个空间(比如时间)中的相应函数变成了它自身的移位版本的总和. 事实上,它可能不会完全像那样出现。典型情况是采样时域在频域中创建周期性复制(顺便说一句:这就是人们所说的混叠的原因)。但是,您可以应用傅里叶变换的对偶属性来获得 (*)。

等式(3)现在有意义吗?它是周期性的。越接近 0 越像sinc功能。

因此,您的练习是通过在频率空间中采样和截断并应用逆变换来导出等式(3)。