數(shù)字信號(hào)處理課程設(shè)計(jì)--用matlab對(duì)語(yǔ)音信號(hào)進(jìn)行分析及處理

1、<p><b>　　《數(shù)字信號(hào)處理》</b></p><p><b>　　課程設(shè)計(jì)報(bào)告</b></p><p>　　用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行分析及處理</p><p><b>　　摘要</b></p><p>　　MATLAB 語(yǔ)言是一種數(shù)據(jù)分析和處理功能十分強(qiáng)大

2、的計(jì)算機(jī)應(yīng)用軟件 ,它可以將聲音文件變換為離散的數(shù)據(jù)文件 , 然后利用其強(qiáng)大的矩陣運(yùn)算能力處理數(shù)據(jù) ,如數(shù)字濾波 、傅里葉變換 、時(shí)域和頻域分析 、聲音回放以及各種圖的呈現(xiàn)等 , 信號(hào)處理是MATLAB 重要應(yīng)用的領(lǐng)域之一 。 </p><p>　　本課程設(shè)計(jì)介紹了基于Matlab的對(duì)語(yǔ)音信號(hào)采集、處理及濾波器的設(shè)計(jì)，并使之實(shí)現(xiàn)的過(guò)程。理解與掌握課程中的基本概念、基本原理、基本分析方法，用Matlab進(jìn)行數(shù)字語(yǔ)音

3、信號(hào)處理，并闡述了課程設(shè)計(jì)的具體方法、步驟和內(nèi)容。綜合運(yùn)用本課程的理論知識(shí)進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過(guò)理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為工具進(jìn)行實(shí)現(xiàn)，從而復(fù)習(xí)鞏固課堂所學(xué)的理論知識(shí)，提高對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)對(duì)數(shù)字信號(hào)的處理。</p><p><b>　　目錄</b></p><p><b>　　一、設(shè)計(jì)目的：4<

4、/b></p><p><b>　　二、設(shè)計(jì)內(nèi)容：4</b></p><p><b>　　三、設(shè)計(jì)原理：4</b></p><p>　　1)．直接設(shè)計(jì)數(shù)字濾波器；5</p><p>　　2)．脈沖響應(yīng)不變法；6</p><p>　　3)．雙線性變換；6<

5、/p><p>　　4). 設(shè)計(jì)思路：7</p><p><b>　　四、設(shè)計(jì)步驟：7</b></p><p>　　1). 語(yǔ)音信號(hào)的采集7</p><p>　　2). 語(yǔ)音信號(hào)的頻譜分析8</p><p>　　3). 設(shè)計(jì)數(shù)字濾波器，畫(huà)出其頻率響應(yīng)曲線10</p><p&

6、gt;　　五、心得體會(huì)：13</p><p>　　六、參考文獻(xiàn)：14</p><p><b>　　一、設(shè)計(jì)目的：</b></p><p>　　1)．掌握數(shù)字信號(hào)處理的基本概念、基本理論和基本方法； </p><p>　　2)．掌握MATLAB設(shè)計(jì)FIR和IIR數(shù)字濾波器的方法； </p><p&g

7、t;　　3)．掌握在Windows環(huán)境下語(yǔ)音信號(hào)采集以及時(shí)域、頻域分析；</p><p>　　4)．學(xué)會(huì)MATLAB的使用，掌握MATLAB的程序設(shè)計(jì)方法；</p><p>　　5)．學(xué)會(huì)用MATLAB對(duì)信號(hào)進(jìn)行分析和處理。</p><p><b>　　二、設(shè)計(jì)內(nèi)容：</b></p><p>　　1)．錄制一段自己的語(yǔ)音

8、信號(hào),對(duì)錄制的信號(hào)進(jìn)行采樣；</p><p>　　2)．畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形和頻譜圖;</p><p>　　3)．給定濾波器的性能指標(biāo),采MATLAB設(shè)計(jì)數(shù)字濾波器,并畫(huà)出濾波器的頻率響應(yīng)；</p><p>　　4)．然后用自己設(shè)計(jì)的濾波器對(duì)采集的信號(hào)進(jìn)行濾波,畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜,并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比,分析信號(hào)的變化；</p>

