DSP实验报告DTMF信号产生与检测.docx

DTMF 信号的产生与检测 一、试验目的 1、 了解语音拨号的产生和原理； 2、 理解和把握 DTMF 信号的产生和检测； 3、学会使用 MATLAB 软件进展仿真分析； 二、试验原理 实际的拨号系统中，每一个按键都对应着一个高频信号和一个低频信号。按下一个键之后，就会产生有一个相应的双音频信号。通过一个窗函数就可以将其检测出来。但在实际的过程中要考虑拨号时间之间的间隔。窗的长短要适宜，不能让窗同时滑过两个语音信号，否则就检测不出来。 各按键所对应的语音信号： 1209Hz 1336Hz 1477Hz 1633Hz 697Hz 1 2 3 A 770Hz 4 5 6 B 852Hz 7 8 9 C 941Hz * 0 # D 任意按键产生的信号为： Y[n]=cos(wl*t)+cos(wh*t) 承受 Goertzel 算法进展谱分析： Goertzel 算法本质上是一种非标准 DFT，它是 DTMF 检测器的根底。利用这种方法从输入信号中提取频谱信息有效且快速。该算法实际上是利用双极点 IIR 滤波器来有效地计算 DFT 值。一般的 DFT 在开头处理之前都需要有肯定的数据量，而 Goertzel 算法的逐点递归构造就有明显的优点，可以有效地提高系统的响应速度。Goertzel 滤波器在经过肯定数量的样点 N(等效于 DFT 块的大小)后，输出便聚合到一个 Go-ertzel 算法的内部变量 Vii(n)上，并由它得出幅度的平方。简洁的数学描述如下： 递归计算：vdn)=2cos(2kn／N)·vdn．1)．vk(n．2)+x(n)其中 vd．1)=0，v。(-2)=0，x(n)=input，n 0·N。 每 N 次循环计算一次：Ix(k)12=yk(N)yk+(N)=V，(N)+Vk2(N．1)．2cos(2kn ／N)W(N)Vd(N．1) 由上面的分析可以看出，Goertzel 算法只需要对行频和列频及其相应的二次谐波，共计 8 个频点进展能量计算。由于只需要少数的频率值并且只对这些值进展滤波，所以 Goertzel 算法比 DFT 快得多。其中的二次谐波能量能够把 DTMF 信号与声音信号区分开来。Goertezl 算法的系统函数和流图如以下图所示。 1 ? cosw* 1z1 1 ? cosw* 1z 1 ? z cos w * 1z ? z ?2 计算 y k? 计算 y k 在 k={8 个值}： k 的准确值 k 取整 误差 697 17.861 18 0.139 770 19.731 20 0.269 852 21.833 22 0.167 741 24.113 24 0.113 1209 30.981 31 0.019 1336 34.235 35 0.235 1477 37.848 38 0.152 1633 41.816 42 0.154 三、 试验内容及结果 1、 产生一个双音频信号并参加噪声tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; N=205; K=[18,20,22,24,31,34,38,42]; f1=[697,770,852,941]; f2=[1209,1336,1477,1633]; TN=input(”键入8位 号码=”) TNr=0; A=zeros(1,400); x=zeros(8,400); X=0; Y=1*rand(1,400*16+1); for l=1:8; d=fix(TN/10^(8-l)); TN=TN-d*10^(8-l); for p=1:4; for q=1:4; if tm(p,q)==abs(d);break,end end end if tm(p,q)==abs(d);break,end end n=1:400; x(l,:)=2*sin(2*pi*n*f1(p)/8000)+2*sin(2*pi*n*f2(q)/8000); X=[X,x(l,:),A]; X=X+Y; for i=0:7; X1=goertzel(X(1+800*i:205+800*i),K+1); val=abs(X1); for s=5:8; if val(s)80,break,end end for r=1:4; if val(r)80,break,end end end TNr=TNr+tm(r,s-4)*10^(7-i); disp(”承受端检测到得号码为：”) disp(TNr) 运行结果： 键入 8 位 号码TN =接收端检测到的号码为：产生的信号如下： 2、 对给定的信号进展检测 加滑动平均窗，对信号进展采样： L = length(z1); f