語(yǔ)音數(shù)字信號(hào)處理與分析及Matlab實(shí)現(xiàn)數(shù)字信號(hào)課程設(shè)計(jì)

1、吉林工程技術(shù)師范學(xué)院信息工程學(xué)院數(shù)字信號(hào)處理課程設(shè)計(jì)報(bào)告題 目：語(yǔ)音數(shù)字信號(hào)處理與分析及matlab實(shí)現(xiàn)專(zhuān) 業(yè)： 電子信息工程 班 級(jí)： 姓 名： 學(xué) 號(hào)： 指導(dǎo)教師： 時(shí) 間： 2013/11/4-2013/11/22 目 錄第一章 緒 論11.1 設(shè)計(jì)的目的及意義11.2 設(shè)計(jì)要求11.3 設(shè)計(jì)內(nèi)容21.3.1 理論依據(jù)及信號(hào)采集21.3.2 構(gòu)造受干擾信號(hào)并對(duì)其進(jìn)行fft頻譜分析21.3.3 數(shù)字濾波器的設(shè)計(jì)及信號(hào)處理2第二章 信號(hào)頻譜分析42.1原始信號(hào)及頻譜分析42.2加入干擾噪聲后的信號(hào)及頻譜分析5第三章 數(shù)字濾波器的設(shè)計(jì)73.1 高通濾波器的設(shè)計(jì)73.2 低通濾波器的設(shè)計(jì)73.3

2、 帶通濾波器的設(shè)計(jì)83.4 帶阻濾波器的設(shè)計(jì)9第四章 gui用戶(hù)圖形界面設(shè)計(jì)114.1 gui的基本設(shè)計(jì)流程114.2 數(shù)字濾波器的實(shí)現(xiàn)164.2.1 高通濾波器的實(shí)現(xiàn)164.2.2 低通濾波器的實(shí)現(xiàn)174.2.3 帶通濾波器的實(shí)現(xiàn)174.2.4 帶阻濾波器的實(shí)現(xiàn)18第五章 課程設(shè)計(jì)總結(jié)22參考文獻(xiàn)i附錄i 第一章 緒 論1.1 設(shè)計(jì)的目的及意義目的:通過(guò)數(shù)字信號(hào)處理課程設(shè)計(jì),使學(xué)生深入了解和掌握噪聲信號(hào)的頻譜分析,進(jìn)一步明確了數(shù)字濾波器設(shè)計(jì)方法和實(shí)際應(yīng)用,為下一步學(xué)習(xí)和工作打下良好的基礎(chǔ)。意義:語(yǔ)音信號(hào)處理是一門(mén)新興的邊緣學(xué)科，它和認(rèn)知科學(xué)、心理學(xué)、語(yǔ)言學(xué)、計(jì)算機(jī)科學(xué)、模式識(shí)別和人工智能等學(xué)

3、科有著緊密的聯(lián)系。語(yǔ)音信號(hào)處理的發(fā)展依賴(lài)于這些學(xué)科的發(fā)展，而語(yǔ)音信號(hào)處理技術(shù)的進(jìn)步也會(huì)促進(jìn)這些領(lǐng)域的進(jìn)步。通過(guò)matlab可以對(duì)數(shù)字化的語(yǔ)音信號(hào)進(jìn)行處理。由于matlab是一種面向科學(xué)和工程計(jì)算的高級(jí)語(yǔ)言，允許用數(shù)字形式的語(yǔ)言編程，又有大量的庫(kù)函數(shù)，所以編程簡(jiǎn)單、編程效率高、易學(xué)易懂。我們可以對(duì)信號(hào)進(jìn)行加噪和去噪、濾波、截取語(yǔ)音等。1.2 設(shè)計(jì)要求 要求利用matlab對(duì)語(yǔ)音信號(hào)進(jìn)行數(shù)字信號(hào)處理和分析，要求學(xué)生采集語(yǔ)音信號(hào)后，在matlab軟件平臺(tái)進(jìn)行頻譜分析；并對(duì)所采集語(yǔ)音信號(hào)加入干擾噪聲，對(duì)加入噪聲的信號(hào)進(jìn)行頻譜分析，設(shè)計(jì)合適的濾波器濾除噪聲，恢復(fù)原信號(hào)。 1、對(duì)采集的語(yǔ)音信號(hào)進(jìn)行頻譜分析

