語(yǔ)音信號(hào)處理與濾波課程設(shè)計(jì)

1、<p><b>　　《數(shù)字信號(hào)處理》</b></p><p><b>　　課程設(shè)計(jì)報(bào)告</b></p><p>　　題 目： 語(yǔ)音信號(hào)處理與濾波 </p><p>　　專 業(yè)： 信息與計(jì)算科學(xué) </p><p>　　二〇一〇 年 一 月 八 日</

2、p><p><b>　　一、課程設(shè)計(jì)目的</b></p><p>　　綜合運(yùn)用數(shù)字信號(hào)處理的理論知識(shí)進(jìn)行頻譜分析和濾波器設(shè)計(jì)，通過(guò)理論推導(dǎo)得出相應(yīng)結(jié)論，再利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)，從而加深對(duì)所學(xué)知識(shí)的理解。</p><p>　　二、課程設(shè)計(jì)基本要求及分組</p><p>　　1、對(duì)知識(shí)點(diǎn)的掌握要求：<

3、/p><p>　?。?）熟悉離散信號(hào)和系統(tǒng)的時(shí)域特性。</p><p>　?。?）掌握序列快速傅里葉變換FFT方法。</p><p>　?。?）學(xué)會(huì)MATLAB的使用，掌握MATLAB的程序設(shè)計(jì)方法。</p><p>　　（4）利用MATLAB對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。</p><p>　　（5）掌握MATLAB設(shè)計(jì)FIR和

4、IIR數(shù)字濾波器的方法。</p><p><b>　　2、分組情況：</b></p><p><b>　　組長(zhǎng)：張鳳</b></p><p>　　組員：張鳳、張志廣、李云、黃生涯、曾勝、沈呈潔</p><p><b>　　分工情況：</b></p><p&g

5、t;　　張鳳：設(shè)計(jì)全過(guò)程的監(jiān)督及協(xié)助、部分源程序代碼的編寫和整個(gè)源程序代碼的整理。</p><p>　　張志廣：語(yǔ)音信號(hào)的錄制以及對(duì)語(yǔ)音信號(hào)的采樣，F(xiàn)IR低通濾波器的設(shè)計(jì)。</p><p>　　李云：FIR高通濾波器、FIR帶通濾波器的設(shè)計(jì)以及用此兩種濾波器對(duì)語(yǔ)音信號(hào)進(jìn)行濾波。</p><p>　　黃生涯：IIR低通濾波器、IIR高通濾波器的設(shè)計(jì)。</p>

6、;<p>　　曾勝：IIR帶通濾波器的設(shè)計(jì)和用此濾波器對(duì)語(yǔ)音信號(hào)進(jìn)行濾波。</p><p>　　沈呈潔：系統(tǒng)界面設(shè)置和用FIR低通濾波器對(duì)語(yǔ)音信號(hào)進(jìn)行濾波。</p><p><b>　　三、詳細(xì)設(shè)計(jì)過(guò)程</b></p><p>　　1、問(wèn)題描述：錄制一段自己的語(yǔ)音信號(hào)，取不同的數(shù)據(jù)點(diǎn)對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析；對(duì)所有數(shù)據(jù)進(jìn)行插值和抽取

7、處理，改變抽樣率再對(duì)信號(hào)進(jìn)行頻譜分析；設(shè)計(jì)FIR和IIR數(shù)字濾波器，并對(duì)被抽樣后的語(yǔ)音信號(hào)進(jìn)行濾波，分析濾波后信號(hào)的時(shí)域和頻域特征，回放語(yǔ)音信號(hào)。</p><p>　　2、詳細(xì)操作步驟和部分運(yùn)行結(jié)果</p><p>　　（1）采集語(yǔ)音信號(hào)：</p><p>　　利用Windows下的錄音機(jī)，錄制一段自己的話音，時(shí)間控制在1s左右；然后在MATLAB軟件平臺(tái)下，利用函

