1周学FFT——第1天离散傅里叶变换

作者: 理耳兔子 | 来源:发表于2020-04-13 00:27 被阅读0次

离散傅里叶变换（Discrete Fourier Transform，DFT）常用来处理采样得到的电压、电流信号，以便获得其频域特性。

原始信号由 $x(t)$ 描述。

采样动作由单位脉冲序列描述：

$\delta (n) = \sum_{m=-\infty}^{\infty}{\delta(n-m\Delta t)} \quad m=0, \pm 1, \pm 2, ...$

式中 $\Delta t$ 为采样间隔。

采样之后得到离散序列[从玉良，2009，p12]：

$x_{\infty}(n) = x(t)\delta(n) = \sum_{m=-\infty}^{\infty}{x(m\Delta t)\delta(n-m\Delta t)} \quad m=0, \pm 1, \pm 2, ...$

连续信号、单位采样序列、离散序列示意图

截断之后得到有限长（N点）离散序列：

$x(n) = x(t)\delta(n) \quad n = 0, 1, 2, ..., N-1$

截断后的离散序列（10点）

$x(n)$ 的离散傅里叶变换是一个长为N的频域有限长序列（ $0\le k \le N-1$ ），其正变换为：

$X(k) = \text{DFT}[x(n)] = \sum_{n=0}^{N-1}{x(n)e^{-j\frac{2\pi}{N}nk}} = \sum_{n=0}^{N-1}{x(n)W_N^{nk}} \quad k=0,1, ..., N-1$

其反变换为：

$x(n) = \text{IDFT}[X(k)] = \frac{1}{N}\sum_{k=0}^{N-1}{X(k)e^{j\frac{2\pi}{N}nk}} = \sum_{k=0}^{N-1}{X(k)W_N^{-nk}} \quad k=0,1, ..., N-1$

其中 $W_N = e^{-j\frac{2\pi}{N}}$ 。

习题：

假设一个连续信号为 $x(t) = \text{sin}(2\pi t) + 2\text{sin}(4\pi t)$ 。尝试编写matlab脚本程序，对 $x(t)$ 进行离散化采样，采样频率8Hz，输出序列为 $x(n)$ ，并绘制 $x(t)$ 和 $x(n)$ 。

快速傅里叶变换FFT(Fast Fourier Transfor
快速傅里叶变换(Fast Fourier Transform,FFT)用来计算离散傅里叶变换（Discrete F...
Python科学计算——复杂信号FFT
FFT (Fast Fourier Transform, 快速傅里叶变换) 是离散傅里叶变换的快速算法，也是数字信...
NumPy API（十四）——离散傅里叶变换
离散傅里叶变换标准的 FFTs fft(a[, n, axis, norm]) 计算一维离散傅立叶变换。 iff...
快速傅里叶变换和离散傅里叶变换
快速傅里叶变换（FFT）离散傅里叶变换（DFT）基础理论是傅里叶变换的分离形式，和采样定理（香菜定理）采样定...
基2-FFT简记
快速傅里叶变换（Fast Fourier Transform，FFT）是一种可在时间内完成的离散傅里叶变换（Dis...
离散傅里叶变换和FFT
1周学FFT——第1天离散傅里叶变换
离散傅里叶变换（Discrete Fourier Transform，DFT）常用来处理采样得到的电压、电流信号，...
2018-07-25 傅里叶变化和高通滤波器学习笔记
1 图片的傅里叶变换： f = np.fft.fft2(norm_image) fshift = np.fft...
离散傅里叶变换 DFT
离散傅里叶变换 DFT 周期离散信号（离散时间傅里叶变换：非周期，离散；傅里叶变换：非周期，连续；傅里叶级数：...
跃马扬鞭，看量子计算机如何加速傅里叶变换
快速傅里叶变换快速傅里叶变换（Fast Fourier Transform, FFT）是现代生活中的幕后数字主力...