聲音的延時(shí)和混響仿真課程設(shè)計(jì)任務(wù)書(shū)

1、 PAGE11 / NUMPAGES36 課程設(shè)計(jì)任務(wù)書(shū)學(xué)生：專(zhuān)業(yè)班級(jí)： 通信0906班 指導(dǎo)教師： 黃錚 工作單位： 信息工程學(xué)院題 目:聲音的延時(shí)和混響仿真初始條件：MATLAB軟件，電腦要求完成的主要任務(wù):錄制一段自己的語(yǔ)音信號(hào)，時(shí)間控制在1s左右，并對(duì)錄制的信號(hào)進(jìn)行采樣。語(yǔ)音信號(hào)的頻譜分析，畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形和頻譜圖。將信號(hào)加入延時(shí)和混響，并與原始信號(hào)頻譜進(jìn)行比較。設(shè)計(jì)單回聲濾波器，多重回聲濾波器，無(wú)限個(gè)回聲濾波器，全通結(jié)構(gòu)的混響器，并畫(huà)出濾波器的頻域響應(yīng)。用濾波器對(duì)語(yǔ)音信號(hào)進(jìn)行濾波。畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比，分析信號(hào)的變化?；胤耪Z(yǔ)音信號(hào)。時(shí)

2、間安排：序號(hào)階段容所需時(shí)間1查閱資料2天2各部分編程3天4調(diào)試程序2天5GUI設(shè)計(jì)3天答辯1天合 計(jì)11天指導(dǎo)教師簽名： 2012 年 月 日系主任（或責(zé)任教師）簽名： 2012 年 月 日 摘要數(shù)字信號(hào)處理(Digital Signal Processing)技術(shù)，從20世紀(jì)60年代以來(lái)，隨著計(jì)算機(jī)科學(xué)和信息科學(xué)發(fā)展，數(shù)字處理技術(shù)應(yīng)運(yùn)而生并得以快速發(fā)展。 語(yǔ)言是人們進(jìn)行信息溝通的主要方式之一，它具有直接、自然、方便等優(yōu)點(diǎn)。語(yǔ)音則是語(yǔ)言的物理層表達(dá)方式。語(yǔ)音處理主要是對(duì)語(yǔ)音進(jìn)行機(jī)器處理，以達(dá)到傳輸、自動(dòng)識(shí)別、機(jī)器理解等目的。進(jìn)行了語(yǔ)音處理過(guò)程的濾波、采樣、傅立葉變換和譜包絡(luò)提取的算法實(shí)現(xiàn)研究，

3、討論了在算法的DSP實(shí)現(xiàn)方法，Matlab 語(yǔ)言是一種廣泛應(yīng)用于工程計(jì)算與數(shù)值分析領(lǐng)域的新型高級(jí)語(yǔ)言, Matlab 功能強(qiáng)大、簡(jiǎn)單易學(xué)、編程效率高。特別是Matlab 還具有信號(hào)分析工具箱, 不需具備很強(qiáng)的編程能力, 就可以很方便地進(jìn)行信號(hào)分析、處理和設(shè)計(jì)關(guān)鍵詞 語(yǔ)音信號(hào)處理 MATLAB語(yǔ)言 濾波器 頻譜分析 GUIAbstractDigital Signal Processing (Digital Signal Processing) technology from the 1960s,along with the development of computer science and

4、information science, digital processing techniques have emerged and to rapid development. Language is the communication of people to one of the main, it has a direct, natural, and easy. Voice is the language of the physical layer of expression. Voice processing machines mainly for voice processing,

5、in order to achieve transmission, automatic recognition, machine understanding of other purposes. For voice processing filtering, sampling, Fourier transform and spectral envelope extraction algorithm research, discusses the implementation of DSP algorithms, Matlab language is a widely used engineer

6、ing calculation and numerical analysis of the field of new high-level language, Matlab powerful, easy to learn, programming and high efficiency. Matlab also has a particular signal analysis toolbox, need not have strong programming skills, you can easily carry out signal analysis, processing and des

7、ign Keywords Speech signal processing Spectral analysis of MATLAB Filter Spectrum analysis Graphical User Interface目 錄TOC o 1-3 h u HYPERLINK l _Toc17827 摘要 PAGEREF _Toc17827 I HYPERLINK l _Toc30728 Abstract PAGEREF _Toc30728 II HYPERLINK l _Toc25851 1.引言 PAGEREF _Toc25851 5 HYPERLINK l _Toc22636 1.

8、1簡(jiǎn)介 PAGEREF _Toc22636 5 HYPERLINK l _Toc15147 1.2設(shè)計(jì)任務(wù)與要求 PAGEREF _Toc15147 6 HYPERLINK l _Toc22505 2 設(shè)計(jì)原理 PAGEREF _Toc22505 7 HYPERLINK l _Toc32300 2.1語(yǔ)音信號(hào)的采樣 PAGEREF _Toc32300 7 HYPERLINK l _Toc29601 2.2時(shí)域分析 PAGEREF _Toc29601 7 HYPERLINK l _Toc378 2.3頻域分析 PAGEREF _Toc378 8 HYPERLINK l _Toc22643 2.4

9、 延時(shí) PAGEREF _Toc22643 8 HYPERLINK l _Toc12393 2.5 混響 PAGEREF _Toc12393 9 HYPERLINK l _Toc13970 2.6 濾波器的設(shè)計(jì) PAGEREF _Toc13970 10 HYPERLINK l _Toc19428 2.7 圖形用戶界面的設(shè)計(jì) PAGEREF _Toc19428 10 HYPERLINK l _Toc24249 3 設(shè)計(jì)步驟 PAGEREF _Toc24249 12 HYPERLINK l _Toc3917 3.1 總體設(shè)計(jì)流程圖 PAGEREF _Toc3917 12 HYPERLINK l _

10、Toc24071 3.2 詳細(xì)設(shè)計(jì) PAGEREF _Toc24071 13 HYPERLINK l _Toc15719 3.2.1 語(yǔ)音的錄制 PAGEREF _Toc15719 13 HYPERLINK l _Toc13777 3.2.2 原始信號(hào) PAGEREF _Toc13777 13 HYPERLINK l _Toc23907 3.2.3 延時(shí) PAGEREF _Toc23907 14 HYPERLINK l _Toc30442 3.2.4 混響 PAGEREF _Toc30442 15 HYPERLINK l _Toc29070 3.2.5 單回聲濾波器 PAGEREF _Toc2

