【文章內容簡介】 %求出語音信號的長度noise=*randn(n,2)。 %隨機函數(shù)產生噪聲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)。 %隨機函數(shù)產生噪聲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)。 %隨機函數(shù)產生噪聲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)。 %隨機函數(shù)產生噪聲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)。 %隨機函數(shù)產生噪聲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)。 %隨機函數(shù)產生噪聲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高通濾波 經過以上的加噪聲處理后，可在Matlab中用函數(shù)sound對聲音進行回放。其調用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)?？梢圆煊X濾波前后的聲音有明顯的變化。七、設計系統(tǒng)界面 系統(tǒng)界面設計工具—GUI概述圖形用戶界面（graphical user interfaces ，GUI）則是由窗口、光標、按鍵、菜單、文字說明等對象（objects）構成的一個用戶界面。用戶通過一定的方法（如鼠標或鍵盤）選擇、激活這些圖形對象，使計算機產生某種動作或變化，比如實現(xiàn)計算、繪圖等。在MATLAB中GUI是一中包含多種對象的圖形窗口，并為GUI開發(fā)提供一個方便高效的集成開發(fā)環(huán)境GUIDE。GUIDE主要是一個界面設計工具集，MAYLAB將所有GUI支持度控件都集成在這個環(huán)境中，并提供界面外觀、屬性和行為響應方式的設置方法。GUIDE將設計好的GUI保存在一個FIG文件中，同時生成M文件框架。 FIG文件：包括GUI圖形窗口及其所有后裔的完全描述，包括所有對象屬性的屬性值。它是一個二進制文件調用hsave課保存圖形窗口時將生車該文件。M文件包括GUI設計、控件函數(shù)以及定義為子函數(shù)的用戶控件回調函數(shù)，主要用于控制GUI展開時的各種特征。 GUI創(chuàng)建包括界面設計和控件編程兩部分，主要步驟如下。第一步：通過設置GUIDE應用程序的選項來運行GUIDE；第二步：使用界面設計編輯器進行面設計；第三步：編寫控件行為響應控制（即回調函數(shù)）代碼。首先我們新建一個GUI文件：File/New/GUI 如下圖所示：圖15 GUI創(chuàng)建界面選擇Blank GUI(Default)。布置如下圖所示；圖16 GUI開發(fā)環(huán)境布置好各控件以后，我們就可以來為這些控件編寫程序來實現(xiàn)兩數(shù)相加的功能了。最后，使用說明。單擊相應的按鈕，便跳出對應的界面，如按下按鈕“低通濾波器”，就得到如圖所示結果。圖17 GUI布局（1）選擇要分析的信號 圖18 GUI分析信號選擇沒有選擇信號回報如下錯誤??? Data must have one or two columns.Error in == sound at 41playsnd(y,fs,bits)。Error in ==