主题
Search

离散傅里叶变换


连续傅里叶变换定义为

f(nu)=F_t[f(t)](nu)
(1)
=int_(-infty)^inftyf(t)e^(-2piinut)dt.
(2)

现在考虑推广到离散函数的情况,f(t)->f(t_k) 通过令 f_k=f(t_k),其中 t_k=kDelta,其中 k=0, ..., N-1。写出来得到离散傅里叶变换 F_n=F_k[{f_k}_(k=0)^(N-1)](n)

 F_n=sum_(k=0)^(N-1)f_ke^(-2piink/N).
(3)

逆变换 f_k=F_n^(-1)[{F_n}_(n=0)^(N-1)](k) 然后是

 f_k=1/Nsum_(n=0)^(N-1)F_ne^(2piikn/N).
(4)

离散傅里叶变换 (DFT) 非常有用,因为它们揭示了输入数据中的周期性以及任何周期性成分的相对强度。然而,在解释离散傅里叶变换时存在一些微妙之处。一般来说,实数序列的离散傅里叶变换将是相同长度的复数序列。特别地,如果 f_k 是实数,那么 F_(N-n)F_n 通过下式相关

 F_(N-n)=F^__n,
(5)

对于 n=0, 1, ..., N-1,其中 z^_ 表示复共轭。这意味着对于实数数据,分量 F_0 始终是实数。

由于上述关系,周期函数将在两个位置而不是一个位置包含变换后的峰值。发生这种情况是因为输入数据的周期被分成“正”和“负”频率复数分量。

快速傅里叶变换是一种特别高效的算法,用于执行包含特定点数的样本的离散傅里叶变换。

可能影响离散傅里叶变换的两种主要误差类型是:混叠泄漏

DiscreteFourierTransform

上面的图显示了函数 f(x)=sinx (左图)和 f(x)=sinx+sin(3x)/2 (右图)在两个周期内采样 50 次的离散傅里叶变换的实部(红色)、虚部(蓝色)和复模(绿色)。在左图中,左侧和右侧的对称尖峰是单个正弦波的“正”和“负”频率分量。类似地,在右图中,有两对尖峰,较大的绿色尖峰对应于较低频率的较强分量 sinx,较小的绿色尖峰对应于较高频率的较弱分量。离散傅里叶变换的复模的适当缩放图通常称为功率谱

Wolfram 语言复数列表 l 实现了离散傅里叶变换,如下所示Fourier[list]。

离散傅里叶变换是Z 变换的一个特例。

可以使用快速傅里叶变换有效地计算离散傅里叶变换。

在离散傅里叶变换的指数中添加一个额外的因子 b 得到所谓的(线性)分数傅里叶变换

DiscreteFourierTransform2D

离散傅里叶变换也可以推广到二维和更多维度。例如,上面的图显示了函数 sin(x+y) 的二维离散傅里叶变换的复模


另请参阅

混叠, 快速傅里叶变换, 傅里叶变换, 分数傅里叶变换, 哈特利变换, 泄漏, 功率谱, Winograd 变换, Z 变换

使用 Wolfram|Alpha 探索

参考文献

Arfken, G. "Discrete Orthogonality--Discrete Fourier Transform." §14.6 in Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 787-792, 1985.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. "Fourier Transform of Discretely Sampled Data." §12.1 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 494-498, 1989.Roberts, S. Lecture 7-The Discrete Fourier Transform. pp. 82-96. http://www.robots.ox.ac.uk/~sjrob/Teaching/SP/l7.pdf.

在 Wolfram|Alpha 上引用

离散傅里叶变换

请引用本文为

Weisstein, Eric W. "离散傅里叶变换。" 来自 MathWorld--一个 Wolfram Web 资源。 https://mathworld.net.cn/DiscreteFourierTransform.html

主题分类