雙音多頻(DTMF)信號(hào)的MATLAB仿真

1、數(shù)字信號(hào)處理課程設(shè)計(jì)雙音多頻（DTMF）信號(hào)的MATLAB仿真院 系：閩江學(xué)院計(jì)算機(jī)科學(xué)系專 業(yè)：通信工程組 長(zhǎng): 蔡小紅（3121102133）組 員：王垚天（3121102119） 何杰坤（3121102131） 王曉清（3121102134） 蘇月琴（3121102136） 羅 燁（3121102216）指導(dǎo)老師：方 薈 時(shí) 間：2015年1月5日-2015年1月9日 1雙音多頻（DTMF）信號(hào)的MATLAB仿真任務(wù)分配安排表任務(wù)分配負(fù)責(zé)人資料查找：電話中雙頻撥號(hào)的頻率分配，以及雙音多頻信號(hào)的產(chǎn)生方法。蘇月琴、羅 燁課程設(shè)計(jì)方案、思路和具體步驟何杰坤、蔡小紅基本代碼編寫王垚天、何杰坤gu

2、i編寫、界面設(shè)計(jì)王垚天、王曉清代碼調(diào)試蔡小紅、羅 燁實(shí)驗(yàn)報(bào)告編寫王曉清、蘇月琴實(shí)驗(yàn)答辯王垚天、何杰坤、蔡小紅王曉清、蘇月琴、羅 燁目 錄一、摘要1二、引言1三、設(shè)計(jì)要求13.1設(shè)計(jì)要求13.2設(shè)計(jì)步驟13.3設(shè)計(jì)程序23.4設(shè)計(jì)GUI23.5原理33.5.1電話中的DTMF信號(hào)的組成33.5.2電話中的雙音多頻（DTMF）信號(hào)的產(chǎn)生43.5.3 Goerztel函數(shù)與DTMF信號(hào)的產(chǎn)生5四、設(shè)計(jì)過(guò)程5圖4.1設(shè)計(jì)流程圖6五、調(diào)試與結(jié)果75.1運(yùn)行結(jié)果如下圖所示75.2調(diào)試中遇到的問(wèn)題及解決方法7六、分析和總結(jié)76.1小組總結(jié)76.2成員總結(jié)8七、參考文獻(xiàn)9八、附錄：主要程序10一、摘要DTMF

3、信號(hào)系統(tǒng)是一個(gè)典型的小型信號(hào)處理系統(tǒng)，它要用數(shù)字方法產(chǎn)生模擬信號(hào)并進(jìn)行傳輸，應(yīng)用D/A變換器，采用快速的頻率檢測(cè)算法戈澤爾（GEORTZEL）算法的原理，以及該算法在雙音多頻撥號(hào)檢測(cè)系統(tǒng)中的應(yīng)用。對(duì)讀入的電話號(hào)碼，根據(jù)鍵入的號(hào)碼數(shù)產(chǎn)生時(shí)域離散雙音多頻（DTMF）信號(hào)，并連續(xù)發(fā)出號(hào)碼對(duì)應(yīng)的雙音頻聲音。同時(shí)進(jìn)行時(shí)域離散DTMF信號(hào)進(jìn)行檢測(cè)，畫出幅度譜，并用MATLAB在計(jì)算機(jī)上對(duì)雙音多頻的信號(hào)產(chǎn)生及檢測(cè)接收系統(tǒng)進(jìn)行了仿真測(cè)試。 關(guān)鍵詞 雙音多頻戈澤爾算法DFT二、引言DTMF（DoubleToneMulitiFrequency,雙音多頻）作為實(shí)現(xiàn)電話號(hào)碼快速可靠傳輸?shù)囊环N技術(shù)，具有很強(qiáng)的抗干擾能力

4、和較高的傳輸速度，廣泛用于電話通信系統(tǒng)中。另外，它也可以在通信系統(tǒng)中廣泛用來(lái)實(shí)現(xiàn)各種數(shù)據(jù)和語(yǔ)音等信息的遠(yuǎn)程傳輸，技術(shù)非常成熟、專業(yè)芯片種類繁多，在民用、金融、工業(yè)等領(lǐng)域的應(yīng)用非常普及，甚至有通過(guò)軟件來(lái)產(chǎn)生和檢測(cè)DTMF信令方面的應(yīng)用。但隨著現(xiàn)代工業(yè)的持續(xù)發(fā)展，利用DTMF信令作為通信系統(tǒng)中的控制信號(hào)的通信設(shè)備正遭受著前所未有的各種日益復(fù)雜的干擾，尤其是與DTMF信令的同頻干擾，使通信設(shè)備頻頻出現(xiàn)誤動(dòng)作等嚴(yán)重事故，所以消除與DTMF同頻干擾已是利用該控制方式控制的通信系統(tǒng)不得不解決的一個(gè)問(wèn)題。三、設(shè)計(jì)要求3.1設(shè)計(jì)要求（1）按一個(gè)數(shù)字鍵如“1”，則產(chǎn)生頻率為697Hz和1209Hz的兩個(gè)正弦波，

