【通达】数字信号处理实验二

实验三 快速变换及其应用

一、实验目的

(1)在理论学习的基础上，通过本实验，加深对 FFT 的理解，熟悉 MATLAB 中的有关函

数。

(2)应用 FFT 对典型信号进行频谱分析。

(3)了解应用 FFT 进行信号频谱分析过程中可能出现的问题以便在实际中正确应用 FFT。

(4)应用 FFT 实现序列的线性卷积和相关。

二、实验内容

一、

实验题目：

观察高斯序列的时域和幅频特性，固定信号 x(n)中参数 p=8，改变 q 的值,使 q 分别等

于 2、4、8,观察它们的时域和幅频特性，了解当 q 取不同值时，对信号序列的时域和幅频

特性的影响;固定 q=8,改变 p,使 p 分别等于 8、13、14,观察参数 p 变化对信号序列的时域及

幅频特性的影响,注意 p 等于多少时，会发生明显的泄漏现象，混叠是否也随之出现?记录

实验中观 察到的现象,绘出相应的时域序列和幅频特性曲线。

实验代码：

n=0:1:15;

p=8;q=2;

xa1=exp(-((n-p).^2)/q);

subplot(5,2,1);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('xa(n)');

title('p=8 q=2')

xk1=abs(fft(xa1));

subplot(5,2,2);stem(n,xk1)

xlabel('k');

ylabel('Xa(k)');

title('p=8 q=2')

p=8;q=4;

xa1=exp(-((n-p).^2)/q);

subplot(5,2,3);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('xa(n)');

title('p=8 q=4')

xk1=abs(fft(xa1));

subplot(5,2,4);

stem(n,xk1) xlabel('k');

ylabel('Xa(k)');

title('p=8 q=4')

p=8;q=8;

xa1=exp(-((n-p).^2)/q);

subplot(5,2,5);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('xa(n)');

xk1=abs(fft(xa1));

title('p=8 q=8')

subplot(5,2,6);

stem(n,xk1)

xlabel('k');

ylabel('Xa(k)');

title('p=8 q=8')

p=8;q=8;

xa1=exp(-((n-p).^2)/q);

subplot(5,2,5);

plot(n,xa1,'-*');

xlabel('t/T');ylabel('xa(n)');

xk1=abs(fft(xa1));

title('p=8 q=8') </