DTMF电话中的拨号信号发生与检测系统仿真设计

1、燕山大学课 程 设 计 说 明 书题目：dtmf电话中的拨号信号发生与检测系统仿真设计学院（系）： 电气工程学院 年级专业： 09级仪表（2）班 学 号： 090103020074 学生姓名： 吴季哲 指导教师： 谢平 吴晓光 教师职称： 教授 讲师 燕山大学课程设计（论文）任务书院（系）：电气工程学院 基层教学单位：自动化仪表系学 号090103020074学生姓名吴季哲专业（班级）09级仪表2班设计题目dtmf电话中的拨号信号发生与检测系统仿真设计设计技术参数采样频率是8000hz；每125ms输出一个样本；所用频率：697hz、770hz、852hz、942hz、1209hz、1336h

2、z、1477hz、633hz设计要求1、设计一个电话拨号模拟系统2、设计一个检测系统3、绘制接收到信号的频谱图4、撰写设计报告工作计划周一，设计一个电话拨号模拟系统；周二，设计一个检测系统；周三、四，绘制接收到信号的频谱图；周五，撰写设计报告。参考资料数字信号处理（第三版） 高西全，丁玉美编著数字信号处理教程 第三版 程佩青 编著数字信号处理原理及应用 谢平 王娜 林洪彬 编著指导教师签字谢平 吴晓光基层教学单位主任签字谢平说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。2012年6 月20 日 燕山大学课程设计评审意见表指导教师评语：成绩： 指导教师： 年 月 日答辩小组评语：

3、成绩： 组长： 年 月 日课程设计总成绩：答辩小组成员签字：年 月 日摘 要双音多频（dual tone multi frequency, dtmf）信号是音频电话中的拨号信号，由美国at&t贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动检测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其他通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送dtmf信号选择语音菜单进行操作。 dtmf信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了d/a

4、变换器；在接收端用a/d变换器将其转换成数字信号，并进行数字信号处理。本文主要介绍用matlab在计算机上对双音多频的信令产生及检测接收系统进行了仿真测试。关键词：双音多频，数字信号处理，傅里叶变换；目 录摘要21.绪论41.1 引言41.2 数字信号简介41.3 数字滤波器61.4 频率分析7 1.4.1 傅里叶级数及傅里叶变换 7 1.4.2 z变换及系统函数 8 1.4.3 离散系统傅里叶变换（dft）92.双音多频（dtmf）通信系统6 2.1 双音多频（dtmf）信号简介10 2.2 双音多频（dtmf）信号的产生10 2.3 双音多频（dtmf）信号的检测115.双音多频（dtmf

5、）系统的matlab实现 165.1 matlab简介 165.2 gui界面介绍175.3 双音多频（dtmf）系统的仿真19第4章 学习心得20参考文献 21程序清单 221 绪论1.1 引言双音多频（dual tone multi frequency，dtmf）信号，由美国at&t贝尔实验室研制，用于音频电话网络中的拨号信号。一方面这种信号具有非常高的拨号速度，另一方面它便于自动检测识别及电话业务的拓展，所以很快代替了原有的脉冲计数方式的拨号形式。自1876年美国贝尔发明电话以来，电话交换技术发展主要经历了三个阶段：人工交换、机电交换、电子交换。早在1878年就出现了人工供电制交换机，它

6、借助话务员进行话务连接。15年后步进制的交换机问世，这种交换机属于“直接控制”方式，用户通过话机拨号脉冲直接控制步进接续器做升降和旋转动作。从而自动完成用户间的接续。这种交换机虽然实现了自动接续，但存在着速度慢、效率低、杂音大与机械磨损严重等缺点。直到1938年发明了纵横制交换机由直接控制过渡到间接控制方式，随后又出现了电子交换机及程控交换机等。用户的拨号脉冲不再直接控制接线器动作，而先由记发器接收、存储，然后通过标志器驱动接线器，以完成用户间接续。由于标志器记取的是标识信号，所以可以不采用不稳定的脉冲计数方式，逐步发展到数字号码标识，所以快速、稳定灵活的dtmf音频拨号方式取代了脉冲计数方式

