【正文】 程序：[y,fs,nbits]=wavread(39。G:\Windows 39。)。figure(1)subplot(2,1,1)。plot(abs(y))。xlabel(39。n39。)。ylabel(39。y39。)。title(39。原語音信號時域波形39。)Y=fft(y,32768)。f=22050*(0:16384)/32768。subplot(2,1,2)。plot(f,abs(Y(1:16385)))。xlabel(39。f39。)。ylabel(39。Y39。)。title(39。原語音信號頻譜圖39。)t=0::。x1=*sin(2*pi*100*t)+*sin(2*pi*220*t)+*sin(2*pi*6000*t)。x=[rot90(x1,3),rot90(x1,3)]。y=y(1:20191,:)z=x+y。figure(2)subplot(2,1,1)。plot(z)。xlabel(39。n39。)。ylabel(39。z39。)。title(39。被污染語音信號時域波形39。)Z=fft(z,32768)。f=22050*(0:16384)/32768。subplot(2,1,2)。plot(f,abs(Z(1:16385)))。xlabel(39。f39。)。ylabel(39。Z39。)。title(39。被污染語音信號頻譜圖39。)wavplay(y,fs,39。sync39。)。wavplay(z,fs,39。sync39。)。wpz=[2*700/fs,2*4500/fs]。wsz=[2*500/fs,2*4700/fs]。wp=2*fs*tan(wpz*pi/2)。ws=2*fs*tan(wsz*pi/2)。rp=20。rs=50。[N,wc]=cheb2ord(wp,ws,rp,rs,39。s39。)。[b,a]=cheby2(N,rs,wc,39。s39。)。[bz,az]=bilinear(b,a,fs)。figure(3)Hk=freqz(bz,az,16384,fs)。f=22050*(0:16383)/32768。plot(f,20*log10(abs(Hk)))xlabel(39。f39。)。ylabel(39。Hk39。)。title(39。系統(tǒng)函數幅頻特性39。)figure(4)H=rot90(Hk)。Hk=[H fliplr(H)]。H=Hk39。Hk=[H H]。Yk=Z.*Hk。f=22050*(0:16382)/32768。subplot(2,1,2)。plot(f,abs(Yk(1:16383)))。xlabel(39。f39。)。ylabel(39。Yk39。)。title(39。濾波后語音信號頻譜圖39。)yk=ifft(Yk,32768)。subplot(2,1,1)。plot(abs(yk))。axis([0,25000,0,])。xlabel(39。n39。)。ylabel(39。yk39。)。title(39。濾波后語音信號時域波形39。)wavplay(abs(yk),fs)。 (4) hamming濾波器完整程序：[x1,fs,bits]=wavread(39。G:\Windows 39。)。subplot(2,2,1)。plot(x1)。title(39。原始聲音波形39。)。X1=fft(x1)。mag1=abs(X1)。subplot(2,2,2)。plot(mag1)。title(39。聲音信號頻譜39。)。t=1/fs。x2=x1+*cos(2*pi*10000*t/22050000)。sound(x2,fs,bits)。subplot(2,2,3)。plot(x2)。title(39。加噪聲聲音波形39。)。X2=fft(x2)。mag2=abs(X2)。subplot(2,2,4)plot(mag2)。title(39。加噪聲信號頻譜39。)。figure(2)。wpg=pi/2。wsg=pi/4。Btg=wpgwsg。N0g=ceil(*pi/Btg)。Ng=N0g+mod(N0g+1,2)。wcg=(wpg+wsg)/2/pi。hng=fir1(Ng1,wcg,39。high39。,hanning(Ng))。[magg,radg]=freqz(hng)。subplot(2,2,1)。plot(radg,20*log10(abs(magg)))。xlabel(39。Normalized Frequency(rad)39。)。ylabel(39。Normaliaed Magnitude(dB)39。)。title(39。FIR高通濾波39。)。yg=fftfilt(hng,x2)。Yg=fft(yg)。subplot(2,2,2)。plot(abs(Yg))。title(39。濾波后信號頻譜39。)。 xlabel(39。Frequency（hz）39。)。ylabel(39。Amplitude39。)。subplot(2,2,3)。plot(x2)。title(39。加噪聲聲音原始波形39。)。subplot(2,2,4)。plot(yg)。title(39。加噪聲聲音濾波后波形39。)。sound(yg,fs,bits)。%低通figure(3)。wpd=pi*。wsd=pi*。Btd=wsdwpd。N0d=ceil(*pi/Btd)。Nd=N0d+mod(N0d+1,2)。wcd=(wpd+wsd)/2/pi。hnd=fir1(Nd,wcd,hanning(Nd+1))。[magd,radd]=freqz(hnd)。subplot(2,2,1)。plot(radd,20*log10(abs(magd)))。xlabel(39。Normalized Frequency(rad)39。)。ylabel(39。Normaliaed Magnitude(dB)39。)。title(39。FIR低通濾波39。)。yd=fftfilt(hnd,x2)。Yd=fft(yd)。subplot(2,2,2)。plot(abs(Yd))。title(39。濾波后信號頻譜39。)。 xlabel(39。Frequency（hz）39。)。ylabel(39。Amplitude39。)。subplot(2,2,3)。plot(x2)。title(39。加噪聲聲音原始波形39。)。subplot(2,2,4)。plot(yd)。title(39。加噪聲聲音濾波后波形39。)。sound(yd,fs,bits)。%帶通figure(4)wpdai1=pi*。wpdai2=pi*。wsdai1=pi*。wsdai2=pi*。Btdai=wpdai1wsdai1。Mdai=ceil(*pi/Btdai)。wcdai=[(wpdai1+wsdai1)/2/pi,(wpdai2+wsdai2)/2/pi]。hndai=fir1(Mdai,wcdai,hanning(Mdai+1))。[magdai,raddai]=freqz(hndai)。subplot(2,2,1)。plot(raddai,20*log10(abs(magdai)))。xlabel(39。Normalized Frequency(rad)39。)。ylabel(39。Normaliaed Magnitude(dB)39。)。title(39。FIR帶通濾波39。)。ydai=fftfilt(hnd,x2)。Ydai=fft(ydai)。subplot(2,2,2)。plot(abs(Ydai))。title(39。濾波后信號頻譜39。)。 xlabel(39。Frequency（hz）39。)。ylabel(39。Amplitude39。)。subplot(2,2,3)。plot(x2)。title(39。加噪聲聲音原始波形39。)。subplot(2,2,4)。plot(ydai)。title(39。加噪聲聲音濾波后波形39。)。sound(ydai,fs,bits)。畢業(yè)設計（論文）原創(chuàng)性聲明和使用授權說明原創(chuàng)性聲明本人鄭重承諾：所呈交的畢業(yè)設計（論文），是我個人在指導教師的指導下進行的研究工作及取得的成果。盡我所知，除文中特別加以標注和致謝的地方外，不包含其他人或組織已經發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機構的學位或學歷而使用過的材料。對本研究提供過幫助和做出過貢獻的個人或集體，均已在文中作了明確的說明并表示了謝意。作 者 簽 名：　　 　 　　日　 期：　　 　 　　 指導教師簽名：　　 　 　　 日　　期：　　 　 　　 使用授權說明本人完全了解 大學關于收集、保存、使用畢業(yè)設計（論文）的規(guī)定，即：按照學校要求提交畢業(yè)設計（論文）的印刷本和電子版本；學校有權保存畢業(yè)設計（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務；學?？梢圆捎糜坝　⒖s印、數字化或其它復制手段保存論文；在不以贏利為目的前提下，學校可以公布論文的部分或全部內容。作者簽名：　　 　 　　 日　 期：　　 　 　　 學位論文原創(chuàng)性聲明本人鄭重聲明：所呈交的論文是本人在導師的指導下獨立進行研究所取得的研究成果。除了文中特別加以標注引用的內容外，本論文不包含任何其他個人或集體已經發(fā)表或撰寫的成果作品。對本文的研究做出重要貢獻的個人和集體，均已在文中以明確方式標明。本人完全意識到本聲明的法律后果由本人承擔。作者簽名： 日期： 年 月 日學位論文版權使用授權書本學位論文作者完全了解學校有關保留、使用學位論文的規(guī)定，同意學校保留并向國家有關部門或機構送交論文的復印件和電子版，允許論文被查閱和借閱。本人授權　　 　 　大學可以將本學位論文的全部或部分內容編入有關數據庫進行檢索，可以采用影印、縮印或掃描等復制手段保存和匯編本學位論文。涉密論文按學校規(guī)定處理。作者簽名： 日期： 年 月 日導師簽名： 日期： 年 月 日注 意 事 項（論文）的內容包括：1）封面（按教務處制定的標準封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關鍵詞4）外文摘要、關鍵詞 5）目次頁（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結論7）參考文獻8）致謝9）附錄（對論文支持必要時）：理工類設計（論文）正文字數不少于1萬字（不包括圖紙、程序清單等）。：任務書、開題報告、外文譯文、譯文原文（復印件）。、圖表要求：1）文字通順，語言流暢，書寫字跡工整，打印字體及大小符合要求，無錯別字，不準請他人代寫2）工程設計類題目的圖紙，要求部分用尺規(guī)繪制，部分用計算機繪制，所有圖紙應符合國家技術標準規(guī)范。圖表整潔，布局合理，文字注釋必須使用工程字書寫，不準用徒手畫3）畢業(yè)論文須用A4單面打印，論文50頁以上的雙面打印4）圖表應繪制于無格子的頁面上5）軟件工程類課題應有程序清單，并提供電子文檔1）設計（論文）2）附件：按照任務書、開題報告、外文譯文、譯文原文（復印件）次序裝訂3）其它