8、數(shù)wavread對(duì)語(yǔ)音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。通過(guò)使用wavread函數(shù)，理解采樣頻率、采樣位數(shù)等概念。</p><p>　　wavread函數(shù)調(diào)用格式如下：</p><p>　　y=wavread(file)，讀取file所規(guī)定的wav文件，返回采樣值放在向量y中。</p><p>　　[y,fs,nbits]=wavread(file)，采樣值放在

9、向量y中，fs表示采樣頻率(Hz)，nbits表示采樣位數(shù)。</p><p>　　y=wavread(file,N)，讀取前N點(diǎn)的采樣值放在向量y中。</p><p>　　y=wavread(file,[N1，N2])，讀取從N1點(diǎn)到N2點(diǎn)的采樣值放在向量y中。</p><p>　　采集語(yǔ)音的程序代碼：</p><p>　　[x1,fs,nb

10、its]=wavread(‘e:\yy.wav’);</p><p>　　sound(x1,fs,bits);</p><p>　　y=fft(x1,1024);</p><p>　　subplot(2,1,1); plot(x1); title(‘原始采樣后時(shí)域信號(hào)’);</p><p>　　原始語(yǔ)音采樣后的的時(shí)域信號(hào)波形圖：</p&

11、gt;<p>　?。?）對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析：</p><p>　　在MATLAB中，利用函數(shù)fft對(duì)信號(hào)進(jìn)行快速傅里葉變換，得到信號(hào)的頻譜特性。首先畫(huà)出語(yǔ)音信號(hào)的時(shí)域波形，然后對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析。</p><p><b>　　程序代碼：</b></p><p>　　[x1,fs,bits]=wavread(‘D:\yy.w

12、av’);</p><p>　　y1=fft(x1,1024);</p><p>　　f=fs*(0:511)/1024;</p><p><b>　　figure(1)</b></p><p>　　subplot(2,1,1); plot(f,abs(1:512))); title(’原始語(yǔ)音信號(hào)頻譜’);</p

13、><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplot(2,1,2); plot(abs(y1(1:1024))); title(‘原始語(yǔ)音信號(hào)FFT頻譜’);</p><p>　　xlabel('點(diǎn)數(shù)N');</p&

14、gt;<p>　　ylabel('幅值');</p><p>　　原始語(yǔ)音信號(hào)的頻率響應(yīng)圖：</p><p>　　原始語(yǔ)音信號(hào)頻譜和FFT轉(zhuǎn)換后的頻譜圖</p><p>　?。?）設(shè)計(jì)數(shù)字濾波器和畫(huà)出頻率響應(yīng)：</p><p>　　根據(jù)語(yǔ)音信號(hào)的特點(diǎn)給出有關(guān)濾波器的性能指標(biāo)：</p><p&g

15、t;　　低通濾波器性能指標(biāo)：fp=1000Hz，fc=1200Hz，As=100dB，Ap=1dB；</p><p>　　高通濾波器性能指標(biāo)：fc=4800Hz，fp=5000Hz，As=100dB，Ap=1dB；</p><p>　　帶通濾波器性能指標(biāo)：fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，fc2=3200Hz，As=100dB，Ap=1dB。</p&g

16、t;<p>　　首先用窗函數(shù)法（矩形窗(Rectangular window)、三角窗(Triangular window)、漢寧窗(Hanning window)、海明窗(Hamming window)、布拉克曼窗(Blackman window)、切比雪夫窗(Chebyshev window)、巴特里特窗(Bartlett window)及凱塞窗(Kaiser window)。）設(shè)計(jì)上面要求的三種濾波器。在MATLAB

17、中，利用函數(shù)fir1設(shè)計(jì)FIR濾波器；然后再用雙線性變換法設(shè)計(jì)上面要求的三種濾波器，利用函數(shù)butte、cheby1和ellip設(shè)計(jì)IIR濾波器；最后，利用MATLAB中的函數(shù)freqz畫(huà)出各種濾波器的頻率響應(yīng)。</p><p>　　（4）低通濾波器的設(shè)計(jì)過(guò)程</p><p>　　用窗函數(shù)法設(shè)計(jì)低通濾波器：</p><p><b>　　程序代碼：</