7、。本文主要阐述goertzel算法在双音拨号系统中的应用，并利用matlab的强大功能，用软件实现dtmf（双音多频）信号的检测。1.2 数字信号简介信号是反映消息的物理量，是消息的表现形式。由于非电的物理量可以通过各种传感器较容易地转换成电信号，而电信号又容易传送和控制，所以使其成为应用最广的信号。电信号是指随着时间而变化的电压或电流，因此在数学描述上可将它表示为时间的函数，并可画出其波形。信息可以通过电信号进行传送、交换、存储、提取等。电信号的形式是多种多样的，可以从不同的角度进行分类。根据信号的机性可以分为确定信号和随机信号；根据信号的周期性可分为周期信号和非周期信号；根据信号的连续性可

8、以分为连续时间信号和离散信号；在电子线路中将信号分为模拟信号和数字信号。我们每天遇到的信号，大部分都是模拟信号。离散时间信号只被定义在一种特定的时间瞬间集合上，表现为具有连续数值范围的数值序列。而数字信号在时间上和幅值上都有离散的数字。因此在实际的操作中，数字信号更容易处理。随着数字电子技术的发展，在语音处理、通信控制、图象处理等领域dsp(digital signal processor 数字信号处理器 )也得到了越来越广泛的应用。数字通信抗干扰能力强、无噪声积累，由于数字信号的幅值为有限个离散值(通常取两个幅值)，在传输过程中虽然也受到噪声的干扰，但当信噪比恶化到一定程度时，即在适当的距离

9、采用判决再生的方法，再生成没有噪声干扰的和原发送端一样的数字信号，所以可实现长距离高质量的传输。较之传统的方法，数字信号处理系统自身有着无可比拟的优点：（1）便于加密处理信息传输的安全性和保密性越来越重要，数字通信的加密处理的比模拟通信容易得多，以话音信号为例，经过数字变换后的信号可用简单的数字逻辑运算进行加密、解密处理。（2）便于存储、处理和交换数字通信的信号形式和计算机所用信号一致，都是二进制代码，因此便于与计算机联网，也便于用计算机对数字信号进行存储、处理和交换，可使通信网的管理、维护实现自动化、智能化。（3）设备便于集成化、微型数字通信采用时分多路复用，不需要体积较大的滤波器。设备中大

10、部分电路是数字电路，可用大规模和超大规模集成电路实现，因此体积小、功耗低。（4）便于构成综合数字网和综合业务数字网采用数字传输方式，可以通过程控数字交换设备进行数字交换，以实现传输和交换的综合。另外，电话业务和各种非话业务都可以实现数字化，构成综合业务数字网。（5）占用信道频带较宽一路模拟电话的频带为4khz带宽，一路数字电话约占64khz，这是模拟通信目前仍有生命力的主要原因。随着宽频带信道(光缆、数字微波)的大量利用以及数字信号处理技术的发展，数字电话的带宽问题已不是主要问题了。1.3 数字滤波器数字滤波器由数字加法器、单位延时和常数乘法器组成的一种算法或装置。数字滤波器的功能是对输入离散

11、信号的数字代码进行运算处理，即把输入序列通过一定的运算变换成输出序列，以达到改变信号频谱的目的。数学公示表示如下： （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滤波器函数

12、。无限长单位冲激响应（iir）滤波器的单位冲激响应h（n）是无限长的；系统函数h(z)在有限z平面（0z）上有极点存在；结构上存在着输出到输入的反馈，也就是结构上是递归型的。1.4 频率分析 任何给定信号的频率分析，都包含着时域信号向其频率分量的转换。之所以要在频域内对信号进行描述，因为在实现信号处理时，通常都利用以频率响应形式描述的系统。而这些频域信号处理概念时傅立叶变换是必不可少的。同时z变换对进行数字系统分析和实现也极其重要。1.4.1 傅立叶级数及傅立叶变换任何周期信号x(t)均可以表示为无穷多个谐波调谐相关的正弦波和复指数之和。具有周期t0(秒)周期信号x(t)，其基本的数学表达式为