4、。 2、對(duì)原采集的語(yǔ)音信號(hào)加入干擾噪聲并對(duì)其進(jìn)行頻譜分析。 3、設(shè)計(jì)合適的濾波器。1.3 設(shè)計(jì)內(nèi)容1.3.1 理論依據(jù)及信號(hào)采集根據(jù)設(shè)計(jì)要求分析系統(tǒng)功能，掌握設(shè)計(jì)中所需理論（采樣頻率、采樣定理；時(shí)域信號(hào)的fft分析；數(shù)字濾波器設(shè)計(jì)原理和方法，各種不同類(lèi)型濾波器的性能比較），闡明設(shè)計(jì)原理。采集語(yǔ)音信號(hào)，把語(yǔ)音信號(hào)保存為.wav文件，長(zhǎng)度小于3秒，并對(duì)語(yǔ)音信號(hào)進(jìn)行采樣；錄制的軟件可以使用windows自帶的錄音機(jī)，或者也可以使用其他專(zhuān)業(yè)的錄音軟件，錄制時(shí)需要配備錄音硬件，為了方面比較，需要在安靜、無(wú)噪音、干擾小的環(huán)境下錄，并對(duì)其進(jìn)行fft頻譜分析，畫(huà)出信號(hào)的時(shí)域波形圖和頻譜圖。1.3.2 構(gòu)造受干

5、擾信號(hào)并對(duì)其進(jìn)行fft頻譜分析對(duì)所采集的語(yǔ)音信號(hào)加入干擾噪聲，對(duì)語(yǔ)音信號(hào)進(jìn)行回放，感覺(jué)加噪前后聲音的變化，分析原因，得出結(jié)論。并對(duì)其進(jìn)行fft頻譜分析，比較加噪前后語(yǔ)音信號(hào)的波形及頻譜，對(duì)所得結(jié)果進(jìn)行分析，闡明原因，得出結(jié)論。1.3.3 數(shù)字濾波器的設(shè)計(jì)及信號(hào)處理 根據(jù)待處理信號(hào)特點(diǎn)，設(shè)計(jì)合適數(shù)字濾波器，繪制所設(shè)計(jì)濾波器的幅頻和相頻特性。 用所設(shè)計(jì)的濾波器對(duì)含噪語(yǔ)音信號(hào)進(jìn)行濾波。對(duì)鋁箔后的語(yǔ)音信號(hào)進(jìn)行fft頻譜分析。畫(huà)出處理過(guò)程中所得各種波形及頻譜圖。對(duì)語(yǔ)音信號(hào)進(jìn)行回放，感覺(jué)濾波前后聲音的變化。比較濾波前后語(yǔ)音信號(hào)的波形及頻譜，對(duì)所得結(jié)果的濾波器性能進(jìn)行頻譜分析，闡明原因，得出結(jié)論。設(shè)計(jì)處理系

6、統(tǒng)的用戶(hù)界面，在所設(shè)計(jì)界面上可以選擇濾波器的參數(shù)，顯示濾波器的頻率響應(yīng)，選擇信號(hào)等。 第二章 信號(hào)頻譜分析2.1原始信號(hào)及頻譜分析下面的一段程序是語(yǔ)音信號(hào)在matlab中的最簡(jiǎn)單表現(xiàn)，它實(shí)現(xiàn)了語(yǔ)音的讀入打開(kāi)，以及繪出了語(yǔ)音信號(hào)的波形頻譜圖。編程詳見(jiàn)附錄結(jié)果如圖2-1-1 圖2-1-1 原始信號(hào)的圖形已知一個(gè)語(yǔ)音信號(hào)，數(shù)據(jù)采樣頻率為100hz，試分別繪制n=128點(diǎn)dft的幅頻圖和n=1024點(diǎn)dft幅頻圖。編程詳見(jiàn)附錄結(jié)果如圖2-1-2 圖2-1-2 n=128點(diǎn)dft的幅頻圖和n=1024點(diǎn)dft幅頻圖 上圖（a）、(b)為n=128點(diǎn)幅頻譜圖，(c)、(d)為n=1024點(diǎn)幅頻譜圖。由于采

