【文章內(nèi)容簡介】 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲S=fft(s)。 %傅里葉變換z11=filter(num11,den11,s)。 %IIR濾波器濾波sound(z11)。m11=fft(z11)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。b39。)。title(39。濾波前信號的頻譜39。)。grid。subplot(2,2,2)。plot(abs(m11),39。r39。)。title(39。濾波后信號的頻譜39。)。grid。subplot(2,2,3)。plot(s,39。b39。)。title(39。濾波前信號的波形39。)。grid。subplot(2,2,4)。plot(z11,39。r39。)。title(39。濾波后的信號波形39。)。圖9 雙線性法IIR低通濾波%**********************IIR帶通濾波器*****************************[y,fs,nbits]=wavread (39。luyin139。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲S=fft(s)。 %傅里葉變換z12=filter(num12,den12,s)。 %IIR帶通濾波 sound(z12)。m12=fft(z12)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。subplot(2,2,2)。plot(abs(m12),39。r39。)。title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(s)。title(39。濾波前信號的波形39。)。subplot(2,2,4)。plot(z12)。title(39。濾波后的信號波形39。)。圖10 雙線性法IIR帶通濾波%**********************IIR高通濾波器***************************** [y,fs,nbits]=wavread (39。luyin139。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲S=fft(s)。 %傅里葉變換z13=filter(num13,den13,s)。 %IIR高通濾波sound(z13)。m13=fft(z13)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。subplot(2,2,2)。plot(abs(m13),39。r39。)。title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(s)。 title(39。濾波前信號的波形39。)。subplot(2,2,4)。plot(z13)。title(39。濾波后的信號波形39。)。圖11 雙線性法IIR高通濾波%========================窗函數(shù)法==========================%**************************FIR低通濾波器*************************[y,fs,nbits]=wavread (39。luyin139。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲 S=fft(s)。 %傅里葉變換z21=fftfilt(b21,s)。 %FIR低通濾波 sound(z21)。m21=fft(z21)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。subplot(2,2,2)。plot(abs(m21),39。r39。)。 title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(s)。 title(39。濾波前信號的波形39。)。subplot(2,2,4)。plot(z21)。 title(39。濾波后的信號波形39。)。圖12 窗函數(shù)法FIR低通濾波%**************************FIR帶通濾波器*************************[y,fs,nbits]=wavread (39。luyin139。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲 S=fft(s)。 %傅里葉變換z22=fftfilt(b22,s)。 %FIR帶通濾波 sound(z22)。m22=fft(z22)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。subplot(2,2,2)。plot(abs(m22),39。r39。)。title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(s)。title(39。濾波前信號的波形39。)。subplot(2,2,4)。plot(z22)。title(39。濾波后的信號波形39。)。圖13 窗函數(shù)法FIR帶通濾波%**************************FIR高通濾波器*************************[y,fs,nbits]=wavread (39。luyin139。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲 S=fft(s)。 %傅里葉變換z23=fftfilt(b23,s)。 %FIR高通濾波sound(z23)。m23=fft(z23)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。subplot(2,2,2)。plot(abs(m23),39。r39。)。title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(s)。title(39。濾波前信號的波形39。)。subplot(2,2,4)。plot(z23)。title(39。濾波后的信號波形39。)。 圖14 窗函數(shù)法FIR高通濾波 經(jīng)過以上的加噪聲處理后，可在Matlab中用函數(shù)sound對聲音進(jìn)行回放。其調(diào)用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)。可以察覺濾波前后的聲音有明顯的變化。七、設(shè)計(jì)系統(tǒng)界面 系統(tǒng)界面設(shè)計(jì)工具—GUI概述圖形用戶界面（graphical user interfaces ，GUI）則是由窗口、光標(biāo)、按鍵、菜單、文字說明等對象（objects）構(gòu)成的一個用戶界面。用戶通過一定的方法（如鼠標(biāo)或鍵盤）選擇、激活這些圖形對象，使計(jì)算機(jī)產(chǎn)生某種動作或變化，比如實(shí)現(xiàn)計(jì)算、繪圖等。在MATLAB中GUI是一中包含多種對象的圖形窗口，并為GUI開發(fā)提供一個方便高效的集成開發(fā)環(huán)境GUIDE。GUIDE主要是一個界面設(shè)計(jì)工具集，MAYLAB將所有GUI支持度控件都集成在這個環(huán)境中，并提供界面外觀、屬性和行為響應(yīng)方式的設(shè)置方法。GUIDE將設(shè)計(jì)好的GUI保存在一個FIG文件中，同時(shí)生成M文件框架。 FIG文件：包括GUI圖形窗口及其所有后裔的完全描述，包括所有對象屬性的屬性值。它是一個二進(jìn)制文件調(diào)用hsave課保存圖形窗口時(shí)將生車該文件。M文件包括GUI設(shè)計(jì)、控件函數(shù)以及定義為子函數(shù)的用戶控件回調(diào)函數(shù)，主要用于控制GUI展開時(shí)的各種特征。 GUI創(chuàng)建包括界面設(shè)計(jì)和控件編程兩部分，主要步驟如下。第一步：通過設(shè)置GUIDE應(yīng)用程序的選項(xiàng)來運(yùn)行GUIDE；第二步：使用界面設(shè)計(jì)編輯器進(jìn)行面設(shè)計(jì)；第三步：編寫控件行為響應(yīng)控制（即回調(diào)函數(shù)）代碼。首先我們新建一個GUI文件：File/New/GUI 如下圖所示：圖15 GUI創(chuàng)建界面選擇Blank GUI(Default)。布置如下圖所示；圖16 GUI開發(fā)環(huán)境布置好各控件以后，我們就可以來為這些控件編寫程序來實(shí)現(xiàn)兩數(shù)相加的功能了。最后，使用說明。單擊相應(yīng)的按鈕，便跳出對應(yīng)的界面，如按下按鈕“低通濾波器”，就得到如圖所示結(jié)果。圖17 GUI布局（1）選擇要分析的信號 圖18 GUI分析信號選擇沒有選擇信號回報(bào)如下錯誤??? Data must have one or two columns.Error in == sound at 41playsnd(y,fs,bits)。Error in ==