5、并相加；（2）可用查表法求數(shù)字鍵對(duì)應(yīng)的頻率，為此先建立撥號(hào)數(shù)字表矩陣；（3）電話音頻信號(hào)在數(shù)字信號(hào)處理時(shí)，取樣頻率為8KHz，每個(gè)數(shù)字信號(hào)持續(xù)時(shí)間為100ms，后面加上100ms的間隔（用0表示）；（4）設(shè)計(jì)好各按鍵提示音，并能播放演示；（5）編制GUI用戶界面；3.2設(shè)計(jì)步驟（1）總體設(shè)計(jì)搭建框架，確定人機(jī)對(duì)話界面，確定輸入輸出函數(shù)接口；（2）編寫主要程序，包括戈澤爾算法和N點(diǎn)DFT；（3）編寫實(shí)現(xiàn)輸入輸出功能的程序；（4）編寫波形輸出程序；（5）用matlab的相應(yīng)函數(shù)進(jìn)行檢驗(yàn)；（6）編寫GUI，并綁定各功能，加以實(shí)現(xiàn)。3.3設(shè)計(jì)程序讀入電話號(hào)碼，根據(jù)鍵入的電話號(hào)碼產(chǎn)生時(shí)域離散雙音多頻（D

6、TMF）信號(hào)，并連續(xù)發(fā)出號(hào)碼對(duì)應(yīng)的雙音頻聲音。對(duì)時(shí)域離散DTMF信號(hào)進(jìn)行頻率檢測(cè)，畫出幅度譜。根據(jù)幅度譜的兩個(gè)峰值，分別查找并確定輸入電話號(hào)碼，并將其撥號(hào)并響應(yīng)。3.4設(shè)計(jì)GUI實(shí)現(xiàn)更加友好的人機(jī)對(duì)話的界面，加入了界面窗口，各種數(shù)據(jù)可以從窗口直接輸入，并且設(shè)置開始按鈕和清除按鈕，可以從窗口選擇六位或者八位電話號(hào)碼的輸入，在窗口中可以直接顯示出電話號(hào)碼基頻的幅值，簡(jiǎn)潔、方便。 圖3.4.1具體的界面窗口。3.5原理 雙音多頻（Dual Tone Multi Frequency, DTMF）信號(hào)是音頻電話中的撥號(hào)信號(hào)，由美國(guó)AT&T貝爾公司實(shí)驗(yàn)室研制，并用于電話網(wǎng)絡(luò)中。這種信號(hào)制式具有很高的撥號(hào)速

7、度，且容易自動(dòng)監(jiān)測(cè)識(shí)別，很快就代替了原有的用脈沖計(jì)數(shù)方式的撥號(hào)制式。這種雙音多頻信號(hào)制式不僅用在電話網(wǎng)絡(luò)中，還可以用于傳輸十進(jìn)制數(shù)據(jù)的其它通信系統(tǒng)中，用于電子郵件和銀行系統(tǒng)中。這些系統(tǒng)中用戶可以用電話發(fā)送DTMF信號(hào)選擇語(yǔ)音菜單進(jìn)行操作。DTMF信號(hào)是一個(gè)典型的小型信號(hào)處理系統(tǒng)，它要用數(shù)字方法產(chǎn)生模擬信號(hào)并進(jìn)行傳輸，用到了D/A變換器。為了系統(tǒng)的檢測(cè)速度并降低成本，還開發(fā)一種特殊的DFT算法，稱為戈澤爾(Goertzel)算法，這種算法既可以用硬件（專用芯片）實(shí)現(xiàn)，也可以用軟件實(shí)現(xiàn)。下面首先介紹雙音多頻信號(hào)的產(chǎn)生方法和檢測(cè)方法，包括戈澤爾算法，最后進(jìn)行模擬實(shí)驗(yàn)。下面先介紹電話中的DTMF信號(hào)的