7、樣頻率f=100hz，故nyquist頻率為50hz，（a）、(c)是0100hz頻譜圖，(b) 、(d)是050hz頻譜圖。由（a）或(c)可見(jiàn)，整個(gè)頻譜圖是以nyquist頻率為軸對(duì)稱(chēng)的。因此利用fft對(duì)信號(hào)做頻譜分析，只要考察0nyquist頻率（采樣頻率一半）范圍的幅頻特性。比較（a）和(c)或(b)和 (d)可見(jiàn)，幅值大小與fft選用點(diǎn)數(shù)n有關(guān)，但只要點(diǎn)數(shù)n足夠不影響研究結(jié)果。從上圖幅頻譜可見(jiàn)，信號(hào)中包括15hz和40hz的正弦分量。2.2加入干擾噪聲后的信號(hào)及頻譜分析若信號(hào)長(zhǎng)度t=25.6s，即抽樣后x(n)點(diǎn)數(shù)為t/ts=256，所得頻率分辨率為hz，以此觀察數(shù)據(jù)長(zhǎng)度n的變化對(duì)d

8、tft分辨率的影響：編程詳見(jiàn)附錄結(jié)果如圖2-2-1 圖2-2-1 加入干擾噪聲后的信號(hào)及頻譜分析第三章 數(shù)字濾波器的設(shè)計(jì)3.1 高通濾波器的設(shè)計(jì) n階高通濾波器的設(shè)計(jì)，高通濾波器是一個(gè)使高頻率比較容易通過(guò)而阻止低頻率通過(guò)的系統(tǒng)。它去掉了信號(hào)中不必要的低頻成分或者說(shuō)去掉了低頻干擾。其特性在時(shí)域及頻域中可分別用沖激響應(yīng)及頻率響應(yīng)描述。后者是用以頻率為自變量的函數(shù)表示，一般情況下它是一個(gè)以復(fù)變量為自變量的的復(fù)變函數(shù)，以表示。它的模和幅角為角頻率的函數(shù)，分別稱(chēng)為系統(tǒng)的“幅頻響應(yīng)”和“相頻響應(yīng)”，它分別代表激勵(lì)源中不同頻率的信號(hào)成分通過(guò)該系統(tǒng)時(shí)所遇到的幅度變化和相位變化?？梢宰C明，系統(tǒng)的“頻率響應(yīng)”就是

9、該系統(tǒng)“沖激響應(yīng)”的傅里葉變換。當(dāng)線(xiàn)性無(wú)源系統(tǒng)可以用一個(gè)n階線(xiàn)性微分方程表示時(shí)，頻率響應(yīng)為一個(gè)有理分式，它的分子和分母分別與微分方程的右邊和左邊相對(duì)應(yīng)。3.2 低通濾波器的設(shè)計(jì) n階低通濾波器的設(shè)計(jì)，低通濾波器是容許低于截止頻率的信號(hào)通過(guò)， 但高于截止頻率的信號(hào)不能通過(guò)的電子濾波裝置。對(duì)于不同濾波器而言，每個(gè)頻率的信號(hào)的強(qiáng)弱程度不同。當(dāng)使用在音頻應(yīng)用時(shí)，它有時(shí)被稱(chēng)為高頻剪切濾波器, 或高音消除濾波器。低通濾波器概念有許多不同的形式，其中包括電子線(xiàn)路（如音頻設(shè)備中使用的hiss 濾波器）、平滑數(shù)據(jù)的數(shù)字算法、音障（acoustic barriers）、圖像模糊處理等等，這兩個(gè)工具都通過(guò)剔除短期波

