数字信号处理实验报告(1)

是周期的整倍数，这种截断效应也会小一些。

３、实验步骤及内容

（１）复习DFT的定义、性质和用DFT作谱分析的有关内容。

（２）复习FFT算法原理与编程思想。 （３）编制信号产生程序，产生以下典型信号供谱分析用： x(n)?R(n)

14?n?1,?x2(n)??8?n,?0,??4?n,?x3(n)??n?3,?0,?

0?n?34?n?7其他n0?n?34?n?7其他n

x4(n)?cosn4? 式中频率f自己选择；

x(t)?cos(8?t)?cos(16?t)?cos(20?t)

x5(t)?sin(2?ft??/8)6（４）分别以变换区间N?8，16，32，对x(n)进行

1FFT，画出相应的幅频特性曲线。

（５）分别以变换区间N?8，16，对x(n)，x(n)进

23行FFT，画出相应的幅频特性曲线。

（６）分别以变换区间N?4，8，16，对x(n)进行

4FFT，画出相应的幅频特性曲线。

（７）分别对模拟信号x(t)和x(t)选择采样频率和采

56 13

样点数。

对x(t)?sin(2?ft??/8)，周期T?1/f，频率f自己选

5择，采样频率f样点数用T对

fs?64Hzps?4f，观测时间Tp?0.5T，T，2T，采

?fs计算。

，选择采样频率

x6(t)?cos(8?t)?cos(16?t)?cos(20?t)，采样点数为16，32，64。

56（８）分别将模拟信号x(t)和x(t)转换成序列，用

x5(n)，x(n)表示，再分别对它们进行FFT，并画出

6相应的幅频特性曲线。

4、实验用MATLAB函数介绍

fft(); figure(); plot(); stem();

abs();title(); xlabel(); ylabel(); text(); hold on; axis(); grid on; subplot(); sin(); cos(); 等。

５、思考题

x(n)和x(n)的幅频特性会相同吗? （１）在N=8时，

23为什么? N=16呢?

（２）如果周期信号的周期预先不知道，如何用FFT进行谱分析?

６、实验报告要求

（１）简述实验目的及实验原理。

（２）编程实现各实验内容，列出实验清单及说

14

明。

（３）将实验结果和理论分析结果进行比较，分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。并总结实验所得的主要结论。

（４）简要回答思考题。

15

思考题

在N=8时，x2(n)和x3(n)的幅频特性会相同吗？为什么？N=16呢？

答：Ｎ＝８时一样，Ｎ＝１６时不一样。因为DFT变换可以看成是将该序列进行周期延拓后的傅里叶级数变换的主值序列。当N=8时，两序列进行周期延拓后序列相同，所以其傅里叶级数变换的主值序列也相同，进而DFT变换

16