北交大DSP性学习报告频谱计算2014试题.docVIP

《数字信号处理》课程研究性学习报告姓名学号同组成员指导教师胡健时间2014年4月8日星期二DFT近似计算信号频谱专题研讨 【目的】 1 掌握利用DFT近似计算不同类型信号频谱的原理和方法。 2 理解误差产生的原因及减小误差的方法。 3 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 1．D大调和弦由频率为293.66, 369.99, 440Hz的正弦信号组成，即xt cos2?f1t+ cos2?f2t+ cos2?f3t的频谱，其中f1293.66Hz，f2369.99Hz，f3440Hz。 （1）确定合适的DFT参数。 （2）利用DFT分析其频谱，并比较不同窗函数对谱分析结果的影响。 （3）若乐曲全音符的持续时间为0.2s，该和弦为16分音符，利用DFT分析其频谱会出现什么问题？ 【题目分析】 首先确定合适的DFT参数：抽样频率 fs；数据长度：矩形窗N，hamming窗N；L，利用fft函数分析频谱。 【仿真结果】 【结果分析】 对实验结果进行比较，总结出选择合适DFT参数的原则。 （1）抽样频率 fs900Hz；数据长度：矩形窗N30，hamming窗N60；L512， （2）利用矩形窗分析信号dft时，增加窗口的长度N可以减少WNejΩ的主瓣宽度，改善频率分辨率。但增加N并不能减少频率泄漏。利用Hamming窗分析信号dft时，可减小旁瓣引起的频率泄漏。 （3）乐曲全音符的持续时间为0.2s，该和弦为16分音符，则该和旋持续时间为0.2/160.0125s,如果仍以抽样频率 fs900Hz抽样，N最大只能取到11，引起频率分辨率下降。 【自主学习内容】 相关乐理知识 【阅读文献】 《数字信号处理（第二版）》 【发现问题】 专题研讨或相关知识点学习中发现的问题： 时域截短时对频域造成的影响。 【问题探究】 时域截短会造成的频域分辨率下降。 【仿真程序】 N100;L512; f1293.66;f2368.99;f3440; fs1000;ws2*pi*fs; T1/fs; t0:N-1*T; xcos2*pi*f1*t+cos2*pi*f2*t+cos2*pi*f3*t; subplot1,2,1; Xfftshiftfftx,L; w-ws/2+0:L-1*ws/L/2*pi; plotw,absX; ylabel矩形窗截短后幅度谱; subplot1,2,2; whhammingN;xx.*wh; Xfftshiftfftx,L; w-ws/2+0:L-1*ws/L/2*pi; plotw,absX; ylabel汉明窗截短后幅度谱;【研讨题目】 基本题 2．已知某离散序列为 x[k]cos?0k+0.75cos?1k, 0 k 63 其中?00.4?, ?1?0+?/64 1 对x[k]做64点FFT, 画出此时信号的频谱。 2 如果1中显示的谱不能分辨两个谱峰，是否可对1中的64点信号补零而分辨出两个谱峰。通过编程进行证实，M2-4） 【温磬提示】 在计算离散非周期序列频谱时常用?/?作为横坐标，称?/?为归一化频率? 1 2补零后 3通过延长波形增加分辨率： 【结果分析】 【自主学习内容】 【阅读文献】 【发现问题】 专题研讨或相关知识点学习中发现的问题： 【问题探究】 【仿真程序】 1不补零clear all; N64; W00.4*pi; W10.4*pi+pi/64; k0:N-1; xcosW0*k+0.75*cosW1*k; LN; Xfftshiftfftx,L; w-pi+0:L-1*2*pi/L; plotw,absX; axis[-pi;pi;0;40]; ylabel幅度谱; 2补零后 clear all; N64; W00.4*pi; W10.4*pi+pi/64; k0:N-1; xcosW0*k+0.75*cosW1*k; L4*N;%补3N个零 Xfftshiftfftx,L; w-pi+0:L-1*2*pi/L; plotw,absX; axis[-pi;pi;0;40]; ylabel幅度谱; 2延长序列： clear all; N128; W00.4*pi; W10.4*pi+pi/64; k0:N-1; xcosW0*k+0.75*cosW1*k; LN; Xfftshiftfftx,L; w-pi+0:L-1*2*pi/L; plotw,absX; axis[-pi;pi;0;50]; ylabel幅度谱;【研讨题目】 基本题 3．已知一离散序列为 x [k]Acos?0k+Bcos0+??k。用长度N64的哈明窗对信号截短后近似计算其频谱。试用不同的