數(shù)字信號處理課程設計

1、<p>　　《數(shù)字信號處理》課程設計指導書</p><p><b>　　常波 編</b></p><p><b>　　淮陰工學院</b></p><p>　　計算機工程學院通信系</p><p><b>　　2009.10</b></p><p&g

2、t;<b>　　一、課程設計目的</b></p><p>　　《數(shù)字信號處理》課程是一門理論和技術發(fā)展十分迅速、應用非常廣泛的前沿性學科， 它的理論性和實踐性都很強，學生在學習這門課程時， 普遍感到數(shù)字信號處理的概念抽象， 對其中的分析方法與基本理論不能很好地理解與掌握。因此， 如何幫助學生理解與掌握課程中的基本概念、基本原理、基本分析方法以及綜合應用所學知識解決實際問題的能力， 是本課程教

3、學中所要解決的關鍵問題。</p><p>　　為了鞏固所學的數(shù)字信號處理理論知識， 使學生對信號的采集、處理、傳輸、顯示和存儲等有一個系統(tǒng)的掌握和理解， 安排了該課程的課程設計。 通過課程設計，應能加強學生如下能力的培養(yǎng)：</p><p>　　獨立工作能力和創(chuàng)造力；</p><p>　　綜合運用專業(yè)及基礎知識，解決實際工程技術問題的能力；</p>&l

4、t;p>　　查閱圖書資料、產(chǎn)品手冊和各種工具書的能力；</p><p>　　編寫技術報告和編制技術資料的能力。</p><p><b>　　二、一般設計要求</b></p><p><b>　　獨立完成設計任務。</b></p><p>　　按照設計任務寫出設計工作小結。對設計過程所進行的有關

5、步驟作出說明，并對完成的設計作出評價，總結自己整個設計工作中的經(jīng)驗教訓，收獲和今后研究方向。</p><p>　　編寫課程設計說明書，說明書必須按統(tǒng)一格式打印，裝訂成冊，字數(shù)一般不少于三千字。</p><p><b>　　三、設計課題</b></p><p>　　題目一 基于MATLAB 的語音處理系統(tǒng)設計</p><p&

6、gt;<b>　?。ㄒ唬﹥?nèi)容</b></p><p>　　錄制一段個人自己的語音信號，并對錄制的信號進行采樣；畫出采樣后語音信號的時域波形和頻譜圖；給定濾波器的性能指標， 采用窗函數(shù)法和雙線性變換設計濾波器， 并畫出濾波器的頻率響應；然后用自己設計的濾波器對采集的信號進行濾波，畫出濾波后信號的時域波形和頻譜， 并對濾波前后的信號進行對比，分析信號的變化；回放語音信號；最后， 設計一個信號處理

7、系統(tǒng)界面。</p><p><b>　　（二）步驟</b></p><p><b>　　1．語音信號的采集</b></p><p>　　利用Windows 下的錄音機， 錄制一段自己的話音， 時間在1 s 內(nèi)。然后在Matlab 軟件平臺下， 利用函數(shù)wavread 對語音信號進行采樣， 記住采樣頻率和采樣點數(shù)。</

8、p><p>　　2．語音信號的頻譜分析</p><p>　　利用函數(shù)fft 對信號進行快速傅里葉變換， 得到信號的頻譜特性， 要求首先畫出語音信號的時域波形；然后對語音信號進行頻譜分析。</p><p>　　3．設計數(shù)字濾波器和畫出其頻率響應</p><p>　　給出各濾波器的性能指標:</p><p>　　(1) 低通濾

9、波器性能指標　fb= 1000Hz， fc=1200 Hz， As=100dB， Ap= 1 dB。</p><p>　　(2) 高通濾波器性能指標　fc= 4800 Hz， fb=5 000 Hz As= 100 dB， Ap= 1 dB。</p><p>　　(3) 帶通濾波器性能指標　fb1=1200 Hz， fb2=3 000 Hz， fc1= 1 000 Hz， fc2= 3 2