10、動(dòng)、保留長(zhǎng)期發(fā)展趨勢(shì)提供了信號(hào)的平滑形式。低通濾波器在信號(hào)處理中的作用等同于其它領(lǐng)域如金融領(lǐng)域中移動(dòng)平均數(shù)（moving average)所起的作用；低通濾波器有很多種，其中，最通用的就是巴特沃斯濾波器和切比雪夫?yàn)V波器。低通濾波器允許從直流到某個(gè)截止頻率 的信號(hào)通過(guò)。將通用濾波器二階傳遞函數(shù)的高通和帶通系數(shù)均設(shè)為零，即得到一個(gè)二階低通濾波器傳遞公式。對(duì)于高于的頻率，信號(hào)按該頻率平方的速率下降。在頻率處，阻尼值使輸出信號(hào)衰減。您可以級(jí)聯(lián)多個(gè)這樣的濾波器部分來(lái)得到一個(gè)更高階的（更陡峭的轉(zhuǎn)降）濾波器。假定設(shè)計(jì)要求一個(gè)截止頻率為10khz的四階貝塞爾(bessel) 低通濾波器。根據(jù)參考文獻(xiàn)1，每部分

11、的轉(zhuǎn)降頻率分別為16.13及18.19 khz，阻尼值分別為1.775及0.821，并且這兩個(gè)濾波器分區(qū)的高通、帶通和低通系數(shù)分別為0、0與1。可以使用這兩個(gè)帶有上述參數(shù)的濾波器部分來(lái)實(shí)現(xiàn)所要求的濾波器。截止頻率為輸出信號(hào)衰減3 db的頻率點(diǎn)。3.3 帶通濾波器的設(shè)計(jì) 2n階帶通濾波器設(shè)計(jì)，帶通濾波器（band-pass filter）是一個(gè)允許特定頻段的波通過(guò)同時(shí)屏蔽其他頻段的設(shè)備。比如rlc振蕩回路就是一個(gè)模擬帶通濾波器。帶通濾波器是指能通過(guò)某一頻率范圍內(nèi)的頻率分量、但將其他范圍的頻率分量衰減到極低水平的濾波器，與帶阻濾波器的概念相對(duì)。一個(gè)模擬帶通濾波器的例子是電阻-電感-電容電路(rlc

12、 circuit)。這些濾波器也可以用低通濾波器同高通濾波器組合來(lái)產(chǎn)生。一個(gè)理想的帶通濾波器應(yīng)該有一個(gè)完全平坦的通帶，在通帶內(nèi)沒(méi)有放大或者衰減， 有源帶通濾波器電路并且在通帶之外所有頻率都被完全衰減掉，另外，通帶外的轉(zhuǎn)換在極小的頻率范圍完成。實(shí)際上，并不存在理想的帶通濾波器。濾波器并不能夠?qū)⑵谕l率范圍外的所有頻率完全衰減掉，尤其是在所要的通帶外還有一個(gè)被衰減但是沒(méi)有被隔離的范圍。這通常稱(chēng)為濾波器的滾降現(xiàn)象，并且使用每十倍頻的衰減幅度的db數(shù)來(lái)表示。通常，濾波器的設(shè)計(jì)盡量保證滾降范圍越窄越好，這樣濾波器的性能就與設(shè)計(jì)更加接近。然而，隨著滾降范圍越來(lái)越小，通帶就變得不再平坦，開(kāi)始出現(xiàn)“波紋”。這

