數(shù)字信號(hào)處理實(shí)驗(yàn)第二次實(shí)驗(yàn)matlab系統(tǒng)的頻率響應(yīng)(編輯修改稿)
2024-11-13 18:01 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 t(2,1,2)。stem(k,abs(xkk),39。g39。)。grid on。xlabel(39。k39。)。ylabel(39。2*NFFT39。)。title(39。Matlab2N點(diǎn)FFT39。)。七、問題思考，結(jié)果一致嗎？說明原因。、運(yùn)算量上的比較。八、心得體會(huì)第四篇：數(shù)字信號(hào)處理實(shí)驗(yàn)二北京信息科技大學(xué)數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告題 目：數(shù)字信號(hào)處理課程設(shè)計(jì)實(shí)驗(yàn)學(xué) 院： 信息與通信工程學(xué)院 專 業(yè)： 通信工程專業(yè) 姓 名： 班 級(jí)：學(xué) 號(hào)： 指導(dǎo)老師：實(shí)驗(yàn)?zāi)康氖煜IR數(shù)字濾波器的設(shè)計(jì)原理與方法。掌握數(shù)字濾波器的計(jì)算機(jī)軟件實(shí)現(xiàn)方法。通過觀察對(duì)實(shí)際心電圖信號(hào)的濾波作用，學(xué)習(xí)數(shù)字濾波器在實(shí)際中的應(yīng)用。實(shí)驗(yàn)儀器及材料計(jì)算機(jī)，MATLAB軟件實(shí)驗(yàn)內(nèi)容及要求1．設(shè)計(jì)巴特沃斯低通數(shù)字濾波器對(duì)人體心電信號(hào)進(jìn)行濾波(1)人體心電圖信號(hào)在測(cè)量過程中會(huì)受到工業(yè)高頻干擾，所以必須經(jīng)過低通濾波處理，才能作為判斷心臟功能的有用信息。以下為一個(gè)實(shí)際心電圖信號(hào)采樣序列x(n)，其中存在高頻干擾，抽樣周期Ts=1秒。在實(shí)驗(yàn)中，以x(n)作為輸入序列，濾除其中干擾成分。x(n)=[4,2,0,4,6,4,2,4,6,6,4,4,6,6,2,6,12,8,0,16,38,60,84,90,66,32,4,2,4,8,12,12,10,6,6,6,4,0,0,0,0,0,2,4,0,0,0,2,2,0,0,2,2,2,2,0] 對(duì)序列x(n)用FFT做頻譜分析，生成x(n)的頻譜圖。(2)設(shè)計(jì)一個(gè)巴特沃斯低通IIR數(shù)字濾波器H(z)。設(shè)計(jì)指標(biāo)參數(shù)為：，最大衰減小于1dB； 在阻帶內(nèi) [, π]頻率區(qū)間上，最小衰減大于15dB。jw|H(e)|。寫出數(shù)字濾波器H(z)的表達(dá)式，畫出濾波器的幅頻響應(yīng)曲線(3)用所設(shè)計(jì)的濾波器對(duì)實(shí)際心電圖信號(hào)采樣序列x(n)進(jìn)行濾波處理，編寫程序，求濾波后的序列y(n)，并分別畫出濾波前后的心電圖信號(hào)波形圖和頻譜圖。y(n)= [0,0,0,0, 0,0,0,0, ,, , ,, ,,,,,, ,, ,]源程序clear all，clcx=[4,2,0,4,6,4,2,4,6,6,4,4,6,6,2,6,12,8,0,16,38,60,84,90,66,32,4,2,4,8,12,12,10,6,6,6,4,0,0,0,0,0,2,4,0,0,0,2,2,0,0,2,2,2,2,0]。%未經(jīng)濾波的心電圖信號(hào) L=length(x)。l=0:L1。y=fft(x,L)。Wp=*pi。Ws=*pi。Rp=1。Rs=15。[N,Wn] = buttord(Wp,Ws,Rp,Rs,39。s39。)。[b,a] = butter(N,Wn,39。s39。)。[numa,dena]=impinvar(b,a,1)。w=linspace(0,pi,1024)。h=freqz(numa,dena,w)。norm=max(abs(h))。numa=numa/norm。[z,p]=tf2zp(b,a)。figure(1)plot(w,20*log10(abs(h)/norm))。grid。xlabel(39。數(shù)字頻率39。)。ylabel(39。幅度響應(yīng)dB39。)。figure(2)plot(w,abs(h))。grid。xlabel(39。數(shù)字頻率39。)。ylabel(39。幅度響應(yīng)|H(e^(jw))|39。)。figure(3)zplane(z,p)。xx=filter(b,a,x)。yy=fft(xx,L)。figure(4)subplot(2,1,1)stem(l,x)。title(39。未經(jīng)濾波的心電圖信號(hào)39。)。xlabel(39。n39。)。subplot(2,1,2)stem(l,xx)。title(39。經(jīng)濾波之后的心電圖信號(hào)39。)。xlabel(39。n39。)。figure(5)subplot(2,1,1)plot(l,abs(y))。title(39。未經(jīng)濾波的心電圖信號(hào)的頻譜39。)。subplot(2,1,2)plot(l,abs(yy))。title(39。經(jīng)濾波處理的心電圖信號(hào)的頻譜39。)。2．，了解調(diào)用格式。編程設(shè)計(jì)教材習(xí)題62，求模擬濾波器Ha(s)的表達(dá)式。源程序close all clear all clcWp=2*pi*3000。Rp=2。Ws=2*pi*12000。Rs=50。Fs=24000。w=linspace(0,pi,1024)。[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,39。s39。)。e=sqrt(10^(Rp/10)1)。[b,a]=cheby1(N,e,Wn,39。s39。)[numa,dena]=impinvar(b,a,Fs)。h=freqz(numa,dena,w)。norm=max(abs(h))。plot(w*Fs/pi,20*log10(abs(h)/norm))title(39。幅度響應(yīng)39。)xlabel(39。頻率(Hz)39。)ylabel(39。幅度(dB)39。)grid3．模擬濾波器的數(shù)字化用內(nèi)部函數(shù)impinvar及bilinear實(shí)現(xiàn)教材習(xí)題65，求數(shù)字濾波器H(z)的表達(dá)式。源程序close all clear all clc b1=[0 0 1]。a1=[1 1 1]。b2=[0 0 1]。a2=[2 3 1]。[numa1,dena1]=impinvar(b1,a1,)[numa2,dena2]=bilinear(b1,a1,)[numa3,dena3]=impinvar(b2,a2,)[numa4,dena4]=bilinear(b2,a2,)本實(shí)驗(yàn)所用的部分MATLAB函數(shù)178。 L=length(x)：求序列x長(zhǎng)度。178。 y=fft(x,L)：將序列x(n)做L點(diǎn)
數(shù)學(xué)相關(guān)推薦
鄂ICP備17016276號(hào)-1