11、9070 17 HYPERLINK l _Toc1067 3.2.6 多回聲濾波器 PAGEREF _Toc1067 18 HYPERLINK l _Toc31790 3.2.7 無(wú)限回聲濾波器 PAGEREF _Toc31790 19 HYPERLINK l _Toc7747 3.2.8 全通濾波器 PAGEREF _Toc7747 21 HYPERLINK l _Toc3220 4圖形用戶界面的設(shè)計(jì) PAGEREF _Toc3220 23 HYPERLINK l _Toc11040 4.1 GUI布局 PAGEREF _Toc11040 23 HYPERLINK l _Toc7533 4.

12、2 設(shè)置屬性 PAGEREF _Toc7533 23 HYPERLINK l _Toc32103 4.3 編程 PAGEREF _Toc32103 24 HYPERLINK l _Toc5837 4.4 運(yùn)行結(jié)果 PAGEREF _Toc5837 24 HYPERLINK l _Toc17614 5總結(jié) PAGEREF _Toc17614 2 HYPERLINK l _Toc6556 參考文獻(xiàn) PAGEREF _Toc6556 3 HYPERLINK l _Toc22123 附錄 PAGEREF _Toc22123 41.引言1.1簡(jiǎn)介數(shù)字信號(hào)處理主要是研究用數(shù)字或符號(hào)序列表示和處理信號(hào)。處理

13、的目的可以是削弱信號(hào)中的多余容，濾除混雜的噪聲和干擾，或者是將信號(hào)變換為容易分析和識(shí)別的形式，便于估計(jì)和選擇它的特征參數(shù)。語(yǔ)音信號(hào)分析是語(yǔ)音信號(hào)處理的前提和基礎(chǔ)，只有分析出可表示語(yǔ)音信號(hào)本質(zhì)特征的參數(shù)，才有可能利用這些參數(shù)進(jìn)行高效的語(yǔ)音通信、語(yǔ)音合成和語(yǔ)音識(shí)別等處理。而且，語(yǔ)音合成的音質(zhì)好壞，語(yǔ)音識(shí)別率的高低，也都取決于對(duì)語(yǔ)音信號(hào)分橋的準(zhǔn)確性和精確性。因此語(yǔ)音信號(hào)分析在語(yǔ)音信號(hào)處理應(yīng)用中具有舉足輕重的地位。根據(jù)所分析出的參數(shù)的性質(zhì)的不同，可將語(yǔ)音信號(hào)分析分為時(shí)域分析、頻域分析、倒領(lǐng)域分析等；時(shí)域分析方法具有簡(jiǎn)單、計(jì)算量小、物理意義明確等優(yōu)點(diǎn)，但由于語(yǔ)音信號(hào)最重要的感知特性反映在功率譜中，而相位

14、變化只起著很小的作用，所以相對(duì)于時(shí)域分析來(lái)說(shuō)頻域分析更為重要。本文將簡(jiǎn)要介紹時(shí)域分析、頻域分析。聲音信號(hào)是一維連續(xù)信號(hào),而計(jì)算機(jī)只能處理離散信號(hào)。為了從離散信號(hào)還原連續(xù)信號(hào),根據(jù)采樣定理,可以確定采樣頻率的最小值。wav文件是一種數(shù)字聲音文件格式,本課程設(shè)計(jì)基于Matlab分析了wav聲音文件頻譜與聲音的關(guān)系。通過(guò)采集個(gè)人的一段聲音進(jìn)行頻譜分析等處理，然后設(shè)計(jì)數(shù)字濾波器處理這個(gè)原始聲音的wav文件,并比較濾波以后輸出聲音信號(hào)與原聲音信號(hào)的異同。1.2設(shè)計(jì)任務(wù)與要求設(shè)計(jì)任務(wù)：(1)利用Windows下的錄音機(jī)或其他軟件，錄制一段自己的語(yǔ)音信號(hào)，時(shí)間控制在1s左右，并對(duì)錄制的信號(hào)進(jìn)行采樣。(2)語(yǔ)

15、音信號(hào)的頻譜分析，畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形和頻譜圖。(3)將信號(hào)加入延時(shí)和混響，再分析其頻譜，并與原始信號(hào)頻譜進(jìn)行比較。(4)設(shè)計(jì)幾種特殊類(lèi)型的濾波器：?jiǎn)位芈暈V波器，多重回聲濾波器，無(wú)限個(gè)回聲濾波器，全通結(jié)構(gòu)的混響器，并畫(huà)出濾波器的頻域響應(yīng)。(5)用自己設(shè)計(jì)的濾波器對(duì)采集的語(yǔ)音信號(hào)進(jìn)行濾波。(6)分析得到信號(hào)的頻譜，畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比，分析信號(hào)的變化。(7)回放語(yǔ)音信號(hào)。設(shè)計(jì)要求(1)熟悉離散信號(hào)和系統(tǒng)的時(shí)域特性。(2)熟悉語(yǔ)音信號(hào)的特點(diǎn)。(3)掌握數(shù)字信號(hào)處理的基本概念，基本理論和基本方法。(4)掌握序列快速傅里葉變換方法。(5)學(xué)會(huì)MATLAB的

16、使用，掌握MATLAB的程序設(shè)計(jì)方法。(6)掌握MATLAB設(shè)計(jì)各種數(shù)字濾波器的方法和對(duì)信號(hào)進(jìn)行濾波的方法2 設(shè)計(jì)原理2.1語(yǔ)音信號(hào)的采樣采樣頻率采樣頻率是指計(jì)算機(jī)每秒鐘采集多少個(gè)聲音樣本，是描述聲音文件的音質(zhì)、音調(diào)，衡量聲卡、聲音文件的質(zhì)量標(biāo)準(zhǔn)。采樣頻率越高，即采樣的間隔時(shí)間越短，則在單位時(shí)間計(jì)算機(jī)得到的聲音樣本數(shù)據(jù)就越多，對(duì)聲音波形的表示也越精確。這就是說(shuō)采樣頻率是衡量聲卡采集、記錄和還原聲音文件的質(zhì)量標(biāo)準(zhǔn)。(2)采樣位數(shù)采樣位數(shù)即采樣值或取樣值，用來(lái)衡量聲音波動(dòng)變化的參數(shù)，是指聲卡在采集和播放聲音文件時(shí)所使用數(shù)字聲音信號(hào)的二進(jìn)制位數(shù)。采樣頻率是指錄音設(shè)備在一秒鐘對(duì)聲音信號(hào)的采樣次數(shù)，采樣