13、種現(xiàn)象在通帶的邊緣處尤其明顯，這種效應(yīng)稱(chēng)為吉布斯現(xiàn)象。除了電子學(xué)和信號(hào)處理領(lǐng)域之外，帶通濾波器應(yīng)用的一個(gè)例子是在大氣科學(xué)領(lǐng)域，很常見(jiàn)的例子是使用帶通濾波器過(guò)濾最近3到10天時(shí)間范圍內(nèi)的天氣數(shù)據(jù)，這樣在數(shù)據(jù)域中就只保留了作為擾動(dòng)的氣旋。在頻帶較低的剪切頻率和較高的剪切頻率之間是共振頻率，這里濾波器的增益最大，濾波器的帶寬就是和之間的差值。3.4 帶阻濾波器的設(shè)計(jì) 2n階帶阻濾波器的設(shè)計(jì)，阻濾波器是指能通過(guò)大多數(shù)頻率分量、但將某些范圍的頻率分量衰減到極低水平的濾波器，與帶通濾波器的概念相對(duì)。其中點(diǎn)阻濾波器（notch filter）是一種特殊的帶阻濾波器，它的阻帶范圍極小，有著很高的q值（q fa

14、ctor）。帶阻濾波器是指能通過(guò)大多數(shù)頻率分量、但將某些范圍的頻率分量衰減到極低水平的濾波器，與帶通濾波器說(shuō)明的概念相對(duì)。其中點(diǎn)阻濾波器（notch filter）是一種特殊的帶阻濾波器，它的阻帶范圍極小，有著很高的q值（q factor）。 將輸入電壓同時(shí)作用于低通濾波器和高通濾波器，再將兩個(gè)電路的輸出電壓求和，就可以得到帶阻濾波器，如下圖所示。其中低通濾波器的截止頻率 應(yīng)小于高通濾波器的截止頻率 ，因此，電路的阻帶為( - )。帶阻濾波器分為腔體帶阻濾波器、lc帶阻濾波器和有源帶阻濾波電路實(shí)用電路常利用無(wú)源低通濾波器和高通濾波器并聯(lián)構(gòu)成無(wú)源帶阻濾波電路，然后接同相比例運(yùn)算電路，從而得到有源

15、帶阻濾波電路。 以上幾種濾波，我們都可以從信號(hào)濾波前后的波形圖以及頻譜圖上看出變化。當(dāng)然，也可以用sound()函數(shù)來(lái)播放濾波后的語(yǔ)言，從聽(tīng)覺(jué)上直接感受語(yǔ)音信號(hào)的變化，但由于人耳聽(tīng)力的限制，有些情況下我們是很難聽(tīng)出異同的。同樣，通過(guò)函數(shù)的調(diào)用，也可以將信號(hào)的頻譜進(jìn)行“分離觀察”，如顯出信號(hào)的幅值或相位。下面，通過(guò)改變系統(tǒng)函數(shù)的分子與分母系數(shù)比，來(lái)觀察信號(hào)濾波前后的幅值與相位。并且使結(jié)果更加明顯，使人耳得以很容易的辨聽(tīng)。第四章 gui用戶(hù)圖形界面設(shè)計(jì)4.1 gui的基本設(shè)計(jì)流程打開(kāi)matlab軟件,文件-新建-gui,如圖4-1-1,即可出現(xiàn)圖4-1-2 圖4-1-1 新建gui 圖4-1-2

16、guide quick start點(diǎn)擊ok，即可出現(xiàn)圖4-1-3 gui界面 圖4-1-3 gui界面點(diǎn)擊filepreferences,如圖4-1-4，即可出現(xiàn)圖4-1-5 圖4-1-4 改變圖形界面將show names in component palette 勾選即可出現(xiàn)出現(xiàn)圖4-1-5，點(diǎn)擊ok，即可出現(xiàn)圖4-1-6 圖4-1-5 gui preferences 圖4-1-6 gui繪制界面用push button按鍵及axes按鍵繪制如圖4-1-7圖 圖4-1-7 整體界面右擊push button按鍵，選擇view callbackscallback，如圖4-1-8 圖4-1-8

