DTMF信号系统的Matlab仿真毕业设计.doc

dtmf信号系统的matlab仿真摘 要双音多频（dual tone multi frequency, dtmf）信号是音频电话中的拨号信号，由美国at&t贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动检测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其他通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送dtmf信号选择语音菜单进行操作。 dtmf信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了d/a变换器；在接收端用a/d变换器将其转换成数字信号，并进行数字信号处理。本文着重介绍了一种快速的频率检测算法geortzel算法的原理，以及该算法在双音多频拨号检测系统中的应用。并根据此算法设计出了相应的滤波器，利用matlab在计算机上对双音多频的信令产生及检测接收系统进行了仿真测试。关键词：双音多频，goertzel算法，数字信号处理，傅里叶变换；abstract dual-tone multi-frequency signaling (dtmf) is tone dialing signals , developed by at&t , and used in telephone networks. this signal system has a high rate of speed dial and easy automatic identification, soon took the place of the original dial-up form of pulse counting method used. this dual-tone multi-frequency signaling system is not only used in the telephone network , can also be used in other communication systems to transfer decimal data, e-mail and banking system. these system users can select send dtmf signals by telephone voice menu for operation.dtmf signaling system is a typical small signal processing system , it is to use digital methods to produce the analog signal and transmission , which also used the d / a converter; at the receiving end with d / a converter to convert it into a ditigal signal, and digital signal processing. this article specially introduced the theory of a fast check algorithm -geortzel algorithm, and its application in detecting the double tone multi frequency system. we designed the corresponding filter according to the algorithm the same time，and we imitate the whole dialing process of the system through the matlab in the computer at last.keywords： dual-tone multi-frequency (dtmf); geortzel algorithm; digital signal processing(dsp); fourier transfer;目 录摘要 abstract 1.绪论11.1 引言11.2 数字信号简介11.3 数字滤波器31.4 频率分析4 1.4.1 傅里叶级数及傅里叶变换 4 1.4.2 z变换及系统函数 4 1.4.3 离散系统傅里叶变换（dft）52.双音多频（dtmf）通信系统6 2.1 双音多频（dtmf）信号简介6 2.2 双音多频（dtmf）信号的产生6 2.3 双音多频（dtmf）信号的检测73.戈泽尔(goertzel)算法 113.1 戈泽尔(goertzel)算法原理 113.2 戈泽尔(goertzel)算法程序 134.检测双音多频信号（dtmf）的dft参数 144.1 频谱分析的分辨率 144.2 频谱分析的频率范围 144.3 检测频率的准确性 145.双音多频（dtmf）系统的matlab实现 16 5.1 matlab简介 165.2 双音多频（dtmf）系统的仿真176.结束语 21致谢 22参考文献 231 绪论1.1 引言双音多频（dual tone multi frequency，dtmf）信号，由美国at&t贝尔实验室研制，用于音频电话网络中的拨号信号。一方面这种信号具有非常高的拨号速度，另一方面它便于自动检测识别及电话业务的拓展，所以很快代替了原有的脉冲计数方式的拨号形式。自1876年美国贝尔发明电话以来，电话交换技术发展主要经历了三个阶段：人工交换、机电交换、电子交换。早在1878年就出现了人工供电制交换机，它借助话务员进行话务连接。15年后步进制的交换机问世，这种交换机属于“直接控制”方式，用户通过话机拨号脉冲直接控制步进接续器做升降和旋转动作。从而自动完成用户间的接续。这种交换机虽然实现了自动接续，但存在着速度慢、效率低、杂音大与机械磨损严重等缺点。直到1938年发明了纵横制交换机由直接控制过渡到间接控制方式，随后又出现了电子交换机及程控交换机等。用户的拨号脉冲不再直接控制接线器动作，而先由记发器接收、存储，然后通过标志器驱动接线器，以完成用户间接续。由于标志器记取的是标识信号，所以可以不采用不稳定的脉冲计数方式，逐步发展到数字号码标识，所以快速、稳定灵活的dtmf音频拨号方式取代了脉冲计数方式。本文主要阐述goertzel算法在双音拨号系统中的应用，并利用matlab的强大功能，用软件实现dtmf（双音多频）信号的检测。1.2 数字信号简介信号是反映消息的物理量，是消息的表现形式。由于非电的物理量可以通过各种传感器较容易地转换成电信号，而电信号又容易传送和控制，所以使其成为应用最广的信号。电信号是指随着时间而变化的电压或电流，因此在数学描述上可将它表示为时间的函数，并可画出其波形。信息可以通过电信号进行传送、交换、存储、提取等。电信号的形式是多种多样的，可以从不同的角度进行分类。根据信号的随机性可以分为确定信号和随机信号；根据信号的周期性可分为周期信号和非周期信1号；根据信号的连续性可以分为连续时间信号和离散信号；在电子线路中将信号分为模拟信号和数字信号。我们每天遇到的信号，大部分都是模拟信号。离散时间信号只被定义在一种特定的时间瞬间集合上，表现为具有连续数值范围的数值序列。而数字信号在时间上和幅值上都有离散的数字。因此在实际的操作中，数字信号更容易处理。随着数字电子技术的发展，在语音处理、通信控制、图象处理等领域dsp(digital signal processor 数字信号处理器 )也得到了越来越广泛的应用。数字通信抗干扰能力强、无噪声积累，由于数字信号的幅值为有限个离散值(通常取两个幅值)，在传输过程中虽然也受到噪声的干扰，但当信噪比恶化到一定程度时，即在适当的距离采用判决再生的方法，再生成没有噪声干扰的和原发送端一样的数字信号，所以可实现长距离高质量的传输。较之传统的方法，数字信号处理系统自身有着无可比拟的优点：（1）便于加密处理信息传输的安全性和保密性越来越重要，数字通信的加密处理的比模拟通信容易得多，以话音信号为例，经过数字变换后的信号可用简单的数字逻辑运算进行加密、解密处理。（2）便于存储、处理和交换数字通信的信号形式和计算机所用信号一致，都是二进制代码，因此便于与计算机联网，也便于用计算机对数字信号进行存储、处理和交换，可使通信网的管理、维护实现自动化、智能化。（3）设备便于集成化、微型数字通信采用时分多路复用，不需要体积较大的滤波器。设备中大部分电路是数字电路，可用大规模和超大规模集成电路实现，因此体积小、功耗低。（4）便于构成综合数字网和综合业务数字网采用数字传输方式，可以通过程控数字交换设备进行数字交换，以实现传输和交换的综合。另外，电话业务和各种非话业务都可以实现数字化，构成综合业务数字网。2（5）占用信道频带较宽一路模拟电话的频带为4khz带宽，一路数字电话约占64khz，这是模拟通信目前仍有生命力的主要原因。随着宽频带信道(光缆、数字微波)的大量利用以及数字信号处理技术的发展，数字电话的带宽问题已不是主要问题了。1.3 数字滤波器数字滤波器由数字加法器、单位延时和常数乘法器组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，即把输入序列通过一定的运算变换成输出序列，以达到改变信号频谱的目的。数学公示表示如下： （1.3-1)通常可以用两种方法实现数字滤波器：一种方法是把滤波器所完成的运算编成程序并让计算机执行，也就是采用就算机软件实现；另一种方法是设计专用的数字硬件、专用的数字信号处理器或采用通用的数字信号处理器来实现。存在着两种数字滤波器的的基本类型：无限长单位冲击响应（iir）滤波器和有限长单位冲击响应（fir）滤波器。它们的系统函数分别为： (1.3-2) (1.3-3)式（1.3-2）中的h(z)称为n阶iir滤波器函数式，（1.3-3）中的h(z)称为n-1阶fir滤波器函数。无限长单位冲激响应（iir）滤波器的单位冲激响应h（n）是无限长的；系统函数h(z)在有限z平面（0z）上有极点存在；结构上存在着输出到输入的反馈，也就是结构上是递归型的。1.4 频率分析 3任何给定信号的频率分析，都包含着时域信号向其频率分量的转换。之所以要在频域内对信号进行描述，因为在实现信号处理时，通常都利用以频率响应形式描述的系统。而这些频域信号处理概念时傅立叶变换是必不可少的。同时z变换对进行数字系统分析和实现也极其重要。1.4.1 傅立叶级数及傅立叶变换任何周期信号x(t)均可以表示为无穷多个谐波调谐相关的正弦波和复指数之和。具有周期t0(秒)周期信号x(t)，其基本的数学表达式为具有系列定义的傅立叶级数： (1.4.1-1)式中ck为傅立叶级系数，第k个傅立叶系数ck可以表示为 (1.4.1-2)用周期信号的傅立叶级数通过极限的方法导出的非周期信号的频谱表示式，称为傅立叶变换：傅立叶逆变换： (1.4.1-3)傅立叶正变换： = (1.4.1-4)1.4.2 z变换及系统函数对于离散时间系统来说，与拉普拉斯对应的变换是是z变换，z变换提供了离散时间信号和系统的频域描述，并为数字滤波器的设计和实现，提供了有力的工具。序列x(n)的z变换定义为 ： (1.4.2-1)式中z是一个复变量，定义中对n求和是在和之间求和，称为双边z变换。还有一种称为单边z变换的定义，如下式： 4 (1.4.2-2)将系统的单位脉冲响应h（n）进行z变换，得到h（z）,一般称h（z）为系统的系统函数，它表征了系统的频域特性。对n阶差分方程进行z变换，得到系统函数的一般表示式 (1.4.2-3)令z=ej即得系统的传输函数h(ej)。1.4.3 离散系统傅立叶变换（dft）离散傅立叶变换是针对有限长序列或周期序列从存在的，相当把序列的连续傅立叶变换加以抽样，频率的离散化造成时间函数也呈周期，故级数应限制在一个周期内。令，则，我们即可得到离散傅立叶变换对为： (1.4.31) (1.4.32)其中表示有限长序列的抽样点数，或周期序列一个周期的抽样点数。52 双音多频（dtmf）通信系统2.1 双音多频（dtmf）信号过去的电话拨号时靠脉冲计数确定09这10个数字的，拨号速度慢，也不能扩展电话上网其他服务功能。现在均采用双音拨号。每一位号码由两个不同的单音频组成，所有的频率可分成高频带和低频带两组。每一位号码由一个高频信号和一个低频信号叠加形成，dtmf信号有16个编码。一般情况下，声音信号很难造成对dtmf接收器的错误触发。电话中的双音多频信号有两个作用：用拨号信号去控制交换机接通被叫的用户电话机；控制电话机的各种动作，如播放留言、语音信箱等。2.2 双音多频（dtmf）信号的产生dtmf编码器基于两个二阶数字正弦波振荡器，一个用于产生行频，一个用于产生列频。典型的dtmf信号频率范围是7001700hz，选取8000hz作为采样频率，dtmf双音频信号由两个二阶数字正弦振荡器产生，一个用来产生行音频信号，另一个产生列音频信号。在电话键盘上的每次按键，会产生式（2.2-1）表示的两个音调之和。规定用8khz对dtmf信号进行采样，采样后得到时域离散信号为 (2.2-1)两个频率f1,f2唯一确定了被按压的键。具体数字对应的两个频率如表2.2.1所示。6表2-1 dtmf拨号频率（第四列常被省略）频低频高1209hz1336hz1477hz633hz697hz123a770hz456b852hz789c942hz*0#d形成上面序列的方法有两种，即计算法和查表法。用计算法求正弦波的序列值容易，但实际中要占用一些计算时间，影响运行速度。查表法是预先将正弦波的各序列值计算出来，寄存在存储器中，运行时只要按顺序和一定的速度取出便可。这种方法要占用一定的存储空间，但是速度快。因为采样频率是8000hz，因此要求每125ms输出一个样本，得到的序列再送到d/a变换器和平滑滤波器，输出便是连续时间的dtmf信号。dtmf信号通过电话线路送到交换机。2.3双音多频（dtmf）信号的检测在接收端，要对收到的将收到的双音多频信号进行检测，即检测两个正弦波的频率，以判断其对应的十进制数字或者符号。用数字方法进行检测，需要将接收到的模拟音频信号进行a/d变换，恢复为数字信号，然后检测其中的音频频谱来确定所发送的数字。检测方法有两种，一种是用一组滤波器来提取所需频率，判断对应的数字或符号；另一种是用快速傅立叶变换（fft）算法的dft对双音多频信号进行频谱分析，有信号的幅度谱，判断信号的两个频率，最后确定对应的数字或符号。以下章节对两种方法进行分析并对goertzel算法进行详细的推导。由上面可以知道一个dtmf信号是由两个不通频率f1和f2的正弦波组成，它可以用下式表示 (2.3-1)7我们可以考虑通过离散傅立叶变换进行信号的频谱分析来检测离散的双音多频码。这是因为只有时域及频域都是离散的情况下，才能适合于在计算机上运算；也就是周期的离散时间信号与离散频率间的变换对。周期性离散时间信号x（n）的离散傅立叶变换周期性离散频率函数x（k）两相邻谱线分量之间的角频率增量与周期之间的关系可表示为， 代表信号的基频。取样频率与取样周期t的关系是 取样角频率 t为时域取样间隔，在一个周期内取样点数为n。在自变量为t及f的情况下，在一个时域中对函数进行取样，两取样点间增量的倒数，必是另一个域中函数的周期。现序列的周期为nt，所以对频谱取样的间距是。以数字频率表示时，则频谱间距是 k0，1，n1令并称之为因子。离散傅立叶变换（dft）可写成如下形式：， 0 kn1 (2.3-2)dft的表达式可以直接写成：， k0，1，n1 ( 2.3-3) 一般来说，x（n）和都是复数，x（k）也是复数，每计算一个x（k）值，需要n次复数乘法(x（n）与相乘)以及(n1)次复数加法。而x（k）一共有n个点（k从0取到n1），所以完成整个dft运算总共需要次8复数乘法及n（n1）次复数加法。复数运算实际上是由实数运算来完成的，因此上式可写成： (2.3-4)由此可见，一次复数乘法需用四次实数乘法和两次实数加法；一次复数加法则需两次复数加法。因而每运算一个需要4n次复数乘法和2n2（n-1）=2(2n-1)次实数加法。 所以，直接计算dft，乘法次数和加法次数都是和成正比的，当n很大时，运算量是很可观的，例如，当n8时，dft需64次复数乘法，而当n1024时，dft所需复数乘法为1，048，576次，即一百多万次复数乘法运算，这对实时性很强的信号处理来说，对计算速度的要求太高了。因而需要改进对dft的计算方法，以大大减少运算次数。仔细观察dft的运算量可以看出，利用系数的以下固有特性，可以减少dft的运算量。(1)的共轭对称性： (2)的周期性： (3) 的可约性：由此得出，这样，（1）利用这些特性，使dft运算中有些项可以合并；（2）利用的周期性和对称性，可以将长序列的dft分解为短序列的dft。而前面已经说到，dft的运算量使与成正比的，所以n越小越有利，因而小点数的dft比大9点数的dft的运算量要小。基于这样的思路，快速傅立叶变换（fft）发展起来了。如果按时间抽取的fft算法，系统的运算量是，运算量明显比dft运算量减少。但它有一个缺点，需要全部抽样数据到达时，才能进行计算；不适合双音多频信号频谱的计算，满足不了双音拨号系统实时性的要求。因此开发出了适合于dmft的特殊算法goertzel算法。103 戈泽尔(goertzel)算法3.1 戈泽尔(goertzel)算法原理戈泽尔算法利用dft中的旋转因子的周期性，将dft的运算转换成一种线性滤波运算。下面推导戈泽尔算法的计算公式和实现结构。假设长度为n的序列x（n）的n点dft用x（k）表示，因为 （3.1-1）因此 k=0,1,12,n-1 (3.1-2)按照上式定义序列 （3.1-3）这个方程可以解释为有限长序列x(n)，与序列的卷积，因此表示为 （3.1-4）令 （3.1-5）则 （3.1-6）由上式，将看成是序列x（n）通过单位脉冲响应为的滤波器的输出，根据式（3.1-2）和式（3.1-3），得到 (3.1-7)那么，x（n）的dft的第k点就是序列x(n)通过滤波器输出的第n=n点样值。那么n点dft就是这n个滤波器分别对序列x(n)的响应序列的第n点输出。对式（3.1-5）进行z变换，得到滤波器的系统函数11 （3.1-8）该滤波器是一个一阶纯极点滤波器，极点为，极点频率为。该滤波器的结构图如图3-1所示，戈泽尔算法的原理方框图如图3-2所示。图3-1单极点滤波器图3-2戈泽尔算法的滤波器组在图3-1中存在一次复乘算法，为了避免复乘算法，将一阶纯极点滤波器变为二阶滤波器，推导如下 （3.1-9）按照图3-2所示。再按照该结构图，可以用两个差分方程表示该二阶滤波器，即 （3.1-10）= （3.1-11）对于音调检测，只有的幅值平方是需要的，通过实现这种需求，可以使算法作进一步简化，由（3.1-11）式，的平方幅值可以计算为：12 （3.1-12）这样输入信号是实序列，用式（3.1-10）计算中间变量，用式（3.1-12）计算输出信号的幅度，避免了复数的运算。由此得到=。因为有8种音频要检测，所以需要8个式（3.1-8）表示的滤波器，或者8个式（3.1-9）表示的滤波器。8个滤波器的中心频率分别对应8种音频。3.2 goertzel算法程序下面是按（3.1-12）式编写的goertzel算法子程序。它根据输入序列x（n）和指定的dft样本的序号k计算待求的dft样本x。其核心语句只有四句，具体如下：function x=gfft(x,k)n=length(x);x1=x,0;d1=2*cos(2*pi*k/n)v=filter(1,1,-d1,1,x1);w=exp(-i*2*pi*k/n);x=v(n+1)-w*v(n);134 检测双音多频信号（dtmf）的dft参数 用dft检测模拟dtmf信号所含有的两个音频频率，即为用dft对模拟信号进行频谱分析的问题，要确定三个参数：采样频率fs，dtf的变换点数n，需要对时间的观察长度tp。4.1 频谱分析的分辨率 观察表2-1要检测的8个频率，相邻间隔最小的是第一个和第二个频率，间隔是73hz,要求dft至少能够分辨相隔73hz的两个信号，即要求fmin=73hz。dft的分辨率和对信号的观察时间tp有关，tpmin=1/f=1/73=13.7ms。考虑到可靠性，应留有富余量，要求按键的时间在40ms以上。4.2 频谱分析的频率范围 要检测信号的频率范为6971633hz，但考虑存在语音干扰，除了检测这8个频率外，还要检测它们的二次倍频的幅度大小。波形正常且干扰小的正弦波的二次倍频是很小的。如果发现二次谐波很大，则认为不是dtmf信号。这样频谱分析的范围为6973266hz。按照采样定理，信号的最高频率不能超过折叠频率，即0.5fs3266hz，由此要求最小采样频率应为6.53khz。总系统已经规定fs=8khz，因此一定满足对频谱分析的要求。按照tpmin=13.7ms,fs=8khz,算出对信号最少的采样点数为nmin=tpmin*fs110。4.3 检测频率的准确性这是一个用dft检测正弦波频率是否准确的问题。序列的n点dft是对序列频谱函数在0区间的n点等间隔采样，如果是一个周期序列，截取周期序列的整数倍周期，进行dft，其采样点刚好在周期信号的频率上，dft的幅度最大处就是信号的准确频率。分析这些dtmf信号，不可能经过采样得到周期序列，因此存在检测频率的准确性问题。 dft的频率采样点频率为（k=0,1,2,-,n-1），相应的模拟域采样14点频率为（k=0,1,2,-,n-1），希望选择一个合适的n，使用该公式算出的能接近要检测的频率，或者用8个频率中的任一个频率代入公式中时，得到的k值最接近整数值，这样虽然用幅度最大点检测的频率有误差，但可以准确判断所对应的dtmf频率，即可以准确判断所对应的数字或符号。经过分析研究认为n205是最好的。按照8khz，n205，算出8个频率及其二次谐波对应k值，和k取整数时的频率误差见表2。表4-1频率误差表8个基频hz最近的整数k值dft的k值绝对误差二次谐波hz对应的k值最近的整数k值绝对误差69717.861180.139139435.024350.02477019.531200.269154038.692390.30885221.833220.167170442.813430.18794124.113240.113188247.285470.285120930.981310.019241860.752610.248133634.235340.235267267.134670.134147737.848380.152295474.219740.219163341.846420.154326682.058820.058通过以上分析，确定8khz，n205，。155 双音多频（dtmf）系统的matlab实现5.1 matlab简介matlab是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如c、fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。matlab 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真数字图像处理技术数字信号处理技术 通讯系统设计与仿真 财务与金融工程matlab 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 matlab 函数集）扩展了 matlab 环境，以解决这些应用领域内特定类型的问题。它有如下的特点：（1）功能强的数值运算 在matlab环境中，有超过500种数学、统计、科学及工程方面的函数可使用，函数的标示自然，使得问 题和解答像数学式子一般简单明了，让使用者可全力发挥在解题方面，而非浪费在电脑操作上。（2）先进的资料视觉化功能 matlab的物件导向图形架构让使用者可执行视觉数据分，并制作高品质的图形，完成科学性或工程 性图文并茂的文章。16（3）高阶但简单的程式环境 作为一种直译式的程式语言，matlab容许使用者在短时间内写完程式，所花的时间约为用 fortran 或 c 的几分之一，而且不需要编译(compile)及联结 (link) 即能执行，同时包含了更多及更容易使用的内建功能。（4）开放及可延伸的架构 matlab容许使用者接触它大多数的数学原始码，检视运算法，更改现存函数，甚至加入自己的函数使 matlab成为使用者所须要的环境。（5）丰富的程式工具箱 matlab的程式工具箱融合了套装前软体的优点，与一个灵活的开放但容易操作之环境，这些工具箱提 供了使用者在特别应用领域所需之许多函数。现有工具箱有：符号运算（利用maple v的计算核心执行 ）、影像处理、统计分析、讯号处理、神经网路、模拟分析、控制系统、即时控制、系统确认、强建控 制、弧线分析、最佳化、模糊逻辑、mu分析及合成、化学计量分析。因为有了这些工具箱，大量繁杂的数学运算和分析可通过调用matlab函数直接求解，大大提高了编程效率，其程序编译和执行速度远远超过了basic、fortran和c等语言。5.2 双音多频（dtmf）系统的仿真dtmf信号的产生与识别仿真实验在matlab环境下进行，编写仿真程序，运行程序，送入8位电话号码，程序自动产生每一位号码数字相应的dtmf信号，并送出双频声音，再用dft进行谱分析，显示每一位号码数字的dtmf信号的dft幅度谱，安照幅度谱的最大值确定对应的频率，再安照频率确定每一位对应的号码数字，最后输出8位电话号码，程序名为exp8。下面用matlab程序对该系统进行模拟，程序分四段：第一段（第27行）设置参数，并读入8位电话号码；第二段（第920行）根据输入的8位电话号码产生时域离散dtmf信号，并连续发出8位号码对应的双音频声音；第三段（第2225行）对时域离散dtmf信号进行频率检测，画出幅度谱；第四段（第2633行）根据幅度谱的两个峰值，分别查找并确定输入的8位电话号码。程序如下：17% dtmf双频拨号信号的生成和检测程序：%clear all;clc;tm=1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68; % dtmf信号代表的16个数n=205;k=18,20,22,24,31,34,38,42;%8个基频对应的8个k值f1=697,770,852,941; % 行频率向量f2=1209,1336,1477,1633; % 列频率向量tn=input(键入8位电话号码= ); % 输入8位数字tnr=0; %接收端电话号码初值为零for l=1:8; %分别对每位号码数字处理：产生信号，发声，检测 d=fix(tn/10(8-l); %计算出第l位号码数字 tn=tn-d*10(8-l); for p=1:4; for q=1:4; if tm(p,q)=abs(d); break,end % 检测与l位号码相符的列号q endif tm(p,q)=abs(d); break,end % 检测与l位号码相符的行号p end n=0:1023; % 为了发声，加长序列 x = sin(2*pi*n*f1(p)/8000) + sin(2*pi*n*f2(q)/8000);% 构成双音频信号 sound(x,8000); % 发出声音 pause(0.1) %相邻信号响声之间加0.1秒停顿 % 接收检测端的程序 x=goertzel(x(1:n),k+1); % 用goertzel算法计算8点dft样本 val = abs(x); % 列出八点dft的摸 subplot(4,2,l); stem(k,val,.);grid;xlabel(k);ylabel(|x(k)|) % 画出8点dft(k)的幅度 axis(10 50 0 120) limit = 80; %基频检测门限为80 for s=5:8; if val(s) limit, break, end % 查找列号 end for r