17、頻率越高聲音的還原就越真實(shí)越自然。采樣位數(shù)和采樣率對(duì)于音頻接口來(lái)說(shuō)是最為重要的兩個(gè)指標(biāo)，也是選擇音頻接口的兩個(gè)重要標(biāo)準(zhǔn)。無(wú)論采樣頻率如何，理論上來(lái)說(shuō)采樣的位數(shù)決定了音頻數(shù)據(jù)最大的力度圍。每增加一個(gè)采樣位數(shù)相當(dāng)于力度圍增加了6dB。采樣位數(shù)越多則捕捉到的信號(hào)越精確。對(duì)于采樣率來(lái)說(shuō)你可以想象它類(lèi)似于一個(gè)照相機(jī)，44.1kHz意味著音頻流進(jìn)入計(jì)算機(jī)時(shí)計(jì)算機(jī)每秒會(huì)對(duì)其拍照達(dá)441000次。顯然采樣率越高，計(jì)算機(jī)攝取的圖片越多，對(duì)于原始音頻的還原也越加精確。(3)采樣定理在進(jìn)行模擬/數(shù)字信號(hào)的轉(zhuǎn)換過(guò)程中，當(dāng)采樣頻率fsmax大于信號(hào)最高頻率fmax的2倍時(shí)，即：fsmax=2fmax,則采樣之后的數(shù)字信

18、號(hào)完整地保留了原始信號(hào)中的信息，一般實(shí)際應(yīng)用中保證采樣頻率為信號(hào)最高頻率的510倍，采樣定理又稱(chēng)奈奎斯特定理。1924年奈奎斯特(Nyquist)就推導(dǎo)出：在理想低通信道的最大碼元傳輸速率=2W*log2 N (其中W是理想低通信道的帶寬,N是電平強(qiáng)度)。2.2時(shí)域分析語(yǔ)音信號(hào)的時(shí)域分析就是分析和提取語(yǔ)音信號(hào)的時(shí)域參數(shù)。進(jìn)行語(yǔ)音分析時(shí)，最先接觸到并且也是最直觀的是它的時(shí)域波形。語(yǔ)音信號(hào)本身就是時(shí)域信號(hào)，因而時(shí)域分析是最早使用，也是應(yīng)用最廣泛的一種分析方法，這種方法直接利用語(yǔ)音信號(hào)的時(shí)域波形。時(shí)域分析通常用于最基本的參數(shù)分析與應(yīng)用，如語(yǔ)音的分割、預(yù)處理、大分類(lèi)等。這種分析方法的特點(diǎn)是： = 1

19、* GB3 表示語(yǔ)音信號(hào)比較直觀、物理意義明確。 = 2 * GB3 實(shí)現(xiàn)起來(lái)比較簡(jiǎn)單、運(yùn)算且少。 = 3 * GB3 可以得到語(yǔ)音的一些重要的參數(shù)。只使用示波器等通用設(shè)備，使用較為簡(jiǎn)單等。2.3頻域分析在MATLAB的信號(hào)處理工具箱中函數(shù)FFT和IFFT用于快速傅立葉變換和逆變換。下面介紹這些函數(shù)：（1）函數(shù)FFT用于序列快速傅立葉變換。函數(shù)的一種調(diào)用格式為：y=fft(x)其中，x是序列，y是序列的FFT，x可以為一向量或矩陣，若x為一向量，y是x的FFT。且和x一樣長(zhǎng)度。若x為一矩陣，則y是對(duì)矩陣的每一列向量進(jìn)行FFT。如果x長(zhǎng)度是2的冪次方，函數(shù)fft執(zhí)行高速基2FFT算法；否則fft

20、執(zhí)行一種混合基的離散傅立葉變換算法，計(jì)算速度較慢。（2）函數(shù)FFT的另一種調(diào)用格式為：y=fft(x,N)式中，x，y意義同前，N為正整數(shù)。函數(shù)執(zhí)行N點(diǎn)的FFT。若x為向量且長(zhǎng)度小于N，則函數(shù)將x補(bǔ)零至長(zhǎng)度N。若向量x的長(zhǎng)度大于N，則函數(shù)截短x使之長(zhǎng)度為N。若x 為矩陣，按一樣方法對(duì)x進(jìn)行處理。（3）經(jīng)函數(shù)fft求得的序列y一般是復(fù)序列，通常要求其幅值和相位。MATLAB提供求復(fù)數(shù)的幅值和相位函數(shù)：abs，angle，這些函數(shù)一般和FFT同時(shí)使用。函數(shù)abs(x)用于計(jì)算復(fù)向量x的幅值，函數(shù)angle(x)用于計(jì)算復(fù)向量的相角，以弧度表示。函數(shù)unwrap(p)用于展開(kāi)弧度相位角p 。（4）f

21、ftshift作用：將零頻點(diǎn)移到頻譜的中間用法：Y=fftshift(X)Y=fftshift(X，dim)描述：fftshift移動(dòng)零頻點(diǎn)到頻譜中間，重新排列fft,fft2和fftn的輸出結(jié)果。將零頻點(diǎn)放到頻譜的中間對(duì)于觀察傅立葉變換是有用的。2.4 延時(shí) 延時(shí)就是將音源延遲一段時(shí)間后，再欲播放的效果處理。依其延遲時(shí)間的不同，可分別產(chǎn)生合唱、鑲邊、回音等效果。當(dāng)延遲時(shí)間在335ms之間時(shí)人耳感覺(jué)不到滯后音的存在，并且他與原音源疊加后，會(huì)因其相位干涉而產(chǎn)生梳狀濾波效應(yīng)，這就是鑲邊效果。如果延遲時(shí)間在50ms以上時(shí)，其延遲音就清晰可辨，此時(shí)的處理效果才是回音?；匾籼幚硪话愣际怯糜诋a(chǎn)生簡(jiǎn)單的混響