13、具有系列定义的傅立叶级数： (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变换的定义，如下式： (

14、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)其中表示有限长序列的抽样点数，或周期序列一个周期的抽样点数。2 双音多频（dtmf）通信系统2.1 双

15、音多频（dtmf）信号过去的电话拨号时靠脉冲计数确定09这10个数字的，拨号速度慢，也不能扩展电话上网其他服务功能。现在均采用双音拨号。每一位号码由两个不同的单音频组成，所有的频率可分成高频带和低频带两组。每一位号码由一个高频信号和一个低频信号叠加形成，dtmf信号有16个编码。一般情况下，声音信号很难造成对dtmf接收器的错误触发。电话中的双音多频信号有两个作用：用拨号信号去控制交换机接通被叫的用户电话机；控制电话机的各种动作，如播放留言、语音信箱等。2.2 双音多频（dtmf）信号的产生dtmf编码器基于两个二阶数字正弦波振荡器，一个用于产生行频，一个用于产生列频。典型的dtmf信号频率范

16、围是7001700hz，选取8000hz作为采样频率，dtmf双音频信号由两个二阶数字正弦振荡器产生，一个用来产生行音频信号，另一个产生列音频信号。在电话键盘上的每次按键，会产生式（2.2-1）表示的两个音调之和。规定用8khz对dtmf信号进行采样，采样后得到时域离散信号为 (2.2-1)两个频率f1,f2唯一确定了被按压的键。具体数字对应的两个频率如表2.2.1所示。表2-1 dtmf拨号频率（第四列常被省略）频低频高1209hz1336hz1477hz633hz697hz123a770hz456b852hz789c942hz*0#d形成上面序列的方法有两种，即计算法和查表法。用计算法求正

17、弦波的序列值容易，但实际中要占用一些计算时间，影响运行速度。查表法是预先将正弦波的各序列值计算出来，寄存在存储器中，运行时只要按顺序和一定的速度取出便可。这种方法要占用一定的存储空间，但是速度快。因为采样频率是8000hz，因此要求每125ms输出一个样本，得到的序列再送到d/a变换器和平滑滤波器，输出便是连续时间的dtmf信号。dtmf信号通过电话线路送到交换机。2.3双音多频（dtmf）信号的检测在接收端，要对收到的将收到的双音多频信号进行检测，即检测两个正弦波的频率，以判断其对应的十进制数字或者符号。用数字方法进行检测，需要将接收到的模拟音频信号进行a/d变换，恢复为数字信号，然后检测其

18、中的音频频谱来确定所发送的数字。检测方法有两种，一种是用一组滤波器来提取所需频率，判断对应的数字或符号；另一种是用快速傅立叶变换（fft）算法的dft对双音多频信号进行频谱分析，有信号的幅度谱，判断信号的两个频率，最后确定对应的数字或符号。还有一种goertzel算法更合适，但本文没有做详细介绍，有兴趣的可以查阅资料。由上面可以知道一个dtmf信号是由两个不通频率f1和f2的正弦波组成，它可以用下式表示 (2.3-1)我们可以考虑通过离散傅立叶变换进行信号的频谱分析来检测离散的双音多频码。这是因为只有时域及频域都是离散的情况下，才能适合于在计算机上运算；也就是周期的离散时间信号与离散频率间的变

19、换对。周期性离散时间信号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）和都

20、是复数，x（k）也是复数，每计算一个x（k）值，需要n次复数乘法(x（n）与相乘)以及(n1)次复数加法。而x（k）一共有n个点（k从0取到n1），所以完成整个dft运算总共需要次复数乘法及n（n1）次复数加法。复数运算实际上是由实数运算来完成的，因此上式可写成： (2.3-4)由此可见，一次复数乘法需用四次实数乘法和两次实数加法；一次复数加法则需两次复数加法。因而每运算一个需要4n次复数乘法和2n2（n-1）=2(2n-1)次实数加法。 所以，直接计算dft，乘法次数和加法次数都是和成正比的，当n很大时，运算量是很可观的，例如，当n8时，dft需64次复数乘法，而当n1024时，dft所需复

21、数乘法为1，048，576次，即一百多万次复数乘法运算，这对实时性很强的信号处理来说，对计算速度的要求太高了。因而需要改进对dft的计算方法，以大大减少运算次数。仔细观察dft的运算量可以看出，利用系数的以下固有特性，可以减少dft的运算量。(1)的共轭对称性： (2)的周期性： (3) 的可约性：由此得出，这样，（1）利用这些特性，使dft运算中有些项可以合并；（2）利用的周期性和对称性，可以将长序列的dft分解为短序列的dft。而前面已经说到，dft的运算量使与成正比的，所以n越小越有利，因而小点数的dft比大点数的dft的运算量要小。基于这样的思路，快速傅立叶变换（fft）发展起来了。如

22、果按时间抽取的fft算法，系统的运算量是，运算量明显比dft运算量减少。但它有一个缺点，需要全部抽样数据到达时，才能进行计算；不太适合双音多频信号频谱的计算，满足不了双音拨号系统实时性的要求。有一种dmft的特殊算法goertzel算法，更适合适合双音多频信号频谱的计算。这篇文章介绍的是用快速傅里叶算法实现的。5 双音多频（dtmf）系统的matlab实现5.1 matlab简介matlab是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境

23、中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如c、fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。matlab 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真数字图像处理技术数字信号处理技术 通讯系统设计与仿真 财务与金融工程matlab 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 matlab 函数集）扩展了 matlab 环境，以解决这

24、些应用领域内特定类型的问题。它有如下的特点：（1）功能强的数值运算 在matlab环境中，有超过500种数学、统计、科学及工程方面的函数可使用，函数的标示自然，使得问 题和解答像数学式子一般简单明了，让使用者可全力发挥在解题方面，而非浪费在电脑操作上。（2）先进的资料视觉化功能 matlab的物件导向图形架构让使用者可执行视觉数据分，并制作高品质的图形，完成科学性或工程 性图文并茂的文章。（3）高阶但简单的程式环境 作为一种直译式的程式语言，matlab容许使用者在短时间内写完程式，所花的时间约为用 fortran 或 c 的几分之一，而且不需要编译(compile)及联结 (link) 即能

25、执行，同时包含了更多及更容易使用的内建功能。（4）开放及可延伸的架构 matlab容许使用者接触它大多数的数学原始码，检视运算法，更改现存函数，甚至加入自己的函数使 matlab成为使用者所须要的环境。（5）丰富的程式工具箱 matlab的程式工具箱融合了套装前软体的优点，与一个灵活的开放但容易操作之环境，这些工具箱提 供了使用者在特别应用领域所需之许多函数。现有工具箱有：符号运算（利用maple v的计算核心执行 ）、影像处理、统计分析、讯号处理、神经网路、模拟分析、控制系统、即时控制、系统确认、强建控 制、弧线分析、最佳化、模糊逻辑、mu分析及合成、化学计量分析。因为有了这些工具箱，大量繁

26、杂的数学运算和分析可通过调用matlab函数直接求解，大大提高了编程效率，其程序编译和执行速度远远超过了basic、fortran和c等语言。5.2 gui界面介绍 gui的广泛应用是当今计算机发展的重大成就之一，他极大地方便了非专业用户的使用。人们从此不再需要死记硬背大量的命令，取而代之的是可以通过窗口、菜单、按键等方式来方便地进行操作。而嵌入式gui具有下面几个方面的基本要求：轻型、占用资源少、高性能、高可靠性、便于移植、可配置等特点。 matlab中拥有gui界面的功能，可以将想要实现的效果通过gui界面展现出来，在matlab主界面，通过路径file-new-gui打开gui界面的编辑

27、窗口，将需要的按钮拖出，双击可进行编辑。右击找到callback,在该函数下编写该按钮要实现的功能，可实现该按钮的功能。5.3 双音多频（dtmf）系统的仿真dtmf信号的产生与识别仿真实验在matlab环境下进行，编写仿真程序，运行程序，送入电话号码，程序自动产生每一位号码数字相应的dtmf信号，产生频率普图，并送出双频声音，再用dft进行谱分析，在gui界面上显示每一位号码数字的dtmf信号的dft幅度谱。（如下图）6 学习心得本文采用快速傅里叶（fft）算法，实现了对双音多频信号的检测，在matlab环境下，对双音多频信号的解码检测过程进行了模拟仿真。本文只是完成了双音多频拨号系统的软件

28、仿真，如何实现高性价比的硬件系统是下一步需要研究和完成的工作。在做matlab界面，做好以后同时会产一个m文件，只要把代码编写在对应的控件中即可以方便的实现仿真。试验结果表明基于matlab 的dtmf 编解码具有速度快、精度高、稳定性好以及便于灵活应用等特点。而利用matlab语言，可以较好地、快速地进行算法的仿真。本文是在老师的悉心指导下完成的。每当有什么不解之处，老师总是给予详细的解答，而老师对于我们未来的生活更是给予了耐心的指导，并提出了宝贵的建议，使我终身受教。在整个设计过程中，我学习到了许多东西，同时也锻炼了独立思考和设计能力，使自己明白了学习是一个长期积累的过程，在以后的工作和生

29、活中都应该不断学习，提高自身的知识水平和综合素质。相信此次课程设计让我所体会到的人生哲理，必将成为今后生活和学习中的巨大财富，使我终生受益。在此深切的感谢这些帮助，鼓励和关心过我的老师们!参考文献数字信号处理（第三版） 高西全，丁玉美编著数字信号处理教程 第三版 程佩青 编著数字信号处理原理及应用 谢平 王娜 林洪彬 编著程序清单function varargout = untitled(varargin)gui_singleton = 1;gui_state = struct(gui_name, mfilename, . gui_singleton, gui_singleton, . gui

30、_openingfcn, untitled_openingfcn, . gui_outputfcn, untitled_outputfcn, . gui_layoutfcn, , . gui_callback, );if nargin & ischar(varargin1) gui_state.gui_callback = str2func(varargin1);end if nargout varargout1:nargout = gui_mainfcn(gui_state, varargin:);else gui_mainfcn(gui_state, varargin:);endfunct

31、ion untitled_openingfcn(hobject, eventdata, handles, varargin)handles.output = hobject; guidata(hobject, handles);function varargout = untitled_outputfcn(hobject, eventdata, handles) varargout1 = handles.output;function pushbutton1_callback(hobject, eventdata, handles)%按键4对应功能f1=697,770,852,941; f2=

32、1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(2)/8000) + sin(2*pi*n*f2(1)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,4) num=ge