9、<p>　　5)．回放語(yǔ)音信號(hào)。</p><p><b>　　三、設(shè)計(jì)原理：</b></p><p>　　濾波器的設(shè)計(jì)原理 設(shè)計(jì)數(shù)字濾波器的任務(wù)就是尋求一個(gè)因果穩(wěn)定的線性時(shí)不變系統(tǒng)，并使系統(tǒng)函數(shù)H(z)具有指定的頻率特性。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時(shí)域特性，可分為兩種，即無(wú)限長(zhǎng)沖激響應(yīng)（IIR）濾波器和有限長(zhǎng)沖激響應(yīng)（FIR）濾波器。IIR濾波器的特

10、征是，具有無(wú)限持續(xù)時(shí)間沖激響應(yīng)。這種濾波器一般需要用遞歸模型來(lái)實(shí)現(xiàn)，因而有時(shí)也稱(chēng)之為遞歸濾波器。FIR濾波器的沖激響應(yīng)只能延續(xù)一定時(shí)間，在工程實(shí)際中可以采用遞歸的方式實(shí)現(xiàn)，也可以采用非遞歸的方式實(shí)現(xiàn)。</p><p>　　濾波器的性能指標(biāo) 在設(shè)計(jì)一個(gè)濾波器之前，必須首先根據(jù)工程實(shí)際的需要確定濾波器的技術(shù)指標(biāo)。在很多實(shí)際應(yīng)用中，數(shù)字濾波器常常被用來(lái)實(shí)現(xiàn)選頻操作。因此，指標(biāo)的形式一般在頻域中給出幅度和相位響應(yīng)。幅度

11、指標(biāo)主要以兩種方式給出。第一種是絕對(duì)指標(biāo)。它提供對(duì)幅度響應(yīng)函數(shù)的要求，一般應(yīng)用于FIR濾波器的設(shè)計(jì)。第二種指標(biāo)是相對(duì)指標(biāo)。它以分貝值的形式給出要求。在工程實(shí)際中，這種指標(biāo)最受歡迎。對(duì)于相位響應(yīng)指標(biāo)形式，通常希望系統(tǒng)在通頻帶中人有線性相位。運(yùn)用線性相位響應(yīng)指標(biāo)進(jìn)行濾波器設(shè)計(jì)具有如下優(yōu)點(diǎn)：①只包含實(shí)數(shù)算法，不涉及復(fù)數(shù)運(yùn)算；②不存在延遲失真，只有固定數(shù)量的延遲；③長(zhǎng)度為N的濾波器（階數(shù)為N-1），計(jì)算量為N/2數(shù)量級(jí)。下面給出了有關(guān)濾波器的性

12、能指標(biāo) 1、 低通濾波器的性能指標(biāo)：Fp=1000Hz，F(xiàn)s=1200Hz，As=50dB，Ap=1dB； 2、 高通濾波器的性能指標(biāo)：Fp=4000Hz，F(xiàn)s=3500Hz，As=50dB，Ap=1dB； 3、 帶通濾波器的性能指標(biāo)：Fp1=1200Hz，F(xiàn)p2=3000HZ，F(xiàn)s1=1000Hz， Fs</p><p>　　1)．直接設(shè)計(jì)數(shù)字濾波器；&

13、lt;/p><p>　　數(shù)字濾波器是數(shù)字信號(hào)處理的一個(gè)重要組成部分。數(shù)字濾波實(shí)際上是一種運(yùn)算過(guò)程，其功能是將一組輸入的數(shù)字序列通過(guò)一定的運(yùn)算后轉(zhuǎn)變?yōu)榱硪唤M輸出的數(shù)字序列。</p><p>　　數(shù)字濾波器一般可以用兩種方法來(lái)實(shí)現(xiàn)：一種是根據(jù)描述數(shù)字濾波器的數(shù)學(xué)模型或信號(hào)流圖，用數(shù)字硬件裝配成一臺(tái)專(zhuān)門(mén)的設(shè)備，構(gòu)成專(zhuān)用的信號(hào)處理機(jī)；另一種方法就是直接利用通用計(jì)算機(jī)，將所需要的運(yùn)算編成程序讓計(jì)算機(jī)來(lái)執(zhí)