17、 設(shè)置回調(diào)函數(shù)保存，如圖4-1-9 圖4-1-9 保存文件將程序輸入進(jìn)去，如圖4-1-10所示 圖4-1-10 輸入回調(diào)函數(shù)保存并且運(yùn)行即可出現(xiàn)圖形4.2 數(shù)字濾波器的實(shí)現(xiàn)4.2.1 高通濾波器的實(shí)現(xiàn)高通濾波器是一個(gè)使高頻率比較容易通過(guò)而阻止低頻率通過(guò)的系統(tǒng)。它去掉了信號(hào)中不必要的低頻成分或者說(shuō)去掉了低頻干擾。其特性在時(shí)域及頻域中可分別用沖激響應(yīng)及頻率響應(yīng)描述。一般情況下它是一個(gè)以復(fù)變量為自變量的的復(fù)變函數(shù)，以表示。它的模和幅角為角頻率的函數(shù)，分別稱(chēng)為系統(tǒng)的“幅頻響應(yīng)”和“相頻響應(yīng)”，gui界面成圖如圖4-2-1 圖4-2-1 高通濾波器圖形4.2.2 低通濾波器的實(shí)現(xiàn)低通濾波器是容許低于截止

18、頻率的信號(hào)通過(guò)， 但高于截止頻率的信號(hào)不能通過(guò)的電子濾波裝置對(duì)于不同濾波器而言，每個(gè)頻率的信號(hào)的強(qiáng)弱程度不同。當(dāng)使用在音頻應(yīng)用時(shí)，它有時(shí)被稱(chēng)為高頻剪切濾波器, 或高音消除濾波器。gui界面成圖如圖4-2-2 圖4-2-2 低通濾波器圖形4.2.3 帶通濾波器的實(shí)現(xiàn)帶通濾波器（band-pass filter）是一個(gè)允許特定頻段的波通過(guò)同時(shí)屏蔽其他頻段的設(shè)備。比如rlc振蕩回路就是一個(gè)模擬帶通濾波器。帶通濾波器是指能通過(guò)某一頻率范圍內(nèi)的頻率分量、但將其他范圍的頻率分量衰減到極低水平的濾波器，與帶阻濾波器的概念相對(duì)。一個(gè)模擬帶通濾波器的例子是電阻-電感-電容電路(rlc circuit)。這些濾波

19、器也可以用低通濾波器同高通濾波器組合來(lái)產(chǎn)生。gui界面成圖如圖4-2-3 圖4-2-3 帶通濾波器圖形4.2.4 帶阻濾波器的實(shí)現(xiàn)阻濾波器是指能通過(guò)大多數(shù)頻率分量、但將某些范圍的頻率分量衰減到極低水平的濾波器，與帶通濾波器的概念相對(duì)。因此，電路的阻帶為( - )。帶阻濾波器分為腔體帶阻濾波器、lc帶阻濾波器和有源帶阻濾波電路，實(shí)用電路常利用無(wú)源低通濾波器和高通濾波器并聯(lián)構(gòu)成無(wú)源帶阻濾波電路，然后接同相比例運(yùn)算電路，從而得到有源帶阻濾波電路。 圖4-2-4 帶阻濾波器圖形 編程詳見(jiàn)附錄結(jié)果如圖4-2-5 圖4-2-5 改變系數(shù)比之后的圖形 濾波后，編程詳見(jiàn)附錄結(jié)果如圖4-2-6 圖4-2-6 濾