10、00 Hz， As=100 dB，Ap= 1dB。</p><p>　　要求用窗函數(shù)法和雙線性變換法設計上面要求的3種濾波器。在Matlab中，可以利用函數(shù)fir1設計FIR濾波器，可以利用函數(shù)butte， cheby1和ellip設計IIR濾波器；利用Matlab中的函數(shù)freqz畫出各濾波器的頻率響應。</p><p>　　4．用濾波器對信號進行濾波</p><p&

11、gt;　　要求用自己設計的各濾波器分別對采集的信號進行濾波， 在M at lab 中， FIR 濾波器利用函數(shù)fftfilt 對信號進行濾波。</p><p>　　5．比較濾波前后語音信號的波形及頻譜</p><p>　　要求在一個窗口同時畫出濾波前后的波形及頻譜。</p><p><b>　　6．回放語音信號</b></p>&

12、lt;p>　　在Matlab 中， 函數(shù)sound 可以對聲音進行回放。其調(diào)用格式: sound (x， fs， bit s) ； 可以感覺濾波前后的聲音有變化。</p><p><b>　　7．設計系統(tǒng)界面</b></p><p>　　為了使編制的程序操作方便， 要求有能力的學生， 設計處理系統(tǒng)的用戶界面。在所設計的系統(tǒng)界面上可以選擇濾波器的類型， 輸入濾波器

13、的參數(shù)， 顯示濾波器的頻率響應， 選擇信號等。</p><p><b>　?。ㄈ┱f明書要求</b></p><p>　　1．研究語音信號的特點；</p><p>　　2．探討語音分析與識別、處理的基本理論基礎、關鍵技術；</p><p>　　3．寫出各步驟的MATLAB的程序代碼；</p><p&g

14、t;　　4．分析采集的語音信號的時域波形與頻譜；分析濾波前后語音信號的波形與頻譜；</p><p>　　5．濾波器的設計可采用圖形化設計工具FDATool和圖形化信號處理設計工具SPTool；</p><p>　　題目二 基于MATLAB 有噪聲語音信號的處理</p><p><b>　　（一）內(nèi)容</b></p><p&

15、gt;　　選擇一個語音信號作為分析的對象，或錄制一段各人自己的語音信號，對其進行頻譜分析；利用MATLAB中的隨機函數(shù)產(chǎn)生噪聲加入到語音信號中，模仿語音信號被污染，并對其進行頻譜分析；設計F IR和IIR數(shù)字濾波器， 并對被噪聲污染的語音信號進行濾波， 分析濾波后信號的時域和頻域特征，回放語音信號。最后，設計一個信號處理系統(tǒng)界面。</p><p><b>　?。ǘ┎襟E與要求</b><

16、/p><p><b>　　1．語音信號的采集</b></p><p>　　利用Windows 下的錄音機， 錄制一段自己的話音， 時間在1 s 內(nèi)。然后在Matlab 軟件平臺下， 利用函數(shù)wavread 對語音信號進行采樣， 記住采樣頻率和采樣點數(shù)。</p><p>　　2．語音信號的頻譜分析</p><p>　　在Mat

17、lab中，可以利用函數(shù)fft 對信號進行快速傅里葉變換， 得到信號的頻譜特性， 要求學生首先畫出語音信號的時域波形； 然后對語音信號進行頻譜分析。</p><p>　　3．利用MATLAB中的隨機函數(shù)產(chǎn)生噪聲加入到語音信號中，模仿語音信號被污染，并對其進行頻譜分析。</p><p>　　4．設計數(shù)字濾波器和畫出其頻率響應給出濾波器的性能指標。</p><p>　　確

