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）
