20xx用雙線性變換法設(shè)計原型低通為巴特沃茲型的數(shù)字iir高通濾波器(參考版)
2025-01-17 01:18本頁面
【正文】 資料共分享，我們負(fù)責(zé)傳遞知識。此次課程設(shè)計對自己在許多方面都有了很大的幫助。同時，用MATLAB計算有關(guān)數(shù)字濾波器的設(shè)計參數(shù)，如H(z)、h(n)等，對于數(shù)字濾波器的硬件實現(xiàn)也提供了一條簡單而準(zhǔn)確的途徑和依據(jù)。與其他高級語言的程序設(shè)計相比，MATLAB環(huán)境下可以更方便、快捷地設(shè)計出具有嚴(yán)格線性相位的FIR濾波器，節(jié)省大量的編程時間，提高編程效率，且參數(shù)的修改也十分方便，還可以進(jìn)一步進(jìn)行優(yōu)化設(shè)計。因此在設(shè)計過程中，我是看了遍數(shù)字信號處理的教科書并且特別注重看了關(guān)于數(shù)字濾波器的設(shè)計的內(nèi)容，再通過利用參考文獻(xiàn)與網(wǎng)絡(luò)，完成了用Matlab進(jìn)行數(shù)字信號處理課程設(shè)計。 title(＇濾波器輸出y的波形＇) 6 調(diào)試分析 z = + + + + + + + + p = + + + + + + + + k = 分子系數(shù)是 Columns 1 through 6 Columns 7 through 12 Columns 13 through 18 分母系數(shù)是 Columns 1 through 6 Columns 7 through 12 Columns 13 through 18 f1 = 300 7 測試結(jié)果 8 實驗總結(jié) 通過此次課程設(shè)計我對數(shù)字信號處理有了進(jìn)一步的認(rèn)識，但是由于自己數(shù)字信號處理的理論課程沒有達(dá)到很深的了解，對數(shù)字濾波器的設(shè)計缺少認(rèn)識，因此在此次課程設(shè)計中我還是遇到了很多問題。 xlabel(＇t＇)。 grid on。%數(shù)字濾波器輸出 subplot(2,2,4)。ylabel(＇x(t)＇)。 axis([0,50*pi,3,3])。 subplot(2,2,3)%繪制輸入x的波形 plot(x)。ylabel(＇x2(t)＇)。 axis([0,50*pi,3,3])。 subplot(2,2,2)%繪制x1的波形 plot(x2)。ylabel(＇x1(t)＇)。 axis([0,50*pi,3,3])。 subplot(2,2,1)%繪制x1的波形 plot(x1)。 x=x1+x2。 x1=sin(2*pi*f1*t)。 n=0:600。%繪制傳輸函數(shù)零極點 title(＇傳輸函數(shù)的零極點＇)。 figure。 ylabel(＇增益/ dB＇)。grid %繪制巴特沃茲高通濾波器幅頻特性 axis([0 1 60 5])。 % 繪制增益響應(yīng) figure。 h = freqz(num,den,w)。disp(den)。disp(num)。 [num,den]=bilinear(num2,den2,)。 %設(shè)計濾波器 [num1,den1] = butter(N,Wn,＇s＇)。 Lemta_S_EquivalentLowPass=1/Eta_S。 %歸一化 Eta_S=OmegaS1/OmegaP1。 % 通帶頻率非線性變換 OmegaS1=tan(Ws1/2)。 % 通帶頻率模數(shù)轉(zhuǎn)換 Ws1=2*pi*OmegaS1_1/Fp。 % 阻帶邊界頻率 Fp=2000。 %阻帶衰減/dB OmegaP1_1=500。clear all Rp = 1。
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號-1