濾波器實(shí)驗(yàn)5基于matlab的數(shù)字濾波器設(shè)計(jì)-文庫(kù)吧
2024-08-10 03:49 本頁(yè)面
【正文】 解 在MATLAB中可以用remezord 和remez兩個(gè)函數(shù)設(shè)計(jì)，其結(jié)果如圖2，MATLAB程序如下： fedge=[800 1000]。 mval=[1 0]。 dev=[ ]。 fs=4000。 [N,fpts,mag,wt]=remezord(fedge,mval,dev,fs)。 b=remez(N,fpts,mag,wt)。 [h,w]=freqz(b,1,256)。 plot(w*2000/pi,20*log10(abs(h)))。 grid。 xlabel(39。頻率/Hz39。) 。 ylabel(39。幅度/dB39。)。 一、實(shí)驗(yàn)內(nèi)容： 利用MATLAB編程設(shè)計(jì)一個(gè)數(shù)字帶通濾波器，指標(biāo)要求如下：通帶邊緣頻率：，通帶峰值起伏：。阻帶邊緣頻率：，最小阻帶衰減： 。 分別用窗函數(shù)法和等波紋濾波器法設(shè)計(jì)兩種FIR數(shù)字濾波器。實(shí)驗(yàn)要求：給出FIR數(shù)字濾波器的沖激響應(yīng)，繪出它們的幅度和相位頻響曲線(xiàn)，討論它們各自的實(shí)現(xiàn)形式和特點(diǎn)。11）用窗函數(shù)法實(shí)現(xiàn)：調(diào)用函數(shù)[n,wn,bta,ftype]=kaiserord（f，a，dev，fs）參數(shù)：f=[ ]為對(duì)應(yīng)數(shù)字頻率，， a=[0 1 0]為由f指定的各個(gè)頻帶上的幅值向量，一般只有0和1表示；和f長(zhǎng)度關(guān)系為（2*a的長(zhǎng)度）—2=（f的長(zhǎng)度）devs=[ ]用于指定各個(gè)頻帶輸出濾波器的頻率響應(yīng)與其期望幅值之間的最大輸出誤差或偏差，長(zhǎng)度與a相等，計(jì)算公式：阻帶衰減誤差=，通帶波動(dòng)衰減誤差=fs缺省值為2HZ [n,wn,bta,ftype]=kaiserord([ ],[0 1 0],[ ])。%用kaiserord函數(shù)估計(jì)出濾波器階數(shù)n和beta參數(shù) h1=fir1(n,wn,ftype,kaiser(n+1,bta),39。noscale39。)。 [hh1,w1]=freqz(h1,1,256)。 figure(1) subplot(2,1,1) plot(w1/pi,20*log10(abs(hh1))) grid xlabel(39。歸一化頻率w39。)。ylabel(39。幅度/db39。)。 subplot(2,1,2) plot(w1/pi,angle(hh1)) grid xlabel(39。歸一化頻率w39。)。ylabel(39。相位/rad39。)。h1 = Columns 1 through 8
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1