我一直在阅读信号处理方面的一些论文,对问题的标题感到非常困惑。考虑时间$ t $,$ f(t)$的连续函数,我在不均匀的时间$ t_k $采样,其中$ k = 1,2,...,N $。对我而言,采样函数是有意义的:
$$ f_s(t)= \ sum_ {k = 1} ^ N \ delta_ {t,t_k} f(t),\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \(1)$$
其中$ \ delta_ {t,t_k} $是Kronecker的增量(当$ t = t_k $时等于$ 1 $,在其他位置为零)。但是,在本文中,作者将采样信号定义为:
$$ f_s(t)= \ frac {1} {N} \ sum_ {k = 1} ^ {N} f(t)\ delta (t-t_k),\ \ \(2)$$
其中$ \ delta(t-t_k)$是Dirac的delta函数,我真的不明白为什么$ 1 / N $出现在这里(作者声称采样函数实际上是增量函数的加权和
$$ s(t)= C \ frac {\ sum_ {k = 1} ^ {N} w_k \ delta(t-t_k)} {\ sum_ {k = 1} ^ {N} w_k},$$
,他在这里选择$ C = w_k = 1 $。我真的不明白为什么)。最后一条语句对我来说意义不大:采样信号在$ t = t_k $处具有无限的振幅!

尽管如此,在第二种情况下(等式$(2)$)定义$ f_s(t)$的傅立叶变换要容易得多,因为这只是窗口函数(狄拉克梳子的FT)和连续信号$ f(t)$的FT,而在等式$(1)$上,FT则稍微复杂一些,因为我们有一个整数函数(Kronecker的增量)相乘通过连续函数($ f(t)$)。有什么亮点吗?

#1 楼

在我的经验中,通过将连续时间信号乘以一列狄拉克脉冲来对采样过程进行建模是最常见的解释。如果您深入研究它,您会发现对这种方法的数学精度有一些分歧*,但我不会担心。这只是该过程的便捷模型。手机的ADC内没有脉冲发生器,可以产生周期性的闪电信号来增加其模拟输入。

您已经注意到,您无法计算Kronecker delta函数的连续时间傅立叶变换,因为它的域不是连续的(仅限于整数)。相比之下,狄拉克增量函数具有简单的傅立叶变换,并且由于其筛选特性而容易表现出将信号与一系列狄拉克脉冲相乘的效果。

*:例如,如果您要精确地进行数学计算,您会说Dirac增量根本不是函数,而是分布。但是,从工程学的角度来看,这些问题实际上只是语义。

编辑:我将在下面解决这个问题。您给出了采样过程的心理模型:

$$
f_s(t)= \ sum_ {k = 1} ^ N \ int_ {t_k- \ epsilon_k} ^ {t_k + \ epsilon_k} f(t)\ delta(t-t_k)dt。
$$$

这种解释的问题是典型的理想采样模型没有内置该积分。相反,它是输入信号与狄拉克脉冲串的纯乘积。如果您更仔细地看一下为$ f_s(t)$显示的等式,您会发现右侧实际上没有自变量; $ t $是积分的虚拟变量。对于上面的任何$ \ epsilon_k> 0 $,根据Dirac脉冲的筛选属性,您将得到:

$$
f_s(t)= \ sum_ {k = 1} ^ N f(t_k),
$$

这是不正确的。取而代之的是,采样信号的模型为:

$$
f_s(t)= \ sum_ {k =-\ infty} ^ {\ infty} f(t)\ delta(t-kT)
$$

与除了对沿时间轴的无限长的脉冲序列进行一般化,并假设在时刻$ t_k = kT $均匀采样数据外,其他内容与上述相同。所得信号的傅立叶变换为:

$$
\ begin {align}
F_s(\ omega)&= \ int _ {-\ infty} ^ {\ infty } f_s(t)e ^ {-j \ omega t} dt \\
&= \ int _ {-\ infty} ^ {\ infty} \ sum_ {k =-\ infty} ^ {\ infty} f (t)\ delta(t-kT)e ^ {-j \ omega t} dt \\
&= \ sum_ {k =-\ infty} ^ {\ infty} \ int _ {-\ infty} ^ {\ infty} f(t)\ delta(t-kT)e ^ {-j \ omega t} dt \\
&= \ sum_ {k =-\ infty} ^ {\ infty} f(kT )e ^ {-j \ omega kT}
\ end {align}
$$

如果定义信号$ f(t)$的离散采样版本等于$ x [n] = f(nT)$,那么您将剩下:

$$
$$

,这正是离散时间傅立叶变换的定义。

评论


$ \ begingroup $
您将如何解决振幅为“无限”的事实?我通常认为,实际上您并不是在离散时间$ t_k $上对信号进行“采样”,而是在给定时间$ \ Delta t_k $内对信号进行积分。但是,由于与Kronecker delta相同的原因,这种解释将违反任何形式的傅立叶变换计算。另外...为什么我给出的链接中的论文作者将狄拉克梳子除以$ N $?这对我来说毫无意义。
$ \ endgroup $
–Néstor
2012年7月23日下午0:35

$ \ begingroup $
实际上,您是对的。由于任何ADC的模拟前端的带宽有限,因此模拟信号总会有一些有效的“积分时间”。但是,理论构架不受此类担忧的限制。粗略地讲,脉冲的“无限高度”通过其“零宽度”来平衡,从而使其整合为一体。如果在这种情况下应用短时积分解释(乘以一个脉冲,在无限短时间内积分),然后通过筛选属性,您将得到$ x [n] = x(nT)$通常会出现。
$ \ endgroup $
–Jason R
2012年7月23日在1:25

$ \ begingroup $
是的,但是我关心的不是解释,而是傅里叶变换。假设我写的采样过程为:$$ f_ {s}(t)= \ sum_ {k = 1} ^ N \ int_ {t_k- \ epsilon_k} ^ {t_k + \ epsilon_k} f(t) \ delta(t-t_k)dt,$$您将如何进行傅立叶变换?我知道把$ \ epsilon_k \设为0 $时的窍门,但这对我来说没有多大意义(而且做FT更难!)。即使假设是这样,我也没有得到与Roberts等人的论文相同的窗口函数。我引用的。而且我坚持... $ 1 / N $对我来说没有任何意义。
$ \ endgroup $
–Néstor
2012年7月23日在2:17



$ \ begingroup $
好的,您的编辑和上面的评论中有我的错误。但是,由于Dirac的Delta,即$ f(t = t_k)\ to \ infty $,而在实际中,脉冲仍然在$ t = t_k $处具有无限幅度的事实,我仍然无法下定决心。观察(并要建模)是$ f(t = t_k)= f(t_k)$
$ \ endgroup $
–Néstor
2012年7月25日在4:48