数字信号《用FFT对信号作频谱分析》实验完整模版.docx

《数字信号《用FFT对信号作频谱分析》实验完整模版.docx》由会员分享，可在线阅读，更多相关《数字信号《用FFT对信号作频谱分析》实验完整模版.docx（11页珍藏版）》请在第一文库网上搜索。

1、计算机科学与工程学院数字信号处理实验报告3专业班级实验时间2012 年 5 月 27 0学生学号实验地点学生姓名指导教师实验工程数字信号处理用FFT对信号作频谱分析实验类别实验学时实验目的及要求D掌握线性卷积的计算机编程方法，利用卷积的方法观察、分析系统 响应的时域特性。2）验证卷积定理。3）掌握循环卷积的计算机编程方法，并比拟与线性卷积的差异。利用 循环卷积的方法观察、分析系统响应的时域特性。成绩评定表类 另IJ评分标准分值得分合计上机表现按时出勤、遵守纪律 认真完成各项实验内容30分报告质量程序代码标准、功能正确 填写内容完整、表达收获70分说明：评阅教师：日期：2010年月日实验内容一.

2、实验目的学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析， 误差及其原因，以便正确应用FFT。二.实验内容及步骤(1)对以下序列进行谱分析。M ()= &()n +1, 0 3x2(n) = 8-，4 H 70 ,其它n4- /?, 0 3/?-3, 4 ylabeC 幅度)；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);%

3、计算 6nT 的 32 点 DFTX6k32=fftshift(X6k32);%将零频率移到频谱中心Tp=N*TjF=lp; %频率分辨率Fk=-N/2:N/2-l;fk=k*F;%产生16点DFT对应的采样点频率以零频率为中心subplot(3,l,2)jstem(fk,abs(X6k32).,)jbox on %绘制 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 = 16x6n

4、T=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); % 对6(t)64点采样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); box On%绘制 8 点 DFT 的幅频特性图title(,(6a) 64 点IDFTx_6(nT)|);XlabelCf(HZ)ylab

5、el(幅度)；axis(-N*F2-lzN*F2-l,0,1.2*max(abs(X6k64)、思考题(1)对于周期序列，如果周期不知道，如何用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均是其周期的整数倍，得到正