8、組成。3.5.1電話中的DTMF信號(hào)的組成過(guò)去的電話撥號(hào)時(shí)靠脈沖計(jì)數(shù)確定0-9這10個(gè)數(shù)字的，不僅撥號(hào)速度慢，也不能擴(kuò)展電話上網(wǎng)其他服務(wù)功能?，F(xiàn)在均采用雙音撥號(hào)。在電話中，數(shù)字0-9中每一個(gè)都用兩個(gè)不同的單音頻傳輸，所用的8個(gè)頻率分成高頻帶和低頻帶兩組，低頻帶有四個(gè)頻率：679Hz,770Hz,852Hz和941Hz；高頻帶也有四個(gè)頻率：1209Hz,1336Hz,1477Hz和1633Hz.。每一個(gè)數(shù)字均由高、低頻帶中各一個(gè)頻率構(gòu)成，例如1用697Hz和1209Hz兩個(gè)頻率，信號(hào)用表示，其中，。這樣8個(gè)頻率形成16種不同的雙頻信號(hào)。具體號(hào)碼以及符號(hào)對(duì)應(yīng)的頻率如表1所示。表中最后一列在電話中暫

9、時(shí)未用。表1雙頻撥號(hào)的頻率分配列行1209Hz1336Hz1477Hz1633Hz697Hz123A770Hz456B852Hz789C941Hz*0#DDTMF信號(hào)在電話中有兩種作用，一個(gè)是用撥號(hào)信號(hào)去控制交換機(jī)接通被叫的用戶電話機(jī)，另一個(gè)作用是控制電話機(jī)的各種動(dòng)作，如播放留言、語(yǔ)音信箱等。3.5.2電話中的雙音多頻（DTMF）信號(hào)的產(chǎn)生DTMF編碼器基于兩個(gè)二階數(shù)字正弦波振蕩器，一個(gè)用于產(chǎn)生行頻，一個(gè)用于產(chǎn)生列頻。典型的DTMF信號(hào)頻率范圍是7001700Hz，選取8000Hz作為采樣頻率，DTMF雙音頻信號(hào)由兩個(gè)二階數(shù)字正弦振蕩器產(chǎn)生，一個(gè)用來(lái)產(chǎn)生行音頻信號(hào)，另一個(gè)產(chǎn)生列音頻信號(hào)。在電話

10、鍵盤上的每次按鍵會(huì)產(chǎn)生式2.1表示的兩個(gè)音調(diào)之和。規(guī)定用8KHz對(duì)DTMF信號(hào)進(jìn)行采樣，采樣后得到時(shí)域離散信號(hào)為 :兩個(gè)頻率f1,f2唯一確定了被按壓的鍵。具體數(shù)字對(duì)應(yīng)的兩個(gè)頻率如上表1所示。形成上面序列的方法有兩種，即計(jì)算法和查表法。用計(jì)算法求正弦波的序列值容易，但實(shí)際中要占用一些計(jì)算時(shí)間，影響運(yùn)行速度。本實(shí)驗(yàn)采用查表法即預(yù)先將正弦波的各序列值計(jì)算出來(lái)，寄存在存儲(chǔ)器中，運(yùn)行時(shí)只要按順序和一定的速度取出便可。這種方法要占用一定的存儲(chǔ)空間，但是速度快。因?yàn)椴蓸宇l率是8000Hz，因此要求每125ms輸出一個(gè)樣本，得到的序列再送到D/A變換器和平滑濾波器，輸出便是連續(xù)時(shí)間的DTMF信號(hào)。DTMF信

11、號(hào)通過(guò)電話線路送到交換機(jī)。3.5.3 Goerztel函數(shù)與DTMF信號(hào)的產(chǎn)生下面先介紹MATLAB工具箱函數(shù)goertzel，然后介紹DTMF信號(hào)的產(chǎn)生與識(shí)別仿真實(shí)驗(yàn)程序。Goerztel函數(shù)的調(diào)用格式為：Xgk=goertzel(xn,K)xn是被變換的時(shí)域序列，用于DTMF信號(hào)檢測(cè)時(shí)，xn就是DTMF信號(hào)的205個(gè)采樣值。K是要求計(jì)算的DFTxn的頻點(diǎn)序號(hào)向量，用N表示xn的長(zhǎng)度，則要求1KN。由表2可知，如果只計(jì)算DTMF信號(hào)8個(gè)基頻時(shí),K=18，20，22，24，31，34，38，42，如果同時(shí)計(jì)算8個(gè)基頻及其二次諧波時(shí)，K=18，20，22，24，31，34，35，38，39，42

