1、计算机科学与工程学院数字信号处理实验报告3专业班级实验时间2012年5月270学生学号实验地点学生姓名指导教师实验工程数字信号处理用FFT对信号作频谱分析实验类别实验学时实验目的及要求D掌握线性卷积的计算机编程方法，利用卷积的方法观察、分析系统响应的时域特性。2）验证卷积定理。3）掌握循环卷积的计算机编程方法，并比拟与线性卷积的差异。利用循环卷积的方法观察、分析系统响应的时域特性。成绩评定表类别评分标准分值得分合计上机表现按时出勤、遵守纪律认真完成各项实验内容30分报告质量程序代码标准、功能正确填写内容完整、表达收获70分说明：评阅教师：日期：2010年月日实验内容一.实验目的学习用FFT对连
2、续信号和时域离散信号进行谱分析的方法，了解可能出现的分析，误差及其原因，以便正确应用FFT。二.实验内容及步骤(1)对以下序列进行谱分析。M()=&()n+1,03x2(n)=8-，4H70,其它n4-/?,03/?-3,4ylabeC幅度)；axis(-N*F2-lzN*F2-lz0z1.2*max(abs(X6kl6)N=32;n=0:N-l;%FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);%对6(t)32点采样X6k32=fft(x6nT);%计算6nT的32点DFTX6k32=fftshift(X6k32);%
3、将零频率移到频谱中心Tp=N*TjF=lp;%频率分辨率Fk=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率以零频率为中心subplot(3,l,2)jstem(fk,abs(X6k32).,)jboxon%绘制8点DFT的幅频特性图title(,(6b)32点IDFT区6(nT);XIabelCf(HZ)TyIabeIC幅度)；axis(-N*F2-lzN*F2-l,0z1.2*max(abs(X6k32)N=64;n=0:N-l;%FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T);%对6(t)6
4、4点采样X6k64=fft(6nT);%计算6nT的64点DFTX6k64=fftshift(X6k64);%将零频率移到频谱中心TP=N*T;F=1/Tp;%频率分辨率Fk=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率以零频率为中心SUbPlOt(3,1,3);Stem(fk,abs(X6k64),Y);boxOn%绘制8点DFT的幅频特性图title(,(6a)64点IDFTx_6(nT)|);XlabelCf(HZ)ylabel(幅度)；axis(-N*F2-lzN*F2-l,0,1.2*max(abs(X6k64)、思考题(1)对于周期序列，如果周期不知道，如
5、何用FFT进行谱分析？周期信号的周期预先不知道时,可先截取M点进行DFT,再将截取长度扩大1倍截取,比拟结果,如果二者的差异满足分析误差要求,那么可以近似表示该信号的频谱,如果不满足误差要求就继续将截取长度加倍,重复比拟,直到结果满足要求(2)如何选择FFT的变换区间？(包括非周期信号和周期信号)一、对于非周期信号：有频谱分辨率F,而频谱分辨率直接和FFT的变换区间有关，因为FFT能够实现的频率分辨率是2/N.因此有最小的N2/Fo就可以根据此式选择FFT的变换区间。二、对于周期信号，周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT,得到的离散谱才能代表周期信号的频谱。(3)当N=8时，七5)和13()的幅频特性会相同吗？为什么？N=16呢？%()=COSl的周期为8,所以N=8和N=16均是其周期的整数倍，得到正确的单一频率正弦波的频谱，仅在0.25处有1根单一谱线。如下图。5()=。(%(刀72/4)+以%(71/8)的周期为16,所以N=8不是其周期的整数倍，得到的频谱不正确，如下图。N=16是其一个周期，得到正确的频谱，仅在0.25元和0.125n处有2根单一一谱线，如下图。实验总结通过这次实验，学习了线性卷积的计算机编程方法，利用卷积的方法观察、分析用FFT对信号作频谱分析和计算的方法。对matlab有了更多的了解。
