4DFT分析信号频谱ch2数字信号处理.ppt

数字信号处理(Digital Signal Processing) 利用DFT分析信号频谱 问题的提出 四种信号频谱之间的关系 利用DFT分析连续非周期信号频谱 混叠现象、泄漏现象、栅栏现象 DFT参数选取 求x(t)=e-t u(t)的幅度谱 求x(t)=e-t u(t)的幅度谱 小结 1.掌握X[m]与X(j?)的关系 利用DFT进行确定信号谱分析作业 1.利用DFT分析x(t)=Acos(2pf1t)+Bcos(2pf2t)的频谱，其中f1=100Hz，f2=120Hz。 (1)A=B=1; (2)A=1,B=0.2 要求选择不同的DFT参数及窗函数，并对实验结果进行比较,总结出选择合适DFT参数的原则. 2.P133 M2-1(1), M2-2, M2-3,M2-4,M2-5 3.分析矩形窗、 Hann窗、 Hamming窗、 Blackman窗、 Kaiser窗的频谱，并进行比较。 4.利用DFT分析x(t)=e -|t| u(t)的频谱。 5.利用DFT分析计算x(t )=cos2pt的频谱,设w0=p / 2 思考题 (1)既然可以直接计算CTFT，为什么利用DFT分析连续信号谱? (2)在利用DFT分析连续信号频谱时，会出现哪些误差？如何克服或减弱？ (3)在利用DFT分析连续信号频谱时，如何选择窗函数？ (4)若截短信号造成泄漏而导致频谱分辨率下降，可否通过在截短序列后补零得到改善？ (5)窗函数对频谱分辨率有何影响吗？怎样提高频谱分辨率？ (6) 如何选择取样频率？ (7) 既然频谱分辨率与信号采集时间成反比，是否意味着在实际中频谱分辨率可以很容易实现？ (8) 讨论补零对计算结果的影响。 要求 分组进行，4-5名同学一组。各组选择一位组长，本周五组长提交组员名单。 考核 以小组为单位进行考核 1.各小组交纳一份实验报告。 2.小组汇报。各组随机抽出一名组员进行汇报（做PPT），汇报内容当场随机抽取。 3.成绩10分，自评和同学互评。 考核 1.如何利用DFT分析连续非周期信号频谱？在进行谱分析时，会出现哪些误差？如何克服或减弱？ 2.如何利用DFT分析连续周期信号频谱？在进行谱分析时，会出现哪些误差？如何克服或减弱？ 3.如何利用DFT分析离散非信号频谱？在进行谱分析时，会出现哪些误差？如何克服或减弱？ 4.常用窗函数的特性，在利用DFT进行谱分析时，加窗对谱分析有何影响？如何选择窗函数？ 5. 在利用DFT进行谱分析时，增加序列长度（窗函数的长度）、序列补零对计算结果有何影响？两者的作用相同吗？ N=20 例: 频谱泄漏 不能分辨 两个谱峰 加窗 抽样 DFT 如何提高频率分辨率? 减小窗函数频谱主瓣宽度 即增加窗函数长度N N=30 提高频率分辨率 利用矩形窗计算有限长余弦信号频谱 N=30; %数据的长度 L=512; %DFT的点数 f1=100; f2=120; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+cos(2*pi*f2*t); F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); hd=plot(w,abs(F)); ylabel(幅度谱) 矩形窗 N=25 矩形窗 N=50 例: 如何解决? 由于泄漏使得信号中幅度小的频率分量难以检测 选择旁瓣幅度小的窗函数 哈明窗 哈明窗N=50 改用哈明窗截短 利用Hamming窗计算有限长余弦信号频谱 N=50; %数据的长度 L=512; %DFT的点数 f1=100;f2=150; fs=600; %抽样频率 T=1/fs; %抽样间隔 ws=2*pi*fs; t=(0:N-1)*T; f=cos(2*pi*f1*t)+0.15*cos(2*pi*f2*t); wh=(hamming(N)); f=f.*wh; F=fftshift(fft(f,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(F)); ylabel(幅度谱) (3 ) 栅栏效应： 对信号补零，抽样率不变，频域抽样点数增加，可以更多地显示出频谱中的细节 DFT的两个相邻点对应的连续信号的频谱间隔为 混叠现象、泄漏现象、栅栏现象 内插公式 从信息表达的的角度，给定了N个序列值 ，就可以 唯一地确定其频谱函数 的所有信息。同样，只需要 中的N个等间隔样点值 （m=0,1,…,N-1）就可 以唯一地确定 其中， 内插公式 式中，