傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（二）

由于本文公式所占用的字符比较多，无法在一篇博客中完整发布，所以将其分为两篇博客。本篇主要介绍了离散傅里叶变换 DFT 的内容，以及相关的总结。对于前置内容，包括傅里叶级数 FS，连续时间傅里叶变换 CTFT，以及离散时间傅里叶变换 DTFT 等，请移步以下链接。

傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（一）
https://blog.csdn.net/qq_33552519/article/details/130255215

对于离散傅里叶变换所衍生出的离散余弦变换和离散正弦变换，在我的另一篇文章中有所介绍，欢迎阅读与提出相关修改意见。

• DFT：https://blog.csdn.net/qq_33552519/article/details/124917126
• DCT：https://blog.csdn.net/qq_33552519/article/details/124917473
• DST：https://blog.csdn.net/qq_33552519/article/details/124962042

4 离散傅里叶变换

在前面的章节中，我们从周期函数的傅里叶级数扩展到了非周期函数的傅里叶变换，再从非周期函数的傅里叶变换扩展到了离散序列的离散时间傅里叶变换。然而对比三者可以发现，只有离散时间傅里叶变换的时域是离散的，而只有傅里叶级数的频域是离散的，并不存在一种变换的时域和频域都是离散的，这意味着它们都不适合数字计算机的处理。为了克服这种问题，一种自然而然的想法就是综合傅里叶级数的频域与离散时间傅里叶变换的时域，这就引出了我们接下来所要讨论的离散傅里叶变换。

对于一个定义在有限区间 $t\in \left(0,T\right)$ 上的连续时间函数 $f_a\left(t\right)$，对其以 $T$ 为周期进行周期性延拓，得到周期函数 ${\tilde{f}}_a\left(t\right)$，其中

$${\tilde{f}}_a\left(t+nT\right)=f_a\left(t\right),-\infty <n<\infty .\text{\hspace{1em}(4.1)}$$

那么可得 $f_a\left(t\right)$ 的连续时间傅里叶变换 $F_a\left(j\Omega \right)$ 与 ${\tilde{f}}_a\left(t\right)$ 的傅里叶级数系数 ${\tilde{F}}_a\left[n\right]$ 分别为

$$\begin{array}{c}{F}_a\left(j\Omega \right)=\underset{0}{\overset{T}{\int }}{f}_a\left(t\right)e^{-j\Omega t}dt,\\ {\tilde{F}}_a\left[n\right]=\frac{1}{T}\underset{0}{\overset{T}{\int }}{\tilde{f}}_a\left(t\right)e^{-jn{\Omega }_{T}t}dt=\frac{1}{T}{F}_a\left(jn{\Omega }_T\right),\\ where{\Omega }_T=2\pi /T.\end{array}\text{\hspace{1em}(4.2)}$$

也就是说，周期函数的傅里叶级数系数实际上是该周期函数在一个周期内的傅里叶变换以 ${\Omega }_T=2\pi /T$ 为间隔均匀采样与幅度缩放的结果。在讨论离散时间傅里叶变换时就可以发现，在时域中对连续时间函数以 $T$ 为间隔均匀采样会导致该函数的傅里叶变换以 ${\Omega }_T$ 为周期进行延拓与叠加。那以相反的视角来看式 (4.2)，是否是因为 $f_a\left(t\right)$ 的傅里叶变换 $F_a\left(j\Omega \right)$ 以 ${\Omega }_T$ 为间隔均匀采样，导致了 $f_a\left(t\right)$ 在时域以 $T$ 为周期进行延拓与叠加从而获得了周期函数 ${\tilde{f}}_a\left(t\right)$ 呢？而且因为 $f_a\left(t\right)$ 只在 $t\in \left(0,T\right)$ 上的区间定义，所以以 $T$ 为周期进行周期延拓所获得的周期函数 ${\tilde{f}}_a\left(t\right)$ 可以完整地恢复 $f_a\left(t\right)$，不会造成信息失真。

由于一般的连续时间函数的傅里叶变换通常不具有周期性，且在整个频率数轴 $\Omega \in \left(-\infty ,+\infty \right)$ 上都有定义，所以讨论其频域的采样对于计算机的实用性并不具有太大意义。然而，对于连续时间函数经过均匀采样所获得的离散时域序列，其离散时间傅里叶变换则以 $2\pi$ 为周期，因此在 $2\pi$ 的范围内对其进行采样，即可获得有限长度的频域离散序列。如果限制时域序列的长度，即可满足计算机对时域与频域转换的可实现性。为此，接下来我们主要基于离散时间序列及其离散时间傅里叶变换进行讨论。

对于离散序列 $f\left[n\right]$，其离散时间傅里叶变换为