33、t(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton1,string);set(handles.edit7,string,nu);endfunction pushbutton7_callback(hobject, eventdata, handles) %按键#对应功能f1=697,770,852,941; f2=1209,1336,1477

34、,1633; n=0:1023; x = sin(2*pi*n*f1(4)/8000) + sin(2*pi*n*f2(3)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000;plot(t,fftx) set(handles.text5,string,#) num=get(handles.edit7,string);n=length(num);if n=11

35、; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton7,string);set(handles.edit7,string,nu);end function pushbutton2_callback(hobject, eventdata, handles)%按键7对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(3)/8000) +

36、 sin(2*pi*n*f2(1)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,7) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框

37、已满！);set(handles.edit7,string,);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton2,string);set(handles.edit7,string,nu);endfunction pushbutton3_callback(hobject, eventdata, handles)%按键6对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(2)/8000) + sin(2*pi*n*f2(

38、3)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,6) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handle

39、s.edit7,string,);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton3,string);set(handles.edit7,string,nu);endfunction pushbutton4_callback(hobject, eventdata, handles)%按键5对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(2)/8000) + sin(2*pi*n*f2(2)/8000); sound

40、(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,5) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,

41、);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton4,string);set(handles.edit7,string,nu);endfunction pushbutton5_callback(hobject, eventdata, handles)%按键8对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(3)/8000) + sin(2*pi*n*f2(2)/8000); sound(x,8000); axes(

42、handles.axes1); plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,8) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strca

