MATLAB的一个FFT程序

FFT信号流图：

 

程序实现是这样：

 

程序流程如下图：

 

 

首先进行位逆转，其实很简单，就是把二进制的位逆转过来：

Matlab的位逆转程序：

function a=bitreverse(Nbit, num)

%Nbit = 4;

%num = 8;

a = 0;

b = bitshift(1,Nbit-1);

for i = 1:Nbit;

if((bitand(num,1)) == 1)

a = bitor(a,b);

end

num = bitshift(num,-1);

b = bitshift(b,-1);

end;

 

说明：Nbit是逆转位是几位，num是逆转的数即变量。

 

 

三个循环，第一个循环是进行N阶的FFT运算

第二个循环其实就是，每一阶FFT的时候，有多少组DFT对象，拿8点来说，第一阶的时候，有4组DFT对象，到了第二阶，就有