实验四 离散信号的频域分析
一、 实验目的
1.
掌握序列的傅里叶变换、离散傅里叶级数、离散傅里叶变换、快速傅里叶变换的Matlab 实现; 2.
学习用FFT 对连续信号和离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确应用FFT 。
二、 实验内容及步骤
1. 计算序列的DTFT 和DFT ,观察栅栏效应
设)()(4n R n x =,要求用MATLAB 实现:
(1)计算)(n x 的傅里叶变换)(ωj e X ,并绘出其幅度谱;
(2)分别计算)(n x 的4点DFT 和8点DFT ,绘出其幅度谱。并说明它们和)(ωj e X 的关系。
(提示:DFT 变换可用MA TLAB 提供的函数fft 实现,也可以自己用C 语言或matlab 编写)
2.计算序列的FFT ,观察频谱泄漏
已知周期为16的信号)16
12cos()1610cos()(n n n x ππ+=。 (1) 截取一个周期长度M=16点,计算其16点FFT ,并绘出其幅度谱;
(2) 截取序列长度M=10点,计算其16点FFT ,绘出其幅度谱,并与(1)的结果进行比
较,观察频谱泄漏现象,说明产生频谱泄漏的原因。
三、 实验报告要求
1.
结合实验中所得给定典型序列幅频特性曲线,与理论结果比较,并分析说明误差产生的原因以及用FFT 作谱分析时有关参数的选择方法。 2. 总结实验所得主要结论。
1. 计算序列的DTFT 和DFT ,观察栅栏效应
设)()(4n R n x =,要求用MATLAB 实现:
(1)计算)(n x 的傅里叶变换)(ωj e X ,并绘出其幅度谱;
(2)分别计算)(n x 的4点DFT 和8点DFT ,绘出其幅度谱。并说明它们和)(ωj e
X 的关
系。
(1)代码:
n=0:3;
M=10;