應(yīng)用matlab對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析及濾波

數(shù)字信號(hào)處理—綜合實(shí)驗(yàn)報(bào)告綜合實(shí)驗(yàn)名稱(chēng)：應(yīng)用MatLab對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析及濾波系：學(xué)生姓名：班級(jí)：學(xué)號(hào)：成績(jī)：指導(dǎo)教師：開(kāi)課時(shí)間學(xué)年學(xué)期目錄TOC\o"1-3"\h\u29788一．綜合實(shí)驗(yàn)題目114094二、綜合實(shí)驗(yàn)?zāi)康暮鸵饬x1243752.1綜合實(shí)驗(yàn)?zāi)康?310232.2綜合實(shí)驗(yàn)的意義12654三．綜合實(shí)驗(yàn)的主要內(nèi)容和要求1138833.2綜合實(shí)驗(yàn)的要求：224151四．實(shí)驗(yàn)的原理2176904.1數(shù)字濾波器的概念233234.2數(shù)字濾波器的分類(lèi)213453〔1〕根據(jù)單位沖激響應(yīng)h(n)的時(shí)間特性分類(lèi)29514五．實(shí)驗(yàn)的步驟311914下面對(duì)各步驟加以具體說(shuō)明。5.1語(yǔ)音信號(hào)的采集333765.2語(yǔ)音信號(hào)的頻譜分析;370655.3設(shè)計(jì)數(shù)字濾波器和畫(huà)出其頻率響應(yīng)5248925.3.1設(shè)計(jì)數(shù)字濾波器的性能指標(biāo)：5117535.3.2用Matlab設(shè)計(jì)數(shù)字濾波器6242115.6設(shè)計(jì)系統(tǒng)界面1931282六、心得體會(huì)202412參考文獻(xiàn)：21一．綜合實(shí)驗(yàn)題目應(yīng)用MatLab對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析及濾波二、綜合實(shí)驗(yàn)?zāi)康暮鸵饬x2.1綜合實(shí)驗(yàn)?zāi)康臑榱朔€(wěn)固所學(xué)的數(shù)字信號(hào)處理理論知識(shí)，使學(xué)生對(duì)信號(hào)的采集、處理、傳輸、顯示和存儲(chǔ)等有一個(gè)系統(tǒng)的掌握和理解，再者，加強(qiáng)學(xué)生對(duì)Matlab軟件在信號(hào)分析和處理的運(yùn)用綜合運(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í)的理解，建立概念。2.2綜合實(shí)驗(yàn)的意義語(yǔ)言是我們?nèi)祟?lèi)所特有的功能，它是傳承和記載人類(lèi)幾千年文明史，沒(méi)有語(yǔ)言就沒(méi)有我們今天人類(lèi)的文明。語(yǔ)音是語(yǔ)言最根本的表現(xiàn)形式，是相互傳遞信息最重要的手段，是人類(lèi)最重要、最有效、最常用和最方便的交換信息的形式。語(yǔ)音信號(hào)處理屬于信息科學(xué)的一個(gè)重要分支，大規(guī)模集成技術(shù)的高度開(kāi)展和計(jì)算機(jī)技術(shù)的飛速前進(jìn)，推動(dòng)了這一技術(shù)的開(kāi)展；它是研究用數(shù)字信號(hào)處理技術(shù)對(duì)語(yǔ)音信號(hào)進(jìn)行處理的一門(mén)新興學(xué)科，同時(shí)又是綜合性的多學(xué)科領(lǐng)域和涉及面很廣的交叉學(xué)科，因此我們進(jìn)行語(yǔ)言信號(hào)處理具有時(shí)代的意義。三．綜合實(shí)驗(yàn)的主要內(nèi)容和要求3.1綜合實(shí)驗(yàn)的主要內(nèi)容：錄制一段個(gè)人自己的語(yǔ)音信號(hào)，并對(duì)錄制的信號(hào)進(jìn)行采樣；畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形和頻譜圖；給定濾波器的性能指標(biāo)，采用窗函數(shù)法和雙線性變換設(shè)計(jì)濾波器，并畫(huà)出濾波器的頻率響應(yīng)；然后用自己設(shè)計(jì)的濾波器對(duì)采集的信號(hào)進(jìn)行濾波，畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行比照，分析信號(hào)的變化；回放語(yǔ)音信號(hào)；綜合實(shí)驗(yàn)應(yīng)完成的工作：(1〕語(yǔ)音信號(hào)的采集;(2〕語(yǔ)音信號(hào)的頻譜分析;(3〕數(shù)字濾波器的設(shè)計(jì);(4〕對(duì)語(yǔ)音信號(hào)進(jìn)行濾波處理;(5〕對(duì)濾波前后的語(yǔ)音信號(hào)頻譜進(jìn)行比照分析;3.2綜合實(shí)驗(yàn)的要求：(1〕學(xué)生能夠根據(jù)設(shè)計(jì)內(nèi)容積極主動(dòng)查找相關(guān)資料;(2〕濾波器的性能指標(biāo)可以根據(jù)實(shí)際情況作調(diào)整;(3〕對(duì)設(shè)計(jì)結(jié)果進(jìn)行獨(dú)立思考和分析;(4〕設(shè)計(jì)完成后,要提交相關(guān)的文檔;①綜合實(shí)驗(yàn)報(bào)告書(shū)(紙質(zhì)和電子版各一份,具體格式參照學(xué)校課程設(shè)計(jì)管理規(guī)定),報(bào)告內(nèi)容要涵蓋設(shè)計(jì)題目、設(shè)計(jì)任務(wù)、詳細(xì)的設(shè)計(jì)過(guò)程、原理說(shuō)明、、頻譜圖的分析、調(diào)試總結(jié)、心得體會(huì)、參考文獻(xiàn)〔在報(bào)告中參考文獻(xiàn)要做標(biāo)注，不少于5篇〕。②可運(yùn)行的源程序代碼(電子版)(5〕在根本要求的根底上,學(xué)生可以根據(jù)個(gè)人對(duì)該綜合實(shí)驗(yàn)的理解,添加一些新的內(nèi)容。四．實(shí)驗(yàn)的原理4.1數(shù)字濾波器的概念：假設(shè)濾波器的輸入、輸出都是離散時(shí)間信號(hào)，那么該濾波器的單位沖激響應(yīng)h(n)也必然是離散的，這種濾波器稱(chēng)為數(shù)字濾波器〔DF〕?！?】數(shù)字濾波器實(shí)質(zhì)上是一種運(yùn)算過(guò)程——用來(lái)描述離散系統(tǒng)輸入與輸出關(guān)系的差分方程的計(jì)算或卷積計(jì)算。數(shù)字濾波器的實(shí)質(zhì)是用一有限精度算法實(shí)現(xiàn)的離散時(shí)間線性時(shí)不變系統(tǒng)，以完成對(duì)信號(hào)進(jìn)行濾波處理的過(guò)程。它是數(shù)字信號(hào)處理的一個(gè)重要分支，具有穩(wěn)定性好、精度高、靈活性強(qiáng)、體積小、質(zhì)量輕等諸多優(yōu)點(diǎn)?！?】4.2數(shù)字濾波器的分類(lèi)：數(shù)字濾波器根據(jù)不同的分類(lèi)標(biāo)準(zhǔn)可以將濾波器分成不同的類(lèi)別。〔1〕根據(jù)單位沖激響應(yīng)h(n)的時(shí)間特性分類(lèi)無(wú)限沖激響應(yīng)〔IIR〕數(shù)字濾波器有限沖激響應(yīng)〔FIR〕數(shù)字濾波器[1]〔2〕根據(jù)實(shí)現(xiàn)方法和形式分類(lèi)遞歸型數(shù)字濾波器非遞歸型數(shù)字濾波器快速卷積型〔3〕根據(jù)頻率特性分類(lèi)低通數(shù)字濾波器、高通數(shù)字濾波器、帶通數(shù)字波器4.3濾波器的作用：濾波器是一種用來(lái)消除干擾雜訊的器件，將輸入或輸出經(jīng)過(guò)過(guò)濾而得到純潔的直流電【3】。對(duì)特定頻率的頻點(diǎn)或該頻點(diǎn)以外的頻率進(jìn)行有效濾除的電路，就是濾波器，其功能就是得到一個(gè)特定頻率或消除一個(gè)特定頻率。五．實(shí)驗(yàn)的步驟錄制一段個(gè)人自己的語(yǔ)音信號(hào)，并對(duì)錄制的信號(hào)進(jìn)行采樣；畫(huà)出采樣后語(yǔ)音信號(hào)的時(shí)域波形和頻譜圖；給定濾波器的性能指標(biāo)，采用窗函數(shù)法和雙線性變換設(shè)計(jì)濾波器，并畫(huà)出濾波器的頻率響應(yīng)；然后用自己設(shè)計(jì)的濾波器對(duì)采集的信號(hào)進(jìn)行濾波，畫(huà)出濾波后信號(hào)的時(shí)域波形和頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行比照，分析信號(hào)的變化；回放語(yǔ)音信號(hào)；最后，設(shè)計(jì)一個(gè)信號(hào)處理系統(tǒng)界面。下面對(duì)各步驟加以具體說(shuō)明。

5.1語(yǔ)音信號(hào)的采集：利用PC機(jī)上的聲卡和Windows操作系統(tǒng)可以進(jìn)行數(shù)字信號(hào)的采集。將話筒插入計(jì)算機(jī)的語(yǔ)音輸入插口上,啟動(dòng)錄音機(jī)。按下錄音按鈕，對(duì)話筒讀一段話,讀完后停止錄音,屏幕左側(cè)將顯示所錄聲音的長(zhǎng)度?！?】5.2語(yǔ)音信號(hào)的頻譜分析：(1)首先畫(huà)出由wavread函數(shù)采樣后的信號(hào)的時(shí)域波形和它的頻域響應(yīng)波形；語(yǔ)音信號(hào)的程序如下:[x1,fs,bits]=wavread('c:\wudan2.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1figure(20)plot(x1)%做原始語(yǔ)音信號(hào)以44.1k采樣后的時(shí)域圖形title('原始語(yǔ)音采樣后時(shí)域信號(hào)');xlabel('時(shí)間軸n');ylabel('幅值A(chǔ)');figure(21)freqz(x1)%繪制原始語(yǔ)音信號(hào)采樣后的頻率響應(yīng)圖title('原始語(yǔ)音信號(hào)采樣后頻率響應(yīng)圖');運(yùn)行結(jié)果:圖1圖2〔2)對(duì)語(yǔ)音信號(hào)進(jìn)行FFT變換，得到FFT頻譜特性曲線，與原語(yǔ)音信號(hào)的頻譜特性曲線進(jìn)行比擬。語(yǔ)音信號(hào)FFT的程序如下:[x1,fs,bits]=wavread('c:\wudan2.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1y1=fft(x1,1024);%對(duì)采樣后信號(hào)做1024點(diǎn)FFT變換f=fs*(0:511)/1024;figure(19)subplot(2,1,1);plot(f,abs(y1(1:512)));title('原始語(yǔ)音信號(hào)頻譜')xlabel('頻率/Hz');ylabel('幅值')subplot(2,1,2);plot(abs(y1(1:1024)))%采樣后信號(hào)的FFT頻譜圖title('原始語(yǔ)音信號(hào)FFT頻譜')xlabel('點(diǎn)數(shù)N');ylabel('幅值');運(yùn)行結(jié)果：圖35.3設(shè)計(jì)數(shù)字濾波器和畫(huà)出其頻率響應(yīng)設(shè)計(jì)數(shù)字濾波器的性能指標(biāo)：①低通濾波器性能指標(biāo)fb=1000Hz，fc＝1200Hz，As＝100db，Ap＝1db。②高通濾波器性能指標(biāo)fs＝4800Hz，fb＝5000Hz，As＝100db，Ap＝1db。③帶通濾波器性能指標(biāo)fb1＝1200Hz，fb2＝3000Hz，fc1＝1000Hz，fc2＝3200Hz，As＝100db，Ap＝1db。用Matlab設(shè)計(jì)數(shù)字濾波器①用雙線性變換法設(shè)計(jì)上面要求的3種濾波器。在Matlab中，可以利用函數(shù)butter，cheby1和ellip設(shè)計(jì)IIR濾波器；利用Matlab中的函數(shù)frepz畫(huà)出各濾波器的頻率響應(yīng)?！?】②用窗函數(shù)法設(shè)計(jì)上面要求的3種濾波器。在Matlab中，可以利用函數(shù)fir1設(shè)計(jì)FIR濾波器；利用Matlab中的函數(shù)freqz畫(huà)出各濾波器的頻率響應(yīng)。

5.4對(duì)語(yǔ)音信號(hào)進(jìn)行濾波處理用自己設(shè)計(jì)的各濾波器分別對(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)行濾波。并比擬濾波前后語(yǔ)音信號(hào)的波形及頻譜?！?】IIR低通濾波器程序如下:fs=22050;[x2]=wavread('c:\wudan2.wav');Ts=1/fs;R1=10;wp=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rl=100;wp1=2/Ts*tan(wp/2);%將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)ws1=2/Ts*tan(ws/2);[N,Wn]=buttord(wp1,ws1,Rp,R1,'s');%選擇濾波器的最小階數(shù)[Z,P,K]=buttap(N);%創(chuàng)立butterworth模擬濾波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2lp(Bap,Aap,Wn);[bz,az]=bilinear(b,a,fs);%用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換[H,W]=freqz(bz,az);%繪制頻率響應(yīng)曲線figure(16)plot(W*fs/(2*pi),abs(H))gridxlabel('頻率／Hz')ylabel('頻率響應(yīng)幅度')title('IIR低通濾波器')f1=filter(bz,az,x2);figure(17)subplot(2,1,1)plot(x2)%畫(huà)出濾波前的時(shí)域圖title('IIR低通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f1);%畫(huà)出濾波后的時(shí)域圖title('IIR低通濾波器濾波后的時(shí)域波形');sound(f1,44100);%播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;figure(18)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));%畫(huà)出濾波前的頻譜圖title('IIR低通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)F1=plot(f,abs(F0(1:512)));%畫(huà)出濾波后的頻譜圖title('IIR低通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');【5】運(yùn)行結(jié)果：圖4圖5圖6IIR高通濾波器程序如下:fs=22050;x1=wavread('c:\wudan2.wav');Fs=22050;Ts=1/Fs;R1=50;Wp=2*pi*5000/fs;Ws=2*pi*4800/fs;Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2);%將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)Ws1=2/Ts*tan(Ws/2);[N,Wn]=cheb2ord(Wp1,Ws1,Rp,Rl,'s');%選擇濾波器的最小階數(shù)[Z,P,K]=cheb2ap(N,Rl);%創(chuàng)立切比雪夫模擬濾波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2hp(Bap,Aap,Wn);[bz,az]=bilinear(b,a,Fs);%用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換[H,W]=freqz(bz,az);%繪制頻率響應(yīng)曲線figure(13)plot(W*fs/(2*pi),abs(H))gridxlabel('頻率／Hz')ylabel('頻率響應(yīng)幅度')title('IIR高通濾波器')f1=filter(bz,az,x1);figure(14)subplot(2,1,1)plot(x1)%畫(huà)出濾波前的時(shí)域圖title('IIR高通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f1);%畫(huà)出濾波后的時(shí)域圖title('IIR高通濾波器濾波后的時(shí)域波形');sound(f1,44100);%播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;figure(15)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));%畫(huà)出濾波前的頻譜圖title('IIR高通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)plot(f,abs(F0(1:512)));%畫(huà)出濾波后的頻譜圖title('IIR高通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');運(yùn)行結(jié)果：圖7圖8圖9IIR帶通濾波器程序如下:fs=22050;x1=wavread('c:\wudan2.wav');Fs=22050;Ts=1/Fs;R1=30;fb1=1200;fb2=3000;fc1=1000;fc2=3200;fs=22050;W1=2*fb1*pi/fs;W2=2*fc1*pi/fs;W3=2*fb2*pi/fs;W4=2*fc2*pi/fs;Wp=[W1,W3];Ws=[W2,W4];Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2);%將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)Ws1=2/Ts*tan(Ws/2);[N,Wn]=cheb2ord(Wp1,Ws1,Rp,R1,'s');%選擇濾波器的最小階數(shù)[Z,P,K]=cheb2ap(N,Rl);%創(chuàng)立切比雪夫模擬濾波器[Bap,Aap]=zp2tf(Z,P,K);[b,a]=lp2bp(Bap,Aap,2100*2*pi,1800*2*pi);[bz,az]=bilinear(b,a,Fs);%用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換[H,W]=freqz(bz,az);%繪制頻率響應(yīng)曲線figure(10)plot(W*fs/(2*pi),abs(H))gridxlabel('頻率／Hz')ylabel('頻率響應(yīng)幅度')title('IIR帶通濾波器')f1=filter(bz,az,x1);figure(11)subplot(2,1,1)plot(x1)%畫(huà)出濾波前的時(shí)域圖title('IIR帶通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f1);%畫(huà)出濾波后的時(shí)域圖title('IIR帶通濾波器濾波后的時(shí)域波形');sound(f1,44100);%播放濾波后的信號(hào)F0=fft(f1,1024);f=fs*(0:511)/1024;figure(12)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));%畫(huà)出濾波前的頻譜圖title('IIR帶通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)plot(f,abs(F0(1:512)));%畫(huà)出濾波后的頻譜圖title('IIR帶通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');【6】運(yùn)行結(jié)果：圖10圖11圖12FIR低通濾波器程序如下:fs=10000;x1=wavread('c:\wudan2.wav');wp=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rs=100;wdelta=ws-wp;N=ceil(8*pi/wdelta);%取整wn=(wp+ws)/2;[b,a]=fir1(N,wn/pi,hamming(N+1));%選擇窗函數(shù)，并歸一化截止頻率figure(7)freqz(b,a,512);title('FIR低通濾波器');f2=filter(b,a,x1);figure(8)subplot(2,1,1)plot(x1)title('FIR低通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f2);title('FIR低通濾波器濾波后的時(shí)域波形');sound(f2,44100);%播放濾波后的語(yǔ)音信號(hào)F0=fft(f2,1024);f=fs*(0:511)/1024;figure(9)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR低通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)F2=plot(f,abs(F0(1:512)));title('FIR低通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');【7】運(yùn)行結(jié)果：圖13圖14圖15FIR高通濾波器程序如下:fs=22050;x1=wavread('c:\wudan2.wav');wp=2*pi*5000/fs;ws=2*pi*4800/fs;Rp=1;Rs=100;wdelta=wp-ws;N=ceil(8*pi/wdelta);%取整wn=(wp+ws)/2;[b,a]=fir1(N,wn/pi,'high');figure(4)freqz(b,a,512);title('FIR高通濾波器');f2=filter(b,a,x1);figure(5)subplot(2,1,1)plot(x1)title('FIR高通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f2);title('FIR高通濾波器濾波后的時(shí)域波形');sound(f2,44100);%播放濾波后的語(yǔ)音信號(hào)F0=fft(f2,1024);f=fs*(0:511)/1024;figure(6)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR高通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)F2=plot(f,abs(F0(1:512)));title('FIR高通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');【9】運(yùn)行結(jié)果:圖16圖17圖18FIR帶通濾波器程序如下:fs=22050;x1=wavread('c:\wudan2.wav');wp1=2*pi*1200/fs;wp2=2*pi*3000/fs;ws1=2*pi*1000/fs;ws2=2*pi*3200/fs;Rp=1;Rs=100;wp=(wp1+ws1)/2;ws=(wp2+ws2)/2;wdelta=wp1-ws1;N=ceil(8*pi/wdelta);%取整wn=[wpws];[b,a]=fir1(N,wn/pi,'bandpass');figure(1)freqz(b,a,512);title('FIR帶通濾波器');f2=filter(b,a,x1);figure(2)subplot(2,1,1)plot(x1)title('FIR帶通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f2);title('FIR帶通濾波器濾波后的時(shí)域波形');sound(f2,44100);%播放濾波后的語(yǔ)音信號(hào)F0=fft(f2,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR帶通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)plot(f,abs(F0(1:512)));title('FIR帶通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值');運(yùn)行結(jié)果：圖19圖20圖215.5回放語(yǔ)音信號(hào)

在Matlab中，函數(shù)sound可以對(duì)聲音進(jìn)行回放。其調(diào)用格式：sound(x，fs，bits)；可以感覺(jué)濾波前后的聲音有變化?；胤耪Z(yǔ)音程序：[x1,fs,bits]=wavread('c:\wudan2.wav');%讀取語(yǔ)音信號(hào)的數(shù)據(jù)，賦給變量x1sound(x1,fs);%播放語(yǔ)音信號(hào)【8】5.6設(shè)計(jì)系統(tǒng)界面本系統(tǒng)界面主要采用if…else語(yǔ)句和menu菜單函數(shù)進(jìn)行設(shè)計(jì)。程序如下：b=menu('請(qǐng)選擇選項(xiàng)','原始信號(hào)采樣后時(shí)域圖和頻譜圖','FIR濾波器','IIR濾波器','退出');ifb==4b==0;endifb==1temp=menu('請(qǐng)選擇選項(xiàng)','播放原始語(yǔ)音','原始語(yǔ)音時(shí)域圖和頻率響應(yīng)圖','原始語(yǔ)音頻譜圖及做FFT變換后頻譜圖','返回');iftemp==1voice;main;elseiftemp==2time;main;elseiftemp==3pingpu;main;elsemain;endelseifb==2temp=menu('請(qǐng)選擇選項(xiàng)','FIR低通濾波器','FIR高通濾波器','FIR帶通濾波器','返回');iftemp==1FIR_LP;main;elseiftemp==2FIR_HP;main;elseiftemp==3FIR_BP;main;elsemain;endelseifb==3temp=menu('請(qǐng)選擇選項(xiàng)','IIR低通濾波器','IIR高通濾波器','IIR帶通濾波器','返回');iftemp==1IIR_LP;main;elseiftemp==2IIR_HP;main;elseiftemp==3IIR_BP;main;else