数字信号处理(第三版)上机实验.ppt

　　2. 电话中的双音多频（DTMF）信号的产生与检测 　　(1) 双音多频信号的产生。 假设时间连续的 DTMF信号用x(t)=sin(2πf1t)+sin(2πf2t)表示， 式中f1和f2是按照表8.6.1选择的两个频率， f1代表低频带中的一个频率， f2代表高频带中的一个频率。 显然， 采用数字方法产生DTMF信号， 优点是方便而且体积小。 下面介绍采用数字方法产生DTMF信号。 规定用8 kHz对DTMF信号进行采样， 采样后得到时域离散信号为 形成上面序列的方法有两种， 即计算法和查表法。 用计算法求正弦波的序列值容易， 但实际中要占用一些计算时间， 影响运行速度。 查表法是预先将正弦波的各序列值计算出来， 寄存在存储器中， 运行时只要按顺序和一定的速度取出便可。 这种方法要占用一定的存储空间， 但是速度快。  因为采样频率是8000 Hz， 因此要求每125 ms输出一个样本， 得到的序列再送到D/A变换器和平滑滤波器， 输出便是连续时间的DTMF信号。 DTMF信号通过电话线路送到交换 机。 　　（2） 双音多频信号的检测。 在接收端， 要对收到的双音多频信号进行检测， 检测两个正弦波的频率是多少， 以判 断所对应的十进制数字或者符号。 显然， 这里仍然要用数字方法进行检测， 因此要将收到的时间连续DTMF信号经过A/D变换， 变成数字信号进行检测。 检测的方法有两种， 一种是用一组滤波器提取所关心的频率， 根据有输出信号的两个滤波器判断相应的数字或符号。 另一种是用DFT（FFT）对双音多频信号进行频谱分析， 由信号的幅度谱， 判断信号的两个频率， 最后确定相应的数字或符号。 　　当检测的音频数目较少时， 用滤波器组实现更合适。 FFT是DFT的快速算法， 但当DFT的变换区间较小时， FFT快速算法的效果并不明显， 而且还要占用很多内存， 因此不如直接用DFT合适。 下面介绍Goertzel算法， 这种算法的实质是直接计算DFT的一种线性滤波方法。 这里略去Goertzel算法的介绍（请参考文献［19］）， 可以直接调用MATLAB信号处理工具箱中戈泽尔算法的函数Goertzel， 计算N点DFT的几个感兴趣的频点的值。 　　3. 检测DTMF信号的DFT参数选择 　　用DFT检测模拟DTMF信号所含有的两个音频频率， 是一个用DFT对模拟信号进行频谱分析的问题。 根据第3章用DFT对模拟信号进行谱分析的理论， 确定三个参数： ① 采样频率Fs； ② DFT的变换点数N； ③ 需要对信号的观察时间的长度Tp。 这三个参数不能随意选取， 要根据对信号频谱分析的要求进行确定。 这里对信号频谱分析也有三个要求： ① 频谱分析的分辨率； ② 频谱分析的频率范围； ③ 检测频率的准确性。 　　(1) 频谱分析的分辨率。 观察要检测的8个频率， 相邻间隔最小的是第一和第二个频率， 间隔是73 Hz， 要求 DFT最少能够分辨相隔73 Hz的两个频率， 即要求Fmin=73 Hz。 DFT的分辨率与对信号的观察时间Tp有关， Tp　min=1/F=1/73 =13.7 ms。 考虑到可靠性， 留有富裕量， 要求按键的时间大于40 ms。 　　(2) 频谱分析的频率范围。 要检测的信号频率范围是697～1633 Hz， 但考虑到存在语音干扰， 除了检测这8个频率外， 还要检测它们的二次倍频的幅度大小， 波形正常且干扰小的正弦波的二次倍频是很小的， 如果发现二次谐波很大， 则不能确定这是DTMF信号。 这样频谱分析的频率范围为697～3266 Hz。 按照采样定理， 最高频率不能超过折叠频率， 即0.5Fs≥3622 Hz， 由此要求最小的采样频率应为7.24 kHz。 因为数字电话总系统已经规定Fs＝8 kHz， 因此对频谱分析范围的要求是一定满足的。按照Tp min=13.7 ms， Fs＝8 kHz， 算出对信号最少的采样点数为Nmin=Tp min · Fs≈110。 　　(3) 检测频率的准确性。 这是一个用DFT检测正弦波频率是否准确的问题。 序列的N点DFT是对序列频谱函数在0～2π区间的N点等间隔采样。 如果是一个周期序列， 截取周期序列的整数倍周期， 进行DFT， 其采样点刚好在周期信号的频率上， DFT的幅度最大处就是信号的准确频率。 分析这些DTMF信号， 发现不可能经过采样而得到周期序列， 因此存在检测频率的准确性问题。 　　DFT的频率采样点频率为ωk=2πk（k=0， 1， 2， …， N－1）， 相应