【正文】 )。加躁信號頻域波形subplot(3,1,1)。39。plot(y)。%plot(w*21000*,abs(h))。%求模擬的低通濾波器階數(shù)和截止頻率[b13,a13]=cheby1(n13,1,wn13,39。%隨機(jī)函數(shù)產(chǎn)生噪聲si=y+Noise。%求出語音信號的長度%sound(y,fs,nbits)。[y,fs,nbits]=wavread1.首先根據(jù)通帶截頻確定wc wc=wp 。=1/[1+e2C2N(w/wc)]這里，ε是紋波系數(shù)，用來描述波動得大小。如圖11所示：在通帶內(nèi)：title(39。subplot(3,1,2)。title(39。subplot(3,1,2)。sound(f1,45000)。[z,p,k]=buttap(M)。ws=2*pi*fs/Fs。fc=20000。c:\39。lg[（）/() ]/2lg(ws/wp)可以得到濾波器的階數(shù)N由公式wp/()1/2N163。H（jw）189。 信號提?。和ㄟ^圖形用戶界面上的菜單功能按鍵采集電腦上的一段音頻信號，完成音頻信號的頻率，幅度等信息的提取，并得到該語音信號的波形圖。 信號調(diào)整：在設(shè)計(jì)的用戶圖形界面下對輸入的音頻信號進(jìn)行各種變化，如變化幅度、改變頻率等操作，以實(shí)現(xiàn)對語音信號的調(diào)整。2=1/[1+(w/wc)2N] 其中，wc是3dB截止頻率，N為濾波器階數(shù)。wc163。)。%載波頻率%Noise=modulate(y,fc,fs,39。wap=2*Fs*tan(wp/2)。[bp,ap]=zp2tf(z,p,k)。plot(si)。濾波信號時域波形39。plot(abs(SI))。濾波信號頻域波形39。在阻帶中：其中 為通帶截止頻率,為阻帶截止頻率，為通帶最大衰減, 為阻帶最大衰減。ε一般在0～1之間，ε=0表示沒有波動，ε=1表示波動達(dá)到半功率點(diǎn)。由通帶衰減Ap確定e e=()求階數(shù)N，根據(jù)公式N179。(39。%原信號聲音Y=fft(y)。%高通到低通濾波器參數(shù)轉(zhuǎn)換ws1=tan(pi*Fs/Ft)。s39。title(39。F1=fft(f1)。title(39。)。plot(abs(Y))。39。plot(w*21000*,abs(h))。切比雪夫?yàn)V波器是使通帶內(nèi)誤差分布均勻的濾波器。(39。%原信號聲音Y=fft(y)。s39。subplot(2,1,1)。|H(z)|39。)。sound(f1,fs)。plot(si)。title(39。)。subplot(3,1,3)。用MATLAB設(shè)計(jì)橢圓濾波器可以大大簡化設(shè)計(jì)過程。c:\39。)。Ap=1。b=fir1(N,ws,wn)。title(39。subplot(3,1,3)。原信號頻域波形39。plot(abs(X))。六、幾類濾波器的比較分析比較結(jié)果：巴特沃斯濾波器的頻率特性曲線,無論在通帶還是阻帶都是頻率的單調(diào)增函數(shù).,因此,當(dāng)通帶邊界處滿足指標(biāo)要求時,可以得出,切比雪夫?yàn)V波器比巴特沃斯濾波器有較窄的過渡帶。為了設(shè)計(jì)簡單方便，通常選擇為具有片段常數(shù)特性的理想濾波器。在課程設(shè)計(jì)剛剛開始的時候，覺得很無助，于是請教同學(xué)、上網(wǎng)搜資料、去圖書館查找等等，終于功夫不負(fù)有心人，在不斷的資料搜尋當(dāng)中我漸漸的明白了我們將要做的設(shè)計(jì)的基本原理和相關(guān)知識。通過實(shí)際問題分析了離散經(jīng)典數(shù)字濾波器分類和設(shè)計(jì)方法。八、 參考文獻(xiàn) [1] 程佩青．?dāng)?shù)字信號處理教程〔M〕．北京：清華大學(xué)出版社，2002． [2] 樓順天、李博菡. 基于MATLAB的系統(tǒng)分析與設(shè)計(jì)－：西安電子［3］ Sanjit 著，孫洪，余翔宇 譯，《數(shù)字信號處理實(shí)驗(yàn)指導(dǎo)書》，電子工業(yè)出版社，2005 年。其設(shè)計(jì)過程都是由模擬濾波器的系統(tǒng)函數(shù)去變換出相應(yīng)的數(shù)字濾波器的系統(tǒng)函數(shù)。窗函數(shù)設(shè)計(jì)法就是截取為有限長的一段因果序列，并用合適的窗函數(shù)進(jìn)行加權(quán)作為FIRDF的單位脈沖響應(yīng)h(n)。 切比雪夫?yàn)V波器的 在通帶范圍內(nèi)是等幅起伏的，所以同樣的通帶衰減，其階數(shù)較巴特沃茲濾波器要小。濾波信號頻域波形39。subplot(3,1,2)。title(39。)39。freqz(b,1)。 wp=2*pi*fp/fs。 %隨機(jī)函數(shù)產(chǎn)生噪聲si=y+Noise。N = length (y) 。因此，系統(tǒng)函數(shù)為 ()由式（）可得濾波器的結(jié)構(gòu)見圖1圖1 FIR濾波器結(jié)構(gòu)示意圖由于FIR濾波器的單位脈沖響應(yīng)是一個有限長序列，是的次多項(xiàng)式，它在平面上有個