利用fft计算时域卷积重叠法保留_离散傅立叶变换DTFT、DFT和FFT在工程与数学结合的通俗理解...

1.离散时间傅里叶变换DTFT

何为DTFT？就是对连续时间非周期信号进行抽样（乘积），得到的离散时间非周期信号再求傅里叶变换的过程就是DTFT。其实等同于信号频谱与脉冲信号频谱的卷积，这样得到的就是连续且周期DTFT频谱（把频谱进行了周期延拓），如下图所示。

2.离散傅里叶变换DFT

DFT其实是在DTFT的频谱上进行采样，使得频谱离散化（频谱和抽样信号进行乘积），这也等同于将时域信号进行了卷积，与采样信号卷积就是周期延拓，得到周期离散时间信号。

那么在数字信号角度分析，首先从DFT公式入手：

首先x(n)的最大周期我们可以确定下来，也就是数字处理时所截取的那一段周期，由上面公式可知其采样点的数目为N，最大周期T=1（角频率w=2π，T=w/2π），即可得知x(n)的最小频率为1（除零外），且频率最小间隔也为1；

为了简单化理解，我们把指数部分通过欧拉变回正余弦表示，且幅度相位都不要了，那么该公式就可以解释为x(n)的N个离散值与sin(2πkn)被分成了N份的值进行相乘：

当k=1时，有

此时n取0-N，就是把2π分成了N份，每一个都对应的一个正弦值，再将每个值分别与x(n)对应值相乘后再累加，最后得到的就是这个频率2π/2π=1时的分量！

同理可得，X(2)、X(3).....、X(N)都是按照这个方式进行计算，这样就可以得到每个频率的分量大小，这也就是工程上DFT变换的本质。

看了这么多，似乎有点费脑子，要不休息一会，打把王者？？

我们继续，经过DFT实例计算后我们会发现，频谱是中间对称的，这又是咋回事呢？

其实这个不难理解，把每一个k对应的正弦在坐标轴里面画出来后，会发现其实就是在圆上取点，如下情况：

k=2时，θ=0、4π/N、...... 、4π（N-1）/N

.....

k=N-2时，θ=0、（2πN-4π）/N、（4πN-8π）/N、...... 、(2πN*(N-1)-4π*(N-1))/N

k=N-2进行三角变换后：

θ=0、4π/N、8π/N、...... 、4π（N-1）/N

这下就可以发现，k=N-2时与k=2时相等，这是因为他们在圆上取点是相同的，因此有中心对称的值都是相等的。既然所有的角度相等了，那么正弦取值也相等，再与x(n)相乘求和必定也相等，因此频谱是中间对称的，因此进行傅里叶变换后只采用其一半的频谱进行计算。

简单总结一下：

1.在DFT中，是将2pi、4pi、6pi、...、2kpi分别分成了N份，每一个都是对应k的正余弦的最小变化步长。

2.当被采集信号的时间长度被确定时，其最小频率和最小频率间隔也被确定，即为1/T；那么每一个三角函数的间隔也被确定，即原有最小频率的整数倍k；

3.最小角频率以及最小角频率的整数倍角频率的三角函数可以表示出待求函数x(n)；（注意频率和角频率的区别）

4.离散傅里叶变换得到的频谱是中心对称的频谱，只取前N/2分析即可

那么为何不同频率的三角函数与原函数相乘求和就可以得到每一个频率的分量，怎么区别开每一个频率的，不同频率间不会互相干扰吗？

这就涉及到了傅里叶级数完备正交基的概念，可以参考之前写的文章：

舞曜：傅里叶级数和傅里叶变换的巧妙之处

3.快速傅里叶变换FFT

简单来说FFT是DFT的优化算法，是在计算机出来之后，减少DFT的计算复杂度的算法，他的核心精髓是利用系数

的特性对DFT进行简化：

（1） 对称性

（2）周期性

（3）可约性

具体FFT推导咱不说，咱们关注的是结果，是工程什么的应用，故在进行一个N个采样点时，所需计算量如下所示：

1) DFT常规计算时

计算一个X(K)的值需要运算量

复数乘法次数：N ；复数加法次数：N-1

DFT计算全部X(K)的值需要运算量

复数乘法次数：N^2 ；复数加法次数：N*(N-1)

2) 在使用FFT计算时

DFT计算全部X(K)的值需要运算量

复数乘法次数：

；复数加法次数：

要注意的是，在计算机中，一次复数乘法需要通过4次实数乘法和2次实数加法实现。因此FFT在很大程度上减少了计算量，其实无论DFT就是把数据以某种方程变成另外一种数据，而FFT就是把这个方程的计算量进行了优化简化，仅此而已。

在FFT中，值得借鉴的是FFT优化DFT的方法思路，利用了系数

的特征，深究其数学特性的意义，对以后的算法设计优化有很好的参考价值，也对未来的科研道路有所帮助。

以上乃是本人小生对离散傅里叶变换在工程上使用和数学公式所结合的通俗理解，希望能够和大家共同探讨，也愿能够对大家的学习和科研有所帮助，与知友、智友共勉之。