12、，43，47，61，67，74，82。Xgk是變換結(jié)果向量，其中存放的是由K指定的頻率點(diǎn)的DFTx(n)的值。設(shè)X(k)= DFTx(n)，則 DTMF信號(hào)的產(chǎn)生與識(shí)別仿真實(shí)驗(yàn)在MATLAB環(huán)境下進(jìn)行，編寫仿真程序，運(yùn)行程序，送入電話號(hào)碼，程序自動(dòng)產(chǎn)生每一位號(hào)碼數(shù)字相應(yīng)的DTMF信號(hào)，并送出雙頻聲音，撥打電話號(hào)碼。對(duì)不同號(hào)碼發(fā)生模擬響應(yīng)。四、設(shè)計(jì)過(guò)程1. 設(shè)置參數(shù)并讀入電話號(hào)碼。設(shè)置每位號(hào)碼所對(duì)應(yīng)的高頻分量和低頻分量。2. 根據(jù)鍵入的電話號(hào)碼產(chǎn)生相應(yīng)的時(shí)域離散DTMF信號(hào)，添加撥號(hào)程序（Call.m），連續(xù)發(fā)出電話號(hào)碼對(duì)應(yīng)的雙頻聲音，并模擬撥號(hào)，程序根據(jù)所撥號(hào)碼響應(yīng)。3. 鍵入按鍵時(shí)對(duì)時(shí)域離散

13、DTMF信號(hào)進(jìn)行檢測(cè)，畫出幅度譜發(fā)生響應(yīng)。4. 設(shè)計(jì)GUI界面將程序可視化。用Goertzel算法計(jì)算八點(diǎn)DFT樣本構(gòu)成雙音頻信號(hào)檢測(cè)與m位號(hào)碼相符的高頻帶檢測(cè)與m位號(hào)碼相符的低頻帶輸入電話號(hào)碼開 始5. 界面優(yōu)化，素材整合。畫出DFT樣本幅度撥打電話號(hào)碼響應(yīng)不同提示音結(jié) 束圖4.1設(shè)計(jì)流程圖五、調(diào)試與結(jié)果5.1運(yùn)行結(jié)果如下圖所示圖5.1.1 撥果 圖5.1.2撥打10000結(jié)果鍵入電話號(hào)碼: 10000，接收端接檢測(cè)的號(hào)碼為：10000，同時(shí)播放語(yǔ)音，對(duì)時(shí)域離散DTMF信號(hào)進(jìn)行數(shù)據(jù)檢測(cè)，撥其他號(hào)碼時(shí)語(yǔ)音提示所撥打號(hào)碼為空號(hào)。5.2調(diào)試中遇到的問(wèn)題及解決方法（1）音樂(lè)

14、播放函數(shù)報(bào)錯(cuò)，經(jīng)過(guò)查找后改用audiopiayer函數(shù)實(shí)現(xiàn)功能。（2）使用matlab全局變量時(shí)遇到問(wèn)題，通過(guò)使用global定義變量,并在每個(gè)使用到的地方先聲明全局變量來(lái)解決問(wèn)題。 （3）對(duì)gui界面使用不熟悉，通過(guò)網(wǎng)上教學(xué)視頻學(xué)習(xí)gui的界面設(shè)計(jì)。六、分析和總結(jié)6.1小組總結(jié)通過(guò)本次課程設(shè)計(jì)，首先使我們對(duì)于雙音多頻撥號(hào)有了更全面的認(rèn)識(shí)和理解，雙音多頻撥號(hào)是現(xiàn)在最常用的撥號(hào)系統(tǒng)，數(shù)字0-9的中每一個(gè)都用兩個(gè)不同的單音頻傳輸，所用的8個(gè)頻率分成高頻帶和低頻帶兩組，低頻帶有四個(gè)頻率：679Hz、770Hz、852Hz及941Hz；高頻帶也有四個(gè)頻率：1209Hz、1336Hz、1477Hz及16

