基于DFT的信号识别系统实验报告.docx

实验一基于DFT的信号识别系统 专业班级：12电子信息工程 团队成员：1228401141 顾鹏伟 1228401133 陆遥 1209401045 张春辉 一、【实验目的】 （1）通过实验巩固离散福利叶变换DFT的认识和理解。 （2）熟练掌握应用DFT进行频谱分析的方法。 （3）了解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）输入信号1.2sin(0.08πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 （3）输入信号1.5+3cos(0.5πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 （4）输入信号0.7sin(0.14πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 （5）输入信号1.2cos(0.5πn)+ 9.5sin(0.02πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 （6）输入信号cos(0.102πn)，理论计算并画出0≤f≤fs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 三、【程序代码】 MATLAB程序： function [RETURN ] = Dft( A,a,B,b,C ) fs=10000; N=1000; n=0:N-1; Amp=0; max=0; x=A*cos(a*pi*n)+B*sin(b*pi*n)+C; y=x; 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 Amp=abs(y(k)) Amp=abs(y(k)); max=k-1; end end fmax=fs*max/N; belong=type; if fmax=1000 belong=第四类; elseif (fmax=500)(fmax1000) belong=第三类; elseif (fmax=200)(fmax500) belong=第二类; elseif (fmax=0)(fmax200) belong=第一类; end RETURN={max,fmax,Amp,belong}; subplot(2,1,1); n=0:N-1; plot(n,x) xlabel(n); ylabel(x(n)) subplot(2,1,2); k=0:N-1; plot(k,abs(y)) xlabel(k); ylabel(|X(k)|); end 四、【实验分析】 信号的离散傅立叶变换定义如下： ， k=0~N-1 ，所以。 1、输入信号为x(n)=1.2sin(0.08πn)时 理论计算：峰值出现在=400Hz 实验数据： [RETURN ] = Dft(0,0,1.2,0.08,0 ) RETURN = [40] [400] [600.0000] 第二类 返回值表示 K=40 频率f=400 Hz 幅值600 最大峰值频率400Hz 属于第2类 2、输入信号为x(n)=1.5+3cos(0.5πn)时 理论计算： 输入信号有两个频率分量，峰值出现在=0Hz??和f=2500Hz 实验数据： [RETURN ] = Dft(3,0.5,0,0,1.5 ) RETU