14、行，這也是用軟件來(lái)實(shí)現(xiàn)數(shù)字濾波器。</p><p>　　數(shù)字濾波器是離散時(shí)間系統(tǒng)，所處理的信號(hào)是離散時(shí)間信號(hào)。一般時(shí)域離散系統(tǒng)或網(wǎng)絡(luò)可以用差分方程、單位脈沖響應(yīng)以及系統(tǒng)函數(shù)進(jìn)行描述。如果輸入、輸出服從N階差分方程</p><p>　　則其系統(tǒng)函數(shù)，即濾波器的傳遞函數(shù)為</p><p>　　2)．脈沖響應(yīng)不變法；</p><p>　　利用模擬濾

15、波器來(lái)設(shè)計(jì)數(shù)字濾波器，就是從已知的模擬濾波器傳遞函數(shù)Ha(s)設(shè)計(jì)數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。</p><p>　　因此，它歸根結(jié)底是一個(gè)由S平面映射到Z平面的變換，這個(gè)變換通常是復(fù)變函數(shù)的映射變換，為了保證轉(zhuǎn)換后的H(z)穩(wěn)定且滿足技術(shù)要求，這個(gè)映射變換必須滿足以下兩條基本要求：</p><p>　　(1)因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。即 S平面的左半平面必須映

16、射到Z平面單位圓的內(nèi)部。</p><p>　　(2)數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻 響，s平面的虛軸映射為z平面的單位圓，相應(yīng)的頻率之間成線性關(guān)系。</p><p><b>　　核心原理：</b></p><p>　　通過(guò)對(duì)連續(xù)函數(shù)ha(t)等間隔采樣得到離散序列ha(nT)。令h(n)= ha(nT) ，T為采樣間隔。它是一種時(shí)域上

17、的轉(zhuǎn)換方法。</p><p><b>　　轉(zhuǎn)換步驟：</b></p><p><b>　　3)．雙線性變換；</b></p><p>　　脈沖響應(yīng)不變法的主要缺點(diǎn)是頻譜交疊產(chǎn)生的混淆，這是從S平面到Z平面的標(biāo)準(zhǔn)變換的多值對(duì)應(yīng)關(guān)系導(dǎo)致的,為了克服這一缺點(diǎn)，設(shè)想變換分為兩步：</p><p>　　第一步

18、：將整個(gè)S平面壓縮到S1平面的一條橫帶里；</p><p>　　第二步：通過(guò)標(biāo)準(zhǔn)變換關(guān)系將此橫帶變換到整個(gè)Z平面上去。由此建立S平面與Z平面一一對(duì)應(yīng)的單值關(guān)系，消除多值性，也就消除了混淆現(xiàn)象。</p><p>　　與脈沖響應(yīng)不變法相比，雙線性變換的主要優(yōu)點(diǎn)：靠頻率的嚴(yán)重非線性關(guān)系得到S平面與Z平面的單值一一對(duì)應(yīng)關(guān)系，整個(gè)jΩ軸單值對(duì)應(yīng)于單位圓一周。在零頻率附近，Ω～ω接近于線性關(guān)系，Ω進(jìn)一

