南京邮电大学DSP实验报告

10.50-0.5幅度-1-1.5-2-2.50510152025n3035404550

；

实验二 快速傅里叶变换(FFT)及其应用

一、实验目的

(1) 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。 (2) 应用FFT对典型信号进行频谱分析。

(3) 了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。 (4) 应用FFT实现序列的线性卷积和相关。

二、实验内容

实验中用到的信号序列 a) 高斯序列

??(n?p)q?xa(n)??e??020?n?15

其他b) 衰减正弦序列

?e?ansin(2?fn)0?n?15 xb(n)??0其他?c) 三角波序列

0?n?3?n?xc(n)??8?n4?n?7

?0其他?d) 反三角波序列

?4?n0?n?3?xd(n)??n?44?n?7

?0其他?

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

等于2，4，8，观察它们的时域和幅频特性，了解当q取不同值时，对信号序列的时域幅频特性的影响；固定q=8，改变p，使p分别等于8，13，14，观察参数p变化对信号序列的时域及幅频特性的影响，观察p等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。

解：

function gauss(p,q) n=0:1:15; N=length(n); xa=exp(-(n-p).^2/q); M=10000;

w=2*pi/M*(0:1:M-1); Xa=zeros(1,M); for k=1:M

Xa(k)=sum(xa*(exp(-j*w(k)*(0:N-1)'))); end

subplot(2,1,1); stem(n,xa);

xlabel('n'),ylabel('x_a(n)') subplot(2,1,2);

plot(w,abs(Xa))

xlabel('\\omega'),ylabel('幅度谱')

>>gauss〔8，2〕

>>gauss（8，4）

>>gauss（8，8）

>>（13，8）

>>（14，8）