$$F\left(e^{j\omega }\right)=\sum _{n=-\infty }^{\infty }f\left[n\right]e^{-j\omega n}.\text{\hspace{1em}(4.3)}$$

对 $F\left(e^{j\omega }\right)$ 以 $2\pi /M$ 为间隔均匀采样，其中 $M$ 为正整数。 类似于连续时间信号采样，引入周期冲激串

$$p\left(\omega \right)=\sum _{k=-\infty }^{\infty }\delta \left(\omega -\frac{2\pi k}{M}\right).\text{\hspace{1em}(4.4)}$$

则 $F\left(e^{j\omega }\right)$ 采样所得结果可表示为

$$G\left(e^{j\omega }\right)=F\left(e^{j\omega }\right)p\left(\omega \right)=\sum _{k=-\infty }^{\infty }F\left(e^{j\frac{2\pi k}{M}}\right)\delta \left(\omega -\frac{2\pi k}{M}\right).\text{\hspace{1em}(4.5)}$$

因为 $G\left(e^{j\omega }\right)$ 是由狄拉克函数定义的，在推导其离散时间傅里叶反变换时在积分边界上会存在一点问题。因为 $G\left(e^{j\omega }\right)$ 是以 $2\pi$ 为周期的，且反变换的积分周期位置不影响结果，所以我们不妨在 $\omega \in \left(0^{-},2{\pi }^{-}\right)$ 上进行积分，即积分边界是通过左侧极限来定义的，于是有

$$\begin{array}{rl}g\left[m\right]& =\frac{1}{2\pi }\underset{0^{-}}{\overset{2{\pi }^{-}}{\int }}{F}_p\left(e^{j\omega }\right)e^{j\omega m}d\omega \\ & =\frac{1}{2\pi }\underset{0^{-}}{\overset{2{\pi }^{-}}{\int }}\left[\sum _{k=-\infty }^{\infty }F\left(e^{j\frac{2\pi k}{M}}\right)\delta \left(\omega -\frac{2\pi k}{M}\right)\right]e^{j\omega m}d\omega \\ & =\frac{1}{2\pi }\underset{0^{-}}{\overset{2{\pi }^{-}}{\int }}\left[\sum _{k=0}^{M-1}F\left(e^{j\frac{2\pi k}{M}}\right)\delta \left(\omega -\frac{2\pi k}{M}\right)\right]e^{j\omega m}d\omega \\ & =\frac{1}{2\pi }\sum _{k=0}^{M-1}F\left(e^{j\frac{2\pi k}{M}}\right)\underset{0^{-}}{\overset{2{\pi }^{-}}{\int }}\delta \left(\omega -\frac{2\pi k}{M}\right)e^{j\omega m}d\omega \\ & =\frac{1}{2\pi }\sum _{k=0}^{M-1}F\left(e^{j\frac{2\pi k}{M}}\right)e^{j\frac{2\pi km}{M}}.\end{array}\text{\hspace{1em}(4.6)}$$

可以发现，式 (4.6) 所表示的离散时间傅里叶反变换成为了一个只涉及有限项 $M$ 项的求和运算。为了建立起 $G\left(e^{j\omega }\right)$ 的离散傅里叶反变换序列 $g\left[m\right]$ 与原始序列 $f\left[n\right]$ 的联系，将 $F\left(e^{j\frac{2\pi k}{M}}\right)$ 代入式 (4.6) 可得

$$\begin{array}{rl}g\left[m\right]& =\frac{1}{2\pi }\sum _{k=0}^{M-1}F\left(e^{j\frac{2\pi k}{M}}\right)e^{j\frac{2\pi km}{M}}\\ & =\frac{1}{2\pi }\sum _{k=0}^{M-1}\left(\sum _{n=-\infty }^{\infty }f\left[n\right]e^{-j\frac{2\pi k}{M}n}\right)e^{j\frac{2\pi km}{M}}\\ & =\frac{1}{2\pi }\sum _{n=-\infty }^{\infty }f\left[n\right]\left(\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(m-n\right)}\right).\end{array}\text{\hspace{1em}(4.7)}$$

当 $m-n=pM$，其中 $p$ 为整数时，可得

$$\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(m-n\right)}=\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}pM}=M.\text{\hspace{1em}(4.8)}$$

而当 $m-n=pM+q$，其中 $p,q$ 为整数，且 $1⩽q⩽M-1$ 时，可得

$$\begin{array}{c}\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(m-n\right)}=\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(pM+q\right)}=\sum _{k=0}^{M-1}{\left(e^{j\frac{2\pi q}{M}}\right)}^k\\ =\frac{1-{\left(e^{j\frac{2\pi q}{M}}\right)}^M}{1-e^{j\frac{2\pi q}{M}}}=\frac{1-e^{j2\pi q}}{1-e^{j\frac{2\pi q}{M}}}=0.\end{array}\text{\hspace{1em}(4.9)}$$

