【正文】 57答辯情況： 自我陳述、回答問題的正確性、用語準(zhǔn)確性、邏輯思維、是否具有獨(dú)到見解等。103論文選題是否先進(jìn)，是否具有前沿性或前瞻性。當(dāng)然，我也要感謝計(jì)算機(jī)工程學(xué)院，感謝他們給我提供這次實(shí)習(xí)的機(jī)會(huì)。解決不了的問題，我們就像老師和其他同學(xué)虛心請(qǐng)教。指導(dǎo)老師嚴(yán)謹(jǐn)求實(shí)的治學(xué)態(tài)度，踏實(shí)堅(jiān)韌的工作精神，值得我學(xué)習(xí)。設(shè)計(jì)出了一個(gè)簡(jiǎn)易的用戶應(yīng)用界面，可以讓人實(shí)現(xiàn)界面操作，更加方便的進(jìn)行語音的頻譜分析與加噪濾波處理。其實(shí)，這正是干擾所造成的。那么，就可以完全利用數(shù)字信號(hào)處理的知識(shí)來解決語音及加噪處理問題。作為存儲(chǔ)于計(jì)算機(jī)中的語音信號(hào)，其本身就是離散化了的向量，我們只需將這些離散的量提取出來，就可以對(duì)其進(jìn)行處理了。圖17 系統(tǒng)運(yùn)行界面本設(shè)計(jì)圓滿的完成了對(duì)加噪聲語音信號(hào)的讀取與打開，與課題的要求十分相符；也較好的完成了對(duì)原始語音信號(hào)和加噪聲后的語音信號(hào)的頻譜分析，通過fft變換，得出了語音信號(hào)的頻譜圖；在濾波這一部分，課題主要是從巴特沃斯濾波器入手來設(shè)計(jì)濾波器等入手，實(shí)現(xiàn)了預(yù)期的濾波效果；完成了界面的設(shè)計(jì)，但也存在許多的不足，只是很勉強(qiáng)的達(dá)到了打開語音文件、顯示已定濾波前后的波形圖，界面沒有用到其它的一些控件，如列表框、下拉菜單等等，GUI的界面設(shè)計(jì)功能沒能得到充分的發(fā)揮。 界面設(shè)計(jì)及使用說明首先我們新建一個(gè)GUI文件：File/New/GUI 如下圖所示：圖15 GUI創(chuàng)建界面選擇Blank GUI(Default)。它是一個(gè)二進(jìn)制文件調(diào)用hsave課保存圖形窗口時(shí)將生車該文件。在MATLAB中GUI是一中包含多種對(duì)象的圖形窗口，并為GUI開發(fā)提供一個(gè)方便高效的集成開發(fā)環(huán)境GUIDE。其調(diào)用格式：sound(y,Fs)，sound(y)和sound(y,Fs,bits)。title(39。濾波前信號(hào)的波形39。)。r39。g39。sound(z23)。 %求出語音信號(hào)的長度noise=*randn(n,2)。圖13 窗函數(shù)法帶通濾波%*************************高通濾波器*****************************[y,fs,nbits]=wavread (39。plot(z22)。title(39。濾波后信號(hào)的頻譜39。plot(abs(m22),39。title(39。 %求濾波后的信號(hào)subplot(2,2,1)。 %語音信號(hào)加入噪聲 S=fft(s)。)。濾波后的信號(hào)波形39。)。subplot(2,2,3)。)。)。g39。sound(z21)。 %求出語音信號(hào)的長度noise=*randn(n,2)。圖11 雙線性法高通濾波%========================窗函數(shù)法==========================%**************************低通濾波器****************************[y,fs,nbits]=wavread (39。plot(z13)。 title(39。濾波后信號(hào)的頻譜39。plot(abs(m13),39。title(39。 %求濾波后的信號(hào)subplot(2,2,1)。 %語音信號(hào)加入噪聲S=fft(s)。)。濾波后的信號(hào)波形39。)。subplot(2,2,3)。)。)。g39。sound(z12)。 %求出語音信號(hào)的長度noise=*randn(n,2)。圖9 雙線性法低通濾波%**********************帶通濾波器*****************************[y,fs,nbits]=wavread (39。plot(z11)。濾波前信號(hào)的波形39。grid。)。grid。)。m11=fft(z11)。 %隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise。speech39。其中，h是系統(tǒng)單位沖擊響應(yīng)向量；x是輸入序列向量；y是系統(tǒng)的輸出序列向量；M是有用戶選擇的輸入序列的分段長度，缺省時(shí)，默認(rèn)的輸入向量的重長度M=512。axis([3000 6000 0 ])。 %得到頻率響應(yīng)plot(w*12000*,abs(h))。b23=fir1(n23,wn23,39。fpts=[ws wp]。用butter設(shè)計(jì)39。%雙線性變換實(shí)現(xiàn)S域到Z域的轉(zhuǎn)換[h,w]=freqz(num22,den22)。s39。[n22,wn22]=buttord(wp,ws,1,50,39。 w=wp1*wp2/ws2。圖6 FIR低通濾波器%************************FIR帶通濾波器************************wp1=tan(pi*Fp1/Ft)。 %得到頻率響應(yīng)plot(w/pi,abs(h))。dev=[p s]。p=110.^(rp/20)。用cheby1設(shè)計(jì)39。title(39。 %求S域的頻率響應(yīng)的參數(shù)[num,den]=lp2hp(b13,a13,wn13)。)。wp=1。Fp=4000。legend(39。%將S域低通參數(shù)轉(zhuǎn)為帶通的[num12,den12]=bilinear(num2,den2,)。 %求低通濾波器階數(shù)和截止頻率[b12,a12]=butter(n12,wn12,39。ws=(wp1*wp2w.^2)/(bw*w)。ws2=tan(pi*Fs2/Ft)。)。 %雙線性變換實(shí)現(xiàn)S域到Z域的變換 [h,w]=freqz(num11,den11)。 %求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,39。fs=2*Fs*tan(wp/2)。下面我們將給出FIR和IIR數(shù)字濾波器的主要程序。其中可選的窗函數(shù)有Rectangular Barlrtt Hamming Hann Blackman窗，其相應(yīng)的都有實(shí)現(xiàn)函數(shù)。（4）加窗得到設(shè)計(jì)結(jié)果。先按照阻帶衰減選擇窗函數(shù)類型。（4）用雙線性變換法，模擬濾波器系統(tǒng)函數(shù)轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。為了克服之一缺點(diǎn)，可以采用雙線性變換法。（4）設(shè)計(jì)模擬低通濾波器。 對(duì)于數(shù)字高通、帶通濾波器的設(shè)計(jì)，通用方法為雙線性變換法。IIR濾波器設(shè)計(jì)方法有間接法和直接法，間接法是借助于模擬濾波器的設(shè)計(jì)方法進(jìn)行的。 title(39。)。subplot(2,1,1)。N = length (y) 。Randn函數(shù)有兩種基本調(diào)用格式：Randn(n)和Randn（m,n）,前者產(chǎn)生nn服從標(biāo)準(zhǔn)高斯分布的隨機(jī)數(shù)矩陣，后者產(chǎn)生mn的隨機(jī)數(shù)矩陣。我們可以直接應(yīng)用兩個(gè)函數(shù)：一個(gè)是WGN，另一個(gè)是AWGN。plot(abs(Y))。title(39。 %回放語音信號(hào)N= length (y) 。其程序如下：[y,fs,nbits]=wavread (39。 %把語音信號(hào)加載入Matlab 仿真軟件平臺(tái)中然后，畫出語音信號(hào)的時(shí)域波形，再對(duì)語音信號(hào)進(jìn)行頻譜分析。（3）y=wavread（file，N）功能說明：讀取錢N點(diǎn)的采樣值放在向量y中。再對(duì)其進(jìn)行采樣，記住采