% FFT
% 本程序根据《数字信号处理教程》中提供的算法编制;
% 本程序是按时间抽选的基-2FFT 算法,输入倒位序,输出顺序;
% Made by ceng1010。
clear all;
close all;
N=2^8; % 指定序列长度 N,可修改 N
flag=1; % flag=1 时,直接给定 x(n);flag=0 时,
由程序给定
if flag==0
x=[1 2 3 4]; % 给定序列
else
x=0; % 或者通过程序给出序列 x(n)=n-1;用于验
证程序
for i=0:N-1
x(i+1)=i;
end
x
end
NV2=N/2; % 输入倒位序处理
NM1=N-1;
I=0;
J=0;
while I<NM1
if I<J
T=x(J+1);
x(J+1)=x(I+1);
x(I+1)=T;
end
K=NV2;
while K<=J
J=J-K;
K=K/2;
end
J=J+K;
I=I+1;
end
x % 倒位序后的结果
评论0