编考 实验 题数字信号处理实验一1120111260黄昳.pdfVIP

数字信号处理实验姓名：黄昳学号：1120111260班级专业：通信工程

实验1基2-FFT算法实现

一、实验目的

1.掌握基2-FFT的原理及具体实现方法。

2.编程实现基2-FFT算法。

3.加深理解FFT算法的特点。

二、实验设备与环境

计算机、MATLAB软件环境。

三、实验基础理论

FFT是DFT的一种快速算法，能使DFT的计算大大简化，运算时间缩短。FFT利用了

nk

W

N的三个固有特性，即对称性、周期性和可约性，将长序列的DFT分解为短序列的DFT，

合并了DFT运算中的某些项，从而减少了DFT的运算量。

FFT算法基本上可以分为两大类，即按时间抽取法和按频率抽取法。

在实现FFT算法时，要重点考虑两个问题，注意数据的读取和存储：（1）输入输出的

排序；（2）碟形运算的实现。按时间抽取算法中输入反序输出顺序，按频率抽取算法中输

入顺序输出反序；运算过程中的每一级都由N/2个碟形运算构成，每一个碟形运算单元中，

两个节点变量运算后得到的结果为下一列相同位置的节点变量，而和其他节点变量无关，

nk

W

N

可以采用原位运算，节省存储单元。另外，碟形运算中的复系数可以存储为能及时查阅

的系数表，这样可以节约计算量，但是需要N/2个复数存储器。

MATLAB中提供了用于计算FFT的函数fft，可将实验中所得到的结果与利用MATLAB

中fft函数计算的结果相比较，以此验证结果的正确性。

四、实验内容

1.编程实现序列长度为N=8的按时间抽取的基2-FFT算法。给定一个8点序列，采用

编写的程序计算其DFT，并与MATLAB中fft函数计算的结果相比较，以验证结果的正确性。

2.编程实现序列长度为N=8的按频率抽取的基2-FFT算法。给定一个8点序列，采用

编写的程序计算其DFT，并与MATLAB中fft函数计算的结果相比较，以验证结果的正确性。

3.将上述FFT程序推广到序列长度为N=2v的情况，要求利用原位运算。

五、实验结果

1.实验代码：

x=[12345678];

m=3;

N=8;

nxd=bin2dec(fliplr(dec2bin([1:N]-1,m)))+1;

y=x(nxd);

-1-

数字信号处理实验姓名：黄昳学号：1120111260班级专业：通信工程

formm=1:m

Nz=2^mm;u=1;

WN=exp(-i*2*pi/Nz);

forj=1:Nz/2

fork=j:Nz:N

kp=k+Nz/2;

t=y(kp)*u;

y(kp)=y(k)-t;

y(k)=y(k)+t;

end

u=u*WN;

end

end

y1=fft(x)

实验结果：

y1=

Columns1through4

36.0000+0.0000i-4.0000+9.6569i-4.0000+4.0000i-4.0000+1.6569i

Columns5through8

-4.0000+0.0000i