實(shí)驗(yàn)iir濾波器ppt課件-資料下載頁(yè)
2025-04-29 00:58本頁(yè)面
【正文】 0,2,2,2,2,0]。n=0:55。subplot(2,2,1)。 stem(n,x, 39。. 39。)。 %畫(huà)出時(shí)域序列axis([0 60 100 50])。hold on。 %保持原圖n=0:60。 %擴(kuò)大橫軸的范圍 0~60m=zeros(61)。 plot(n,m)。xlabel(39。n39。)。 ylabel(39。x(n)39。)。title(39。心電圖信號(hào)采樣序列 x(n)39。)。X=fft(x,256)。 %求 x的頻譜subplot(2,2,2)k=0:255。plot(2*k/256,abs(X))xlabel(39。k39。)。 ylabel(39。X(k)39。)。title(39。心電圖信號(hào)采樣序列的頻譜 39。)。%用雙線(xiàn)性變換法濾波器設(shè)計(jì) IIRDFT=1。 Fs=1/T。 %采樣頻率Wp=*pi。 Ws=*pi。Rp=1。 Rs=40。Omegap=2*tan(Wp/2)/T。 %模擬波器通帶截止頻率Omegas=2*tan(Ws/2)/T。 %模擬波器通帶截止頻率[N,Wc]=buttord(Omegap,Omegas,Rp,Rs,39。s39。)。 %計(jì)算巴特沃斯模擬濾波器階數(shù) N及 3dB截止頻率 Wc[Z,P,K]=buttap(N)。 %求 butterworth歸一化原型 H(p)，零極點(diǎn)增益形式[b,a]=zp2tf(Z,P,K)。 % 轉(zhuǎn)換成歸化系統(tǒng)函數(shù)形式， b,a是 H(p)的分子分母向量[bt,at]=lp2lp(b,a,Wc)。 %去歸一化轉(zhuǎn)換為實(shí)際低通 Ha(s)% [b1,a1]=butter(N,Wc,39。s39。)。 %可以代替以上三個(gè)步驟 [bb,ab]=bilinear(bt,at,Fs)。 %雙線(xiàn)性變換得 H(z)， Fs為采樣頻率%求數(shù)字濾波器的幅頻特性[H,w]=freqz(bb,ab,100)。 %對(duì)應(yīng)采樣間隔為 *pimag=abs(H)。db=20*log10((mag+eps)/max(mag))。subplot(2,2,3)。plot(w/pi,db)。axis([0,50,10])。grid ontitle(39。濾波器的幅頻響應(yīng)曲線(xiàn) 39。)。%用設(shè)計(jì)濾波器對(duì)輸入信號(hào)進(jìn)行濾波y=filter( )。 %完成對(duì) x(n)的濾波subplot(2,2,4)。n=0:55。stem(n,y,39。.39。)。axis([0 60 100 50])。hold on。n=0:60。m=zeros(61)。plot(n,m)。xlabel(39。n39。)。 ylabel(39。x(n)39。)。title(39。進(jìn)行濾波后的心電圖信號(hào) 39。)。簡(jiǎn)述實(shí)驗(yàn)?zāi)康募霸?。給出上述四個(gè)濾波器的頻響特性，并對(duì)結(jié)果加以分析。根據(jù)對(duì)心電信號(hào)進(jìn)行處理的程序，寫(xiě)出用雙線(xiàn)性變換法設(shè)計(jì)的巴特沃斯低通 IIR數(shù)字濾波器的系統(tǒng)函數(shù)。五、實(shí)驗(yàn)報(bào)告及要求
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1