【正文】 LACKMANHARRIS109dbBOHMANWINCHEBWIN113dbNUTTALLWIN108dbTUKEYWIN22db FIR濾波器　　FIR(Finite Impulse Response)濾波器：有限長單位沖激響應(yīng)濾波器，是數(shù)字信號(hào)處理系統(tǒng)中最基本的元件，它可以在保證任意幅頻特性的同時(shí)具有嚴(yán)格的線性相頻特性，同時(shí)其單位抽樣響應(yīng)是有限長的，因而濾波器是穩(wěn)定的系統(tǒng)。圖 設(shè)置語音屬性 語音信號(hào)的頻譜分析 然后在Matlab軟件平臺(tái)下，利用函數(shù)wavread對(duì)語音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。 Y=Y(1:N/2)。wc1=(fs1+fp1)/2。Y_fil=Y_fil(1:N/2)。在這次課程設(shè)計(jì)中，有許多的知識(shí)都不是很懂，通過課設(shè)中查閱資料等，我拓寬了知識(shí)面，增長了見識(shí)。 % 應(yīng)該可以明顯聽出有尖銳的單頻嘯叫聲X=abs(fft(x))。)。xlabel(39。幅度譜39。 % 將Hz為單位的模擬頻率換算為rad為單位的數(shù)字頻率 ws1=fs1/fs*2*pi。db39。)。 % 計(jì)算頻譜取前一半 subplot(3,2,1)。幅度譜39。plot(f,Y)。幅度39。% 理想低通濾波器計(jì)算% % [hd] = ideal_lp(wc,M)% hd = 0 to M1之間的理想脈沖響應(yīng)% wc = 截止頻率(弧度) % M = 理想濾波器的長度%alpha = (M1)/2。% grd = diff(pha)。db = 20*log10((mag+eps)/max(mag))。title(39。xlabel(39。)。xlabel(39。)。 subplot(2,2,3)。xlabel(39。 % 帶阻濾波器設(shè)計(jì)指標(biāo) df=min((fs1fp1),(fp2fs2))。xlabel(39。) 。幅度39。 % 計(jì)算時(shí)間范圍，樣本數(shù)除以采樣頻率 x=x39。在編程過程中遇到和很多不懂得問題，但是在老師和同學(xué)們的幫助下，問題基本上都得到了解決。 % 調(diào)用自編函數(shù)計(jì)算濾波器的頻率特性： 用濾波器對(duì)信號(hào)進(jìn)行濾波在將濾波器設(shè)計(jì)好后，我們用自己設(shè)計(jì)的帶阻濾波器對(duì)采集的語音信號(hào)進(jìn)行濾波。 Rp=1。 %在原信號(hào)上添加噪聲sound(y,fs,bits)。圖 進(jìn)入錄音機(jī)界面然后點(diǎn)擊最右邊的紅色圓心的標(biāo)志，即可進(jìn)行錄音。為了設(shè)計(jì)簡(jiǎn)單方便，通常選擇為具有片段常數(shù)特性的理想濾波器。數(shù)字信號(hào)處理技術(shù)及設(shè)備具有靈活、精確、抗干擾強(qiáng)、設(shè)備尺寸小、造價(jià)低、速度快等突出優(yōu)點(diǎn)，這些都是模擬信號(hào)處理技術(shù)與設(shè)備所無法比擬的。數(shù)字濾波器的種類很多，但總的來說可以分成兩大類，一類是經(jīng)典濾波器，另一類可稱為現(xiàn)代濾波器。要求：（1）濾波器指標(biāo)必須符合工程實(shí)際。（2）課程設(shè)計(jì)附件（主要是源程序）。與不同信源相同濾波方法的同學(xué)比較各種信源的特點(diǎn)，與相同信源不同濾波方法的同學(xué)比較各種濾波方法性能的優(yōu)劣。而信息是信號(hào)的具體內(nèi)容。因此，F(xiàn)IR濾波器在通信、圖像處理、模式識(shí)別等領(lǐng)域都有著廣泛的應(yīng)用。采集完成后在信號(hào)中加入一個(gè)單頻噪聲。 % 截取前半部分deltaf=fs/N。wc2=(fs2+fp2)/2。 % 計(jì)算頻譜取前一半： 音樂信號(hào)的時(shí)域波形和頻譜圖比較 結(jié)果分析，原信號(hào)與濾波去噪信號(hào)的時(shí)域圖基本相似，只有邊緣部分有點(diǎn)差異；原信號(hào)與濾波去噪信號(hào)的頻譜圖波形也大致相似。最終在和同學(xué)們共同進(jìn)步中圓滿完成課程設(shè)計(jì)任務(wù)。 Y=abs(fft(y))。axis([0 ,300 100])。時(shí)間（單位：s）39。)。ws2=fs2/fs*2*pi。)。ylabel(39。plot(t,x)。)。xlabel(39。)。n = [0:1:(M1)]。% grd = [grd(1) grd]。mag = abs(H)。)。plot(t,y_fil)。幅度39。plot(f,X)。w(n)39。 axis([0 ,0 2])。plot(w/pi,db)。As=15。plot(f,Y)。原始語音信號(hào)幅度譜39。ylabel(39。 % 單頻噪聲頻率，此參數(shù)可自己設(shè)計(jì)t=0:1/fs:(N1)/fs。 語音文件添加好之后就是對(duì)語音添加噪聲，進(jìn)行頻譜分析了，這些基本上都是以前上Matlab課程中接觸到的一些基本知識(shí)，只是很多函數(shù)是老師給的，也給我們省下了不少時(shí)間。 % 用窗口法計(jì)算實(shí)際濾波器脈沖響應(yīng) [db,mag,pha,grd,w]=freqz_m(h,1)。fs2=fn+50。+(fn*2*pi*t)。：開 始對(duì)語音信號(hào)進(jìn)行采樣，并加入一個(gè)單頻噪聲對(duì)語音信號(hào)進(jìn)行頻譜分析，畫出時(shí)域和頻域波形圖 用三角窗設(shè)計(jì)FIR濾波器畫出其頻率響應(yīng)用FIR濾波器對(duì)語音信號(hào)進(jìn)行濾波比較濾波前后語音信號(hào)的波形及頻譜結(jié) 束在左下角點(diǎn)擊開始/程序/附件/娛樂/錄音機(jī)。設(shè)希望逼近的濾波器的頻率響應(yīng)為，其單位脈沖響應(yīng)用表示。2 設(shè)計(jì)原理數(shù)字信號(hào)處理的算法需要利用計(jì)算機(jī)或?qū)Ｓ锰幚碓O(shè)備如數(shù)字信號(hào)處理器（DSP）和專用集成電路（ASIC）等。隨著現(xiàn)代通信的數(shù)字化,數(shù)字濾波器變得更加重要。與不同信源相同濾波方法的同學(xué)比較各種信源的特點(diǎn)，與