快速傅里叶变换FFT的原理及公式_精品文档资料下载.pdf
《快速傅里叶变换FFT的原理及公式_精品文档资料下载.pdf》由会员分享,可在线阅读,更多相关《快速傅里叶变换FFT的原理及公式_精品文档资料下载.pdf(3页珍藏版)》请在冰豆网上搜索。
![快速傅里叶变换FFT的原理及公式_精品文档资料下载.pdf](https://file1.bdocx.com/fileroot1/2022-10/14/94dbc6b7-b64b-4b9b-b5f2-d0087bb071f5/94dbc6b7-b64b-4b9b-b5f2-d0087bb071f51.gif)
因此需要利用离散信号x(nT)来计算信号x(t)的频谱。
有限长离散信号x(n),n=0,1,N-1的DFT定义为:
可以看出,DFT需要计算大约N2次乘法和N2次加法。
当N较大时,这个计算量是很大的。
利用WN的对称性和周期性,将N点DFT分解为两个N2点的DFT,这样两个N2点DFT总的计算量只是原来的一半,即(N2)2+(N2)2=N22,这样可以继续分解下去,将N2再分解为N4点DFT等。
对于N=2m点的DFT都可以分解为2点的DFT,这样其计算量可以减少为(N2)log2N次乘法和Nlog2N次加法。
图1为FFT与DFT-所需运算量与计算点数的关系曲线。
由图可以明显看出FFT算法的优越性。
将x(n)分解为偶数与奇数的两个序列之和,即x1(n)和x2(n)的长度都是N2,x1(n)是偶数序列,x2(n)是奇数序列,则其中X1(k)和X2(k)分别为x1(n)和x2(n)的N2点DFT。
由于X1(k)和X2(k)均以N2为周期,且WNk+N/2=-WNk,所以X(k)又可表示为:
上式的运算可以用图2表示,根据其形状称之为蝶形运算。
依此类推,经过m-1次分解,最后将N点DFT分解为N2个两点DFT。
图3为8点FFT的分解流程。
FFT算法的原理是通过许多小的更加容易进行的变换去实现大规模的变换,降低了运算要求,提高了与运算速度。
FFT不是DFT的近似运算,它们完全是等效的。
关于FFT精度的说明:
因为这个变换采用了浮点运算,因此需要足够的精度,以使在出现舍入误差时,结果中的每个组成部分的准确整数值仍是可辨认的。
为了FFT的舍入误差,应该允许增加几倍log2(log2N)位的二进制。
以256为基数、长度为N字节的数可以产生大到(256)2N阶的卷积分量,所以为了正确存储,需要16+log2N位精度,若数i是浮点尾数的二进制位数,则有条件:
如果i=24,对于任意感兴趣(N256)的N值,单精度是不合适的;
如果i=53,也就是采用双精度,则允许N大于106,相当于几百万十进制位。
所以,用FFT作大数乘法时,向量数组选用双精度类型。