基于单片机的DTMF信号的译码算法

1、 引言

DTMF信号首先用于电话的拨号系统，在频率编码遥控系统及数据编码传输中的应用也很普遍。目前的DTMF译码器中，大多采用通用集成器件（单音译码电路和组合门电路）或专用DTMF信号译码集成电路（如MC145436等）组成译码电路。在很多情况下，DTMF译码器输出的数据仍需送入单片机进行相应的运算及处理，进而控制其它各种设备的动作。因此，如果能找到一种基于单片机的DTMF信号的译码算法，再辅之以简单的整形电路就可以，既可省去成套译码电路，又能达到简化电路降低成本的目的。本文所要介绍的，就是这种构想的初衷，结果通过计算机仿真计算数据论证认为完全能够达到设计要求。

2、DTMF信号频率组成及整形前DTMF信号的幅值密度

在DTMF信号中，16个指令键均由两个单音频率信号组合（见表1）。单音频率有两组，高频组为（1209Hz、1336Hz、1447Hz、1633Hz），低频组为（697Hz、770Hz、852Hz、941Hz），每个指令键，对应的都是一个高频组的频率和一个低频组的频率的组合。以“*”号指令键为例：其DTMF信号是由941Hz的低频组信号和1209Hz的高频组信号组成。图1为其频谱图。

表1 DTMF信号频率组成表

从图1中可以看出，“*”号指令键的DTMF信号在941Hz和1209Hz处有谱线。在理想条件下，可计算出频谱密度函数X（f）在f=941Hz和1209Hz时的模值为|X（941）|，|X（1209）|，即它们模值均不为零。也就是说，要是同时存在频谱密度函数模值为|X（941）|，|X（1209）|，并且它们模值均不为零时，其表征的键号为“*”。以此类推（见表1），可分别计算出低频组信号|X（697）|、|X（770）|、|X（852）|、|X（941）|和高频组信号|X（1209）|、|X（1336）|、|X（1447）|、|X（1633）|的模值。如果高、低频组中均各自有一个X（f）的模值不为零，则再通过f在表1查找出其表征的指令键。

在用单片机进行X（f）运算即离散傅立叶变换（DFT）时，只能对有限长的DTMF信号进行分析与处理，即对有限时间Tp=NT内的N个数据进行离散傅立叶变换（N为采样点数，T为采样时间间隔）。

根据DFT定义式：

同样以“*” 指令键信号为例，在高频率组f=1209 Hz，采样点数N=256，采样时间间隔T=55×10-6 S时，由式K=f×N×T，可得K=17，即|X（17）|为DTMF信号在频率为1209 Hz处的幅值密度其值为121.5。在低频率组f=941 Hz，采样点数N=256，采样时间间隔T=54×10-6 S时，由式K=f×N×T，可得K=13，即|X（13）|为DTMF信号在频率为941 Hz处的幅值密度其值为123.6。同理，可计算出其它15个指令键的幅值密度，见表2（表2为对16个指令键的DTMF信号采用计算机仿真计算后的幅值密度）。

表2 整形前DTMF信号的幅值密度表

从表2看出：由于时域无限长DTMF信号被截断所引起的泄漏效应，如“2”、“3”号键对应的DTMF信号虽然不含有频率为1209 Hz和941 Hz的信号成份，可是｜X（17）|、｜X（13）|不为零，理想时应为零，也就是说存在一定的幅值密度误差。但对于含有f=1209 Hz高频组信号的DTMF信号（如“1”、“4”、“7”、“*”键），其｜X（17）|值远大于不含f=1209 Hz高频组信号的DTMF信号的｜X（17）|值。同样，对于含有f=941 Hz低频组信号的DTMF信号的｜X（13）|值远大于不含f=941Hz低频组信号的DTMF信号的｜X（13）|值，这样就为实际DTMF信号译码识别提供了必要的条件。

因在实际DTMF信号译码应用时，任一键号所对应的DTMF信号的译码过程为：

（1）分别采样DTMF信号计算出谱线为697Hz、770Hz、852Hz、941Hz的幅值密度｜X（k）|。

（2）从中排序找到低频组频率fL对应幅值密度|X（k）|最大的值。

（3）同法计算，谱线为1209Hz、1336Hz、1467Hz、1633Hz的幅值密度｜X（k）|。

（4）从中排序找到高频组频率fH对应幅值密度｜X（k）|最大的值。

（5）根据fL、fH查表1，即可得到其所表征的键号。

2、整形后DTMF信号的幅值密度及其实验数据误差分析

为了把DTMF信号送入单片机进行DTMF信号译码，还必须要对DTMF信号进行整形，见图2所示。DTMF信号经比较限幅，整形为方波后。从DFT变换定义式看出：式中x（nT）的值只能为0或者1，因此计算｜X（k）｜的运算量大大降低，表3即为通过整形后DTMF信号采用计算机仿真计算出的幅值密度。  

图2 DTMF信号比较限幅示意图

表3 整形后DTMF信号的幅值密度表

对比表3和表2可以发现，整形为方波后的DTFM信号泄漏谱线的幅值密度有所增大，而有效谱线的幅值密度也相应变小。例如7号指令键产生的泄漏谱线的幅值密度其值从14.09增大到22.38，而有效幅值密度其值从127.9减小到107.26。造成此类问题的主要原因是：

（1）由于N，T只能选择整数，1/f不可能被N×T所整除，所以这必然会带来栅栏效应，此时计算所得的有效谱线的幅值密度必然 小 于 实 际 值 。同 时 ，由于采样时间Tp=NT有限长而引起的泄漏效应，也必然会导致泄漏谱线的幅值密度增大。

（2）DTMF信号经整形为方波后会产生了十分丰富的谐波干扰，这些干扰信号的频率如果接近泄漏谱线的频率，也会使泄漏谱线的幅值密度增大。

因此在译码过程中，如果有效谱线的幅值密度值变小，而泄漏谱线的幅值密度增大，当泄漏谱线的幅值密度大于有效谱线的幅值密度时，就会引起错译和漏译现象。所以在N×T值选择过小，或者与待测周期的整数倍相差过大，再加上整形为方波后谐波干扰，将有可能引起错译和漏译。

但是从表3中可见泄漏谱线的幅值密度最大值为22.38，而有效谱线的幅值密度最小值为87.92两者相差近4倍，还存在较大的冗余量。例如：对于含有f=1209 Hz高频组信号的DTMF信号（如“1”、“4”、“7”、“*”键），其｜X（17）|值仍远大于不含f=1209 Hz高频组信号的DTMF信号的｜X（17）|值。同样，对于含有f=941 Hz低频组信号的DTMF信号的｜X（13）|值也远大于不含f=941Hz低频组信号的DTMF信号的｜X（13）|值，二者仍可在排序中明显区分出来，所以对DTMF信号整形后产生的一定幅值密度误差增大，完全可以忽略不计，只要译码应用程序设计得当,合理选取N、T值，不会出现错译和漏译现象。

在DTMF信号比较准确或干扰较小的场合，甚至还可以通过减少N，T的值来提高运算速度，减少译码时间。