18、b></p><p>　　fp=1000; fc=1200; As=100; Ap=1; fs=22050;</p><p>　　wc=2*fc/fs; wp=2*fp/fs;</p><p>　　N=ceil（(As-7.95)/(14.36*(wc-wp)/2))+1;</p><p>　　beta=0.1102*(As-8.7);

19、</p><p>　　Win=Kaiser(N+1,beta);</p><p>　　b=fir1(N,wc,Win);</p><p>　　freqz(b,1,512,fs);</p><p><b>　　運(yùn)行結(jié)果如下圖：</b></p><p>　　用雙線性變換法設(shè)計(jì)的低通濾波器:</p

20、><p><b>　　程序代碼：</b></p><p>　　fp=1000; fc=1200; As=100; Ap=1; fs=22050;</p><p>　　wc=2*fc/fs; wp=2*fp/fs;</p><p>　　[n,wn]=ellipord(wp,wc,Ap,As);</p><p

21、>　　[b,a]=ellip(n,Ap,As,wn);</p><p>　　freqz(b,a,512,fs);</p><p>　　程序運(yùn)行結(jié)果如下圖：</p><p>　?。?）高通濾波器的設(shè)計(jì)過(guò)程</p><p><b>　　程序代碼：</b></p><p>　　fp=5000; f

22、c=4800; As=100; Ap=1; fs=22050;</p><p>　　wc=2*fc/fs; wp=2*fp/fs;</p><p>　　[n,wn]=ellipord(wp,wc,Ap,As);</p><p>　　[b,a]=ellip(n,Ap,As,wn);</p><p>　　freqz(b,a,512,fs);<

23、;/p><p>　　程序運(yùn)行結(jié)果如下圖：</p><p>　?。?）帶通濾波器的設(shè)計(jì)全過(guò)程</p><p><b>　　程序代碼:</b></p><p>　　fp1=1200;fp2=3000; fc1=1000; fc2=3200;As=100; Ap=1; fs=22050;</p><p>　

24、　wc=[2*fc1/fs,2*fc2/fs]; wp=[2*fp1/fs,2*fp2/fs];</p><p>　　[n,wn]=ellipord(wp,wc,Ap,As);</p><p>　　[b,a]=ellip(n,Ap,As,wn);</p><p>　　freqz(b,a,512,fs);</p><p><b>　　

25、程序運(yùn)行結(jié)果：</b></p><p>　?。?）IIR濾波器的設(shè)計(jì)</p><p><b>　　IIR低通濾波器：</b></p><p>　　IIR高通濾波器： </p><p><b>　　IIR帶通濾波器：</b></p><p>　?。?）用濾波器對(duì)信號(hào)

26、進(jìn)行濾波</p><p>　　比較兩種濾波器性能，然后用性能好的濾波器對(duì)采集的信號(hào)進(jìn)行濾波。在MATLAB中，F(xiàn)IR濾波器利用函數(shù)fftfilt對(duì)信號(hào)進(jìn)行濾波，IIR濾波器利用函數(shù)filter對(duì)信號(hào)進(jìn)行濾波。</p><p>　　（9）比較濾波前后語(yǔ)音信號(hào)的波形及頻譜，在同一個(gè)窗體中畫(huà)出濾波前后的波形圖:</p><p>　　FIR低通濾波器濾波前后的波形圖：<

27、;/p><p>　　FIR高通濾波器濾波前后的波形圖：</p><p>　　FIR帶通濾波器濾波前后的波形圖：</p><p>　　IIR低通濾波器濾波前后的波形圖:</p><p>　　IIR高通濾波器濾波前后的波形圖：</p><p>　　IIR帶通濾波器濾波前后的波形圖：</p><p>　　

28、（10）回放語(yǔ)音信號(hào)</p><p>　　在MATLAB中，函數(shù)sound可以對(duì)聲音進(jìn)行回放。其調(diào)用格式為</p><p>　　sound(x,fs,bits)，可以感覺(jué)濾波前后的聲音有變化。</p><p>　?。?1）設(shè)計(jì)系統(tǒng)界面</p><p>　　在所設(shè)計(jì)的系統(tǒng)界面上可以選擇濾波器的類型，輸入濾波器的參數(shù)，顯示濾波器的頻率響應(yīng)，選擇信