15、33Hz。每一個(gè)數(shù)字均由高、低頻帶中各一個(gè)頻率構(gòu)成。其次，對(duì)于數(shù)字信號(hào)處理中很多理解不清晰的問(wèn)題也加深了理解，也對(duì)DFT有了更深一層的理解，并且學(xué)習(xí)到計(jì)算DFT的新算法: Goerztel函數(shù)，更加方便了DFT的計(jì)算。若有機(jī)會(huì)接觸到雙音多頻，可針對(duì)音頻信號(hào)的頻率檢測(cè)進(jìn)行深一步的研究雙音頻，如作為控制信號(hào)時(shí)，可以在有線網(wǎng)絡(luò)或無(wú)線網(wǎng)絡(luò)中通過(guò)廣播形式發(fā)出,實(shí)現(xiàn)有效的通信識(shí)別和可靠的遠(yuǎn)程控制,以保證廣播系統(tǒng)運(yùn)行的安全和穩(wěn)定等實(shí)際生活中的運(yùn)用。再次，對(duì)于MATLAB的設(shè)計(jì)也更加清晰明了，matlab用起來(lái)挺順手的，比C語(yǔ)言簡(jiǎn)單。但是用到細(xì)節(jié)處的時(shí)候卻困難重重，因?yàn)楹芏嘀R(shí)都沒(méi)有學(xué)習(xí)，就算知道函數(shù)名，也不

16、知道如何調(diào)用。通過(guò)自己查詢相關(guān)資料，克服了相關(guān)的問(wèn)題，使自己對(duì)MATLAB的認(rèn)識(shí)及應(yīng)用也進(jìn)入了一個(gè)新的層次?？偠灾?，通過(guò)本次試驗(yàn)，加深了對(duì)于Matlab的了解，同時(shí)復(fù)習(xí)了數(shù)字信號(hào)方面的相關(guān)知識(shí)，對(duì)于還不太清楚的相關(guān)知識(shí)，通過(guò)進(jìn)一步的查找得到了很好的解決。了解了雙音多頻信號(hào)的產(chǎn)生、包括對(duì)雙音多頻信號(hào)進(jìn)行DFT時(shí)的參數(shù)選擇等的相關(guān)知識(shí)。在做MATLAB界面，做好以后同時(shí)會(huì)產(chǎn)一個(gè)M文件，只要把代碼編寫在對(duì)應(yīng)的控件中即可以方便的實(shí)現(xiàn)仿真。實(shí)驗(yàn)結(jié)果表明基于MATLAB的DTMF編解碼具有速度快、精度高、穩(wěn)定性好以及便于靈活應(yīng)用等特點(diǎn)。而利用MATLAB語(yǔ)言可以較好地、快速地進(jìn)行算法的仿真。6.2成員總

17、結(jié)何杰坤：通過(guò)這次課程設(shè)計(jì)，對(duì)用GUI做界面有了更深的認(rèn)識(shí)，也對(duì)用MATLAB做音頻處理有了一些基礎(chǔ)的認(rèn)識(shí)；讓我學(xué)到了很多知識(shí)，也學(xué)習(xí)到了很多和matlab相關(guān)的知識(shí)。雖然只是模擬通話，但是我還是希望在硬件實(shí)現(xiàn)上有所提高。蔡小紅：課程設(shè)計(jì)使我熟悉基本的Matlab仿真的使用，還意識(shí)到英文水平的重要性，我們學(xué)習(xí)通信專業(yè)的，還需要了解更多的專業(yè)英文，此外這次課程設(shè)計(jì)還加深了對(duì)于數(shù)字信號(hào)處理原理知識(shí)的理解。王曉清：通過(guò)本次課程設(shè)計(jì)，對(duì)MATLAB軟件有了更深刻的了解和應(yīng)用，把動(dòng)手設(shè)計(jì)和課本知識(shí)緊密的聯(lián)系起來(lái)。明白團(tuán)隊(duì)意識(shí)的重要性，最后讓自己對(duì)設(shè)計(jì)報(bào)告的編寫有了新的認(rèn)識(shí)。蘇月琴：學(xué)習(xí)到關(guān)于matlab

