[信息与通信]数字信号处理资源第4章快速傅里叶变换.pptVIP

第4章　快速傅里叶变换 4.1　引言4.2　基-2时域抽选(DIT)FFT算法4.3　基-2频域抽选(DIF)FFT算法4.4　离散傅里叶反变换(IDFT)的快速算法4.5　任意基FFT算法4.6　调频z变换4.7　其他的快速计算方法 4.1　引言 离散傅里叶变换实现了频域离散化，它可以直接用来分析信号的频谱、计算滤波器的频率响应，以及实现信号通过线性系统的卷积运算等，在数字信号处理中起着极其重要的作用。但是在其付诸实际应用时却遇到了计算量过于庞大的棘手问题。 4.2　基-2时域抽选(DIT)FFT算法 4.2.1　算法原理4.2.2　运算特点4.2.3　矩阵分解表示4.2.4　编程思想4.2.5　硬件实现 4.2.1　算法原理 1.基本原理2.处理过程 1.基本原理 DFT的乘法计算次数与序列的长度有关，复数乘法为N2次，加法为N(N-1)次。如果将N分为两个N2处理，乘法次数将是原来的一半，加法次数也比原来少。 另外由WnkN的周期性和对称性可知,DFT的乘法和加法次数减少就可以提高它的运算速度。 2.处理过程 图4-1　蝶式运算的信号流图 2.处理过程 图4-2　N点DFT一次时域抽选分解 2.处理过程 图4-3　N点DFT第二次时域抽选分解 2.处理过程 图4-4　8点DIT-FFT运算流图 2.处理过程 图4-5　FFT算法与直接计算DFT所需乘法次数的比较曲线 2.处理过程 图4-6　DIT蝶形 4.2.2　运算特点 1. 同址运算2.蝶类、蝶距以及旋转因子的变化规律3.数据重排及倒序二进制数4.存储单元 1. 同址运算 由图4-4可以看出，DIT-FFT的每轮计算都是把N个存储单元中的复数，经N/2个蝶式运算变成另外N个复数，每个蝶式运算如图4-6所示，且可表达为下述的基本迭代运算。 Am+1(i)=Am(i)+Am+1(j)WkN Am+1(j)=Am(i)-Am+1(j)WkN 2.蝶类、蝶距以及旋转因子的变化规律 再次观察图4-4可以发现，第一轮只有一种类型的蝶形，蝶距为1，旋转因子为W0N。第二轮有两种类型的蝶形，蝶距为2，旋转因子为W0N，W2N。第三轮则有4种类型的蝶形，蝶距为4，系数为W0N，W1N，W2N，W3N。推广到2M的一般情况，第一轮只有一种类型的蝶形，蝶距为1，旋转因子也只有一个W0N，以后每轮的蝶类、蝶距和旋转因子个数都比前列增加一倍，到第M轮时有N/2=2M-1个蝶式类型，蝶距也为N/2，旋转因子的个数也是2M-1= N/2，且为W0N，W1N，…，WN/2-1N,第L轮旋转因子的通式为 WPN=Wj2L,j=0,1,…,(2L-1-1) 3.数据重排及倒序二进制数 由图4-4可见，对于同址运算结构，其输出结果X(k)是按自然顺序(正序)排列的，即X(0)，X(1)，X(2)，X(3)，…,X(N-1)。但输入序列x(n)的排列为x(0),x(4),x(2),x(6)，x(1),…，似乎混乱无序，但其实它是很有规律的，谓之倒序排列。 4.存储单元 由于DIT-FFT是同址运算，只需输入序列x(n)(n=0,1,…,N-1)的N个存储单元，加上系数WrN(r=0,1,…,N/2-1)的N/2个存储单元，即一共仅需112N个存储单元即可。 4.2.3　矩阵分解表示 如果把DFT用矩阵来表示，那么FFT就相当于矩阵因式分解处理。 4.2.4　编程思想 1.变址(倒序)运算2.M轮的递推计算 1.变址(倒序)运算 图4-7　雷德算法流程图 2.M轮的递推计算 (1)设N=2M，则共有M轮蝶式运算，每轮有N/2个蝶式。(2)由最后一轮向前每推进一轮，则旋转因子取后轮旋转因子中偶数序号那一半，第三轮的旋转因子为W0,W1,W2,W3，而第二轮的旋转因子为W0，W2。(3)蝶距为2L-1，L为所在的轮数，显然，每当向前推进一轮，蝶距就变成原蝶距的一半。 (3)蝶距为2L-1，L为所在的轮数，显然，每当向前推进一轮，蝶距就变成原蝶距的一半。 表4-1　顺序和倒序二进制数对照表(N=8) 顺　　序 倒　　序 顺序数I 二进制数 二进制数 倒序数J 01234567 0 0 00 0 10 1 00 1 11 0 01 0 11 1 01 1 1 0 0 01 0 00 1 01 1 00 0 11 0 10 1 11 1 1 04261537 (3)蝶距为2L-1，L为所在的轮数，显然，每当向前推进一轮，蝶距就变成原蝶距的一半。 图4-8　基2DIT-FFT流程图 4.2.5　硬件实现 1.顺序处理2.级联处理3.并行迭代处理4.阵列处理 1.顺序处理 (1)只用一个运算单元；(2)输入量、中间量、输