快速福利叶变换原理与程序代码.docVIP

快速傅立叶变换（时域抽取基二fft） 编程思想 根据快速傅立叶变换的信号流图可知，可将整个过程中所有的数据组成一个二维数组data(N,M+1)，数组共有N行，M+1列（傅立叶变换分为M=log2(N)级，再加上第一级倒序数组输入，则共有M+1列）。除第一列单独赋值外，其余列则按照共同的规律来赋值。这里则详细说明其的规律性。 （1）对于第k列（k1）： 可分为2^(M+1-k)个计算单元，各计算单元间相互独立进行离散傅里叶变换。 （2）对于第k列的第Mblock个计算单元 该单元总共含有2^(k-1)个数，数据的起始项是：data( (1+(Mblock-1)*2^(k-1)),k )，结束项是data( ( (Mblock)*2^(k-1)),k ) （3）每个计算单元又可分为2^(k-2)个计算组，即蝶形单元 对于每个组可以运用蝶形运算 则第k列的第Mblock个计算单元的第grop个计算组的两个元素的下表分别为： (Mblock-1)*2^(k-1)+grop与 (Mblock-1)*2^(k-1)+grop+2^(k-2) 碟形单元所出现的系数为： 然后根据蝶形运算便可进行程序的编写。 通过上述方法在MATLAB下编程，大概在35行左右，还包括倒序的程序在里面。 MATLAB程序清单： %DIT(时间抽取的基二FFT算法程序) %程序输入为x(n)，输出为y(n),要求n是2的整数次幂 function y=b2fft(x) N=length(x); L=log2(N); rev_x=zeros(1,N); for m=0:N-1 bit_str=dec2bin(m); %将i变成2进制数，为倒序操作做准备 dec2bin()将十进制数变成二进制数，结果为一字符串 if length(bit_str)L Z=zeros(1,(L-length(bit_str))); %保证倒序操作的成功，要求二进制位数始终是 N Z_str=num2str(Z); %将数字变为字符串 bit_str=[Z_str,bit_str]; end it_str=fliplr(bit_str); %fliplr其作用是将一个数组的行倒置，而列的顺序则不变，倒置的意思是倒序显示 rev_x(m+1)=x(bin2dec(it_str)+1); end data=zeros(N,L+1); %进行第一级赋值 data(:,1)=rev_x; %进行第K级赋值 for k=2:L+1 %多少级（第一级为倒序后的数组） for Mblock=1:2^(L+1-k) %每级多少个块 for grop=1:2^(k-2) %每个块分成多少个组 coeff=cos(2*pi*(grop-1)/(2^(k-1)))-sin(2*pi*(grop-1)/(2^(k-1)))*i; data((Mblock-1)*2^(k-1)+grop,k)=data((Mblock-1)*2^(k-1)+grop,k-1)+coeff*data((Mblock-1)*2^(k-1)+grop+2^(k-2),k-1); data((Mblock-1)*2^(k-1)+grop+2^(k-2),k)=data((Mblock-1)*2^(k-1)+grop,k-1)-coeff*data((Mblock-1)*2^(k-1)+grop+2^(k-2),k-1); end end end %输出L+1级数据，为傅立叶变换后的数据 y=data(:,L+1); 我在时域,你在频域,需要经过傅立叶变换,才能发现你的美丽,我把爱的语言调制到星座,通过伪随机序列,载波到你的频率,并波束赋形到你的接收阵列矩阵,你说我的爱噪声太大,经过层层滤波,原来发现,那是在宇宙开始的时候,我发给你的爱的微波背景辐射衍生出来的傅立叶级数版本傅立叶研究物理提出了三角级数串两百年前粗略的论断催生傅立叶变换不朽的缠绵我在z平面前凝视系统函数的极点祈祷你的收敛域回到我的单位圆卷积采样滤波重建哪怕坐标已变换在抽取插值间守护你不失真的容颜拜读奥本海姆鸿篇以线性系统之名许愿：想像二阶无阻尼振荡到永远当信号只剩下离散的语言DTFT就成了无尽的思念我给你的藏在频谱间在