【正文】 要求，選擇串窗數(shù)類型（矩形窗、三角窗、漢寧窗、哈明窗、凱塞窗等），并估計窗口長度 N。原則是在保證阻帶衰減滿足要求的情況下，盡量選擇主瓣的窗函數(shù)。 （ 3）計算 h(n).。 接下來，我們根據(jù)語音信號的特點給出有關(guān)濾波器的技術(shù)指標(biāo)： ①低通濾波器的性能指標(biāo)： fp=1000Hz， fc=1200Hz， As=100db ,Ap=1dB ②高通濾波器的性能指標(biāo)： 3 fp=3500Hz， fc=4000Hz， As=100dB， Ap=1dB； ③帶通濾波器的性能指標(biāo)： fp1=1200Hz， fp2=3000hZ， fc1=1000Hz， fc2=3200Hz， As=100dB， Ap=1dB 在 Matlab中，可以利用函數(shù) fir1設(shè)計 FIR濾波器，利用函數(shù) butter,cheby1 和 ellip設(shè)計 IIR 濾波器，利用 Matlab 中的函數(shù) freqz 畫出各步步器的頻率響應(yīng)。如果缺省 window 參數(shù)，則fir1 默認(rèn)為哈明窗。 MATLAB 信號處理工具箱函數(shù) buttp buttor butter 是巴特沃斯濾波器設(shè)計函數(shù)，其有 5種調(diào)用格式，本課程設(shè)計中用到的是 [N,wc]=butter(N,wc,Rp,As,’ s’ ),該格式用于計算巴特沃斯模擬濾波器的階數(shù) N和 3dB 截止頻率 wc。我們用到的是 cheeby1 函數(shù)，其調(diào)用格式如下： [B,A]=cheby1(N,Rp,wpo,’ ftypr’ ) [B,A]=cheby1(N,Rp,wpo,’ ftypr’ ,’ s’ ) 函數(shù) butter,cheby1 和 ellip 設(shè)計 IIR 濾波器時都是默認(rèn)的雙線性變換法，所以在設(shè)計濾波器時只需要代入相應(yīng)的實現(xiàn)函數(shù)即可。 %=========================IIR 低通濾波器 ======================= wp=2*pi*Fp/Ft。 fp=2*Ft*tan(wp/2)。 [n11,wn11]=buttord(wp,ws,1,50,39。)。s39。 %求 S域的頻率響應(yīng)的參數(shù) [num11,den11]=bilinear(b11,a11,)。 %根據(jù)參數(shù)求出頻率響應(yīng) plot(w*8000*,abs(h))。用 butter設(shè)計 39。 4 圖 3 IIR 低通濾波器 %======================IIR 帶通 ========================== wp1=tan(pi*Fp1/Ft)。 ws1=tan(pi*Fs1/Ft)。 w=wp1*wp2/ws2。 wp=1。 [n12,wn12]=buttord(wp,ws,1,50,39。)。s39。 %求 S域的頻率響應(yīng)參數(shù) [num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw)。%雙線性變換實現(xiàn) S域到 Z域的轉(zhuǎn)換 [h,w]=freqz(num12,den12)。 axis([0 4000 0 ])。用 butter設(shè)計 39。 5 圖 4 IIR 帶通濾波器 %======================IIR 高通 ========================== Ft=8000。 Fs=3500。 %高通到低通濾波器參數(shù)轉(zhuǎn)換 ws1=tan(pi*Fs/Ft)。 ws=wp1*wp/ws1。s39。 %求模擬的低通濾波器階數(shù)和截止頻率 [b13,a13]=cheby1(n13,1,wn13,39。)。 %將 S域低通參數(shù)轉(zhuǎn)為高通的 [num13,den13]=bilinear(num,den,)。 plot(w*21000*,abs(h))。IIR高通濾波器 39。 legend(39。)。 函數(shù) fftfilt 用 的 是重 疊相 加 法實 現(xiàn) 線 性 卷積 的 計算 。其中， h 是系統(tǒng)單位沖擊響應(yīng)向量； x 是輸入序列向量； y 是系統(tǒng)的輸出序列向量； M是有用戶選擇的輸入序列的分段長度，缺省時，默認(rèn)的輸入向量的重長度 M=512。其中 xn 是輸入信號向量， yn輸出信號向量。再對 其進(jìn)行采樣，記住采樣頻率和采樣點數(shù)。 （ 1） y=wavread（ file） 功能說明：讀取 file 所規(guī)定的 wav 文件，返回采樣值放在向量 y中。 7 （ 3） y=wavread（ file， N） 功能說明：讀取錢 N點的采樣值放在向量 y中。 接下來，對語音信號 進(jìn)行采樣。 %把語音信號加載入 Matlab 仿真軟件平臺中 然后，畫出語音信號的時域波形，再對語音信號進(jìn)行頻譜分析。當(dāng) N 小于 xn的長度時， fft 函數(shù)計算 xn 的前 N個元素，忽略其后面的元素。 原始信號程序如下： [y1,fs,bits]=wavread(39。)。 %播放語音信