电子科大DSP第五次实验副本资料.docx

电子科技大学实 验 报 告学生姓名：学号：指导教师：一、实验室名称：数字信号处理实验室二、实验项目名称：双音多频信号的产生与检测三、实验原理：双音多频信号产生与检测：1. 双音多频信号的工作原理及产生方法双音多频（DTMF, Dual-Tone Multi-Frequency）信号及其产生与检测技术广泛应用于电话信号处理，用来完成拨号、自动重播、自助电话查询等任务。现在所用的电话，每一个数字按键都是由两种频率的单音信号组成的，这两种单音信号被分为高频带和低频带。高低频带各由四个频率组成。每一位号码均由一个低频带频率和一个高频带频率叠加形成。如图1所示是国际标准认可的数字和符号键的频率分配情况。图1 键盘的双音多频方案每个按键对应的DTMF信号为： (1)其中和分别是低频单音和高频单音。一般而言，电话中的双音多频信号有两个作用：(1) 用拨号信号去控制交换机接通被叫的用户电话机；(2) 控制电话机的各种动作，如播放留言等。2. 双音多频信号的常用产生方法如下：直接计算法利用定义式(1)，在MATLAB中用指令直接计算，但是运算量和实现成本较高。查表法该方法的思想是构造一个正/余弦函数查找表，表中所列为正/余弦函数的值，通过将表中的值以不同幅度和不同采样间隔输出，就可以得到任意幅度、任意频率的正弦或余弦波。但此方法运算量虽低，但是对存储量的要求很高。数字振荡器法数字振荡器的本质是，使用一个IIR滤波器，通过把它的极点放在单位圆上来产生振荡。利用正弦波的指数形式，可以得到正弦序列的z变换为 (2)上式在时成立，且根据Z变换的原理和性质，可知对于给定的,可以通过反Z变换，唯一确定。因此，产生正弦波就等价于用上式设计一个IIR滤波器，该IIR滤波器的传递函数为 (3)即如图所示系统。图2 二阶数字振荡器由传递函数可得输入/输出关系为： (4)假设系统的输入为单位冲激函数x[n-1]=δ[n-1],即仅当n=1时，x[n-1]=1,带入上式得 (5)在n≥2以后，y[n]可以用y[n-1]和y[n-2]算出，这是一个递归方程。双音多频信号的检测原理及方法DTMF信号检测的目的，是判断被检测信号中是否含有相应的DTMF频率对。因此完成DTMF检测的基本方法就是对信号进行频谱分析，看对应的频率分量上是否出现能量峰值。双音多频信号的检测方法主要有以下几种：FFT或DFT用FFT(DFT)是对双音多频信号进行频谱分析，由信号的幅度谱，判断信号的两个频率，最后确定相应的数字或符号。FFT是DFT的快速算法，但当DFT的变换区间较小时，FFT快速算法的效果并不明显，而且还要占用很多内存，因此不如直接用DFT合适。DFT的计算式如下所示， (6)2）基于Goertzel算法的DTMF信号检测Goertzel算法的基本思想是将每一个频点的DFT运算转换成一个匹配滤波器，通过滤波来实现DFT，对(6)式进行简单的变换可得， (7)式(7)的右边可以看作是两个序列和的卷积。如果定义，那么 (8)即可以看作是将和进行卷积，在第N个节拍输出的结果。如果将看作是一个系统（或滤波器）的单位冲激响应，那么，X[k]就是经过该系统后的第N个输出点。改变k 值，可以计算不同的频域采样点，k = 0, 1, …, N-1。所以，点DFT可以通过如图3所示的方式实现：将时域序列x[n]输入到个系统中，这些系统的冲激响应分别为，，…,,…, 将所有输出序列的第个输出点收集起来，就得到x[n]的点 DFT。图中的个滤波器构成了一组并行滤波器组，因此,Goertzel算法特别适合DFT的并行实现和模块化实现。图3 Goertzel算法原理如果只需要计算某一个单点DFT，如，那么将输入第个滤波器并输出第个时刻的响应即可。对第个滤波器的冲激响应进行变换，可以得到该滤波器的系统传输函数为， (9)由上式可知，该滤波器有一个极点在处，对应频率为，因此滤波器在频率附近的响应较大。图4 Goertzel算法原理画出该滤波器的幅频响应和相频响应，如图4所示。由图可知实际上可以看作是一个带通滤波器，通带中心为。因此，只有当输入信号的频率在附近时，系统才有较大的输出，从而检测到该频率成分。同理，图3中的每一个滤波器都是一个带通滤波器，通带中心分别为0，。基于式(9)画出滤波器的系统框图，如图5(a)所示由式(9)可得系统输入输出差分方程为 (10)容易发现，式(10)的运算量与标准DFT运算量相同，计算单点DFT仍需要次复数乘法。但是，Goertzel算法不再需要存储旋转因子，节省了大量的存