资源描述：
《用DFT对信号作频谱分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验三用DFT对信号作频谱分析一、实验原理计算机上实现信号的频谱分析及其他方而的处理对信号的要求是：在时域和频域都应该是离散的，而且都应该是有限长的。各种形式的傅里叶级数与变换,只有离散傅里叶级数DFS在时域和频域都是离散的，但是丘⑺)和文伙)都是无限长的周期序列，因此时域频域各取一个周期，即为离散傅里叶变换DFT,是信号离散时间傅里叶变换DTFT某种程度上的近似。频域采样即对离散时间傅里叶变换的连续周期频谱离散化的过程，采样后的周期频谱序列对应时域的周期序列，该吋域序列的周期恰好是频域屮一个周期内的采样点数采样，因此频域采样不失真的条件为：频域采样点数N要大
2、于或等于时域序列氏度Mo二、实验目的(1)学习离散叶变换(即DFT)的计算方法及意义。(2)掌握实数序列的DFT系数的对称特点。(3)利用MATLAB编制DFT/IDFT计算程序的方法。(4)频域采样理论的验证三、实验内容(1)x(n)=R5(n)f求"分别取8,16,32,64时的离散傅里叶变换DFTX伙),最后绘岀图形。程序代码：%产生序列[H]®xl(n)=R5(n)%计算xln的8点DFl
3、%计算xln的16点DFI%计算Kin的32点DFI%计算xln的64点DFIxln=.[ones(1$5)]:Xlk8=fft(xln,8):Xlkl6=fft(
4、xln^16):Xlk32=fft(xln,32):Xlk64=fft(xln,64):%以下绘制幅频特性曲线subplot(4,1,1)：stejTi(Xlk8):%绘制8点DFT的幅频特性图titleC'la)8点DFT[x_1'n)]?);xlabel(?oj/JT-')；ylabelC幅度');axis([0)2,0,1.2*max(abs(Xlk8))])subplot(4,1,2)：stem(Xlkl6);%绘制16点DFT的幅频特性图titleC'lb)16点DFT[x_1'n)]?);xlabel(?oj/JT-')；ylabelC幅度');a
5、xis([0,2^0^1.2*max(abs(Xlkl6))]):subplot(4,1,3)：stem(Xlk32);%绘制16点DFT的幅频特性图titleC'lb)32点DFT[x_l(n)]");xlabelCoj/jt-')；ylabelC幅度');axis([0j2^0,1.2*max(abs(Xlk32))])：subplot(4,1,4);stem(Xlk64);%绘制16点DFT的幅频特性图titleC(lb)64点DFT[x_l(n)]')：xlabelC⑴/心):ylabelCIpw度):axis([0,2,0,1.2*max(abs(X
6、lk64))])：8点DFT[x.[(n)]logJ45w/tt16点DFTIx^n)]10log①仃)广[口'ff810121416W/lT32点DFT[X[(n)]18logIIII15202530w/tt64点DFT[X[(n)]log5~oU0102030m/tt⑵利用如下MATLAB程序生成三角波序列%x=[lz1,1,1,1,1,lz1];M=27;N=32;n=0:M;労产生M长三角波序列x(n)xa=0:floor(M/2);xb=ceil(M/2)-1:-1:0;x=[xa,xb];对该序列分别计算离散吋间傅立叶变换DTFT,8点，16点，
7、32点，64点和128点离散傅立叶变换频谱，并利用反变换求各个频谱对应的是与序列，比较这些频谱和序列。生成的三角波图形：三角波序列图1-1长度为27的三角波其程序代码：M=27：N二32;n=O:M:%产生长度为M的三角波序列x2(n)xa=O:floor(M/2):xb=ceil(M/2)-l:-l:0;x2n=[xa,xb]:%产生长度为M的三角波序列(n)figure(5):stem(n,x2n,・’)：titleC三角波序列');对该序列分别计算离散时间傅立叶变换DTFT,8点,16点,32点,64点和128点离散傅立叶变换频谱。其实验结果为图1-2所
8、示。图1-2三角波计算离散吋间福利叶变换其程序代码:M=27：N=32：n=O:M;%产生长度为M的三角波序列x2(n)xa=0:floor(M/2):xb=ceil(M/2).：-l;x2n=[xa,xb]:%产生长度为M的三角波序列疋〔rOX2k8=fft(x2n,8):%计Mx2n的8点DFTX2kl6=fft(x2n?16):%计算:<2n的16点DFIX2k32=fft(x2n,32):%计算x2n的32点DFTX2k64=fft(x2n,64);%计算x2n的64点DFTX2kl28=fft(x2n,128):%计算戈2n的128点DFTfigur
9、e(2)subplot(5,1,1)：