22、效果。延時(shí)、合唱、鑲邊、回音等效果的可調(diào)參數(shù)都差不多，具體有以下幾項(xiàng)：延時(shí)時(shí)間（Dly），即主延時(shí)電路的延時(shí)時(shí)間調(diào)整。反饋增益（FB Gain），即延時(shí)反饋的增益控制。反饋高頻比（Hi Ratio），即反饋回路上的高頻衰減控制。調(diào)制頻率（Freq），指主延時(shí)的調(diào)頻周期。調(diào)制深度（Depth），指上述調(diào)頻電路的調(diào)制深度。高頻增益（HF），指高頻均衡控制。預(yù)延時(shí)（Ini Dly），指主延時(shí)電路預(yù)延時(shí)時(shí)間調(diào)整。均衡頻率（EQ F），這里的頻率均衡用于音色調(diào)整，此為均衡的中點(diǎn)頻率選擇。由于延時(shí)產(chǎn)生的效果都比較復(fù)雜多變，如果不是效果處理專(zhuān)家，建議使用設(shè)備提供的預(yù)置參數(shù)，因?yàn)檫@些預(yù)置參數(shù)給出的處理效果一般

23、都比較好。2.5 混響 混響效果主要是用于增加音源的融合感。自然音源的延時(shí)聲陣列非常密集、復(fù)雜，所以模擬混響效果的程序也復(fù)雜多變。常見(jiàn)參數(shù)有以下幾種：混響時(shí)間：能逼真的模擬自然混響的數(shù)碼混響器上都有一套復(fù)雜的程序，其上雖然有很多技術(shù)參數(shù)可調(diào)，然而對(duì)這些技術(shù)參數(shù)的調(diào)整都不會(huì)比原有的效果更為自然，尤其是混響時(shí)間。高頻滾降：此項(xiàng)參數(shù)用于模擬自然混響當(dāng)中，空氣對(duì)高頻的吸收效應(yīng)，以產(chǎn)生較為自然的混響效果。一般高頻混降的可調(diào)圍為0.11.0。此值較高時(shí)，混響效果也較接近自然混響；此值較低時(shí)，混響效果則較清澈。擴(kuò)散度：此項(xiàng)參數(shù)可調(diào)整混響聲陣密度的增長(zhǎng)速度，其可調(diào)圍為010，其值較高時(shí)，混響效果比較豐厚、溫暖

24、；其值較低時(shí)，混響效果則較空曠、冷僻。預(yù)延時(shí)：自然混響聲陣的建立都會(huì)延遲一段時(shí)間，預(yù)延時(shí)即為模擬次效應(yīng)而設(shè)置。聲陣密度：此項(xiàng)參數(shù)可調(diào)整聲陣的密度，其值較高時(shí)，混響效果較為溫暖，但有明顯的聲染色；其值較低時(shí)，混響效果較深邃，切聲染色也較弱。頻率調(diào)制：這是一項(xiàng)技術(shù)性的參數(shù)，因?yàn)殡娮踊祉懙穆曣嚸芏缺茸匀换祉懴∈瑁瑸榱耸够祉懙穆曇舯容^平滑、連貫，需要對(duì)混響聲陣列的延時(shí)時(shí)間進(jìn)行調(diào)制。此項(xiàng)技術(shù)可以有效的消除延時(shí)聲陣列的段裂聲，可以增加混響聲的柔和感。調(diào)治深度：指上述調(diào)頻電路的調(diào)治深度。2.6 濾波器的設(shè)計(jì)單回聲濾波器的系統(tǒng)函數(shù)：H(z)= 1 （2-1）無(wú)限個(gè)回聲濾波器的系統(tǒng)函數(shù):H(z)=1 （2-2）

25、全通結(jié)構(gòu)的混響器的系統(tǒng)函數(shù):H(z)=1 （2-3）函數(shù)filter的調(diào)用格式為：y=filter(b,a,x)該格式采用數(shù)字濾波器對(duì)數(shù)據(jù)進(jìn)行濾波，既可以用于IIR濾波器，也可以用于FIR濾波器。其中向量b和a分別表示系統(tǒng)函數(shù)的分子、分母多項(xiàng)式的系數(shù)，若a1，此時(shí)表示FIR濾波器，否則就是IIR濾波器。該函數(shù)是利用給出的向量b和a，對(duì)x中的數(shù)據(jù)進(jìn)行濾波，結(jié)果放入向量y。2.7 圖形用戶界面的設(shè)計(jì)GUI(Graphical User Interface)圖形用戶界面，是在圖形界面下安排顯示與用戶交互的組件元素，用戶可以只通過(guò)鍵盤(pán)、鼠標(biāo)和前臺(tái)界面下的組件發(fā)生交互，而所有的計(jì)算、繪圖等部操作都封裝在

26、部，提高了終端用戶使用MATLAB程序的易用性。設(shè)計(jì)圖形用戶界面應(yīng)遵循以下三個(gè)原則：簡(jiǎn)單性（Simplicity）、一致性（Consistency）與習(xí)常性（Familiarity）。（1）簡(jiǎn)單性設(shè)計(jì)界面時(shí)，應(yīng)力求簡(jiǎn)潔、直接、清晰地體現(xiàn)出界面的功能和特征。設(shè)計(jì)的圖形界面要直觀，為此應(yīng)多采用圖形，而間量避免數(shù)值。設(shè)計(jì)界面應(yīng)間量減少窗口數(shù)目，力避在不同窗口之間進(jìn)行來(lái)回切換。（2）一致性所謂一致性有兩層含義：一是讀者自己開(kāi)發(fā)的界面風(fēng)格要盡量一致；二是新設(shè)計(jì)的界面要與其他已有的界面的風(fēng)格不要截然向左。（3）習(xí)常性設(shè)計(jì)界面時(shí)，應(yīng)盡量使用人們所熟悉的標(biāo)志與符號(hào)。用戶可以根據(jù)熟悉標(biāo)志作出正確猜測(cè)，自學(xué)入門(mén)。

