基4fft原理及matlab实现.docVIP

基4FFT原理及MATLAB实现 一．时域抽取法基4FFT基本原理： 有限长序列的N点DFT为： k=0,1,2…,N-1 设序列长度为N，且满足,M为自然数，可把按下列方法分解为4个长为N/4的子序列： r=0,1,…,N/4-1 r=0,1,…,N/4-1 r=0,1,…,N/4-1 r=0,1,…,N/4-1 则的DFT为X(k)=+ = = 因为： 所以：原式可表示为： X(k)= = k=0,1,…,N/1 其中：分别为的N/4点DFT，即： = = = = 由于均以N/4为周期，且有： ，，， 所以： X(k)= k=0,1,…,N/4-1 X(k+ N/4)= k=0,1,…,N/4-1 X(k+ 2N/4)= k=0,1,…,N/4-1 X(k+ 3N/4)= k=0,1,…,N/4-1 二．运算规律及编程思想： 1.按照上述分解法，再对进行反复分解，直到每个序列的长度等于4为止，这个4点的DFT的数学表达式为： 2.旋转因子与运算级数的关系（参考《数字信号处理（第三版）》西安电子科技大学出版社第1