北交大DSP性学习频谱计算2014程序.doc

《数字信号处理》课程研究性学习报告 姓名 学号 同组成员 指导教师 胡健 时间 2014年4月8日星期二 DFT近似计算信号频谱专题研讨 【目的】 (1) 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 (2) 理解误差产生的原因及减小误差的方法。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 1．D大调和弦由频率为293.66, 369.99, 440Hz的正弦信号组成，即x(t)= cos(2?f1t)+ cos(2?f2t)+ cos(2?f3t)的频谱，其中f1=293.66Hz，f2=369.99Hz，f3=440Hz。 （1）确定合适的DFT参数。 （2）利用DFT分析其频谱，并比较不同窗函数对谱分析结果的影响。 （3）若乐曲全音符的持续时间为0.2s，该和弦为16分音符，利用DFT分析其频谱会出现什么问题？ 【题目分析】 首先确定合适的DFT参数：抽样频率 fs；数据长度：矩形窗N，hamming窗N；L，利用fft函数分析频谱。 【仿真结果】 【结果分析】 对实验结果进行比较，总结出选择合适DFT参数的原则。 （1）抽样频率 fs=900Hz；数据长度：矩形窗N=30，hamming窗N=60；L=512， （2）利用矩形窗分析信号dft时，增加窗口的长度N可以减少WN(ejΩ)的主瓣宽度，改善频率分辨率。但增加N并不能减少频率泄漏。利用Hamming窗分析信号dft时，可减小旁瓣引起的频率泄漏。 （3）乐曲全音符的持续时间为0.2s，该和弦为16分音符，则该和旋持续时间为0.2/16=0.0125s,如果仍以抽样频率 fs=900Hz抽样，N最大只能取到11，引起频率分辨率下降。 【自主学习内容】 相关乐理知识 【阅读文献】 《数字信号处理（第二版）》 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 时域截短时对频域造成的影响。 【问题探究】 时域截短会造成的频域分辨率下降。 【仿真程序】 N=100;L=512; f1=293.66;f2=368.99;f3=440; fs=1000;ws=2*pi*fs; T=1/fs; t=(0:(N-1))*T; x=cos(2*pi*f1*t)+cos(2*pi*f2*t)+cos(2*pi*f3*t); subplot(1,2,1); X=fftshift(fft(x,L)); w=(-ws/2+(0:L-1)*ws/L)/(2*pi); plot(w,abs(X)); ylabel(矩形窗截短后幅度谱); subplot(1,2,2); 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．已知某离散序列为 x[k]=cos(?0k)+0.75cos(?1k), 0( k ( 63 其中?0=0.4?, ?1=?0+?/64 (1) 对x[k]做64点FFT, 画出此时信号的频谱。 (2) 如果(1)中显示的谱不能分辨两个谱峰，是否可对(1)中的64点信号补零而分辨出两个谱峰。通过编程进行证实，M2-4） 【温磬提示】 在计算离散非周期序列频谱时常用?/?作为横坐标，称?/?为归一化频率? (1) (2)补零后 (3)通过延长波形增加分辨率： 【结果分析】 【自主学习内容】 【阅读文献】 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 【问题探究】 【仿真程序】 (1)不补零clear all; N=64; W0=0.4*pi; W1=0.4*pi+pi/64; k=0:N-1; x=cos(W0*k)+0.75*cos(W1*k); L=N; X=fftshift(fft(x,L)); w=-pi+(0:L-1)*2*pi/L; plot(w,abs(X)); axis([-pi;pi;0;40]); yl