29、號(hào)等。</p><p><b>　　界面如下圖所示：</b></p><p><b>　　四、調(diào)試分析</b></p><p>　　1、A．用窗函數(shù)法設(shè)計(jì)低通濾波器中，語(yǔ)句N=ceil(As-7.95)/(14.36*(wc-wp)/2))+1;出錯(cuò)，將其改為N=ceil((As-7.95)/(14.36*(wc-wp)/

30、2))+1;即可。</p><p>　　2、在帶通濾波器的設(shè)計(jì)中程序出現(xiàn)如下錯(cuò)誤：??? Undefined function or variable 's'.</p><p>　　經(jīng)過(guò)調(diào)試分析，觀察后發(fā)現(xiàn)是書(shū)寫錯(cuò)誤，將語(yǔ)句wp=[2*fp1/fs,2*fp2/s];改成wp=[2*fp1/fs,2*fp2/fs];即可</p><p><b&

31、gt;　　五、結(jié)果分析</b></p><p>　　1、題目要求最大衰減Ap=1dB，fp=1000Hz，由圖可看出基本符合，而最小衰減As=100dB，fc=1200Hz，基本符合。</p><p>　　2、濾波前后的頻譜分析：由于所取的采樣點(diǎn)數(shù)比較大，濾波前后的頻譜比較相近，但仔細(xì)對(duì)比仍然可以看出下圖的結(jié)果濾掉了高頻的成分。</p><p>　　3、

32、題目要求最大衰減Ap=1dB，fp=5000Hz，由圖可看出，基本符合，而最小衰減As=100dB，fc=4800Hz，也基本符合。</p><p>　　4、濾波前后的頻譜分析:由圖明顯可以看出，此濾波器將低頻成分濾掉了。</p><p>　　5、題目要求最大衰減Ap=1dB時(shí) ,fp1=1000Hz，fp2=3200Hz，由圖可看出，基本符合，而最小衰減As=100dB時(shí),fc1=100

33、0Hz，fc2=3200Hz，基本符合。</p><p>　　6、濾波前后的頻譜分析：該濾波器將低頻成分和高頻成分都濾掉了，只剩下通帶里的頻率。</p><p>　　7、選用凱澤窗設(shè)計(jì)方法，濾波器的幅度和相位響應(yīng)滿足設(shè)計(jì)的指標(biāo)但濾波器長(zhǎng)度（N=708）太長(zhǎng)，實(shí)現(xiàn)起來(lái)很困難，主要原因是濾波器的指標(biāo)太苛刻，因此，一向不用窗函數(shù)法設(shè)計(jì)這種類型的濾波器。在雙線性變換法的設(shè)計(jì)中選用橢圓函數(shù)設(shè)計(jì)，濾波

34、器的幅度和相位均滿足設(shè)計(jì)的要求。</p><p><b>　　六、心得體會(huì)</b></p><p>　　通過(guò)一學(xué)期的數(shù)字信號(hào)處理的學(xué)習(xí)，在這次做課程設(shè)計(jì)的過(guò)程中，真的發(fā)現(xiàn)自己沒(méi)學(xué)到啥東西，所以導(dǎo)致做的過(guò)程中困難重重，拿著課本只好重頭看，通過(guò)網(wǎng)上查資料，同學(xué)之間的相互合作，最終我們還是戰(zhàn)勝了困難，順利完成了此次的課程設(shè)計(jì)。我深深的體會(huì)到了知識(shí)的重要性，不管是為了應(yīng)付考試

35、還是怎么樣也好，對(duì)于自身的幫助還是很大的，整整兩天坐在電腦面前的功夫沒(méi)有白費(fèi)，我還是收獲了很多，對(duì)于Matlab軟件的使用也更加熟悉了。</p><p>　　一年之計(jì)在于春，而我們一學(xué)期的學(xué)習(xí)就在于這次課程設(shè)計(jì)了，我覺(jué)得這才是真正學(xué)到知識(shí)的，也鍛煉了自己的動(dòng)手實(shí)踐能力，耐心和毅力。學(xué)會(huì)團(tuán)結(jié)，才能把事情做得更好。</p><p><b>　　七、參考文獻(xiàn)和附錄</b>&

