基于DSP的DTMF信令的產(chǎn)生分析與檢測

　　編程的流程如圖3所示，由CCITT的規(guī)定，數(shù)字之間必須有適當(dāng)長度的靜音，因此編碼器有兩個任務(wù)，其一是音頻信號任務(wù)，產(chǎn)生雙音樣本，其二是靜音任務(wù)，產(chǎn)生靜音樣本。每個任務(wù)結(jié)束后，啟動下一個任務(wù)前(音頻信號任務(wù)或靜音任務(wù))，都必須復(fù)位決定其持續(xù)時間的定時器變量。在靜音任務(wù)結(jié)束后，DSP從數(shù)字緩存中調(diào)出下一個數(shù)字,判決該數(shù)字信號所對應(yīng)的行頻和列頻信號，并根據(jù)不同頻率確定其初始化參數(shù)a1=-2cosω0與y(-2)=-Asinω0。

　　該流程圖可采用C語言實現(xiàn)，雙音信號的產(chǎn)生則由54x匯編代碼實現(xiàn)。整個程序作為C54x的多通道緩沖串口(McBsp)的發(fā)射串口中斷服務(wù)子程序，由外部送入的8000Hz串口時鐘觸發(fā)中斷，可實時處理并通過D/A轉(zhuǎn)換器輸出DTMF信令信號。

　　2 DTMF信號的檢測

　　在輸入信號中檢測DTMF信號，并將其轉(zhuǎn)換為實際的數(shù)字，這一解碼過程本質(zhì)是連續(xù)的過程，需要在輸入的數(shù)據(jù)信號流中連續(xù)地搜索DTMF信號頻譜的存在。整個檢測過程分兩步：首先采用Goertzel算法在輸入信號中提取頻譜信息;接著作檢測結(jié)果的有效性檢查。

　　2.1 Goertzel算法

　　DTMF解碼即是在輸入信號中搜索出有效的行頻和列頻。計算數(shù)字信號的頻譜可以采用DFT及其快速算法FFT，而在實現(xiàn)DTMF解碼時，采用Goertzel算法要比FFT更快。通過FFT可以計算得到信號所有譜線，了解信號整個頻域信息，而對于DTMF信號只用關(guān)心其8個行頻/列頻及其二次諧波信息即可(二次諧波的信息用于將DTMF信號與聲音信號區(qū)別開)。此時Goertzel算法能更加快速的在輸入信號中提取頻譜信息。

　　圖 4

　　Goertzel算法實質(zhì)是一個兩極點的IIR濾波器，其算法原理框圖如圖4。由于在DTMF檢測中，輸入的信號是實數(shù)序列，并不需要檢測出8個行頻/列頻的相位，只需要計算出其幅度平方即可。

　　2.2 DTMF檢測器流程

　　檢測流程可參照圖5，把檢測程序作為C54x的McBsp接收中斷服務(wù)子程序，在每一個接收中斷到來時，表明采到一個新樣點。樣點值代入式(2)，迭代計算8個行頻/列頻的中間變量vk(n)(k為8個行頻/列頻分別對應(yīng)的數(shù)字頻率)，直到采到N=125個樣點(在8kHz采樣頻率下，約為15ms)。此時再按式(4)計算8個行頻/列頻的幅度平方|X(k)|2。接下來將|X(k)|2與門限作比較，并作二次諧波檢測，判決出有效的音頻信號。將音頻信號映射為數(shù)字信號后，再與上一個檢測到的數(shù)字信號比較，最終判決出有效的數(shù)字信號。

　　圖 5

　　按圖5所示流程得到DTMF信令檢測程序。整個程序作為C54x的McBsp接收串口中斷服務(wù)子程序，從而可以實時分析來自A/D轉(zhuǎn)換器的DTMF信令信號。

　　3 性能分析

　　基于上述原理與算法代碼，在TI公司的DSP開發(fā)環(huán)境Code Composer Studio(CCS)下，分析上述整個DTMF信令的產(chǎn)生與檢測方案的性能。

　　(1)由CCS給出的如下內(nèi)存印象文件報告，DTMF的產(chǎn)生(gen_dtmf.obj)與DTMF的檢測(de_dtmf.obj)這兩段核心代碼分別占用3e6H和1e0H個字(16bit word)，即約占1K字的存儲器空間，消耗系統(tǒng)資源極低;

　　(2)DTMF信令的產(chǎn)生與檢測程序均放置于C54x的McBSP中斷服務(wù)子程序內(nèi)，由CCS的代碼剖析工具分析代碼執(zhí)行時間，當(dāng) C54x 運行在主頻100MHz時，DTMF產(chǎn)生中斷服務(wù)子程序interrupt transmit()最大消耗283個時鐘周期，即2.83μs，DTMF檢測中斷服務(wù)子程序interrupt receive()最大消耗6148個時鐘周期，約61μs。因此該方案能夠?qū)崟r產(chǎn)生與檢測DTMF信令，還可保證有時間冗余度，與其他程序在用戶系統(tǒng)中并發(fā)執(zhí)行。