20、波后的圖形可以聽(tīng)到聲音明顯變得高亢了。從上面的波形與幅值（即幅頻）圖，也可看出，濾波后的幅值變成了濾波前的20倍。編程詳見(jiàn)附錄 得圖4-2-7 圖4-2-7 濾波后相位譜圖形可以看到相位譜沒(méi)什么變化。 第五章 課程設(shè)計(jì)總結(jié)通過(guò)本學(xué)期對(duì)語(yǔ)音信號(hào)處理這門(mén)課程的學(xué)習(xí)，是我對(duì)信號(hào)處理有了更深的了解，同時(shí)極大的提高了我對(duì)語(yǔ)音信號(hào)的興趣。課堂學(xué)習(xí)主要注重于理論知識(shí)，而我們要將所學(xué)知識(shí)應(yīng)用于實(shí)際，在此階段，課程設(shè)計(jì)便是最好的選擇了，通過(guò)課程設(shè)計(jì)，我們可以溫習(xí)我們所學(xué)的理論知識(shí)，同時(shí)為將理論知識(shí)運(yùn)用于實(shí)際搭建了一個(gè)很好的平臺(tái)，不僅如此，通過(guò)這次的課程設(shè)計(jì)，使我知道了在當(dāng)今的信息技術(shù)如此發(fā)達(dá)的世界中，我們必須運(yùn)

21、用多種渠道，去學(xué)習(xí)研究。并要很好的運(yùn)用計(jì)算機(jī)和一些軟件，只有這樣，我們才能更好地、精確地、快速地解決問(wèn)題。還有就是提高了自主解決問(wèn)題的能力。 這次的課程設(shè)計(jì)也不例外。雖然我們以前學(xué)了一點(diǎn)matlab，不少課本上也有提到過(guò) matlab，但是我們還是不怎么熟練，通過(guò)這次自動(dòng)控制原理課程讓我更好地學(xué)會(huì)了如何去使用這個(gè)軟件。自動(dòng)控制原理的知識(shí)在課堂上掌握的也不是很好，所以這次課程設(shè)計(jì)對(duì)于我來(lái)說(shuō)真是難度不小啊。在很多人眼中為期兩周的課程設(shè)計(jì)或許是一種煎熬，這是可以理解的，在這兩周當(dāng)中，我們不僅要完成這個(gè)課程設(shè)計(jì)，而且還要學(xué)習(xí)其他專(zhuān)業(yè)課。 每一次課程設(shè)計(jì)都會(huì)學(xué)到不少東西，這次當(dāng)然也不例外。不但對(duì)自動(dòng)控制

22、原理的知識(shí)鞏固了，也加強(qiáng)了matlab 這個(gè)強(qiáng)大軟件使用的學(xué)習(xí)，這次課程設(shè)計(jì)終于順利完成了，在設(shè)計(jì)中遇到了很多編程問(wèn)題，最后在自己和同學(xué)相互協(xié)助下，終于迎刃而解了。課程設(shè)計(jì)是我們專(zhuān)業(yè)課程知識(shí)綜合應(yīng)用的實(shí)踐訓(xùn)練，著是我們邁向社會(huì)，從事職業(yè)工作前一個(gè)必不少的過(guò)程”千里之行始于足下”，通過(guò)這次課程設(shè)計(jì)，我深深體會(huì)到這句千古名言的真正含義我今天認(rèn)真的進(jìn)行課程設(shè)計(jì)，學(xué)會(huì)腳踏實(shí)地邁開(kāi)這一步，就是為明天能穩(wěn)健地在社會(huì)大潮中奔跑打下堅(jiān)實(shí)的基礎(chǔ)同樣十分感謝xxx、xxx老師在此次實(shí)驗(yàn)中給我提出的建議與不足，老師嚴(yán)謹(jǐn)細(xì)致、一絲不茍的作風(fēng)一直是我工作、學(xué)習(xí)中的榜樣；老師循循善誘的教導(dǎo)和不拘一格的思路給予我無(wú)盡的啟迪

