关于DFT

DFT的笔记

连续函数傅立叶变换和傅立叶级数

傅立叶级数

在满足一定条件的情况下，周期为$T$的函数可以表示为三角函数张成的空间里的向量，表示为

$$f(t)=\sum _{n=-\infty }^{\infty }{c}_n{e}^{i\frac{2\pi n}{T}t}$$

其中，基函数是正交的，那么，系数可以用$f(t)$与$e^{-i\frac{2\pi n}{T}t}$的内积得到，内积是定义在一个周期上的积分除以周期。

系数为$c_n=\frac{1}{T}{\int }_{-T/2}^{T/2}f(t)e^{-i\frac{2\pi n}{T}t}dt$

$c_n$的值相当于频率为$\frac{2\pi n}{T}$的三角函数的强度。

傅立叶变换

傅立叶级数适用于周期函数，而对于非周期函数来说，是无法适用的，不过如果将周期视为无穷大，可由傅立叶级数推导出傅立叶变换，推导过程比较复杂，这里给出形式。

$F(\mu )={\int }_{-\infty }^{\infty }f(t)e^{-i2\pi \mu t}dt$

$f(t)={\int }_{-\infty }^{\infty }F(\mu )e^{i2\pi \mu t}d\mu$

前者就是函数的频谱，相当于傅立叶级数中的$c_n$，只不过是连续的。

带限函数采样的傅立叶变换

考虑一个带限连续函数经采样后的傅立叶变换。

设原函数$f(t)$，傅立叶变换为$F(\mu )$，采样周期为$\Delta T$，冲击串为$s_{\Delta T}(t)$，傅立叶变换为$S_{\Delta T}(\mu )$

采样后的函数为：

$f(t){\sum }_k\delta (t-k\Delta T)=f(t)s_{\Delta T}(t)$

其傅立叶变换为：

$$\begin{array}{rl}F(\mu )\otimes S_{\Delta T}(\mu )& ={\int }_{-\infty }^{\infty }F(m)S(\mu -m)dm\\ & =\frac{1}{\Delta T}{\int }_{-\infty }^{\infty }F(m)\sum _k\delta (\mu -m-k/\Delta T)dm\\ & =\frac{1}{\Delta T}\sum _k{\int }_{-\infty }^{\infty }F(m)\delta (\mu -m-k/\Delta T)dm\\ & =\frac{1}{\Delta T}\sum _{k}F(\mu -k/\Delta T)\end{array}$$

这说明，采样后的函数傅立叶变换是其傅立叶变换以$1/\Delta T$为周期的延拓。

当然，具体形式也可以求得：

$$\begin{array}{rl}{\int }_{-\infty }^{\infty }\sum _{k}f(t)\delta (t-k\Delta T)e^{-i2\pi \mu t}dt& =\sum _k{\int }_{-\infty }^{\infty }f(t)\delta (t-k\Delta T)e^{-i2\pi \mu t}dt\\ & =\sum _{k}f(k\Delta T)e^{-i2\pi \mu k\Delta T}\end{array}$$

虽然$f(k\Delta T)$是离散的，但关于频率$\mu$，这是一个连续函数，也是周期的。

离散傅立叶变换

表征采样后傅立叶变换在一个周期的值就可以得到完整的信息，在实际中这一个周期的值也是连续的，想表征这个周期的值，同样要对这个频域的周期进行采样。

取0到$1/\Delta T$这个周期，等间距采样M个样本，那么我们的采样值可以表示为$\frac{m}{M\Delta T}$

那么，采样得到的值为

$$\sum _{k}f(k\Delta T)e^{-i2\pi (\frac{m}{M\Delta T}k\Delta T)}=\sum _{k}f(k\Delta T)e^{-i2\pi mk/M}$$

这就是离散傅立叶变换，通常，$f$的样本数是有限的也为M个，这样离散傅立叶变换表示为：

$$F_m=\sum _{k=0}^{M-1}{f}_k{e}^{-i2\pi mk/M}$$

离散傅立叶反变换为：

$$f_k=\frac{1}{M}\sum _{m=0}^{M-1}{F}_m{e}^{i2\pi mk/M}$$

这两个式子不明确依赖$\Delta T$，适用于所有均匀采样。

由此可知，DFT输入为一个M个数据的序列，输出也为一个M个数据的序列，既然变换过程不依赖于采样频率，那怎么表示输入输出的单位呢？实际上，输入的点的横坐标就是$k\Delta T$，输出点的横坐标就是$m/\Delta T$，只需要记得我们的采样频率，就可以给输入输出做出标定。