濾波器實驗5基于matlab的數(shù)字濾波器設(shè)計-在線瀏覽
2025-01-17 17:23本頁面
【正文】 b)為凱塞窗設(shè)計的頻率特性。 nn=[0:1:29]。 hd=sin(*pi*(nnalfa))./(pi*(nnalfa))。 [h1,w1]=freqz(h,1)。 [h1,w1]=freqz(b,1)。 axis([0,1,80,10])。 xlabel(39。) 。幅度 /dB39。 還可以使用 [n,Wn,beta,ftype] = kaiserord(f,a,dev)函數(shù)來估計濾波器階數(shù)等，得到凱塞窗濾波器： fcuts = [ ]。 devs = [ 10^()]。 %計算出凱塞窗 N， beta 的值 hh = fir1(n,Wn,ftype,kaiser(n+1,beta),39。)。 實際中，一般調(diào)用 MATLAB信號處理工具箱函數(shù) remezord 來計算等波紋濾波器階數(shù) N和加權(quán)函數(shù) W(ω )，調(diào)用函數(shù) remez可進行等波紋濾波器的設(shè)計，直接求 出濾波器系數(shù)。 例 2 利用雷米茲交替算法 設(shè)計等波紋濾波器 ，設(shè)計一個線性相位低通 FIR 數(shù)字濾波器，其指標為：通帶邊界頻率 fc=800Hz，阻帶邊界 fr=1000Hz，通帶波動 阻帶最小衰減 At=40dB，采樣頻率 fs=4000Hz。 mval=[1 0]。 fs=4000。 b=remez(N,fpts,mag,wt)。 plot(w*2020/pi,20*log10(abs(h)))。 xlabel(39。) 。幅度 /dB39。 一、 實驗內(nèi)容： 利用 MATLAB 編程設(shè)計一個數(shù)字帶通濾波器，指標要求如下： 通帶邊緣頻率： ? ??P ， ? ??P ，通帶峰值起伏：][1dBp ?? 。 分別用 窗函數(shù)法和等波紋濾波器法設(shè)計 兩種 FIR數(shù)字濾波器。 11）用窗函數(shù)法實現(xiàn)： 調(diào)用函數(shù) [n,wn,bta,ftype]=kaiserord（ f， a， dev， fs） 參數(shù)： f=[ ]為對應(yīng)數(shù)字頻率 ? ??S ， ? ??P ，? ??P ， ? ??S a=[0 1 0]為由 f指定的各個頻帶上的幅值向量，一般只有 0和 1表示；和 f 長度關(guān)系為（ 2*a的長度） — 2=（ f 的長度） devs=[ ]用于指定各個頻帶輸出濾波器的頻率響應(yīng)與其期望幅值之間的最大輸出誤差或偏差，長度與 a相等，計算公式：阻帶衰減誤差 = ，通帶波動衰減誤差 = fs缺省值為 2HZ [n,wn,bta,ftype]=kaiserord([ ],[0 1 0],[ ])。noscale39。 [hh1,w1]=freqz(h1,1,256)。歸一化頻率 w39。ylabel(39。)。歸一化頻率 w39。ylabel(39。)。 畫圖得： 12）用等波紋法設(shè)計： 調(diào)用函數(shù) [n,fpts,mag,wt]=remezord(f， a， dev) f=[ ] a=[0 1 0] dev=[ ] 其含義同 函數(shù) [n,wn,bta,ftype]=kaiserord（ f， a， dev， fs）中的參數(shù)相同。%用 remezord 函數(shù)估算出 remez 函數(shù)要用到的階 n、歸一化頻帶邊緣矢量 fpts、頻帶內(nèi)幅值響應(yīng)矢量 mag 及加權(quán)矢量 w，使 remez函數(shù)設(shè)計出的濾波器滿足 f、 a 及 dev 指定的性能要求。%設(shè)計出等波紋濾波器 [hh2,w2]=freqz(h2,1,256)。歸一化頻率 w39。ylabel(39。)。歸一化頻率 w39。ylabel(39。)。直接得圖： 二、 對課本作業(yè) 。 hd=sin(*pi*(n29))./(pi*(n29))。 h1=2*cos(pi*(n29)/2).*hd.*win。 figure(1) subplot(2,1,1) plot(w1,20*log10(abs(hh1))) grid xlabel(39。)。幅度 /db39。 subplot(2,1,2) plot(w1,angle(hh1)) grid xlabel(39。)。相位 /rad39。
高考資料相關(guān)推薦
鄂ICP備17016276號-1