DSP实验报告_matlab.docVIP

按时间抽取的基 2-FFT算法实现 实验目的 掌握按时间抽取的基 2-FFT 的原理及具体实现方法。 编程实现按时间抽取的基 2-FFT 算法。 加深理解按时间抽取 FFT 算法的特点 实验内容 编程实现序列长度为 N = 8的按时间抽取的基 2-FFT 算法。给定一个 8 点序列，采用编写的程序计算其 DFT，并与 MATLAB 中 fft 函数计算的结果相比较，以验证结果的正确性。 程序代码 disp (请输入一个 N=8的序列 ); for ii=1:8 %用户输入一个 N=8的序列 x(ii)= input([x(,num2str(ii),)=]); end nxd=bin2dec(fliplr(dec2bin([0:7])))+1; %求数列序号的倒序 y=x(nxd); %将 x的倒序序列给 y for m=1:3 %做 3级蝶形运算，每级 2^(3-m)组 Nz=2^m; %蝶形运算间隔序号数 Nz WN=exp(-1i*2*pi/Nz); %定义旋转因子 u=WN^0; fo