【正文】 BackgroundColor39。white39。), get(0,39。,nu)。 if n=11。) fftx=abs(fft(x))。))。string39。)。string39。) set(,39。 plot(n,x)%grid。39。|X(k)|39。 axes()。,39。k39。 n=0:1023。)。 t=n/N*8000。,nu)。 if n=11。) fftx=abs(fft(x))。))。string39。)。string39。) set(,39。 plot(n,x)%grid。39。)。 x = sin(2*pi*n*f1(3)/8000) + sin(2*pi*n*f2(1)/8000)。 set(,39。) fftx=abs(fft(x))。))。string39。)。amp。相信此次課程設(shè)計讓我所體會到的人生哲理，必將成為今后生活和學(xué)習(xí)中的巨大財富，使我終生受益。現(xiàn)有工具箱有：符號運算（利用 Maple V的計算核心執(zhí)行 ）、影像處理、統(tǒng)計分析、訊號處理、神經(jīng)網(wǎng)路、模擬分析、控制系統(tǒng)、即時控制、系統(tǒng)確認、強建控 制、弧線分析、最佳化、模糊邏輯、 mu 分析及合成、化學(xué)計量分析。 (1) nkNW 的 共軛 對稱性： ??? ?? )()( knNknNnNKN WWW (2) nkNW的周期性： )()( NknNkNnNnkN WWW ?? ?? (3) nkNW的可約性： mnk mNmn kmNnkN WWW //?? 由此得出 nkNknNNkNnN WWW ??? ?? )()( ， 12 ?NNW ， kNNkN WW ??? )2( 這樣，（ 1）利用這些特性，使 DFT 運算中有些項可以合并；（ 2）利用 nkNW的周期性和對稱性，可以將長序列 的 DFT 分解為短序列的 DFT。 檢測方法有兩種，一種是 用一組濾波器來提取所需頻率 ，判斷對應(yīng)的數(shù)字或符號；另一種是用 快速傅立葉變換（ FFT）算法的DFT 對雙音多頻信號進行頻譜分析 ， 有信號的幅度譜，判斷信號的兩個頻率，最后確定對應(yīng)的數(shù)字或符號 。每一位號碼由一個高頻信號和一個低頻信號疊加形成， DTMF 信號有 16 個編碼。它們的系統(tǒng)函數(shù)分別為： 燕山大學(xué)課程設(shè)計說明書 第 頁 共 32 頁 7 ???????? NkkkMrrzabzZH101)( () ????? 10)()( NnnznhZH () 式（ ）中的 H(Z)稱為 N 階 IIR 濾波器函數(shù)式 ， （ ）中的 H(Z)稱為 N1 階 FIR 濾波器函數(shù)。而數(shù)字信號在時間上和幅值上都有離散的數(shù)字。早在 1878 年就出現(xiàn)了人工供電制交換機，它借助話務(wù)員進行話務(wù)連接。 2020年 6 月 20 日 燕山大學(xué) 課程設(shè)計評審意見表 指導(dǎo)教師評語： 工 作 態(tài) 度認 真較 認 真不 認 真理 論 分 析正 確 完 善較 為 合 理一 般較 差方 法 設(shè) 計完 善合 理一 般較 差 成績： 指導(dǎo)教師： 年 月 日 答辯小組評語： 工 作 態(tài) 度認 真較 認 真不 認 真理 論 分 析正 確 完 善較 為 合 理一 般較 差方 法 設(shè) 計完 善合 理一 般較 差 成績 ： 組長： 年 月 日 課程設(shè)計總成績： 答辯小組成員簽字： 年 月 日 燕山大學(xué)課程設(shè)計說明書 第 頁 共 32 頁 2 摘 要 雙音多頻（ Dual Tone Multi Frequency, DTMF）信號是音頻電話中的撥號信號，由美國 ATamp。用戶的撥號脈沖不再直接控制接線器動作，而先由記發(fā)器接收、存儲，然后通過標(biāo)志器驅(qū)動接線 器，以完成用戶間接續(xù)。 （ 2）便于存儲、處理和交換 數(shù)字通信的信號形式和 計算機 所用信號一致，都是 二進制代碼 ，因此便于與計算機聯(lián)網(wǎng)，也便于用計算機對數(shù)字信號進行存儲、處理和交換，可使通信網(wǎng)的管理、維護實現(xiàn)自 動化、智能化。同時 Z 變換對進行數(shù)字系統(tǒng)分析和實現(xiàn)也極其重要。 在電話鍵盤上的每次按鍵，會產(chǎn)生 式（ ）表示的兩個音調(diào)之和。 取樣頻率 sf 與取樣周期 T 的關(guān)系是 Tfs 1? 取樣角頻率 T?2?? T 為時域取樣間隔 ，在一個周期內(nèi)取樣點數(shù)為 N。 有一種 DMFT 的特殊算法－ Goertzel 算法 ，更適合 適合雙音多頻信號頻譜的計算 。 MATLAB 中 擁有 GUI 界面的功能，可以將想要實現(xiàn)的效果通過 GUI界面 展現(xiàn)出來，在 MATLAB 主界面，通過路徑 FileNewGUI 打開 GUI界面的編輯窗口，將需要的按鈕拖出，雙擊可進行編輯。gui_Singleton39。 guidata(hObject, handles)。 N=length(x)。數(shù)字框已滿！ 39。 f2=[1209,1336,1477,1633]。string39。 else nu=strcat(get(,39。k39。,7) num=get(,39。string39。|X(k)|39。 n=length(num)。string39。 N=length(x)。數(shù)字框已滿！ 39。 f2=[1209,1336,1477,1633]。xlabel(39。string39。 sound(x,8000)。ylabel(39。)。xlabel(39。string39。),get(,39。ylabel(39。) num=get(,39。string39。|X(k)|39。 n=length(num)。string39。 N=length(x)。數(shù)字框已滿！ 39。 f2=[1209,1336,1477,1633]。xlabel(39。string39。), get(0,39。white39。BackgroundColor39。,39。white39。), get(0,39。)。),get(,39。string39。xlabel(39。)。ylabel(39。 sound(x,8000)。string39。xlabel(39。 f2=[1209,1336,1477,1633]。數(shù)字框已滿！ 39。 plot(t,fftx)%grid。 end function pushbutton9_Callback(hObject, eventdata, handles) %按鍵 *對應(yīng)功能 f1=[697,770,852,941]。 errordlg(39。 axes()。 set(,39。)。ylabel(39。),get(,39。string39。xlabel(39。)。ylabel(39。 sound(x,8000)。string39。xlabel(39。 f2=[1209,1336,1477,1633]。數(shù)字框已滿！ 39。 N=length(x)。string39。 n=length(num)。xlabel(39。)。ylabel(39。 sound(x,8000)。gui_LayoutF39。 試驗結(jié)果表明基于 MATLAB 的 DTMF 編解碼具有速度快、精度高、穩(wěn)定性好以及便于靈活應(yīng)用等特點。 燕山大學(xué)課程設(shè)計說明書 第 頁 共 32 頁 16 它有如下的特點： （ 1） 功能強的數(shù)值運算 ― 在 MATLAB 環(huán)境中，有超過 500 種數(shù)學(xué)、統(tǒng)計、科學(xué)及工程方面的函數(shù)可使用，函數(shù)的標(biāo)示自然，使得問 題和解答像數(shù)學(xué)式子一般簡單明了，讓使用者可全力發(fā)揮在解題方面，而非浪費在電腦操作上。復(fù)數(shù)運算實際上是由實數(shù)運算來完成的，因此上式可寫成： ???? 10 )()( Nn nkNWnxkX }][]Re [}{)]([)](Re [{10 nkNNn nkN WmjWnxmjnx ????? ? ??]}R e [)]([][)](( R e [])[)]([]R e [)]({ ( R e [10nkNnkNNnnkNnkNWnxmWmnxjWmnxmWnx???????? ??? () 由此可見，一次復(fù)數(shù)乘法 需用四次實數(shù)乘法和兩次實數(shù)加法；一次復(fù)數(shù)加法則需兩次復(fù)數(shù)加法。這種方法要占用一定的存儲空間，但是速度快。 燕山大學(xué)課程設(shè)計說明書 第 頁 共 32 頁 9 離散系統(tǒng)傅立葉變換（ DFT） 離散傅立葉變換是針對有限長序列或周期序列從存在的，相當(dāng)把序列的連續(xù)傅立葉變換加以抽樣，頻率的離散化造成時間函數(shù)也呈周期，故級數(shù)應(yīng)限制在一個周期內(nèi)。隨著寬頻帶信道 (光纜 、數(shù)字微波 )的大量利用以及 數(shù)字信號處理 技術(shù)的發(fā)展，數(shù)字電話的帶寬問題已不是主要問題了。信息可以通過電信號進行傳送、交換、存儲、提取等。 本文主要介紹 用 MATLAB 在計算機上對雙音多頻的信令產(chǎn)生及檢測接收系統(tǒng)進行了仿真測試。 DTMF 信號系統(tǒng)是一個典型的小型信號處理系統(tǒng)，它要用數(shù)字方法產(chǎn)生模擬信號并進行傳輸，其中還用到了 D/A 變換器；在接收端用 A/D 變換器將其轉(zhuǎn)換成數(shù)字信號，并進行數(shù)字信號處理。電信號是指隨著時間而變化的電壓或電流 ，因此在數(shù)學(xué)描述上可將它表示為時間的函數(shù)，并可畫出其波形。 （ 5）占用信道頻帶較寬 一路模擬電話的 頻帶 為 4kHz 帶寬，一路數(shù)字電話約占 64kHz，這是模擬通信目前仍有生命力的主要原因。對 N 階差分方程進行 Z 變換，得到系統(tǒng)函數(shù)的一般表示式 ???????? NiiiMiiizazbNzX zYzH00)()()( () 令 Z=ejω即得系統(tǒng)的傳輸函數(shù) H(ejω)。查表法是預(yù)先將正弦波的各序列值計算出來，寄存在存儲器中，運行時只要按順序和一定的速度取出便可。而 X（ k） 一共有 N 個點（ k 從 0 取到 N－ 1），所以完成整個 DFT 運算總共需要 2N 次復(fù)數(shù)乘法及 N（ N－ 1） 次復(fù)數(shù)加法。附加的工具箱（單獨提供的專用 MATLAB 函數(shù)集）擴展了 MATLAB 環(huán)境，以解決這些應(yīng)用領(lǐng)域內(nèi)特定類型的問題。 在 做 MATLAB界面，做好以后同時會產(chǎn)一個 M文件，只要把代碼編寫在對應(yīng)的控件中即可以方便的實現(xiàn)仿真。, untitled_OutputF, ... 39。 x = sin(2*pi*n*f1(2)/8000) + sin(2*pi*n*f2(1)/8000)。)。39。 plot(n,x)%grid。)。 set(,39