【正文】 當(dāng)然，這樣對于我們也是一種警醒和提高，這樣的話，我們就會更加的嚴(yán)格要求自己去完成相應(yīng)的任務(wù)。所以說，這次課程設(shè)計還算比較成功。思路非常清晰，關(guān)鍵是要搞懂濾波器原理，和 MATLAB的使用方法即可。可以說，這次課程設(shè)計，基本上都是自己動手得來的成果。在這次課程設(shè)計中，我所做的題目是 — — 語音信號的采集、分析與處理 ； 我們一共用了兩個星期去完成這次課程設(shè)計。 這次課設(shè)使我認(rèn)識到 matlab軟件功能的強大，可以完成各類信號處理相關(guān)的功能，我的這個課題只是小小的一部分，今后我還需更加努力完善自己，用會，用精 matlab 來分析解決問題，將理論聯(lián)系實際，獲得更大提高。 通過這次課設(shè)，使我對濾波器有了更深的認(rèn)識，特別是濾波器參數(shù)對濾波器性能的影響，因為通帶截止頻率、阻 帶截止頻率、通帶衰減、阻帶衰減都要影響濾波器的階數(shù)，而濾波器的階數(shù)越大，其選頻特性就越好。并改變?yōu)V波器參數(shù)或特性 (低通、高通 )，實現(xiàn)不同的濾波要求。 為以后的工作和 學(xué)習(xí)打下基礎(chǔ)。 通過數(shù)字信號處理的課程設(shè)計，使我對信號的采集，處理，傳輸，顯示，存儲和分析等有一個系統(tǒng)的掌握和理解。)。 title(39。 subplot(2,1,2) y=fftshift(fft(x,datlen))。過濾后語音信號時域波形 39。 figure(3) %畫出濾波后的信號 subplot(2,1,1) 長沙學(xué)院 電子與通信工程系 26 plot(x)。 x=fftfilt(b,x1)。 figure(2)。high39。 wc=(wp+ws)/2/pi。 Bt=wpws。 %采樣頻率必須變大（ N1） ws=2*pi*fc/Fs。Ap=1。fp=5000。 t=0:1/fs:N/fs。原始語音信號頻譜 39。 subplot(2,1,2) plot((datlen/2:datlen/21)/datlen*fs,abs(y1))。原始語音信號時域波形 39。 figure(1) subplot(2,1,1) plot(x1)。 sound(x1,fs,bits)。)。 [x1,fs,bits]=wavread(39。)。 plot((datlen/2:datlen/21)/datlen*fs,abs(y)) title(39。)。 title(39。 %將 信號經(jīng)過濾波器 sound(X,fs,bits)。 %計算通帶截止頻率（關(guān)于 pi歸一化） hn=fir1(N,ws,wn)。 %凱塞窗濾波器階數(shù) N,(h(n)的長度為 N+1) wn= kaiser(N+1,alph)。 %計算過渡 帶寬 alph=*()。 wp=2*pi*fp/Fs。Fs=22050。As=100。 fp=1000。) N=length(x1)1。 title(39。)。 title(39。 y1=fftshift(fft(x1))。 datlen=length(x1)。C:\MATLAB7\work\39。 clear all。)。 title(39。)。 title(39。 %將信號經(jīng)過濾波器 sound(X,fs,bits)。)。 figure(2) plot(w*fs/(2*pi),abs(h)) title(39。 [b,a]=bilinear(bl,al,fs)。 [bap,aap]=zp2tf(z,p,k)。)。 [n,wc]=buttord(wp,ws,rp,as,39。 rp=1。 %巴特沃斯高通濾波器 wp=4800*2*pi。原始語音信號頻譜 39。 subplot(2,1,2) plot((datlen/2:datlen/21)/datlen*fs,abs(y1))。原始語音信號時域波形 39。 figure(1) subplot(2,1,1) plot(x1)。 sound(x1,fs,bits)。)。 [x1,fs,bits]=wavread(39。)。 title(39。)。 title(39。 %將信號經(jīng)過濾波器 sound(X,fs,bits)。)。 figure(2) plot(w*fs/(2*pi),abs(h)) title(39。 [b,a]=bilinear(bl,al,fs)。 [bap,aap]=zp2tf(z,p,k)。)。 [n,wc]=buttord(wp,ws,rp,as,39。 rp=1。) wp=1000*2*pi。 title(39。)。 title(39。 y1=fftshift(fft(x1))。 datlen=length(x1)。C:\MATLAB7\work\39。 clear all。 圖像 長沙學(xué)院 電子與通信工程系 18 用濾波器對信號進(jìn)行濾波并比較前后波形 要求用自己設(shè)計的各種濾波器分別對采集的信號進(jìn)行濾波，在 MATLAB 中， FIR 濾波器利用函數(shù) fftfilt 對信號進(jìn)行濾波， IIR 濾波器利用函數(shù) filter 對信號進(jìn)行濾波。 figure(2)。high39。 wc=(wp+ws)/2/pi。 Bt=wpws。 %采樣頻率必須變大（ N1） ws=2*pi*fc/Fs。Ap=1。fp=5000。 t=0:1/fs:N/fs。 %計算通帶截止頻率（關(guān)于 pi歸一化） hn=fir1(N,ws,wn)。 %凱塞窗濾波器階數(shù) N,(h(n)的長度為 N+1) wn= kaiser(N+1,alph)。 %計算過渡帶寬 alph=*()。 wp=2*pi*fp/Fs。Fs=22050。As=100。 fp=1000。 圖像 長沙學(xué)院 電子與通信工程系 16 有限低通數(shù)字濾波器 程序 N=length(x1)1。IIR高通濾波器頻率響應(yīng) 39。 [h,w]=freqz(b,a)。 [bl,al]=lp2hp(bap,aap,wc)。 [z,p,k]=buttap(n)。s39。 長沙學(xué)院 電子與通信工程系 15 as=100。 ws=10000*2*pi。 圖像 無限高通數(shù)字濾波器 程序 fs=22050。IIR低通濾波