19、步增加時(shí)，ω增長(zhǎng)變得緩慢， (ω終止于折疊頻率處)，所以雙線性變換不會(huì)出現(xiàn)由于高頻部分超過(guò)折疊頻率而混淆到低頻部分去的現(xiàn)象。</p><p>　　雙線性變換法的缺點(diǎn)：Ω與ω的非線性關(guān)系，導(dǎo)致數(shù)字濾波器的幅頻響應(yīng)相對(duì)于模擬濾波器的幅頻響應(yīng)有畸變，(使數(shù)字濾波器與模擬濾波器在響應(yīng)與頻率的對(duì)應(yīng)關(guān)系上發(fā)生畸變)。雙線性變換比脈沖響應(yīng)法的設(shè)計(jì)計(jì)算更直接和簡(jiǎn)單。由于s與z之間的簡(jiǎn)單代數(shù)關(guān)系，所以從模擬傳遞函數(shù)可直接通過(guò)代數(shù)置

20、換得到數(shù)字濾波器的傳遞函數(shù)。</p><p><b>　　4). 設(shè)計(jì)思路：</b></p><p>　　首先將數(shù)字高通濾波器的技術(shù)指標(biāo)轉(zhuǎn)化為模擬濾波器的技術(shù)指標(biāo)，再采用Buttetworth逼近的方法求得其系統(tǒng)函數(shù)H(S)，再利用沖激響應(yīng)不變法將模擬濾波器的系統(tǒng)函數(shù)H(s)變?yōu)橄鄳?yīng)的數(shù)字濾波器的系統(tǒng)函數(shù)H(Z)，然后用MATALB進(jìn)行仿真，然后確定數(shù)字濾波器的結(jié)構(gòu)。

21、</p><p><b>　　四、設(shè)計(jì)步驟：</b></p><p>　　1). 語(yǔ)音信號(hào)的采集</p><p>　　采集一小段語(yǔ)音信號(hào)保存為wav音頻文件格式，放到c目錄下待用。并且在Matlab軟件平臺(tái)下,利用函數(shù)wavread對(duì)語(yǔ)音信號(hào)進(jìn)行采樣,記住采樣頻率和采樣點(diǎn)數(shù)。</p><p>　　圖1 語(yǔ)音采集流程圖&l

22、t;/p><p><b>　　程序?yàn)椋?lt;/b></p><p>　　[y,fs,bite]=wavread('7.wav');</p><p><b>　　sound(y);</b></p><p>　　wavrecord(2*fs,fs);</p><p>&l

23、t;b>　　fft(y);</b></p><p><b>　　Y=fft(y);</b></p><p><b>　　fp=1000;</b></p><p>　　fr=1200;as=40;ap=1;</p><p>　　[N, Wn] = BUTTORD(2*fp/fs, 2

24、*fr/fs, ap, as);</p><p>　　[B,A] = BUTTER(N,Wn);</p><p>　　2). 語(yǔ)音信號(hào)的頻譜分析</p><p>　　要求畫(huà)出語(yǔ)音信號(hào)的時(shí)域波形;然后對(duì)語(yǔ)音號(hào)進(jìn)行快速傅里葉變換,得到信號(hào)的頻譜特性。</p><p><b>　　程序?yàn)椋?lt;/b></p><

25、;p>　　[y,fs,bite]=wavread('7.wav');</p><p><b>　　sound(y);</b></p><p><b>　　Y=fft(y);</b></p><p>　　mag=abs(Y);</p><p>　　phase=angle(Y);&

26、lt;/p><p>　　phase=phase*180/pi;</p><p>　　dB=20*log10(abs(Y)+eps);</p><p>　　t=(0:length(y)-1)/fs;</p><p><b>　　figure(1)</b></p><p>　　subplot(221);p

27、lot(t,y);title('波形圖');xlabel('時(shí)間(s)');ylabel('幅度');</p><p>　　f=(1:length(Y))*fs/length(Y);</p><p>　　subplot(222);plot(f,mag);title('幅度譜');xlabel('頻率(Hz)')

28、;ylabel('頻譜幅度');</p><p>　　subplot(223);plot(f,phase);title('相位譜');xlabel('頻率(Hz)');ylabel('頻譜相位');</p><p>　　subplot(224);plot(f,dB);title('對(duì)數(shù)幅度譜');xlabel(

29、'頻率(Hz)');ylabel('頻譜幅度(dB)');</p><p><b>　　顯示為：</b></p><p>　　3). 設(shè)計(jì)數(shù)字濾波器，畫(huà)出其頻率響應(yīng)曲線</p><p>　　各濾波器的性能指標(biāo)：</p><p>　　1).低通濾波器性能指標(biāo) fp＝1000Hz, fs＝1