18、定濾波器技術指標，分別設計FIR數(shù)字濾波器和IIR數(shù)字濾波器。在Matlab中，可以利用函數(shù)fir1設計FIR濾波器，可以利用函數(shù)butte， cheby1和ellip設計IIR濾波器；利用Matlab中的函數(shù)freqz畫出各濾波器的頻率響應。</p><p>　　5．用濾波器對信號進行濾波</p><p>　　要求用自己設計的各濾波器分別對加噪聲后的信號進行濾波，在Matlab中，F(xiàn)IR

19、濾波器利用函數(shù)fftfilt 對信號進行濾波，IIR濾波器利用函數(shù)filter對信號進行濾波。</p><p>　　5．比較濾波前后語音信號的波形及頻譜</p><p>　　要求在一個窗口同時畫出濾波前后的波形及頻譜。</p><p><b>　　6．回放語音信號</b></p><p>　　在Matlab 中，函數(shù)so

20、und 可以對聲音進行回放。其調(diào)用格式: sound (x，fs， bits) ； 可以感覺濾波前后的聲音有變化。</p><p><b>　　7．設計系統(tǒng)界面</b></p><p>　　為了使編制的程序操作方便， 要求有能力的學生設計處理系統(tǒng)的用戶界面。在所設計的系統(tǒng)界面上可以選擇濾波器的類型， 輸入濾波器的參數(shù)， 顯示濾波器的頻率響應，選擇信號等。</p&

21、gt;<p><b>　?。ㄈ┱f明書要求</b></p><p>　　1．研究語音信號的特點；</p><p>　　2．探討語音分析與識別、處理的基本理論基礎、關鍵技術；</p><p>　　3．寫出各步驟的MATLAB的程序代碼；</p><p>　　4．分析采集的語音信號的時域波形與頻譜；分析濾波前后

22、語音信號的波形與頻譜；</p><p>　　5．濾波器的設計可采用圖形化設計工具FDATool和圖形化信號處理設計工具SPTool；</p><p>　　題目三 基于MATLAB 的語音信號的特技處理</p><p><b>　　（一）內(nèi)容</b></p><p>　　選擇一個語音信號作為分析的對象，或錄制一段各人自己

23、的語音信號，并對其進行頻譜分析； 然后在時域用數(shù)字信號處理方法將信號加入延時和混響， 再分析其頻譜，并與原始信號頻譜進行比較；最后設計一個信號處理系統(tǒng)界面。</p><p><b>　?。ǘ┎襟E</b></p><p>　　1．利用Windows下的錄音機或其他軟件，錄制一段自己的語音信號，時間控制在1s左右，并對錄制的信號進行采樣；</p><

24、p>　　2．語音信號的頻譜分析，畫出采樣后語音信號的時域波形和頻譜圖；</p><p>　　3．將信號加入延時和混響，再分析其頻譜，并與原始信號頻譜進行比較；</p><p>　　4．設計幾種特殊類型的濾波器：單回聲濾波器，多重回聲濾波器，無限個回聲濾波器，全通結構的混響器，并畫出濾波器的頻域響應；</p><p>　　5．用自己設計的濾波器對采集的語音信號進

25、行濾波；</p><p>　　6．分析得到信號的頻譜，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進行對比，分析信號的變化；</p><p><b>　　7．回放語音信號。</b></p><p><b>　?。ㄈ┱f明書要求</b></p><p>　　1．研究語音信號的特點；</p&

26、gt;<p>　　2．探討語音分析與識別、處理的基本理論基礎、關鍵技術；</p><p>　　3．寫出各步驟的MATLAB的程序代碼；</p><p>　　4．分析采集的語音信號的時域波形與頻譜；分析濾波前后語音信號的波形與頻譜；</p><p>　　5．濾波器的設計可采用圖形化設計工具FDATool和圖形化信號處理設計工具SPTool；</p&