27、（4）其他考慮因素注意界面的動(dòng)態(tài)性能。如界面對(duì)用戶操作的響應(yīng)要迅速(Immediate)、連續(xù)（Continous）；對(duì)持續(xù)時(shí)間較長(zhǎng)的運(yùn)算，要給出等待時(shí)間提示，并允許用戶中斷運(yùn)算。圖形用戶界面的制作包括界面設(shè)計(jì)和程序?qū)崿F(xiàn)。具體制作步驟如下： （1）分析界面所要現(xiàn)的主要功能，明確設(shè)計(jì)任務(wù)； （2）在稿紙上繪出界面草圖，并站在使用者的角度來(lái)審查草圖； （3）按構(gòu)思的草圖，上機(jī)制作靜態(tài)界面，并檢查之； （4）編寫(xiě)界面動(dòng)態(tài)功能的程序，對(duì)功能進(jìn)行逐項(xiàng)檢查。 3 設(shè)計(jì)步驟用luyin程序錄制一段1s左右的wav格式的語(yǔ)音，并保存3.1 總體設(shè)計(jì)流程圖用wavread從電腦中導(dǎo)入錄制的語(yǔ)音信號(hào)調(diào)用fft函數(shù)

28、進(jìn)行快速傅里葉變換。畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回放采樣語(yǔ)音信號(hào)對(duì)采樣信號(hào)延時(shí)2000，兩者相加，得到混響。分別畫(huà)出延時(shí)和混響后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖。并回放延時(shí)和混響后語(yǔ)音信號(hào)c.無(wú)限回聲濾波器，畫(huà)出頻域響應(yīng)。用濾波器對(duì)語(yǔ)音信號(hào)濾波,畫(huà)出時(shí)域波形與頻譜圖回放語(yǔ)音信號(hào)D.全通濾波器，畫(huà)出頻域響應(yīng)。用濾波器對(duì)語(yǔ)音信號(hào)濾波,畫(huà)出時(shí)域波形與頻譜圖回放語(yǔ)音信號(hào)b.多聲濾波器，畫(huà)出頻域響應(yīng)。用濾波器對(duì)語(yǔ)音信號(hào)濾波,畫(huà)出時(shí)域波形與頻譜圖回放語(yǔ)音信號(hào)a.單回聲濾波器，畫(huà)出頻域響應(yīng)。用濾波器對(duì)語(yǔ)音信號(hào)濾波,畫(huà)出時(shí)域波形與頻譜圖回放語(yǔ)音信號(hào)分析得到信號(hào)的頻譜，比較回放的聲音，并對(duì)濾波前后的信號(hào)進(jìn)

29、行對(duì)比，分析信號(hào)的變化。圖3-1 總體設(shè)計(jì)流程圖3.2 詳細(xì)設(shè)計(jì)3.2.1 語(yǔ)音的錄制運(yùn)行程序，錄制語(yǔ)音文件，文件名為sound.wav程序：%獲取錄音文件 pausefs=44100; %聲音的采樣頻率為44.1Khz duration=3; %錄音時(shí)間為3s fprintf(按任意鍵開(kāi)始錄音1：n); pausefprintf(錄音中n);sd1=wavrecord(duration*fs,fs); %采樣數(shù)為132300，保存文件名為sd1.wavfprintf(放音中n);wavplay(sd1,fs);fprintf(錄音1播放完畢。n);wavwrite(sd1,fs,sound1

30、.wav); %將錄音文件保存為WAV格式的聲音文件3.2.2 原始信號(hào)讀取聲音,畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回放采樣語(yǔ)音信號(hào)。程序：x1,fs,bits=wavread(D:sound1.wav); %讀取信號(hào),采樣值放在向量x中，fs表示采樣頻率(Hz)，bits表示采樣位數(shù)。x1=x1(:,1); %只取單聲道sound(x1,fs); %對(duì)聲音的回放figure;Subplot(211);plot(x1); %做原始語(yǔ)音信號(hào)的時(shí)域圖形title(原始語(yǔ)音信號(hào)的時(shí)域圖);grid on;xlabel(時(shí)間 n);ylabel(音量 n);y1=fft(x1); %做lengt

