基于DFT的信号识别系统.doc

基于DFT的信号识别系统 姓名__孔亮__ 学号__0928401116__ 专业__电子信息工程__ 一、实验目的 通过实验巩固离散傅立叶变换DFT的认识和理解。 熟练掌握应用DFT进行频谱分析的方法。 理解DFT离散频谱分析的应用价值。 二、实验内容 在语音识别、雷达信号处理、生物医学信号检测与识别等应用领域广泛使用基于离散傅立叶变换的谱分析技术。一个典型的信号识别系统如下图： 短时窗 短时窗 DFT 峰值检测 分类器 x(n) y(n) 设系统的输入信号x(n)是具有单一频谱峰值的正弦信号，短时矩形窗将信号截短为有限长，经过DFT变换得到频谱，频率检测器检测频谱最大峰值的位置，即对应的频率，然后由分类器识别信号的类别。分类器的分类判决规则为： 第一类：最大峰值频率分布范围（Hz）为0≤f≤200。 第二类：最大峰值频率分布范围（Hz）为200≤f≤500。 第三类：最大峰值频率分布范围（Hz）为500≤f≤1000。 第四类：最大峰值频率分布范围（Hz）为f≥1000。 设采样频率fs=10000Hz，短时矩形窗宽度为N=1000，短时加窗信号经过DFT可以得到连续频谱在0≤ω＜2范围内的1000个取样点。 （1）编程实现该系统 （2）输入信号x（n）=1.2sin(0.08πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统的实际识别结果，分析其正确性。 （3）输入信号x（n）=1.5+3cos(0.5πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统的实际识别结果，分析其正确性。 （4）输入信号x（n）=0.7sin(0.14πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统的实际识别结果，分析其正确性。 （5）输入信号x（n）=1.2cos(0.5πn)+ 9.5sin(0.02πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统的实际识别结果，分析其正确性。 （6）输入信号x（n）=cos(0.102πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察系统的实际识别结果，分析其正确性。 三、实验分析 设x(n)是长度为N的有限长信号（注意这个前提），即信号仅仅分布在[0,N-1]区间，其余时间均为0，那么，该信号的离散傅立叶变换定义如下： ， k=0~N-1 f与k的关系为： f= 因为ω=2πf/fs，所以有f=ωfs/2π。以此为依据对输入信号作如下理论分析： (1)???输入信号为x(n)=1.2sin(0.08πn)时，峰值出现在f=ωfs/2π?=400Hz。 (2)???输入信号为x(n)=1.5+3cos(0.5πn)时，峰值出现在f=ωfs/2π=0Hz???或f=2500Hz，而因为直流分量不仅使笫1条谱线（0频率）有很大的值，同时由于泄漏的存在，对低频部分造成影响，所以一般频谱分析要去掉0频率，所以最大峰值出现在f=2500Hz。 (3)???输入信号为x(n)=0.7sin(0.14πn)时，峰值出现在f=ωfs/2/π=700Hz。 (4)???输入信号为x(n)=1.2cos(0.5πn)+ 9.5sin(0.02πn)时，峰值出现在f=ωfs/2π=2500Hz?或f=100Hz，因为幅值9.51.2，所以应该出现在100Hz处。 (5)???输入信号为x(n)=cos(0.102πn)时，峰值出现在f=ωfs/2π=510Hz。 四、实验结果 （1）MATLAB 程序如下： function FS=dft1(A,a,B,b,C) fs=10000; N=1000; n=0:(N-1); x=A*cos(a*pi*n)+B*sin(b*pi*n)+C; y=x; s=0; FS=[0,0,0]; for k=1:N y(k)=0; n=1; while(nN+1) y(k)=x(n)*exp(-j*2*pi*(k-1)*(n-1)/N)+y(k); n=n+1; end if kN/2 if s=abs(y(k)) s=abs(y(k)); m=k-1; end end end f=fs*m/N; FS=[m,f,s]; k=0:N-1; stem(k,abs(y),.); end (2) X(n)=1.2sin(0.08πn) 系统识别结果： 点：k=40 频率：f=400 Hz 幅值600 最大峰值频率为400Hz，属于第2类 与理论值相同，系统识别结果正确。 X(n)= 1.5+3cos(0.5π