通信课程设计——双音多频信号检测.docVIP

通信技术方向课程设计 题目： 双音多频信号检测 物联网工程 学院 专业表示，其中，。这样8个频率形成16种不同的双频信号。具体号码以及符号对应的频率如表1所示。表中最后一列在电话中暂时未用。 表1 双频拨号的频率分配 列 行 1209Hz 1336Hz 1477Hz 633Hz 697Hz 1 2 3 A 770Hz 4 5 6 B 852Hz 7 8 9 C 942Hz * 0 # D DTMF信号在电话中有两种作用，一个是用拨号信号去控制交换机接通被叫的用户电话机，另一个作用是控制电话机的各种动作，如播放留言、语音信箱等。 2、电话中的双音多频（DTMF）信号的产生与检测 （1）双音多频信号的产生 假设时间连续的 DTMF信号用表示，式中是按照表1选择的两个频率，代表低频带中的一个频率，代表高频带中的一个频率。显然采用数字方法产生DTMF信号，方便而且体积小。下面介绍采用数字方法产生DTMF信号。规定用8KHz对DTMF信号进行采样，采样后得到时域离散信号为 形成上面序列的方法有两种，即计算法和查表法。用计算法求正弦波的序列值容易，但实际中要占用一些计算时间，影响运行速度。查表法是预先将正弦波的各序列值计算出来，寄存在存储器中，运行时只要按顺序和一定的速度取出便可。这种方法要占用一定的存储空间，但是速度快。 因为采样频率是8000Hz，因此要求每125ms输出一个样本，得到的序列再送到D/A变换器和平滑滤波器，输出便是连续时间的DTMF信号。DTMF信号通过电话线路送到交换机。 （2）基于Goertzel算法的双音多频信号检测 在接收端，要对收到的双音多频信号进行检测，检测两个正弦波的频率是多少，以判断所对应的十进制数字或者符号。显然这里仍然要用数字方法进行检测，因此要将收到的时间连续 DTMF信号经过A/D变换，变成数字信号进行检测。检测的方法有两种，一种是用一组滤波器提取所关心的频率，根据有输出信号的2个滤波器判断相应的数字或符号。另一种是用DFT（FFT）对双音多频信号进行频谱分析，由信号的幅度谱，判断信号的两个频率，最后确定相应的数字或符号。当检测的音频数目较少时，用滤波器组实现更合适。FFT是DFT的快速算法，但当DFT的变换区间较小时，FFT快速算法的效果并不明显，而且还要占用很多内存，因此不如直接用DFT合适。下面介绍Goertzel算法，这种算法的实质是直接计算DFT的一种线性滤波方法。这里略去Goertzel算法的介绍，可以直接调用MATLAB信号处理工具箱中戈泽尔算法的函数Goertzel，计算N点DFT的几个感兴趣的频点的值。 3、MATLAB工具箱函数goertzel Goerztel函数的调用格式为： Xgk=goertzel(xn,K) xn是被变换的时域序列，用于DTMF信号检测时，xn就是DTMF信号的205个采样值。 K是要求计算的DFT[xn]的频点序号向量，用N表示xn的长度，则要求1≤K≤N。由表2可知，如果只计算DTMF信号8个基频时， K=[18，20，22，24，31，34，38，42]， 如果同时计算8个基频及其二次谐波时， K=[18，20，22，24，31，34，35，38，39，42，43，47，61，67，74，82]。 Xgk是变换结果向量，其中存放的是由K指定的频率点的DFT[x(n)]的值。设X(k)= DFT[x(n)]，则。 4、检测DTMF信号的DFT参数选择 用DFT检测模拟DTMF信号所含有的两个音频频率，是一个用DFT对模拟信号进行频谱分析的问题。根据第三章用DFT对模拟信号进行谱分析的理论，确定三个参数：（1）采样频率，（2）DFT的变换点数N，（3）需要对信号的观察时间的长度。这三个参数不能随意选取，要根据对信号频谱分析的要求进行确定。这里对信号频谱分析也有三个要求： （1）频率分辨率，（2）谱分析的频谱范围，（3）检测频率的准确性。 （1）、频谱分析的分辨率。 观察要检测的8个频率，相邻间隔最小的是第一和第二个频率，间隔是73Hz，要求DFT最少能够分辨相隔73Hz的两个频率，即要求。DFT的分辨率和对信号的观察时间有关， 。考虑到可靠性，留有富裕量，要求按键的时间大于40ms。 （2）、频谱分析的频率范围 要检测的信号频率范围是697～1633Hz，但考虑到存在语音干扰，除了检测这8个频率外，还要检测它们的二次倍频的幅度大小，波形正常且干扰小的正弦波的二次倍频是很小的，如果发现二次谐波很大，则不能确定这是