于是可得恒等式

$$\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(m-n\right)}=\{\begin{array}{cc}M,& m-n=pM,p\in \mathbb{Z},\\ 0,& others.\end{array}\text{\hspace{1em}(4.10)}$$

或者利用冲激序列的定义写为

$$\sum _{k=0}^{M-1}{e}^{j\frac{2\pi k}{M}\left(m-n\right)}=M\sum _{p=-\infty }^{\infty }\delta \left[m-n+pM\right].\text{\hspace{1em}(4.11)}$$

因此，式 (4.7) 可改写为

$$\begin{array}{rl}g\left[m\right]& =\frac{1}{2\pi }\sum _{k=0}^{M-1}F\left(e^{j\frac{2\pi k}{M}}\right)e^{j\frac{2\pi km}{M}}\\ & =\frac{M}{2\pi }\sum _{n=-\infty }^{\infty }f\left[n\right]\sum _{p=-\infty }^{\infty }\delta \left[m-n+pM\right]\\ & =\frac{M}{2\pi }\sum _{p=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }f\left[n\right]\delta \left[m-n+pM\right]\\ & =\frac{M}{2\pi }\sum _{p=-\infty }^{\infty }f\left[m+pM\right].\end{array}\text{\hspace{1em}(4.12)}$$

这样我们就证明了，对序列 $f\left[n\right]$ 的离散时间傅里叶变换 $F\left(e^{j\omega }\right)$ 以 $2\pi /M$ 为间隔进行采样，会导致序列 $f\left[n\right]$ 以 $M$ 为周期进行延拓与叠加，从而获得一个周期序列 $g\left[m\right]$。因此，$f\left[n\right]$ 的有效长度与 $M$ 的选择就成了我们需要讨论的对象。

在现实中我们所采集到的序列 $f\left[n\right]$ 的长度都是有限的，不妨记其为 $N$。由于 $g\left[m\right]$ 是 $f\left[n\right]$ 经过周期延拓与叠加的结果，当 $N⩽M$ 时，$f\left[n\right]$ 的有效数据不会发生重叠，因此可以截取与缩放完整地从 $g\left[m\right]$ 中恢复；而当 $N>M$ 时，$f\left[n\right]$ 部分有限数据发生重叠，因而不能从 $g\left[m\right]$ 中完整地恢复。我们之所以强调要从 $g\left[m\right]$ 中恢复 $f\left[n\right]$，是因为从式 (4.12) 中可知 $g\left[m\right]$ 可以通过有限项的求和运算求得，而无需任何积分运算，这对于计算机处理是至关重要的。总之，由于有限长序列 $f\left[n\right]$ 的长度 $N$ 通常是已知的，不妨令 $M=N$ 即可。综合式 (4.3) 与 (4.12)，即可引出我们所要介绍的离散傅里叶变换。

对于一长度为 $N$ 的有限长序列 $f\left[n\right]$，其中 $n=0,1,...,N-1$，定义其 $N$ 点离散傅里叶变换 (Discrete FT, DFT) 与反变换为

F [ k ] = F { f [ n ] } = ∑ n = 0 N − 1 f [ n ] e − j 2 π k n N , f [ n ] = F − 1 { F [ k ] } = 1 N ∑ k = 0 N − 1 F [ k ] e j 2 π k n N . (4.13) \begin{gathered} F\left[ k \right] = \mathcal{F}\left\{ {f\left[ n \right]} \right\} = \sum\limits_{n = 0}^{N - 1} {f\left[ n \right]{e^{ - j\tfrac{{2\pi kn}}{N}}}} , \\ f\left[ n \right] = {\mathcal{F}^{ - 1}}\left\{ {F\left[ k \right]} \right\} = \frac{1}{N}\sum\limits_{k = 0}^{N - 1} {F\left[ k \right]{e^{j\tfrac{{2\pi kn}}{N}}}} . \\ \end{gathered} \tag{4.13} F[k]=F{f[n]}=n=0∑N−1​f[n]e−jN2πkn​,f[n]=F−1{F[k]}=N1​k=0∑N−1​F[k]ejN2πkn​.​(4.13)

其证明可参考前述讨论，这里不再赘述。因为离散傅里叶变换的正变换与反变换都可以表示为有限项的求和运算，十分适合计算机与更广义的数字电路进行处理，所以其目前的应用非常广泛，不管是语音、图像还是股市、天气等各种数据的处理都可看见其身影。实际上，离散傅里叶变换更让人熟悉的名字是快速傅里叶变换 (Fast FT, FFT)，两者内在原理是相同的，但是FFT利用DFT的性质，合并了部分运算，将 DFT 的运算复杂度从 $O\left(N^2\right)$ 降至 $O\left(N\log N\right)$，而正是这种改进让 DFT 能够真正地在各种需要保证实时性的应用中变得实用。

5 总结

在前面的章节中，我们从上而下地讨论了傅里叶级数 (FS)，连续时间傅里叶变换 (CTFT)，离散时间傅里叶变换 (DTFT) 与离散傅里叶变换 (DFT)。从它们的名字也可以看出，这种区分主要在于从时域与频域周期性，以及时域与频域的连续性。为了更直观地看出它们的区别，这里做个简单的总结，四者关系如图 1 所示。

傅里叶级数 (FS)：

定义：

$$\begin{array}{c}F\left[n\right]=\frac{1}{T}\underset{-T/2}{\overset{T/2}{\int }}f\left(t\right)e^{-jn{\omega }_{0}t}dt,\\ f\left(t\right)=\sum _{n=-\infty }^{\infty }{F}_n{e}^{jn{\omega }_{0}t}.\end{array}\text{\hspace{1em}(5.1)}$$

其中 $f\left(t\right)$ 为以 $T=2\pi /{\omega }_0$ 的周期连续时间函数，$F\left[n\right]$ 为离散序列。

连续时间傅里叶变换 (CTFT)：

定义：

$$\begin{array}{c}F\left(j\omega \right)=\underset{-\infty }{\overset{\infty }{\int }}f\left(t\right)e^{-j\omega t}dt\\ f\left(t\right)=\frac{1}{2\pi }\underset{-\infty }{\overset{\infty }{\int }}F\left(j\omega \right)e^{j\omega t}d\omega .\end{array}\text{\hspace{1em}(5.2)}$$

其中$f\left(t\right)$为非周期连续时间函数，$F\left(j\omega \right)$为非周期连续频率函数。

离散时间傅里叶变换 (DTFT)：

定义：

$$\begin{array}{c}F\left(e^{j\omega }\right)=\sum _{n=-\infty }^{\infty }f\left[n\right]e^{-j\omega n}.\\ f\left[n\right]=\frac{1}{2\pi }\underset{-\pi }{\overset{\pi }{\int }}F\left(e^{j\omega }\right)e^{j\omega n}d\omega .\end{array}\text{\hspace{1em}(5.3)}$$

其中 $f\left[n\right]$ 为无限长离散序列，$F\left(e^{j\omega }\right)$ 为以 $2\pi$ 为周期的连续频率函数。

离散傅里叶变换 (DFT)：

定义：

$$\begin{array}{c}F\left[k\right]=\sum _{n=0}^{N-1}f\left[n\right]e^{-j\frac{2\pi kn}{N}},\\ f\left[n\right]=\frac{1}{N}\sum _{k=0}^{N-1}F\left[k\right]e^{j\frac{2\pi kn}{N}}.\end{array}\text{\hspace{1em}(5.4)}$$

其中 $f\left[n\right]$ 与 $F\left[k\right]$ 为长度为 $N$ 的有限长离散序列，且两者可以 $N$ 为周期进行周期延拓，改变求和区间不影响结果。

图1 四种傅里叶变换形式的联系

由于这四种变换都是在 $e^{jx}$ 即复平面的单位圆上进行分析的，它们是否能够收敛需要满足狄利克雷条件，即绝对可积与绝对可和等，这通常是比较苛刻的。例如信号与系统中最常见的单位阶跃信号明显不满足绝对可积的条件，因此我们无法求得其傅里叶变换。为此，通过加入缩放项，我们可以选择在 $e^{x+jy}$ 即复平面中不限大小的圆上进行讨论，从而引出了对连续时间信号的拉普拉斯变换与对离散时间信号的 z 变换，这两种变换能够提供更大的收敛域，从而能够分析更多类型的信号，包括前面所述的单位阶跃信号，它们是对线性移不变系统的稳定性与延时等性质分析的强有力工具。

实际上，对离散傅里叶变换进一步讨论还可以发展出很多非常实用的变换形式。例如，通过对有限长离散实数序列进行偶对称延拓，可以获得只包含实数余弦函数运算的离散余弦变换。同理，通过对有限长离散实数序列进行奇对称延拓，可以获得只包含实数正弦函数运算的离散正弦变换。另外，离散傅里叶变换主要适用于平稳信号，即信号的频谱分布不随时间而改变，然而不少场景下这种假设是不成立的，例如通信中的调频会使载波的频率随时间而改变，而人说话时声音中的频谱分布也会随着单词的改变而发生变化。除此以外，对于需要持续运行的系统，由于信号的长度是无限的，对信号整体的离散傅里叶变换分析也是不可能的。基于以上原因，我们又引出了短时傅里叶变换与小波变换，这两种变换对于非平稳信号的分析与处理十分有用。