27、gt;<p>　　題目四 基于MATLAB的雙音多頻撥號系統(tǒng)的仿真</p><p><b>　　（一）基本原理</b></p><p><b>　　1、引言</b></p><p>　　雙音多頻（Dual Tone Multi Frequency， DTMF）信號是音頻電話中的撥號信號，由美國AT&

28、T貝爾公司實驗室研制，并用于電話網(wǎng)絡中。這種信號制式具有很高的撥號速度，且容易自動監(jiān)測識別，很快就代替了原有的用脈沖計數(shù)方式的撥號制式。這種雙音多頻信號制式不僅用在電話網(wǎng)絡中，還可以用于傳輸十進制數(shù)據(jù)的其它通信系統(tǒng)中，用于電子郵件和銀行系統(tǒng)中。這些系統(tǒng)中用戶可以用電話發(fā)送DTMF信號選擇語音菜單進行操作。</p><p>　　DTMF信號系統(tǒng)是一個典型的小型信號處理系統(tǒng)，它要用數(shù)字方法產(chǎn)生模擬信號并進行傳輸，其中

29、還用到了D/A變換器；在接收端用A/D變換器將其轉換成數(shù)字信號，并進行數(shù)字信號處理與識別。為了系統(tǒng)的檢測速度并降低成本，還開發(fā)一種特殊的DFT算法，稱為戈澤爾(Goertzel)算法，這種算法既可以用硬件（專用芯片）實現(xiàn)，也可以用軟件實現(xiàn)。下面首先介紹雙音多頻信號的產(chǎn)生方法和檢測方法，包括戈澤爾算法，最后進行模擬實驗。下面先介紹電話中的DTMF信號的組成。</p><p>　　在電話中，數(shù)字0~9的中每一個都用兩

30、個不同的單音頻傳輸，所用的8個頻率分成高頻帶和低頻帶兩組，低頻帶有四個頻率：679Hz，770Hz，852Hz和941Hz；高頻帶也有四個頻率：1209Hz，1336Hz，1477Hz和1633Hz.。每一個數(shù)字均由高、低頻帶中各一個頻率構成，例如1用697Hz和1209Hz兩個頻率，信號用表示，其中，。這樣8個頻率形成16種不同的雙頻信號。具體號碼以及符號對應的頻率如表1所示。表中最后一列在電話中暫時未用。</p>&l

31、t;p>　　表1 雙頻撥號的頻率分配</p><p>　　DTMF信號在電話中有兩種作用，一個是用撥號信號去控制交換機接通被叫的用戶電話機，另一個作用是控制電話機的各種動作，如播放留言、語音信箱等。</p><p>　　2、電話中的雙音多頻（DTMF）信號的產(chǎn)生與檢測</p><p>　　（1）雙音多頻信號的產(chǎn)生</p><p>　

32、　假設時間連續(xù)的 DTMF信號用表示，式中是按照表1選擇的兩個頻率，代表低頻帶中的一個頻率，代表高頻帶中的一個頻率。顯然采用數(shù)字方法產(chǎn)生DTMF信號，方便而且體積小。下面介紹采用數(shù)字方法產(chǎn)生DTMF信號。規(guī)定用8KHz對DTMF信號進行采樣，采樣后得到時域離散信號為</p><p>　　形成上面序列的方法有兩種，即計算法和查表法。用計算法求正弦波的序列值容易，但實際中要占用一些計算時間，影響運行速度。查表法是預先

33、將正弦波的各序列值計算出來，寄存在存儲器中，運行時只要按順序和一定的速度取出便可。這種方法要占用一定的存儲空間，但是速度快。</p><p>　　因為采樣頻率是8000Hz，因此要求每125ms輸出一個樣本，得到的序列再送到D/A變換器和平滑濾波器，輸出便是連續(xù)時間的DTMF信號。DTMF信號通過電話線路送到交換機。</p><p>　?。?）雙音多頻信號的檢測</p>&l