30、200Hz, As＝40dB, Ap＝1dB</p><p>　　2).高通濾波器性能指標(biāo)　fs＝4000Hz, fp＝4300Hz, As＝40dB, Ap＝1dB</p><p>　　3).帶通濾波器性能指標(biāo)　fp1＝1200 Hz, fp2＝3 000 Hz, fs1＝1000 Hz, fs2＝3200Hz, As＝40dB, Ap＝1dB。</p><p>

31、　?。l率變換法設(shè)計(jì)IIR濾波器）</p><p><b>　?、俚屯V波器</b></p><p><b>　　程序?yàn)椋?lt;/b></p><p><b>　　clc</b></p><p>　　[I,Fs,bits]=wavread('7.wav');%讀取語(yǔ)

32、音信號(hào)7.wav</p><p>　　fp=1000; %通帶截止頻率</p><p>　　fr=1200; %阻帶角頻率</p><p>　　ap=1; %通帶紋波</p><p>　　as=40;

33、 %阻帶衰減</p><p>　　[n,fn]= buttord(fp/(Fs/2),fr/(Fs/2),ap,as,'z'); %巴特沃斯濾波器</p><p>　　[b,a]=butter(n,fn); %巴特沃斯IIR濾波器</p><p>　　[h,w]=freqz(b,a); %濾波器

34、的頻率響應(yīng)</p><p>　　I1=filter(b,a,I); %1維數(shù)字濾波器</p><p>　　I2=fftfilt(b,I); %FFT的FIR濾波</p><p>　　subplot(321),plot(I);title('原始語(yǔ)音信號(hào)');</p><p>　　subp

35、lot(322);plot(w*Fs/(2*pi),abs(h));title('巴特沃斯低通濾波');</p><p>　　subplot(323);plot(I1);title('巴特沃斯IIR濾波后語(yǔ)音信號(hào)');</p><p>　　subplot(324);plot(I2);title('巴特沃斯FIR濾波后語(yǔ)音信號(hào)');</p

36、><p>　　subplot(325);plot([h,w]);title('濾波器頻率響應(yīng)');</p><p>　　sound(I1,Fs,bits); %重放語(yǔ)音</p><p><b>　　顯示為：</b></p><p><b>　?、诟咄V波器</b><

37、;/p><p><b>　　程序?yàn)椋?lt;/b></p><p><b>　　clc;</b></p><p>　　[I,Fs,bits]=wavread('7.wav'); %讀取語(yǔ)音信號(hào)7.wav</p><p>　　fp=4300; %通帶截止

38、頻率</p><p>　　fr=4000; %阻帶角頻率</p><p>　　ap=1; %通帶紋波</p><p>　　as=40; %阻帶衰減</p><p>　　[n,fn]=cheb1ord(fp/(Fs/2),

39、fr/(Fs/2),ap,as,'z');%切比雪夫I型濾波器</p><p>　　[b,a]=cheby1(n,ap,fn,'high'); %I型切比雪夫IIR濾波器</p><p>　　[h,w]=freqz(b,a); %濾波器的頻率響應(yīng)</p><p>　　I1=filter(b,a,I);

40、 %1維數(shù)字濾波器</p><p>　　I2=fftfilt(b,I); %FFT的FIR濾波</p><p>　　subplot(321),plot(I);title('原始語(yǔ)音信號(hào)');</p><p>　　subplot(322);plot(w*Fs/(2*pi),abs(h));title(&#

41、39;切比雪夫1高通濾波');</p><p>　　subplot(323);plot(I1);title('切比雪夫1IIR濾波后語(yǔ)音信號(hào)');</p><p>　　subplot(324);plot(I2);title('切比雪夫1FIR濾波后語(yǔ)音信號(hào)');</p><p>　　subplot(325);plot([h,w

42、]);title('濾波器的頻率響應(yīng)');</p><p>　　sound(I1,Fs,bits); %重放語(yǔ)音</p><p><b>　　顯示為：</b></p><p><b>　?、蹘V波器</b></p><p><b>　　程序?yàn)椋?l

