【正文】 (:,1)39。在Matlab中，可以利用矩形窗、三角窗、漢寧窗、漢明窗、布萊克曼窗、凱塞窗等設(shè)計(jì)FIR濾波器。 % 帶阻濾波器設(shè)計(jì)指標(biāo)fcd=(fpd+fsd)/2。不同的濾波器利用不同的函數(shù)對語音信號進(jìn)行濾波，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進(jìn)行濾波。在濾波器頻率特性的圖中可以看到：第一個圖是以db為單位的幅頻特性，第二圖是幅頻特性，第三個圖是濾波器的相頻特性，最后一個圖是濾波器的脈沖響應(yīng)。首先，更加深入理解了濾波器設(shè)計(jì)的各個關(guān)鍵環(huán)節(jié)，包括在什么情況下使用哪種方法設(shè)計(jì)FIR濾波器最好以及在選擇特定的窗函數(shù)進(jìn)行濾波器的設(shè)計(jì)時我們應(yīng)該怎樣確定其性能指標(biāo)；其次，更加深刻的認(rèn)識了語音原始信號與加噪后語音信號的波形及頻譜；再次，較大地提高了綜合運(yùn)用專業(yè)基礎(chǔ)知識及軟件設(shè)計(jì)能力，在一定程度上對自己的動手能力有很大的幫助。 % 單頻噪聲頻率，此參數(shù)可改t=0:1/fs:(N1)/fs。幅度39。subplot(2,2,3)。ylabel(39。wcu=fcu/fs*2*pi。xlabel(39。)。)。t39。subplot(3,2,3)。xlabel(39。)。f39。title(39。f39。plot(f,Y)。)。title(39。subplot(2,2,4)。w/pi39。plot(w/pi,db)。 % 帶阻濾波器設(shè)計(jì)指標(biāo)fcd=(fpd+fsd)/2。plot(f,Y)。)。xlabel(39。)。最后，在MATLAB軟件下編程時最好新建一個File文檔。在第一個圖中：第一幅圖和第二幅圖是原始語音信號的時域圖和頻域圖,第三幅圖和第四幅圖是加入頻率為2100的噪聲。 % 調(diào)用自編函數(shù)計(jì)算理想帶阻濾波器的脈沖響應(yīng)h_bs=w_black39。fsu=1950。 % 計(jì)算頻譜頻率范圍 濾波器設(shè)計(jì)本課程設(shè)計(jì)就是要設(shè)計(jì)一個濾波器慮掉加入的噪聲，使其恢復(fù)原始的語音信號。 % 輸入?yún)?shù)為文件的全路徑和文件名，輸出的第一個參數(shù)是每個樣本的值，fs是生成該波形文件時的采樣率，bits是波形文件每樣本的編碼位數(shù)sound(x,fs,bits)。 布萊克曼窗布萊克曼窗的時域形式可表示為：N(n) (23)它的頻域特性為：WRRRRR (24)其中為矩形窗函數(shù)的幅度頻率特性。確定窗口類型的主要依據(jù)是過渡帶寬和阻帶最小衰耗的指標(biāo)。 FIRDF的設(shè)計(jì)方法主要分為兩類：第一類是基于逼近理想濾波器特性的方法，包括窗函數(shù)法、頻率采樣法和等波紋最佳逼近法；第二類是最優(yōu)設(shè)計(jì)法。MATLAB可以幾乎所有的PC機(jī)和大型計(jì)算機(jī)上運(yùn)行，適用于Window、UNIX等多種系統(tǒng)平臺。正因?yàn)橛羞@些干擾噪聲的存在，接受者接受到的語音已不是原始的純凈語音信號，而是受噪聲干擾污染的帶噪聲語音信號。也可以調(diào)用函數(shù)sound聽濾波前后其語音信號是否帶有噪聲。MATLAB是一個為科學(xué)和工程計(jì)算機(jī)專門設(shè)計(jì)的交互式大型軟件，是一個可以完成各種精確計(jì)算和數(shù)據(jù)處理的、可視化的、強(qiáng)大的計(jì)算工具。就是說，它有（N—1）階極點(diǎn)在z = 0處，有（N—1）個零點(diǎn)位于有限z平面的任何位置因此H(z)是永遠(yuǎn)穩(wěn)定的。這時一般可用離散付里葉變換代替連續(xù)付里葉變換，求得近似值。肩峰值的大小決定了濾波器通帶內(nèi)的平穩(wěn)程度和阻帶內(nèi)的衰減，所以對濾波器的性能有很大的影響。下面是調(diào)用該語言信號以及加入噪聲的程序：[x,fs,bits]=wavread(39。Y=Y(1:N/2)。從畫出的圖中可以清楚的看見濾波器的幅頻和相頻特性。 % 計(jì)算布萊克曼窗設(shè)計(jì)該濾波器時需要的階數(shù)n=0:M1。 % 計(jì)算頻譜取前一半sound (x，fs，bits)。其次，在對加入噪聲后的信號進(jìn)行濾波時沒有使用正確的形式也不能濾去噪聲。設(shè)計(jì)的過程中，我通過查閱大量有關(guān)資料，與同學(xué)交流經(jīng)驗(yàn)和自學(xué)，并向老師請教等方式，使我學(xué)到了不少的東西，雖然有許多的辛酸，但是看到自己課程設(shè)計(jì)完成后心中的那份激動是無法用言語來形容的。 % 計(jì)算頻譜的譜線間隔f=0:deltaf:fs/2deltaf。)。 加入單頻干擾后的語音信號39。fpu=2000。.*hd_bs。濾波器的幅頻特性39。ylabel(39。 %對y進(jìn)行傅里葉變化Y_fil=Y_fil(1:length(Y_fil)/2)。)。y39。加噪聲后的頻域圖39。)。Y_fil39。)。ylabel(39。加噪聲后的時域圖39。x39。)。title(39。ylabel(39。wsu=fsu/fs*2*pi。title(39。時間（單位：s）39。原始語音信號39。sound (y,fs,bits)。在設(shè)置過程中有很多次因?yàn)樵O(shè)置的參數(shù)不合適而導(dǎo)致設(shè)計(jì)的濾波器不能慮出單頻噪聲信號。所以用布萊克曼窗設(shè)計(jì)的濾波器是符合要求的，也就是說該課程設(shè)計(jì)是成功的。最后在同一個圖中畫出原始信號的、加入噪聲的語音信號以及濾波后語音信號的頻域圖和時域圖。 % 計(jì)算上下邊帶中心頻率，和頻率間隔wcd=fcd/fs*2*pi。從六種窗函數(shù)的基本參數(shù)中我們可以得到旁瓣峰值n=5