【正文】 語音信號進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)，采樣率8000Hz，單通道，格式*.wav。用窗函數(shù)wd(n)將hd(n)截斷，并進(jìn)行加權(quán)處理，得到如果要求線性相位特性， 則h(n)還必須滿足：根據(jù)上式中的正、 負(fù)號和長度N的奇偶性又將線性相位FIR濾波器分成四類。脈沖響應(yīng)不變法的主要缺點(diǎn)是產(chǎn)生頻率響應(yīng)的混疊失真。這樣就使S平面與Z平面建立了一一對應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖1所示。首先,把z=ejω，可得 (4)即S平面的虛軸映射到Z平面的單位圓。雙線性變換法優(yōu)缺點(diǎn):雙線性變換法與脈沖響應(yīng)不變法相比，其主要的優(yōu)點(diǎn)是避免了頻率響應(yīng)的混疊現(xiàn)象。由圖2看出，在零頻率附近，模擬角頻率Ω與數(shù)字頻率ω之間的變換關(guān)系接近于線性關(guān)系；但當(dāng)Ω進(jìn)一步增加時，ω增長得越來越慢，最后當(dāng)Ω→∞時，ω終止在折疊頻率ω=π處，因而雙線性變換就不會出現(xiàn)由于高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。圖3雙線性變換法幅度和相位特性的非線性映射對于分段常數(shù)的濾波器，雙線性變換后，仍得到幅頻特性為分段常數(shù)的濾波器，但是各個分段邊緣的臨界頻率點(diǎn)產(chǎn)生了畸變，這種頻率的畸變，可以通過頻率的預(yù)畸來加以校正。對語音信號進(jìn)行快速傅立葉變換，在一個窗口同時畫出信號的時域波形圖和頻譜圖，分析語音信號的頻譜特點(diǎn)[y,fs,nbits]=wavread(39。%回放語音信號 N=length(y)。title(39。plot(abs(Y))。 [x1,Fs,bits]=wavread(39。%設(shè)置頻譜的間隔，分辨率 ，這里保證了x軸的點(diǎn)數(shù)必須和y軸點(diǎn)數(shù)一致fs=Fs。wp1=2*pi*fp1/fs。M1=ceil(()/(*BF1))+1。 %求凱澤窗窗函數(shù)b1=fir1(M1,wc1/pi,Window)。調(diào)用默認(rèn)的N時，其值是512。 x1_low = filter(b1,1, x1)。plot(x1_low)。 subplot(2,1,2)。)。 derta_Fs = Fs/length(x1)。 fs2=2700。wc2=(wp2+ws2)/2。 Window=(kaiser(N2,beta2))。figure(4)。)。subplot(211)。)。信號經(jīng)過FIR高通濾波器（頻域）39。fp3=[1200,3000]。 BF3=wp3(1)ws3(1)。beta3=*()。,Window)。FIR帶通濾波器的頻率響應(yīng)39。figure(7)。信號經(jīng)過FIR帶通濾波器(時域)39。title(39。39。fp1i=1000。 Rp1i=1。Ws1i=2/Ts*tan(ws1i/2)。 %計算模擬濾波器的最小階數(shù) [B1i,A1i]=cheby1(N1i,Rp1i,Wn1i,39。 %運(yùn)用雙線性變換法得到數(shù)字濾波器傳遞函數(shù) figure(10)。)。 subplot(211)。)。信號經(jīng)過IIR低通濾波器（頻域）39。)。%通帶截止頻率fs2i=2700。Ws2i2=2*pi*fs2i/fs。 %歸一化wshi=Wp2i/Ws2i。%設(shè)計濾波器[B2i,A2i]=cheby1(N2i,Rp2i,Wn2i,39。[bz2i,az2i]=bilinear(numi2,deni2,)figure(12)。)。 subplot(211)。)。信號經(jīng)過IIR高通濾波器（頻域）39。)。fdp2i=3000。 Rs3i=100。wds2i=2*pi*fds2i/fs。 Wds2i=tan(wds2i/2)。 W0 = (Wdp1i*Wdp2i)^。 [N3i,Wn3i]=buttord(wp3i,ws3i,Rp3i,Rs3i,39。s39。 %雙線性變換 [H,W] = freqz(bz3i,az3i,512)。grid。 title(39。 plot(W/pi,angle(H))。),ylabel(39。%對信號進(jìn)行帶通濾波 ,Y = filter(B,A,X) ，輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A 提供濾波器系數(shù)，B為分子， A為分母 sound(x1_daitongi,Fs,bits)。title(39。 plot([Fs/2:derta_Fs: Fs/2derta_Fs],abs(fftshift(fft(x1_daitongi))))。通過這兩周的課程設(shè)計，自己能夠更好的使用MATLAB進(jìn)行程序設(shè)計，對濾波器有了進(jìn)一步的了解，初步學(xué)會設(shè)計簡單的濾波器。通過這次課設(shè)，使我對濾波器有了更深的認(rèn)識，特別是濾波器參數(shù)對濾波器性能的影響，因?yàn)橥◣Ы刂诡l率、阻帶截止頻率、通帶衰減、阻帶衰減都要影響濾波器的階數(shù)，而濾波器的階數(shù)越大，其選頻特性就越好。低通濾波器產(chǎn)生的信號相對于原信號變化不大，是因?yàn)樵O(shè)定的選頻范圍是原信號的頻率段中較集中的部分，所以與此同時，高通和帶通濾波器產(chǎn)生的信號相對于原信號的幅度就會大大衰減，其原因也是在于高通和帶通濾波器的選頻范圍是原信號頻率范圍中較少的一部分。這次課設(shè)對我?guī)椭艽?，讓我認(rèn)識到matlab軟件功能的強(qiáng)大，通過實(shí)踐結(jié)合理論，自己對數(shù)字信號處理有更深一步的了解。程佩青. 《數(shù)字信號處理教程》北京：清華大學(xué)出版社出版，2001張威 編. 《MATLAB基礎(chǔ)與編程入門》. 西安：西安電子科技大學(xué)出版社，2005董少平，陳示耕，王洋. 《數(shù)字信號處理基礎(chǔ)》.,2010:1381441.《數(shù)字信號處理》丁玉美，高西全等編著，西安：西安電子科技大學(xué)出版社1