34、t;p>　　在接收端，要對收到的雙音多頻信號進行檢測，檢測兩個正弦波的頻率是多少，以判斷所對應的十進制數(shù)字或者符號。顯然這里仍然要用數(shù)字方法進行檢測，因此要將收到的時間連續(xù) DTMF信號經(jīng)過A/D變換，變成數(shù)字信號進行檢測。檢測的方法有兩種，一種是用一組濾波器提取所關心的頻率，根據(jù)有輸出信號的2個濾波器判斷相應的數(shù)字或符號。另一種是用DFT（FFT）對雙音多頻信號進行頻譜分析，由信號的幅度譜，判斷信號的兩個頻率，最后確定相應的數(shù)字

35、或符號。當檢測的音頻數(shù)目較少時，用濾波器組實現(xiàn)更合適。FFT是DFT的快速算法，但當DFT的變換區(qū)間較小時，F(xiàn)FT快速算法的效果并不明顯，而且還要占用很多內(nèi)存，因此不如直接用DFT合適。下面介紹Goertzel算法，這種算法的實質是直接計算DFT的一種線性濾波方法。這里略去Goertzel算法的介紹（請參考文獻[19]），可以直接調(diào)用MATLAB信號處理工具箱中戈澤爾算法的函數(shù)Goertzel，計算N點DFT的幾個感興趣的頻點的值。&l

36、t;/p><p>　　3、檢測DTMF信號的DFT參數(shù)選擇</p><p>　　用DFT檢測模擬DTMF信號所含有的兩個音頻頻率，是一個用DFT對模擬信號進行頻譜分析的問題。根據(jù)第三章用DFT對模擬信號進行譜分析的理論，確定三個參數(shù)：（1）采樣頻率，（2）DFT的變換點數(shù)N，（3）需要對信號的觀察時間的長度。這三個參數(shù)不能隨意選取，要根據(jù)對信號頻譜分析的要求進行確定。這里對信號頻譜分析也有三個

37、要求： （1）頻率分辨率，（2）譜分析的頻譜范圍，（3）檢測頻率的準確性。</p><p>　?。?） 頻譜分析的分辨率。</p><p>　　觀察要檢測的8個頻率，相鄰間隔最小的是第一和第二個頻率，間隔是73Hz，要求DFT最少能夠分辨相隔73Hz的兩個頻率，即要求。DFT的分辨率和對信號的觀察時間有關， 。考慮到可靠性，留有富裕量，要求按鍵的時間大于40ms。</p>

38、<p>　　（2）頻譜分析的頻率范圍</p><p>　　要檢測的信號頻率范圍是697～1633Hz，但考慮到存在語音干擾，除了檢測這8個頻率外，還要檢測它們的二次倍頻的幅度大小，波形正常且干擾小的正弦波的二次倍頻是很小的，如果發(fā)現(xiàn)二次諧波很大，則不能確定這是DTMF信號。這樣頻譜分析的頻率范圍為697～3266Hz。按照采樣定理，最高頻率不能超過折疊頻率，即，由此要求最小的采樣頻率應為7.24KHz。

39、因為數(shù)字電話總系統(tǒng)已經(jīng)規(guī)定＝8KHz，因此對頻譜分析范圍的要求是一定滿足的。按照，＝8KHz，算出對信號最少的采樣點數(shù)為。</p><p>　?。?）檢測頻率的準確性</p><p>　　這是一個用DFT檢測正弦波頻率是否準確的問題。序列的N點DFT是對序列頻譜函數(shù)在0～區(qū)間的N點等間隔采樣，如果是一個周期序列，截取周期序列的整數(shù)倍周期，進行DFT，其采樣點剛好在周期信號的頻率上，DFT的