23、；這次模具設(shè)計(jì)的每個(gè)實(shí)驗(yàn)細(xì)節(jié)和每個(gè)數(shù)據(jù)，都離不開(kāi)老師您的細(xì)心指導(dǎo)。而您開(kāi)朗的個(gè)性和寬容的態(tài)度，幫助我能夠很順利的完成了這次課程設(shè)計(jì)。24參考文獻(xiàn)1 趙力.語(yǔ)音信號(hào)處理.北京：機(jī)械工業(yè)出版社，2003.2 胡航.語(yǔ)音信號(hào)處理.哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2000.3 姚天任.數(shù)字語(yǔ)音處理.武漢：華中科技大學(xué)出版社，2003.4 楊行峻,遲惠生等.語(yǔ)音信號(hào)處理.北京：電子工業(yè)出版社，19955 劉衛(wèi)國(guó).matlab程序設(shè)計(jì)教程.中國(guó)水利水電出版社，20056 朱衡君.matlab語(yǔ)言及實(shí)踐教程(第二版).清華大學(xué)出版社，2009.87 張靜.matlab語(yǔ)言的應(yīng)用.電子工業(yè)出版社，20078 蔡

24、蓮紅,黃德智,蔡銳.語(yǔ)音技術(shù)基礎(chǔ)與應(yīng)用.北京：清華大學(xué)出版社，2003.9 王炳錫,屈單,彭煊.實(shí)用語(yǔ)音識(shí)別基礎(chǔ).北京：機(jī)械工業(yè)出版社，2005.10 易克初,田斌,付強(qiáng).語(yǔ)音信號(hào)處理.北京：國(guó)防工業(yè)出版社,2000.11 韓紀(jì)慶,張磊,鄭鐵然.語(yǔ)音信號(hào)處理.北京：清華大學(xué)出版社,2004.i附錄x,fs,bits=wavread(jing.wav,1024 5120); sound(x,fs,bits); x=fft(x,4096); magx=abs(x); angx=angle(x); subplot(221);plot(x);title(原始信號(hào)波形); subplot(222);pl

25、ot(x);title(原始信號(hào)頻譜); subplot(223);plot(magx);title(原始信號(hào)幅值); subplot(224);plot(angx);title(原始信號(hào)相位);x=wavread(jing.wav); sound(x);fs=100;n=128;y=fft(x,n);magy=abs(y); f=(0:length(y)-1)*fs/length(y); subplot(221);plot(f,magy); xlabel(頻率(hz); ylabel(幅值); title(n=128(a);grid subplot(222);plot(f(1:n/2),ma

26、gy(1:n/2); xlabel(頻率(hz);ylabel(幅值); title(n=128(b);grid fs=100;n=1024; y=fft(x,n); magy=abs(y); f=(0:length(y)-1)*fs/length(y); subplot(223);plot(f,magy); xlabel(頻率(hz);ylabel(幅值); title(n=1024(c);grid subplot(224);plot(f(1:n/2),magy(1:n/2); xlabel(頻率(hz);ylabel(幅值); title(n=1024(d);gridx,fs,bits=w

27、avread(jing.wav); n=256; f=0:fs/n:fs/2-1/n; x=fft(x);x=abs(x); subplot(211) plot(f(45:60),x(45:60);grid xlabel(hz),ylabel(|h(ejw)|) n=n*4; f=0:fs/n:fs/2-1/n; x=fft(x);x=abs(x); subplot(212) plot(f(45*4:4*60),x(4*45:4*60);grid xlabel(hz),ylabel(|h(ejw)|)x=wavread(jing.wav); sound(x);n=5;wc=0.3;b,a=bu

28、tter(n,wc,high); x=fft(x); subplot(321);plot(x);title(濾波前信號(hào)的波形); subplot(322);plot(x);title(濾波前信號(hào)的頻譜); y=filter(b,a,x);y=fft(y); subplot(323);plot(y);title(iir 濾波后信號(hào)的波形); subplot(324);plot(y);title(iir 濾波后信號(hào)的頻譜);z=fftfilt(b,x); z=fft(z); subplot(325);plot(z);title(fir 濾波后信號(hào)的波形); subplot(326);plot(z);title(fir 濾波后信號(hào)的頻譜);x=wavread(jing.wav);sound(x);n=5;wc=0.3,0.6; b,a=butter(n,wc);x=fft(x); subplot(321);plot(x);title(濾波前信號(hào)的波形); subplot(322);plot(x