43、t(get(handles.edit7,string),get(handles.pushbutton5,string);set(handles.edit7,string,nu);endfunction pushbutton6_callback(hobject, eventdata, handles)%按键9对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(3)/8000) + sin(2*pi*n*f2(3)/8000); sound(x,8000); axes(handles.axes1);

44、 plot(n,x)%grid;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,9) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strcat(get(handles.e

45、dit7,string),get(handles.pushbutton6,string);set(handles.edit7,string,nu);endfunction pushbutton8_callback(hobject, eventdata, handles)%按键0对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(4)/8000) + sin(2*pi*n*f2(2)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid

46、;xlabel(k);ylabel(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,0) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strcat(get(handles.edit7,string),ge

47、t(handles.pushbutton8,string);set(handles.edit7,string,nu);endfunction pushbutton9_callback(hobject, eventdata, handles)%按键*对应功能f1=697,770,852,941; f2=1209,1336,1477,1633; n=0:1023; x = sin(2*pi*n*f1(4)/8000) + sin(2*pi*n*f2(1)/8000); sound(x,8000); axes(handles.axes1); plot(n,x)%grid;xlabel(k);ylab

48、el(|x(k)|) fftx=abs(fft(x); axes(handles.axes2); n=length(x); t=n/n*8000; plot(t,fftx)%grid;xlabel(k);ylabel(|x(k)|) set(handles.text5,string,*) num=get(handles.edit7,string);n=length(num);if n=11; errordlg(数字框已满！);set(handles.edit7,string,);else nu=strcat(get(handles.edit7,string),get(handles.pushbutton9,string);set(handles.edit7,string,nu);endfunction pu