北京理工大学数字信号处理实验二利用DFT分析信号频谱.docxVIP

实验2 利用DFT分析信号频谱 姓名：明眸皓齿王师傅 班级：****** 学号：********* 实验时间：第十二周周三下午第二大节 一．实验目的 加深对DFT原理的理解。 应用DFT分析信号的频谱。 深刻理解利用DFT分析信号频谱的原理，分析实现过程中出现的现象及解决方法。 二．实验设备与环境 计算机、MATLAB软件环境 三．实验原理 DFT与DTFT的关系 有限长序列的离散时间傅里叶变换在频率区间的N个等间隔分布的点上的N个取样值可以由下式表示： 由上式可知，序列的N点DFT，实际上就是序列的DTFT在N个等间隔频率点上样本。 利用DFT求DTFT 方法1：由恢复出的方法如图2.1所示： 由图2.1所示流程可知： 由式2-2可以得到 其中为内插函数 方法2：然而在实际MATLAB计算中，上述插值运算不见得是最好的办法。由于DFT是DTFT的取样值，其相邻两个频率样本点的间距为，所以如果我们增加数据的长度N，使得到的DFT谱线就更加精细，其包络就越接近DTFT的结果，这样就可以利用DFT来近似计算DTFT。如果没有更多数据，可以通过补零来增加数据长度。 利用DFT分析连续时间信号的频谱 采用计算机分析连续时间信号的频谱，第一步就是把连续时间信号离散化，这里需要进行两个操作：一是采样，二是截断。 对于连续时间非周期信号，按采样间隔T进行采样，截取长度M，那么 对进行N点频域采样，得到 因此，可以将利用DFT分析连续非周期信号频谱的步骤归纳如下： 确定时域采样间隔T，得到离散序列； 确定截取长度M，得到M点离散序列，这里w(n)为窗口函数。 确定频域采样点数N，要求。 利用FFT计算离散序列的N点DFT，得到； 根据式(2-6)由计算采样点的近似值。 实验内容 1.已知，完成如下要求： （1）计算其DTFT，并画出区间的波形。 （2）计算4点DFT，并把结果显示在（1）所画图形中。 （3）对补零，计算64点DFT，并显示结果。 （4）根据实验结果，分析是否可以由DFT计算DTFT，如果可以，如何实现。 （1）代码： x=[2,-1,1,1] n=0:3 w=-pi:0.01*pi:pi X=x*exp(-j*n*w) subplot(2,1,1) plot(w,abs(X)) xlabel(\Omega/\pi) title(Magnitude) axis tight subplot(2,1,2) plot(w,angle(X)/pi) xlabel(\Omega/\pi) title(Phase) axis tight （2）代码： x=[2,-1,1,1] n=0:3 w=-pi:0.01*pi:pi X=x*exp(-j*n*w) subplot(2,1,1) plot(w,abs(X)) xlabel(\Omega/\pi) title(Magnitude) axis tight subplot(2,1,2) plot(w,angle(X)/pi) xlabel(\Omega/\pi) title(Phase) axis tight y=fft(x) subplot(2,1,1) hold on stem(0:3,abs(y),fill) title(Magnitude) subplot(2,1,2) hold on stem(0:3,angle(y)/pi,fill) ylabel(/\pi) title(Phase) （3）代码： x=[2,-1,1,1] x=[x,zeros(1,60)] y=fft(x,64) subplot(2,1,1) stem(0:63,abs(y),fill) title(Magnitude) axis tight subplot(2,1,2) stem(0:63,angle(y)/pi,fill) xlabel(\Omega/\pi) ylabel(/\pi) title(Phase) （4）由DFT可以计算出DTFT。从上面几个图中我们可以看出，DFT可以看作是DTFT的取样值，其相邻两个频率样本点的间距为，所以如果我们增加数据的长度N（即减小点与点之间的间距），得到的DFT谱线就更加精细，其包络也就越接近DTFT的结果，这样就可以利用DFT来近似计算DTFT。 2.考察序列 （1），用DFT估计的频谱；将补零加长到长度为100点序列用DFT估计的频谱。要求画出相应波形。 （2）时，用DFT估计的频谱，并画出波形。 （3）根据实验结果，分析怎样提高频谱分辨率。 （1）代码： n=0:10 x=cos(0.48*pi*n)+cos(0.52*pi*n) y=fft(x) subplot(2,1,1) stem(0:10,abs