撒旦范德萨.docVIP

时间抽取（DIT）基2FFT算法的MATLAB实现 ? 目的：利用MATLAB实现时间抽取（DIT）基2FFT算法。 在阅读程序之前，最好首先阅读“时间抽取（DIT）基2FFT算法的讨论”[REF1]，本讨论对算法的分解非常到位，本程序既是基于对此讨论的理解。本程序包括两部分：序列排序和蝶形运算。其中，蝶形运算包括三级循环，分别为：“级”的循环，“组”的循环和“组内”循环。“级”中包含“组”，“组”中包含“组内”。可以留言讨论，下面附程序： %DIT_FFT_2 @ 091211 function [Xk]=DIT_FFT_2(xn,N); %序列排序开始n1=fliplr(dec2bin([0:N-1]));%码位倒置步骤1：将码位转换为二进制，再进行倒序n2=[bin2dec(n1)];%码位倒置步骤2：将码位转换为十进制后翻转for i=1:N??? Xk(i)=xn(n2(i)+1);%根据码位倒置的结果，将xn重新排序,存入Xk中end %蝶形运算开始M=log2(N);%“级”的数量for m=0:M-1 %“级”循环开始??? Num_of_Group=N/2^(m+1);%每一级中组的个数??? Interval_of_Group=2^(m+1);%每一级中组与组之间的间距??? Interval_of_Unit=2^m;%每一组中相关运算单元之间的间距??? Cycle_Count=2^m-1;%每一组内的循环次数??? Wn=exp(-j*2*pi/Interval_of_Group);%旋转因子??? for g=1:Num_of_Group %“组”循环开始??????? Interval_1=(g-1)*Interval_of_Group;%第g组中蝶形运算变量1的偏移量??????? Interval_2=(g-1)*Interval_of_Group+Interval_of_Unit;%第g组中蝶形运算变量2的偏移量??????? for r=0:Cycle_Count;%“组内”循环开始??????????? k=r+1;%“组内”序列的下标??????????? Xk(k+Interval_1)=Xk(k+Interval_1)+Wn^r*Xk(k+Interval_2);%第m级，第g组的蝶形运算式1??????????? Xk(k+Interval_2)=Xk(k+Interval_1)-Wn^r*Xk(k+Interval_2)-Wn^r*Xk(k+Interval_2);%第m级，第g组的蝶形运算式2，注：1和2为同址运算??????? end??? endend