40、幅度最大處就是信號的準確頻率。分析這些DTMF信號，不可能經(jīng)過采樣得到周期序列，因此存在檢測頻率的準確性問題。</p><p>　　DFT的頻率采樣點頻率為（k=0，1，2，---，N-1），相應的模擬域采樣點頻率為（k=0，1，2，---，N-1），希望選擇一個合適的N，使用該公式算出的能接近要檢測的頻率，或者用8個頻率中的任一個頻率代入公式中時，得到的k值最接近整數(shù)值，這樣雖然用幅度最大點檢測的頻率有誤差，但

41、可以準確判斷所對應的DTMF頻率，即可以準確判斷所對應的數(shù)字或符號。經(jīng)過分析研究認為N＝205是最好的。按照＝8KHz，N＝205，算出8個頻率及其二次諧波對應k值，和k取整數(shù)時的頻率誤差見表2。</p><p><b>　　表2 </b></p><p>　　通過以上分析，確定＝8KHz，N＝205，。</p><p>　　4、DTMF信號

42、的產(chǎn)生與識別仿真實驗</p><p>　　下面先介紹MATLAB工具箱函數(shù)goertzel，然后介紹DTMF信號的產(chǎn)生與識別仿真實驗程序。Goerztel函數(shù)的調(diào)用格式額為</p><p>　　Xgk=goertzel(xn，K)</p><p>　　xn是被變換的時域序列，用于DTMF信號檢測時，xn就是DTMF信號的205個采樣值。</p><

43、;p>　　K是要求計算的DFT[xn]的頻點序號向量，用N表示xn的長度，則要求1≤K≤N。由表10.2.2可知，如果只計算DTMF信號8個基頻時，</p><p>　　K=[18，20，22，24，31，34，38，42]，</p><p>　　如果同時計算8個基頻及其二次諧波時，</p><p>　　K=[18，20，22，24，31，34，35，38，3

44、9，42，43，47，61，67，74，82]。</p><p>　　Xgk是變換結果向量，其中存放的是由K指定的頻率點的DFT[x(n)]的值。設X(k)= DFT[x(n)]，則。 </p><p>　　DTMF信號的產(chǎn)生與識別仿真實驗在MATLAB環(huán)境下進行，編寫仿真程序，運行程序，送入6位電話號碼，程序自動產(chǎn)生每一位號碼數(shù)字相應的DTMF信號，并送出雙頻聲音，再用DFT進行譜分析，

45、顯示每一位號碼數(shù)字的DTMF信號的DFT幅度譜，安照幅度譜的最大值確定對應的頻率，再安照頻率確定每一位對應的號碼數(shù)字，最后輸出6位電話號碼，程序名為exp6。</p><p>　　程序分四段：第一段（2—7行）設置參數(shù)，并讀入6位電話號碼；第二段（9—20行）根據(jù)鍵入的6位電話號碼產(chǎn)生時域離散DTMF信號，并連續(xù)發(fā)出6位號碼對應的雙音頻聲音；第三段（22—25行）對時域離散DTMF信號進行頻率檢測，畫出幅度譜；第

46、四段（26—33行）根據(jù)幅度譜的兩個峰值，分別查找并確定輸入6位電話號碼。根據(jù)程序中的注釋很容易分析編程思想和處理算法。程序清單如下：</p><p><b>　　exp6.m</b></p><p>　　% DTMF雙頻撥號信號的生成和檢測程序</p><p>　　%clear all；clc；</p><p>　　t

47、m=[1，2，3，65；4，5，6，66；7，8，9，67；42，0，35，68]； % DTMF信號代表的16個數(shù)</p><p>　　N=205；K=[18，20，22，24，31，34，38，42]；</p><p>　　f1=[697，770，852，941]； % 行頻率向量</p><p>　　f2=[1209，1

48、336，1477，1633]； % 列頻率向量</p><p>　　TN=input('鍵入6位電話號碼= ')； % 輸入6位數(shù)字</p><p>　　TNr=0； %接收端電話號碼初值為零</p><p>　　for l=1:6；<