36、lt;/p><p>　　[1]丁美玉，高西全 數(shù)字信號(hào)處理2版。西安：西安電子科技大學(xué)出版社，2001</p><p>　　[2]懷琛 數(shù)字信號(hào)處理教程——MATLAB釋疑與實(shí)現(xiàn)。北京：電子工業(yè)出版社，2004</p><p>　　[3]王宏 MATLAB6.5及其在信號(hào)處理中的應(yīng)用。北京：清華大學(xué)出版社，2004</p><p>　　[4]劉順

37、蘭，吳杰 數(shù)字信號(hào)處理。西安：西安電子科技大學(xué)出版社，2003</p><p><b>　　附錄：</b></p><p><b>　　源程序代碼：</b></p><p>　　b=menu('請(qǐng)選擇選項(xiàng)','原始信號(hào)采樣后時(shí)域圖和頻譜圖','FIR濾波器','IIR

38、濾波器','退出');</p><p><b>　　if b==4</b></p><p><b>　　b==0;</b></p><p><b>　　end</b></p><p><b>　　if b==1</b></p&

39、gt;<p>　　temp=menu('請(qǐng)選擇選項(xiàng)','播放原始語(yǔ)音','原始語(yǔ)音時(shí)域圖和頻率響應(yīng)圖','原始語(yǔ)音頻譜圖及做FFT變換后頻譜圖','返回');</p><p>　　if temp==1</p><p><b>　　voice; </b></p>

40、<p><b>　　main;</b></p><p>　　else if temp==2</p><p><b>　　pingpu;</b></p><p><b>　　main;</b></p><p>　　else main;</p><p

41、><b>　　end</b></p><p>　　else if b==2</p><p><b>　　temp==1</b></p><p><b>　　FIR_LP;</b></p><p><b>　　main;</b></p>

42、<p>　　else if temp==2</p><p><b>　　FIR_HP;</b></p><p><b>　　main;</b></p><p>　　else if temp==3</p><p><b>　　FIR_BP;</b></p>

43、<p><b>　　main;</b></p><p><b>　　else</b></p><p><b>　　main;</b></p><p><b>　　end</b></p><p>　　else if b==3</p>

44、<p>　　temp=menu('請(qǐng)選擇選項(xiàng)','IIR低通濾波器','IIR高通濾波器','IIR帶通濾波器','返回');</p><p>　　if temp==1</p><p><b>　　IIR_LP;</b></p><p><b>

45、;　　main;</b></p><p>　　else if temp==2</p><p><b>　　IIR_HP;</b></p><p><b>　　main;</b></p><p>　　else if temp==3</p><p><b>

46、　　IIR_BP;</b></p><p><b>　　main;</b></p><p><b>　　else</b></p><p><b>　　main;</b></p><p><b>　　end</b></p><

47、p><b>　　end</b></p><p><b>　　fs=22050;</b></p><p>　　x1=wavread('D:\yy.wav');</p><p>　　wp1=2*pi*1200/fs;wp2=2*pi*3000/fs;</p><p>　　ws1=2

48、*pi*1000/fs;ws2=2*pi*3200/fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rs=100;</b></p><p>　　wn=(wp1+ws1)/2;ws=(wp2+ws2)/2;</p><p>　　wdelta=wp1-ws1;<

49、/p><p>　　N=ceil(8*i/wdelta);</p><p>　　wn=[wp,ws];</p><p>　　[b,a]=firl(N,wn/pi,'bandpass');</p><p><b>　　figure(1)</b></p><p>　　freqz(b,a,51

50、2);</p><p>　　title('FiR帶通濾波器');</p><p>　　f1=filter(b,a,x1);</p><p><b>　　figure(2)</b></p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft(x1,1024

51、);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR帶通濾波器濾波前的頻譜')；</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplot(2,1

52、,2);plot(f,abs(y1(1:512)));title('FIR帶通濾波器濾波后的頻譜')</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p><b>　　fs=22050;</b></p><p>