43、t;/b></p><p>　　[I,Fs,bits]=wavread('7.wav');%讀取語(yǔ)音信號(hào)7.wav</p><p>　　fp1=1200; %通帶截止頻率</p><p>　　fp2=3000; %通帶截止頻率</p><p>　

44、　fp=[fp1,fp2]; %通帶截止頻率</p><p>　　fr1=1000; %阻帶角頻率</p><p>　　fr2=3200; %阻帶角頻率</p><p>　　fr=[fr1,fr2]; %阻帶角頻率</p>

45、;<p>　　ap=1; %通帶紋波</p><p>　　as=40; %阻帶衰減</p><p>　　[n,fn]=ellipord(fp/(Fs/2),fr/(Fs/2),ap,as,'z');%最小橢圓濾波器</p><p>　　[b,a]=

46、ellip(n,ap,as,fn); %橢圓濾波器</p><p>　　[h,w]=freqz(b,a); %波器的頻率響應(yīng)</p><p>　　I1=filter(b,a,I); %1維數(shù)字濾波器</p><p>　　I2=fftfilt(b,I); %FFT的FIR濾波</p><p>　　

47、subplot(321),plot(I);title('原始語(yǔ)音信號(hào)');</p><p>　　subplot(322);plot(w*Fs/(2*pi),abs(h));title('橢圓帶通濾波');</p><p>　　subplot(323);plot(I1);title('橢圓IIR濾波后語(yǔ)音信號(hào)');</p><

48、;p>　　subplot(324);plot(I2);title('橢圓FIR濾波后語(yǔ)音信號(hào)');</p><p>　　subplot(325);plot([h,w]);title('濾波器的頻率響應(yīng)');</p><p>　　sound(I1,Fs,bits); %重放語(yǔ)音</p><p><b>

49、　　顯示為：</b></p><p><b>　　五、心得體會(huì)：</b></p><p>　　以往都是通過(guò)課本來(lái)感性的認(rèn)知語(yǔ)音信號(hào)，通過(guò)本次的課程設(shè)計(jì)，讓我對(duì)語(yǔ)音信號(hào)有了一個(gè)較為實(shí)際的認(rèn)識(shí)。于此同時(shí)，讓我再次把數(shù)字信號(hào)處理及數(shù)字濾波器的設(shè)計(jì)方法重新進(jìn)行了復(fù)習(xí)和學(xué)習(xí)。而對(duì)Matlab，也再次讓我感受到了其功能的強(qiáng)大。最為重要的是，本次課程設(shè)計(jì)讓我重新審視了學(xué)

50、習(xí)的過(guò)程：只去做實(shí)驗(yàn)是不行的，首先還是要思考，遇到了問(wèn)題查書(shū)籍，百度搜索也只是一種手段，更加重要的是想，再理解，只有這樣才能真正的做好實(shí)驗(yàn)。</p><p><b>　　六、參考文獻(xiàn)：</b></p><p>　　1. 陳懷琛.MATLAB及在電子信息課程中的應(yīng)用.北京電子工業(yè)出版社.2008.1 </p><p>　　2. 張文.基于MATL

51、AB的語(yǔ)音信號(hào)的濾波域?qū)崿F(xiàn).山西電子技術(shù).2008.2 </p><p>　　3. 徐靖濤.基于MATLAB的語(yǔ)音信號(hào)分析與處理.重慶科技學(xué)院學(xué)報(bào).2008.1 </p><p>　　4. 程佩青. 數(shù)字信號(hào)處理教程第三版.清華大學(xué)出版社. 2007.2 </p><p>　　5. 陳后金. 數(shù)字信號(hào)處理第2版.高等教育出版社. 2008.11</p>