49、;/p><p>　　d=fix(TN/10^(6-l))；</p><p>　　TN=TN-d*10^(6-l)；</p><p>　　for p=1:4；</p><p>　　for q=1:4；</p><p>　　if tm(p，q)==abs(d)； break，end % 檢測碼相符的列號q</

50、p><p><b>　　end</b></p><p>　　if tm(p，q)==abs(d)； break，end % 檢測碼相符的行號p</p><p><b>　　end</b></p><p>　　n=0:1023； % 為

51、了發(fā)聲，加長序列</p><p>　　x = sin(2*pi*n*f1(p)/8000) + sin(2*pi*n*f2(q)/8000)；% 構成雙頻信號</p><p>　　sound(x，8000)； % 發(fā)出聲音</p><p>　　pause(0.1)</p><p>　

52、　% 接收檢測端的程序</p><p>　　X=goertzel(x(1:205)，K+1)； % 用Goertzel算法計算八點DFT樣本</p><p>　　val = abs(X)； % 列出八點DFT向量</p><p>　　subplot(3，2，l)；</p>&

53、lt;p>　　stem(K，val，'.')；grid；xlabel('k')；ylabel('|X(k)|') % 畫出DFT(k)幅度</p><p>　　axis([10 50 0 120])</p><p>　　limit = 80； %</p><p>　　for s=5

54、:8；</p><p>　　if val(s) > limit， break， end % 查找列號</p><p><b>　　end</b></p><p>　　for r=1:4；</p><p>　　if val(r) > limit， break， end % 查找行號<

55、;/p><p><b>　　end</b></p><p>　　TNr=TNr+tm(r，s-4)*10^(6-l)；</p><p><b>　　end</b></p><p>　　disp('接收端檢測到的號碼為：') % 顯示接收到的字符</p>&l

56、t;p><b>　　disp(TNr)</b></p><p>　　運行程序，根據(jù)提示鍵入6位電話號碼123456，回車后可以聽見6位電話號碼對應的DTMF信號的聲音，并輸出相應的6幅頻譜圖如圖10.10.1所示，左上角的第一個圖在k=18和k=31兩點出現(xiàn)峰值，所以對應第一位號碼數(shù)字1。最后顯示檢測到的電話號碼123456。</p><p>　　圖1 6位電

57、話號碼123456的DTMF信號在8個近似基頻點的DFT幅度</p><p><b>　　（二）設計內(nèi)容</b></p><p>　　任意送入8位電話號碼，打印出相應的幅度譜。觀察程序運行結果，對照表1和表2，判斷程序譜分析的正確性。</p><p><b>　　（三）說明書要求</b></p><p&

58、gt;　　1．研究雙音多頻撥號（DTMF）系統(tǒng)，研究電話中雙音多頻信號的產(chǎn)生與檢測原理；</p><p>　　2．分析程序exp8.m，畫出仿真程序流程圖。</p><p>　　3．打印6位和8位電話號碼DTMF信號的幅度譜。</p><p>　　4．說明DTMF信號的參數(shù)：采樣頻率、DFT的變換點數(shù)以及觀測時間的確定原則。</p><p>

59、<b>　　四、參考文獻</b></p><p>　　1 高西全，丁玉美. 數(shù)字信號處理（第三版）.西安：西安電子科技大學出版社，2008.8</p><p>　　2 王宏. MATLAB6.5及其在信號處理中的應用. 北京：清華大學出版社，2004</p><p>　　3 李正周. MATLAB數(shù)字信號處理與應用. 北京:清華大學出版社，20

60、08</p><p>　　4 徐明遠，劉增力MATLAB仿真在信號處理中的應用.西安:電子科技大學出版社，2007 </p><p>　　5 李顯宏. MATLAB 7.x界面設計與編譯技巧. 北京: 電子工業(yè)出版社， 2006</p><p>　　6 郭仕劍等. MATLAB7.X數(shù)字信號處理. 人民郵電出版社， 2006</p><p>