[工学]09271119_《数字信号处理》课程研究性学习报告_DSP1.doc

《数字信号处理》课程研究性学习报告 姓名 张卓 学号 同组成员 王一明 王欣然 郭鹏飞 马高飞 指导教师 胡健 时间 2011.04.01 DFT近似计算信号频谱专题研讨 【目的】 (1) 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 (2) 理解误差产生的原因及减小误差的方法。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 1．利用DFT分析x(t)=Acos(2?f1t)+Bcos(2?f2t)的频谱，其中f1=100Hz，f2=120Hz。 (1)A=B=1; (2)A=1,B=0.2。要求选择不同的窗函数。 【题目分析】 我们以f=600Hz(fmax(f1,f2))对该信号进行抽样，为了使抽样过程中造成频谱混叠。 【仿真结果】【结果分析】 (1) 通过读书，我们可知由于x[k]为无限长序列，可采用矩形窗对其进行加窗截断。为了能够分辨这两个间隔为△f=f2-f1的相邻谱峰，由N≥f／△f=30,所以我们以N=30，由DFT计算出的频谱如下(第一个图为矩形窗，第二个图为Hamming窗) (2) 由于在信号x（t）中存在一个较弱的频率分量f2，若利用矩形窗函数进行加窗，则由于旁瓣泄露较大，很难检测出x（t）幅度较小的频率分量f2，因而采用Hamming窗函数。当截取信号N=50个样点时，由于所取数据点较少，频率分辨率较低，由DFT计算出的频谱仍不能显示幅度较小的频率分量。当截取信号N=100个样点时，能够清楚地显示出幅度较小的频率分量。 (第一个图为矩形窗，第二个图为Hamming窗) 利用DFT分析连续非周期信号的频谱 混叠现象 泄漏现象 【自主学习内容】 利用DFT分析连续非周期信号的频谱 混叠现象 泄漏现象 【阅读文献】 《数字信号处理》（第三版） 陈后金主编 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 使用矩形窗函数和Hamming窗函数的差别。 【问题探究】 在选择矩形窗函数窗口的长度时，增加窗口的长度N可以减少主瓣宽度，改善频率分辨率，但增加的N不能减少频率泄露。为了减少旁瓣引起的频率泄露，常用的方法是采用非矩形窗。非矩形窗幅度都是逐渐减小，因而可减小矩形窗突然截断而产生的较高的旁瓣分量。此题中利用Hamming窗便可以避免矩形窗出现的问题。 当A，B都很大时，用矩形窗；当A，B都很小或者其中一个很小时，用非矩形窗。 【仿真程序】 (1)用矩形窗 f1=100; f2=120; N=30; L=512; fs=600; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; A=1; B=1; x=A*cos(2*pi.*f1.*t)+B*cos(2*pi.*f2.*t); X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(幅度谱); （1）用Hamming窗 f1=100; f2=120; N=30; L=512; fs=600; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; A=1; B=1; x=A*cos(2*pi.*f1.*t)+B*cos(2*pi.*f2.*t); wh=(hamming(N)); x=x.*wh; X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(幅度谱);? (2)用矩形窗 f1=100; f2=120; N=100; L=512; fs=600; T=1/fs; ws=2*pi*fs; t=(0:N-1)*T; A=1; B=0.2; x=A*c