53、;　　x1=wavread('D:\yy.wav');</p><p>　　wp=2*pi*5000/fs;</p><p>　　ws=2*pi*4800/fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rs=100;</b></p>

54、<p>　　wdelta=wp-ws;</p><p>　　N=ceil(8*pi/wdelta);</p><p>　　wn=(wp+ws)/2;</p><p>　　[b,a]=firl(N,wn/pi,'high');</p><p><b>　　figure(1)</b></p&

55、gt;<p>　　freqz(b,a,512);title('FIR高通濾波器'）;</p><p>　　f2=filter(b,a,x1);</p><p><b>　　figure(2)</b></p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft(x1

56、,1024);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR高通濾波器濾波前的頻譜')；</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplo

57、t(2,1,2);plot(f,abs(y1(1:512)));title('FIR高通濾波器濾波后的頻譜')</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p><b>　　fs=10000;</b></p><

58、;p>　　x1=wavread('D:\yy.wav');</p><p>　　wp=2*pi*1000/fs;</p><p>　　ws=2*pi*1200/fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rs=100;</b></p

59、><p>　　wdelta=ws-wp;</p><p>　　N=ceil(8*pi/wdelta);</p><p>　　wn=(wp+ws)/2;</p><p>　　[b,a]=firl(N,wn/pi,hamming(N+1));</p><p><b>　　figure(1)</b><

60、/p><p>　　freqz(b,a,512);title('FIR低通濾波器'）;</p><p>　　f1=filter(b,a,x1);</p><p><b>　　figure(2)</b></p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft

61、(x1,1024);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR低通濾波器濾波前的頻譜')；</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　sub

62、plot(2,1,2);plot(f,abs(y1(1:512)));title('FIR低通濾波器濾波后的頻譜')</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p><b>　　fs=22050;</b></p>

63、<p>　　x1=wavread('D:\yy.wav');</p><p>　　Ts=1/fs;R1=30;</p><p>　　fb1=1200;fb2=3000;fc1=1000;fc2=3200;fs=22050;</p><p>　　w1=2*fb1*pi/fs;w2=2*fc1*pi/fs;</p><p&g

64、t;　　w3=2*fb2*pi/fs;w4=2*fc2*pi/fs;</p><p>　　wp=[w1,w3];ws=[w2,w4];</p><p><b>　　Rp=1;</b></p><p><b>　　Rl=100;</b></p><p>　　wp1=2/Ts*tan(wp/2);ws1

65、=2/Ts*tan(ws/2);</p><p>　　[N,wn]=cheb2ord(wp1,ws1,Rp,Rl,'s');</p><p>　　[Z,P,K]=cheb2ap(N,R);</p><p>　　[Bap,Aap]=zpwtf(Z,P,K);</p><p>　　[b,a]=Ip2bp(Bap,Aap,2100*

66、2*pi,1800*2*pi);</p><p>　　[bz,az]=bilinear(b,a,fs);</p><p>　　[H,W]=freqz(bz,az);</p><p><b>　　figure(1)</b></p><p>　　plot(w*fs/(2*pi),abs(H))</p><

67、p><b>　　grid</b></p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('頻率響應(yīng)幅度');</p><p>　　title('IIR帶通濾波器')</p><p>　　fl=filter(bz,az,x1);<

68、/p><p><b>　　figure(2)</b></p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft(x1,1024);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1;512)));title('IIR帶通濾波器濾波前的頻譜');&l

