数字信号处理研讨DFT近似计算信号频谱.docVIP

Beijing Jiaotong University 数字信号处理研讨 DFT近似计算信号频谱 学 院： 电子信息工程学院 小组成员： 指导教师: 2013.06.24 DFT近似计算信号频谱专题研讨 【目的】 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 理解误差产生的原因及减小误差的方法。 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】基本题 己知一离散序列为 x[k] = sin(0.2 兀众)，众=0，1，…，31 ⑴用L=32点DFT计算该序列的频谱，求出频谱中谱峰的频率； 对序列进行补零，然后分别用L=64、128、256、512点DFT U?算该序列的频谱，求出频谱中谱峰 的频率； 讨论所获得的结果，给出你的结论。该结论对序列的频谱计算有何指导意义？ 【题目分析】 本题讨论补零对离散序列频谱计算的影响。 【温磬提示】 在计算离散非周期序列频谱吋常用X2/7T作为横坐标，称/2/k为归一化频率(normalized frequency)。 在画频谱时需给出横坐标。每幅图下都需给出简要的文字说明。 由于离散非周期序列频谱是周期的，所以在计算时不必用ffishift函数对fft计算的结果进行重 新排列。 【序列频谱计算的基本方法】 连续信号，通过其抽样的离散信号，和离散信号的DFT变换存在如下关系： N-\ X(ejQ) n =DFT氏㈨?[々]) = X AOTC Q=VZW a-=o 通过如上关系，我们就可以通过DFT来求信号的频谱。 【仿真结果】 □ dd^ k ?、S OS 矚0 fm = 0.9375 fm = 0.9219 fm = 0.9141 fm = 0.9102 fm = 0.9082 【结果分析】 增加DFT的点数可以使频谱更容易观察，即减轻了栅栏效应带来的影响。频谱的横坐标为归一 化频率，所以原信号的峰值第一次应该出现在0.2处，随着DFT点数的增大，频谱表示也越來越精 确。 【自主学习内容】 归一化频率相关知识。 通过matlab计算DFT和matlab的绘图操作。 【阅读文献】 数字信号处理 补零对有限长序列频谱及DFT的影响。 【发现问题】(专题研讨或相关知识点学习屮发现的问题)： DFT点数的增多是否能提高频谱分辨率？ 【问题探究】 虽然DFT点数增加使图像更加细致，但是因为不论DFT点数是多少，抽样点数都是相同的，所 以每个频谱所包含的信息相同，频谱分辨率只与抽样点数相关，与点数无关，频谱分辨率相同。 所以不能通过增大DFT点数而减少信息损失。 DFT点数的增多不能提过频率分辨率。 【仿真程序】 k = 0:31; N = 16; x = sin (0.2*pi*k); for i = 1:5 N = 2*N; X = f ft (x, N); k = 0:N-1; subplot (5, 1, i); plot (k/N,abs(X)); xlabel (▼f/HZ▼}; ylabel (丨Magnitude *); fm = find( X == max(X))/N; fm end 2 已知一离散序列为x [k]=Acosfi)k+Bcos ((Z^)+A/2)Z:)o用长度AMS4的哈明窗对信号截短后近似 计算其频谱。试川不同的和Z?的值(如和打近似相等，＞4和Z?差距较大)，确定用哈明窗 能分辩的最小的谱峰间隔A/2W = 十c的值。 N 【题目分析】 本题讨论用哈明窗H?算序列频谱时的频率分辨率。 【仿真结果】 c=4 A/B=16c=4 AB=420100201002010c=4 VB=12010020100201002010015 20 25c=3 AB=12c=4 AB=122010015 20 25c=22 A/B=1220100152025201002010020100201002010V\15 20 25c=3 A/B=161525201015 20c=22 AB=1625201020101525202515 202515 2025c=2 A?=12c=2 A/B=1Gc=2 AB=415 20c=19 A/B=12201002010025 c=4 A/B=16 c=4 AB=4 20 10 0 20 10 0 20 10 c=4 VB=1 20100 20100 20100 20100 15 20 25 c=3 AB=12 c=4 AB=12 20100 15 20 25 c=22 A/B=12 20100 15 20 25 20100 20100 20100 20100 2010 V\ 15 20 25 c=3 A/B=16 15 25 20 10 15 20 c=22 AB=16 25 2010 2010 15 2