18、的一些新知識(shí)。對(duì)團(tuán)隊(duì)力量有了新的認(rèn)識(shí)。意識(shí)到團(tuán)隊(duì)的重要，和前期策劃的重要。認(rèn)識(shí)到自己還需加強(qiáng)團(tuán)隊(duì)合作能力。羅 燁：在這次課設(shè)中，我不僅學(xué)到了很多編程知識(shí)，更學(xué)到了遇到問(wèn)題時(shí)解決的方法。并且意識(shí)到在一個(gè)項(xiàng)目里團(tuán)隊(duì)合作比任何都來(lái)得重要。王垚天：這次課程設(shè)計(jì)最大的感想是不管學(xué)什么東西，都要自己真真實(shí)實(shí)的動(dòng)手去做。因?yàn)槲易约罕緛?lái)偏向于網(wǎng)絡(luò)硬件方面的，對(duì)軟件一直不是很擅長(zhǎng)，但是在課程設(shè)計(jì)中，自己真的去思考一個(gè)功能的實(shí)現(xiàn)，一個(gè)代碼的編寫時(shí)，才突然發(fā)現(xiàn)其實(shí)并沒(méi)有想象中那么困難。七、參考文獻(xiàn)1.數(shù)字信號(hào)處理課程設(shè)計(jì)指導(dǎo)書2數(shù)字信號(hào)處理教程，程佩青，清華大學(xué)出版社，20073MATLAB數(shù)字信號(hào)處理與應(yīng)用，李正

19、周，北京清華大學(xué)出版社，20084.雙音多頻信號(hào)產(chǎn)生及解碼研究趙霞，微電子學(xué)，2001年，第6期5.數(shù)字信號(hào)處理及MATLAB實(shí)現(xiàn)，余成波、楊菁、楊如民、周登義，清華大學(xué)出版社，2005 八、附錄：主要程序main.mfunction varargout = main(varargin)gui_Singleton = 1;gui_State = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, main_OpeningFcn, . gui_OutputFcn, main_OutputFcn

20、, . 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:);endfunction main_OpeningFcn(hObject, eventdata, handles, varargin)handles.out

21、put = hObject;guidata(hObject, handles);function varargout = main_OutputFcn(hObject, eventdata, handles) varargout1 = handles.output;function pushbutton1_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,1);set(handles.text1,string,num_string);if bool=0;

22、 Music(1);else Signal(1);end function pushbutton2_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,2);set(handles.text1,string,num_string);if bool=0; Music(2);else Signal(2);end function pushbutton3_Callback(hObject, eventdata, handles)global num_string

23、;global bool;num_string=strcat(num_string,3);set(handles.text1,string,num_string);if bool=0; Music(3);else Signal(3);endfunction pushbutton4_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,4);set(handles.text1,string,num_string);if bool=0; Music(4);els

24、e Signal(4);end function pushbutton5_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,5);set(handles.text1,string,num_string);if bool=0; Music(5);else Signal(5);endfunction pushbutton6_Callback(hObject, eventdata, handles)global num_string;global bool;n

25、um_string=strcat(num_string,6);set(handles.text1,string,num_string);if bool=0; Music(6);else Signal(6);endfunction pushbutton7_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,7);set(handles.text1,string,num_string);if bool=0; Music(7);else Signal(7);en

26、d function pushbutton8_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,8);set(handles.text1,string,num_string);if bool=0; Music(8);else Signal(8);endfunction pushbutton9_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strc

27、at(num_string,9);set(handles.text1,string,num_string);if bool=0; Music(9);else Signal(9);end function pushbutton10_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,*);set(handles.text1,string,num_string);if bool=0; Music(42);else Signal(42);end function

28、 pushbutton11_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_string,0);set(handles.text1,string,num_string);if bool=0; Music(0);else Signal(0);end function pushbutton12_Callback(hObject, eventdata, handles)global num_string;global bool;num_string=strcat(num_

29、string,#);set(handles.text1,string,num_string);if bool=0; Music(35);else Signal(35);endfunction pushbuttonCall_Callback(hObject, eventdata, handles)global bool;global num2_string;global num_string;num2_string=num_string;num_string=;Call();if str2num(num2_string)=10000; bool=1; Music(10000);else Musi

30、c(11);endset(handles.text2,string,num2_string);set(handles.text1,string,);function pushbuttonDown_Callback(hObject, eventdata, handles)global num_string;global music;global bool;bool=0;num_string=;set(handles.text1,string,);set(handles.text2,string,);music=0;function text1_CreateFcn(hObject, eventda

31、ta, handles)function text2_CreateFcn(hObject, eventdata, handles)signal.mfunction y=Signal(d);tm=1,2,3,65; 4,5,6,66; 7,8,9,67; 42,48,35,68; N=205; Fs=8000; K=18,20,22,24,31,34,38,42; f1=697,770,852,941; f2=1209,1336,1477,1633; for p=1:4; for q=1:4; if tm(p,q)=abs(d); break;end end if tm(p,q)=abs(d); break;end endn=0:1023; y = sin(2*pi*n*f1(p)/Fs) + sin