69、t;/p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplot(2,1,2);plot(f,abs(y1(1;512)));title('IIR帶通濾波器濾波后的頻譜');</p><p>　　xlabel('頻率/

70、Hz');</p><p>　　ylabel('幅值');</p><p><b>　　fs=22050;</b></p><p>　　x1=wavread('D:\yy.wav');</p><p>　　Ts=1/fs;R1=50;</p><p>　　w

71、p=2*pi*5000/fs;</p><p>　　ws=2*pi*4800/fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rl=100;</b></p><p>　　wp1=2/Ts*tan(wp/2);ws1=2/Ts*tan(ws/2);</p>

72、;<p>　　[N,wn]=chebyord(wp1,ws1,Rp,Rl,'s');</p><p>　　[Z,P,K]=chebyap(N,R);</p><p>　　[Bap,Aap]=zpwtf(Z,P,K);</p><p>　　[b,a]=Ip2bp(Bap,Aap,2100*2*pi,1800*2*pi);</p>

73、;<p>　　[bz,az]=bilinear(b,a,fs);</p><p>　　[H,W]=freqz(bz,az);</p><p><b>　　figure(1)</b></p><p>　　plot(w*fs/(2*pi),abs(H))</p><p><b>　　grid</

74、b></p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('頻率響應(yīng)幅度');</p><p>　　title('IIR高通濾波器')</p><p>　　fl=filter(bz,az,x1);</p><p><b>

75、;　　figure(2)</b></p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft(x1,1024);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1;512)));title('IIR高通濾波器濾波前的頻譜');</p><p>　　xlab

76、el('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplot(2,1,2);plot(f,abs(y1(1;512)));title('IIR高通濾波器濾波后的頻譜');</p><p>　　xlabel('頻率/Hz');</p><

77、p>　　ylabel('幅值');</p><p><b>　　fs=22050;</b></p><p>　　x1=wavread('D:\yy.wav');</p><p>　　Ts=1/fs;R1=10;</p><p>　　wp=2*pi*1000/fs;</p>

78、<p>　　ws=2*pi*1200/fs;</p><p><b>　　Rp=1;</b></p><p><b>　　Rl=100;</b></p><p>　　wp1=2/Ts*tan(wp/2);ws1=2/Ts*tan(ws/2);</p><p>　　[N,wn]=butt

79、ord(wp1,ws1,Rp,Rl,'s');</p><p>　　[Z,P,K]=buttap(N,R);</p><p>　　[Bap,Aap]=zpwtf(Z,P,K);</p><p>　　[b,a]=Ip2bp(Bap,Aap,2100*2*pi,1800*2*pi);</p><p>　　[bz,az]=bilin

80、ear(b,a,fs);</p><p>　　[H,W]=freqz(bz,az);</p><p><b>　　figure(1)</b></p><p>　　plot(w*fs/(2*pi),abs(H))</p><p><b>　　grid</b></p><p>　

81、　xlabel('頻率/Hz');</p><p>　　ylabel('頻率響應(yīng)幅度');</p><p>　　title('IIR低通濾波器')</p><p>　　fl=filter(bz,az,x1);</p><p><b>　　figure(2)</b></

82、p><p>　　y1=fft(f1,1024);</p><p>　　y2=fft(x1,1024);</p><p>　　subplot(2,1,1);plot(f,abs(y2(1;512)));title('IIR低通濾波器濾波前的頻譜');</p><p>　　xlabel('頻率/Hz');</p&

83、gt;<p>　　ylabel('幅值');</p><p>　　subplot(2,1,2);plot(f,abs(y1(1;512)));title('IIR低通濾波器濾波后的頻譜');</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值')

84、;</p><p>　　[x1,fs,bits]=wavread(‘D:\yy.wav’);</p><p>　　y1=fft(x1,1024);</p><p>　　f=fs*(0:511)/1024;</p><p><b>　　figure(1)</b></p><p>　　subplot(

85、2,1,1); plot(f,abs(1:512))); title(’原始語(yǔ)音信號(hào)頻譜’);</p><p>　　xlabel('頻率/Hz');</p><p>　　ylabel('幅值');</p><p>　　subplot(2,1,2); plot(abs(y1(1:1024))); title(‘原始語(yǔ)音信號(hào)FFT頻譜’)

86、;</p><p>　　xlabel('點(diǎn)數(shù)N');</p><p>　　ylabel('幅值');</p><p>　　[x1,fs,bits]=wavread(‘D:\yy.wav’);</p><p>　　sound(x1,fs);</p><p>　　plot(x1);title