31、h(x1)點(diǎn)的FFTy1=fftshift(y1); %平移，是頻率中心為0derta_fs = fs/length(x1); %設(shè)置頻譜的間隔，分辨率Subplot(212);plot(-fs/2:derta_fs:fs/2-derta_fs,abs(y1);%畫(huà)出原始語(yǔ)音信號(hào)的頻譜圖title(原始語(yǔ)音信號(hào)的頻譜);grid on;運(yùn)行結(jié)果：圖3-2 原始語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：采樣點(diǎn)數(shù)為132300，與時(shí)域波形結(jié)果吻合，由頻譜可看出頻率較集中，有一點(diǎn)高頻分量，是由于噪聲的干擾引起的。3.2.3 延時(shí)對(duì)原始語(yǔ)音信號(hào)延時(shí)2000,畫(huà)出延時(shí)后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回放延時(shí)后

32、的語(yǔ)音信號(hào)。程序： %延時(shí)的信號(hào) x1,fs,bits=wavread(D:sound1.wav); x1=x1(:,1); %只取單聲道 z=zeros(2000,1);x1; %對(duì)聲音延時(shí)2000 sound(z,fs); figure(2); Subplot(211); plot(z); title(延時(shí)后的信號(hào)時(shí)域圖);grid on; %畫(huà)出延時(shí)后的信號(hào)時(shí)域圖 z=fft(x1); %做length(x1)點(diǎn)的FFT z=fftshift(z); %平移，是頻率中心為0 derta_fs = fs/length(z); %設(shè)置頻譜的間隔，分辨率 Subplot(212); plot(

33、-fs/2:derta_fs:fs/2-derta_fs,abs(z);%畫(huà)出延時(shí)后的信號(hào)頻譜圖 title(延時(shí)后的信號(hào)頻譜圖);grid on; 運(yùn)行結(jié)果：圖3-3延時(shí)后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，時(shí)域波形向右平移，頻譜無(wú)變化。3.2.4 混響 調(diào)整原始采樣信號(hào)的長(zhǎng)度，將其與延時(shí)信號(hào)相加，就得到混響。畫(huà)出混響后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回放延時(shí)后的語(yǔ)音信號(hào)。程序： %混響后信號(hào) x,fs,bits=wavread(D:sound1.wav); %讀取語(yǔ)音信號(hào) x=x(:,1); %只取單聲道 z=zeros(2000,1);x; %對(duì)語(yǔ)音信號(hào)進(jìn)行延時(shí) x=x

34、;zeros(2000,1); %使語(yǔ)音信號(hào)與延時(shí)后信號(hào)同等長(zhǎng)度 y1=x+z; %信號(hào)的混響 sound(y1,fs); %回放混響后的信號(hào) figure(3);subplot(211);plot(y1); title(混響的時(shí)域圖); %混響時(shí)域圖 Y1=fft(y1); %做length(y1)點(diǎn)的FFT Y1=fftshift(Y1); %平移，是頻率中心為0 derta_fs = fs/length(y1); %設(shè)置頻譜的間隔，分辨率 Subplot(212); plot(-fs/2:derta_fs:fs/2-derta_fs,abs(Y1);%畫(huà)出混響頻譜圖 title(混響的頻

35、譜圖);grid on;運(yùn)行結(jié)果： 圖3-4混響后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，原始信號(hào)的延時(shí)聲陣列非常密集、復(fù)雜，混響后時(shí)域波形音量增大了，各個(gè)脈沖變寬，語(yǔ)音信號(hào)的頻譜展寬了，幅值增大，回放的語(yǔ)音聲音變得融合，能聽(tīng)到回音。3.2.5 單回聲濾波器用單回聲濾波器對(duì)原始語(yǔ)音信號(hào)濾波,畫(huà)出濾波后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回濾波時(shí)后的語(yǔ)音信號(hào)。程序： %單回聲濾波器 x,fs,bits=wavread(D:sound1.wav); %讀取語(yǔ)音信號(hào) x=x(:,1); %只取單聲道 a=0.5; %a取小于等于1 z=zeros(2000,1);x; %對(duì)語(yǔ)音信號(hào)進(jìn)行延時(shí)

36、 x=x;zeros(2000,1); %使語(yǔ)音信號(hào)與延時(shí)后信號(hào)同等長(zhǎng)度 y2=x+z*a; %信號(hào)經(jīng)單回聲濾波 sound(y2,fs); %回放經(jīng)單回聲濾波器后的語(yǔ)音信號(hào) figure;subplot(211);plot(y2);title(單回聲濾波的信號(hào)時(shí)域圖); Y2=fft(y2); %做length(y2)點(diǎn)的FFT Y2=fftshift(Y2); %平移，是頻率中心為0 derta_fs = fs/length(y2); %設(shè)置頻譜的間隔，分辨率 Subplot(212); plot(-fs/2:derta_fs:fs/2-derta_fs,abs(Y2);%畫(huà)出單回聲信號(hào)頻

37、譜圖 title(單回聲濾波的信號(hào)頻譜圖);grid on;運(yùn)行結(jié)果：圖3-5單回聲濾波器濾波后信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，時(shí)域波形的幅值變大，毛刺變多，頻譜的中心頻率兩邊的頻率分量展寬、幅值變大，回放的聲音立體感變強(qiáng)。3.2.6 多回聲濾波器用多回聲濾波器對(duì)原始語(yǔ)音信號(hào)濾波,畫(huà)出濾波后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回濾波時(shí)后的語(yǔ)音信號(hào)。程序：%多回聲濾波器 x,fs,bits=wavread(D:sound1.wav); x1=x(:,1); %變換為單聲道語(yǔ)音信號(hào) delay_t=0.1; %回聲延遲時(shí)間 a=0.2; %衰減系數(shù) N=3; %N階延時(shí) Bz=1;z

38、eros(fs*delay_t,1)*N;aN; Az=1;zeros(fs*delay_t,1);a; y=filter(Bz, Az ,x1 ); sound(y,fs,bits); %回放經(jīng)多回聲濾波器后的語(yǔ)音信號(hào) figure;subplot(211);plot(y);title(多回聲濾波的信號(hào)時(shí)域圖); Y2=fft(y); %做length(y)點(diǎn)的FFT Y2=fftshift(Y2); %平移，是頻率中心為0 derta_fs = fs/length(y); %設(shè)置頻譜的間隔，分辨率 Subplot(212); plot(-fs/2:derta_fs:fs/2-derta_f

39、s,abs(Y2); %畫(huà)出多回聲信號(hào)頻譜圖 title(多回聲濾波的信號(hào)頻譜圖);grid on;運(yùn)行結(jié)果：圖3-6單回聲濾波器濾波后信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，時(shí)域波形的幅值變大，毛刺變得更多，頻譜的中心頻率幅值變大，回放的聲音立體感變的更強(qiáng)。3.2.7 無(wú)限回聲濾波器用無(wú)限回聲濾波器對(duì)原始語(yǔ)音信號(hào)濾波,畫(huà)出濾波后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回濾波時(shí)后的語(yǔ)音信號(hào)。程序： %無(wú)限回聲濾波器 x,fs,bits=wavread(D:sound1.wav); %讀取語(yǔ)音信號(hào) x=x(:,1); %只取單聲道 a=0.5; %a取小于等于1 Bz=0,0,0,0,0,0,

40、0,0,0,0,1; %分子的系數(shù) Az=1,0,0,0,0,0,0,0,0,0,-a; %分母的系數(shù) yy1=filter(Bz,Az,x); %FIR濾波器進(jìn)行濾波 sound(yy1,fs,bits); %回放經(jīng)無(wú)限回聲濾波器后的語(yǔ)音信號(hào) figure(5);subplot(211);plot(yy1); %無(wú)限回聲信號(hào)時(shí)域波形 title(無(wú)限個(gè)回聲濾波的信號(hào)時(shí)域圖); YY1=fft(yy1); %經(jīng)無(wú)限回聲濾波器后的信號(hào)做length(yy1)點(diǎn)的FFT YY1=fftshift(YY1); %平移，是頻率中心為0 derta_fs = fs/length(yy1); %設(shè)置頻譜的

41、間隔，分辨率 Subplot(212);plot(-fs/2:derta_fs:fs/2-derta_fs,abs(YY1);%無(wú)限回聲頻譜圖 title(無(wú)限個(gè)回聲濾波的信號(hào)頻譜圖);grid on; 運(yùn)行結(jié)果：圖3-7無(wú)限回聲濾波器濾波后信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，時(shí)域波形的幅值變大，頻譜的中心頻率不變，中心頻率兩邊的頻率分量變得密集、幅值變小，回放的聲音變得清晰，音量變大。3.2.8 全通濾波器用全通濾波器對(duì)原始語(yǔ)音信號(hào)濾波,畫(huà)出濾波后語(yǔ)音信號(hào)的時(shí)域波形與頻譜圖，并回濾波時(shí)后的語(yǔ)音信號(hào)。程序： %全通結(jié)構(gòu)的混響器 x,fs,bits=wavread(D:sound

42、1.wav); %讀取語(yǔ)音信號(hào) x=x(:,1); %只取單聲道 a=0.5; %a取小于等于1 Bz1=a,0,0,0,0,0,0,0,0,0,1; %分子的系數(shù) Az1=1,0,0,0,0,0,0,0,0,0,a; %分母的系數(shù) yy2=filter(Bz1,Az1,x); %FIR濾波器進(jìn)行濾波 sound(yy2,fs,bits); %回放經(jīng)全通結(jié)構(gòu)的混響器后的語(yǔ)音信號(hào) figure(6); subplot(211);plot(yy2); title(全通濾波信號(hào)的時(shí)域圖) %全通結(jié)構(gòu)的混響器時(shí)域波形 YY2=fft(yy2); %經(jīng)全通濾波器后的信號(hào)做length(yy2)點(diǎn)的FFT

43、 YY2=fftshift(YY2); %平移，是頻率中心為0 derta_fs = fs/length(yy2); %設(shè)置頻譜的間隔，分辨率 Subplot(212); plot(-fs/2:derta_fs:fs/2-derta_fs,abs(YY2); %全通濾波信號(hào)的頻譜圖 title(全通濾波信號(hào)的頻譜圖);grid on;運(yùn)行結(jié)果：圖3-8無(wú)限回聲濾波器濾波后信號(hào)的時(shí)域波形與頻譜圖結(jié)果分析：由運(yùn)行結(jié)果可以看出，時(shí)域波形和頻譜基本不變，回放的聲音跟原語(yǔ)音信號(hào)很接近。4圖形用戶界面的設(shè)計(jì)4.1 GUI布局 在Matlab的主窗口中，選擇File菜單中的New菜單項(xiàng)，再選擇其中的GUI命

44、令，就會(huì)就會(huì)顯示GUI的設(shè)計(jì)模板。新建一個(gè)空白模板jiemian.fig，就會(huì)顯示GUI的設(shè)計(jì)窗口，點(diǎn)擊 和 分別放置按鈕和坐標(biāo)軸，放置結(jié)果如圖4-1。圖4-1GUI布局4.2設(shè)置屬性按鈕(Push Buttons)：執(zhí)行某種預(yù)定的功能或操作；其屬性如圖4-2。Axes屬性為默認(rèn)屬性，如圖4-3。未指定時(shí)將使用系統(tǒng)缺省值，改變了的屬性有：BackgroundColor取值為顏色的預(yù)定義字符或RGB數(shù)值，缺省值為淺灰色；本設(shè)計(jì)為藍(lán)色。String取值為字符串矩陣或塊數(shù)組，定義控件對(duì)象標(biāo)題或選項(xiàng)容；Tag取值為字符串，定義了控件的標(biāo)識(shí)值，在任何程序中都可以通過(guò)這個(gè)標(biāo)識(shí)值控制該控件對(duì)象；FontSi

45、ze取值為數(shù)值；圖4-2按鈕的屬性 圖4-3 坐標(biāo)軸的屬性4.3 編程使用圖形用戶界面開(kāi)發(fā)環(huán)境可方便地創(chuàng)建GUI應(yīng)用程序，它可以根據(jù)設(shè)計(jì)的GUI布局，自動(dòng)生成M文件的框架，使用這一框架編制應(yīng)用程序。將前面編好的各個(gè)功能子程序分別編寫(xiě)到各個(gè)按鈕的callback函數(shù)中，逐個(gè)調(diào)試，運(yùn)行，M文件見(jiàn)附錄。4.4 運(yùn)行結(jié)果圖形用戶界面的運(yùn)行結(jié)果如圖4-4。圖4-4 GUI運(yùn)行結(jié)果5總結(jié)通過(guò)本次課程設(shè)計(jì)，我收獲了很多：一、我對(duì)MATLAB軟件有了較深的了解，并且學(xué)會(huì)了用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行分析和處理。這次課程設(shè)計(jì) ，剛開(kāi)始我無(wú)從下手，不知所措，后來(lái)經(jīng)過(guò)問(wèn)同學(xué)、去圖書(shū)館資查料、上網(wǎng)查資料、看視頻教程，

46、逐漸了解了用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行分析和處理要用到的各種函數(shù)，與濾波器的設(shè)計(jì)與其頻率響應(yīng)的繪制，并在調(diào)試過(guò)程不斷出錯(cuò)的過(guò)程中，不斷加強(qiáng)了對(duì)MATLB編程的掌握程度，能夠用MATLB對(duì)語(yǔ)音信號(hào)進(jìn)行分析和處理，收獲頗深。二、本次課程設(shè)計(jì)鍛煉了我獨(dú)立解決問(wèn)題的能力。剛看到題目時(shí)我一點(diǎn)思路也沒(méi)有，到后來(lái)，自己不斷探索MATLB，不斷摸索各個(gè)函數(shù)的使用方法和格式，最后完成了本次課程設(shè)計(jì)，整個(gè)過(guò)程讓我學(xué)會(huì)了發(fā)現(xiàn)問(wèn)題，并獨(dú)立解決問(wèn)題的方法，也能更好的利用網(wǎng)絡(luò)學(xué)習(xí)知識(shí)。三、本次課程設(shè)計(jì)讓我對(duì)語(yǔ)音信號(hào)的采樣、延時(shí)、混響、濾波器的設(shè)計(jì)與語(yǔ)音信號(hào)的分析有了更深的了解，并且鞏固了數(shù)字信號(hào)處理所學(xué)的容。了解了通過(guò)M

47、ATLAB設(shè)計(jì)單回聲濾波器，多重回聲濾波器，無(wú)限個(gè)回聲濾波器，全通結(jié)構(gòu)的混響器的過(guò)程和方法，使所學(xué)的知識(shí)得到了綜合運(yùn)用與鞏固。本次課程設(shè)計(jì)最重要的是讓我學(xué)會(huì)了獨(dú)立解決問(wèn)題的方法，經(jīng)過(guò)不斷的學(xué)習(xí)、出錯(cuò)、調(diào)試、改正，鍛煉了我的耐心和毅力，讓我體會(huì)到了不斷探索的樂(lè)趣，并體會(huì)到了付出后得到成功的激動(dòng)。此次設(shè)計(jì)也讓我發(fā)現(xiàn)了我的不足：對(duì)MATLB的掌握程度不夠，在原來(lái)學(xué)習(xí)后，我對(duì)MATLAB的許多函數(shù)、屬性的用法還是不很了解，以后應(yīng)當(dāng)多多學(xué)習(xí)軟件，輔助各門(mén)專(zhuān)業(yè)課的學(xué)習(xí)，并激發(fā)對(duì)本專(zhuān)業(yè)的學(xué)習(xí)興趣。參考文獻(xiàn)1 鄒彥，DSP原理與應(yīng)用，第1版，電子工業(yè)，2006；2 啟琮等，DSP技術(shù)的發(fā)展與應(yīng)用，第二版，高等

48、教育，20073 泉，闕大順數(shù)字信號(hào)處理原理與實(shí)現(xiàn)電子工業(yè)20054 雄偉，DSP集成開(kāi)發(fā)與應(yīng)用實(shí)例，第1版，電子工業(yè)，2002；5 洪濤等，數(shù)字信號(hào)處理，第一版，華中科技大學(xué)，20076 述斌等，數(shù)字信號(hào)處理實(shí)踐教程，第一版，華中科技大學(xué)，20077 胡航語(yǔ)音信號(hào)處理工業(yè)大學(xué)20028 威MATLAB基礎(chǔ)與編程入門(mén)電子科技大學(xué)2006附錄圖形用戶界面的M文件：function varargout = jiemian(varargin)% JIEMIAN M-file for jiemian.fig% JIEMIAN, by itself, creates a new JIEMIAN or r

49、aises the existing% singleton*.% H = JIEMIAN returns the handle to a new JIEMIAN or the handle to% the existing singleton*.% JIEMIAN(CALLBACK,hObject,eventData,handles,.) calls the local% function named CALLBACK in JIEMIAN.M with the given input arguments.% JIEMIAN(Property,Value,.) creates a new JI

50、EMIAN or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before jiemian_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to jiemian_OpeningFcn via varargin.% *See GUI

51、Options on GUIDEs Tools menu. Choose GUI allows only one% instance to run (singleton).% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help jiemian% Last Modified by GUIDE v2.5 02-Jul-2012 01:22:36% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_Sta

52、te = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, jiemian_OpeningFcn, . gui_OutputFcn, jiemian_OutputFcn, . gui_LayoutFcn, , . gui_Callback, );if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);endif nargout varargout1:nargout = gui_mainfcn(gu

53、i_State, varargin:);else gui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% Executes just before jiemian is made visible.function jiemian_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdat

54、a reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to jiemian (see VARARGIN)% Choose default command line output for jiemianhandles.output = hObject;% Update handles structureguidata(hObject, handles);

55、% UIWAIT makes jiemian wait for user response (see UIRESUME)% uiwait(handles.figure1);% Outputs from this function are returned to the command line.function varargout = jiemian_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to

56、 figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout1 = handles.output;% Executes on button press in caiyang.function caiyang_Callback(hObject, eventdata,

57、 handles)% hObject handle to caiyang (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)x1,fs,bits=wavread(D:sound1.wav); %讀取信號(hào) x1=x1(:,1); %只取單聲道 sound(x1,fs); %對(duì)聲音的回放 axes(handles.axes1); plot(x1); %做原始語(yǔ)音信號(hào)的時(shí)域圖形 t

58、itle(原始語(yǔ)音信號(hào)的時(shí)域圖,fontname,宋體); y1=fft(x1); %做length(x1)點(diǎn)的FFT y1=fftshift(y1); %平移，是頻率中心為0 derta_fs = fs/length(x1); %設(shè)置頻譜的間隔，分辨率 axes(handles.axes2); plot(-fs/2:derta_fs:fs/2-derta_fs,abs(y1);%畫(huà)出原始語(yǔ)音信號(hào)的頻譜圖 title(原始語(yǔ)音信號(hào)的頻譜,fontname,宋體);% Executes on button press in yanshi.function yanshi_Callback(hObj

59、ect, eventdata, handles)% hObject handle to yanshi (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)x1,fs,bits=wavread(D:sound1.wav); x1=x1(:,1); %只取單聲道 z=zeros(2000,1);x1; %對(duì)聲音延時(shí)2000 sound(z,fs); axes(handles.axe

60、s3); plot(z); title(延時(shí)后的信號(hào)時(shí)域圖,fontname,宋體); %畫(huà)出延時(shí)后的信號(hào)時(shí)域圖 z=fft(x1); %做length(x1)點(diǎn)的FFT z=fftshift(z); %平移，是頻率中心為0 derta_fs = fs/length(z); %設(shè)置頻譜的間隔，分辨率 axes(handles.axes4); plot(-fs/2:derta_fs:fs/2-derta_fs,abs(z);%畫(huà)出延時(shí)后的信號(hào)頻譜圖 title(延時(shí)后的信號(hào)頻譜圖,fontname,宋體);% Executes on button press in hunxiang.functi