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

1、实验01 基于dft的信号识别系统 10级轨道通信 邵锰1042401029许长浩1042401044 蒋屹1042401037 一、实验目的（1）通过实验巩固离散傅利叶变换dft的认识和理解。（2）熟练掌握应用dft进行频谱分析的方法。（3）了解dft离散频谱分析的应用价值。二、实验内容在语音识别、雷达信号处理、生物医学信号检测与识别等领域广泛使用基于离散傅里叶变换的谱分析技术。一个典型的信号识别系统如图1所示：短时窗dft峰值检测分类器x(n)y(n)图1 信号识别系统框图设系统的输入信号x(n)是具有单一频谱峰值的正弦信号，短时矩形窗将信号截短尾有限长，经过dft变换得到频谱，频率检测器

2、检测频谱最大峰值的位置，即对应的频率，然后由分类器识别信号的类别。分类器的分类判决规则为：第一类：最大峰值频率分布范围（hz）为0f200。第二类：最大峰值频率分布范围（hz）为200f500。第三类：最大峰值频率分布范围（hz）为500f1000。第四类：最大峰值频率分布范围（hz）为f1000。设采样频率fs=10000hz，短时矩形窗宽度为n=1000，短时加窗信号经过dft可以得到连续频谱在01.2，所以应该出现在100hz处。(5)输入信号为x(n)=cos(0.102n)时，峰值出现在f=fs/2=510hz。四、实验结果（1）编程实现该系统function fs=dft1(a,a

3、,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; %定义s用于记录最大峰值fs=0,0,0; %将返回值定义为数组用于返回多个数for k=1:n %dft y(k)=0; n=1; while(nn+1) y(k)=y(k)+x(n)*exp(-j*2*pi*(k-1)*(n-1)/n); n=n+1; end if s=0 & fm=200 & fm=500 & fm=1000 title(第四类);end（2）输入信号

4、1.2sin(0.08n)，理论计算并画出0ffs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 dft1(0,0,1.2,0.08,0)ans = 40 400 600（k, f, a）图2理论计算：由=2f/fs得f=0.08*10000/2=400hz 第二类，与matlab显示结果相符。（3）输入信号1.5+3cos(0.5n)，理论计算并画出0ffs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 dft1(3,0.5,0,0,1.5)ans = 0 0 1500（k, f, a）图3理论计算：因为存在的直流成分较大，因此为第一类，与matlab显示结果相符

5、。（4）输入信号0.7sin(0.14n)，理论计算并画出0ffs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 dft1(0,0,0.7,0.14,0)ans = 70 700 350（k, f, a）图4理论计算：由=2f/fs得f=0.14*10000/2=700hz 第三类，与matlab显示结果相符。（5）输入信号1.2cos(0.5n)+ 9.5sin(0.02n)，理论计算并画出0ffs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 dft1(1.2,0.5,9.5,0.02,0)ans = 990 9900 4750（k, f, a）图5理论计算：由

6、=2f/fs得f=0.02*10000/2=100hz 第四类，与matlab显示结果相符。（6）输入信号cos(0.102n)，理论计算并画出0ffs范围的幅度谱，标出峰值频率，观察实际识别结果，分析其正确性。 dft1(1,0.102,0,0,0)ans = 51 510 500（k, f, a）图6理论计算：由=2f/fs得f=0.102*10000/2=510hz 第三类，与matlab显示结果相符。五、思考题1.当矩形窗长度比1000小，例如32，以上实验内容（6）可能出现什么情况？答：时域窗函数不够长，采样点不够多，从而导致频域的频率分辨率不够高，现频谱泄露出现象。实验结果如图7所

7、示：图72. 当输入信号x(n)=cos(0.199n)+0.9sin(0.204n)时，系统能够得到正确的识别结果吗？为什么？答：不能。因为0.199和0.204这两个频率过于接近，时域窗长度仅为n=1000，所以频率分辨率偏小，无法区分开来，因此是错误的识别，为第四类，如下图8。若将时域窗函数加长为n=10000，频域分辨率提高，满足要求，则能正确的识别，为第三类，如下图9.图8 图93.如果输入信号x(n)含有叠加性宽带噪声e(n)会影响识别结果吗？为什么？答：不一定会，噪声信号是宽带信号，在频域内随机分布，相对比较平坦，若幅值不是太大，而输入信号x(n)则是幅值较大的窄带频率单一的信号

8、，因此叠加性宽带噪声不会干扰系统的识别，如图11所示。若噪声信号过大，大到几乎与输入信号幅值接近甚至大于输入信号，则会干扰系统的识别，甚至根本分辨不清输入信号，如图12、13所示。经过粗略计算分析，信噪比至少要大于4db，输入信号才能清晰被识别；而当信噪比小于2db,输入信号将基本不能被识别；当信噪比小于0.3db时，输入信号将完全不能被识别。图10图11图12图13 4.如果系统中的dft要更新为fft，且短时窗不变，则fft计算应做哪些考虑，对识别结果有什么影响？答：fft中要求信号长度是2的整数幂，现在n=1000点，则须将其改为2的整数幂，如1024，则需要进行（1024/2）log21024=5120次复数乘法运算。对识别结果的频谱分析和时间分辨率不会造成影响，可以提高频率分辨率。六、实验总结 本次实验主要是通过matlab编程来实现基于dft的信号识别系统，通过实验来加深我们对dft、频率分辨率、短时窗、时间分辨率等概念的理解，以及对频率、时间分辨率之间的矛盾、信噪比的重要性的认识。由于在做实验前做了较好的准备工作，也使得本次实验的完成相对比较容易。 在这次实验中，我们巩固了离散傅利叶变换dft的认识和理解，熟练掌握应用dft进行频谱分析的方法，了解dft离散频谱分析的应用价值。也让我们了解了matlab这一高级软件，学习了它的编程环