數(shù)字信號處理上機指導(dǎo)

1、數(shù)字信號處理上機實驗指導(dǎo)書光電信息與通信工程學(xué)院電了信息工程專業(yè)北京信息科技大學(xué)2011年11月實驗一離散傅里葉變換的性質(zhì)及應(yīng)用3實驗二iir數(shù)字濾波器設(shè)計5實驗一 離散傅里葉變換的性質(zhì)及應(yīng)用一、實驗?zāi)康?. 了解dft的性質(zhì)及其應(yīng)用2. 熟悉matlab編程特點二、實驗儀器及材料計算機，matlab軟件三、實驗內(nèi)容及要求1驗證dft性質(zhì)。(1) 已知序列兀(力)=45() + 33(n 一 1) + 23(n 一 2) + /( 一 3)，計算序列的6點 dft的幅度、相位、實部及虛部,并用圖形表示相應(yīng)的x(n),|x(k),angle x(k), rex(q)lmx(£)。(2)

2、 編程求解xep (/?)和兀少)，驗證dftlms)卜 rex伙)dftxop(n)= jimx伙)* (3)對序列x(n)=46(ii)4-3(n-1) + 2s(n-2) +-3),自行設(shè)計程序驗證dft時移及頻移性質(zhì)。2. 利用dft實現(xiàn)兩序列的卷積運算，并研究dft點數(shù)與混疊的關(guān)系。(1) 已知兩序列x(n)=/i(n) = 2n = 0,l,2,3(2)用克接法(即用線性卷積的立義計算,結(jié)果，并以圖形方式農(nóng)示結(jié)果;n-1 y(n) = v x(m) x h(n _ m),見下式)計算線性卷積y(n)=x(n)*h(n)的0<n< n-¥ n2-2m=0(3)

3、利用dft計算圓周卷積；分別令圓周卷積的點數(shù)為l=8, 10, 12,以圖形方式 表示結(jié)果。(4) 對比直接法和岡周卷積法所得的結(jié)果。3. 利用dft進行譜分析。(1)對模擬信號進行譜分析x(r) = cos(8 加)+ 2 cos(l 6加)+ 0.5 cos(20 加)請問信號x(r)的最高頻率是多少?(2) 假定抽樣頻率fs = 64hz ,編程實現(xiàn)記錄數(shù)據(jù)點數(shù)為n=16,32,64三種情況下的譜分 析。分別應(yīng)出其郵度譜，討論三種情況幅度譜的區(qū)別？幅度譜與頻率分辨率的對應(yīng)關(guān)系。四、思考題1. 在實驗內(nèi)容1中都驗證了什么dft性質(zhì)？2. 說明為什么在實驗內(nèi)容2中l(wèi)=8, 10, 12時采用

4、fft法對兩序列計算圓周卷積的 結(jié)果不同？與線性卷積y(n)=x(n)*h(n)的結(jié)果關(guān)系如何？3. 說明為什么在實驗內(nèi)容2中為n =16,32,64三種情況下的譜分析結(jié)果不同？周期信 號的譜分析的點數(shù)應(yīng)該如何選??？五、實驗報告要求1. 列出本實驗編寫的所有文件及各項實驗結(jié)果曲線，加注必要的說明；2. 対實驗結(jié)果作理論計算，解釋實驗結(jié)果；3. 冋答思考題；4. 總結(jié)實驗體會及實驗中存在的問題。六、本實驗所用的部分matlab函數(shù)abs(x):対復(fù)數(shù)取模(對實數(shù)収絕對值)。angle(x):対復(fù)數(shù)取相角。real(x):対復(fù)數(shù)取實部。imag(x):對復(fù)數(shù)取虛部。y=conv(xl,x2):將序

5、列xl(n)與序列x2(n)做線性卷積，結(jié)果賦給序列y(n)。y=fft(x,l):將序列門)做1點快速傅立葉變換，結(jié)果賦給序列y(n)oy=iffi(x,l):將序列x(n)做l點快速傅立葉反變換，結(jié)果賦給序列y(n)。stem(x,y):以x為橫坐標(biāo)，y為縱坐標(biāo)畫離散曲線。實驗二iir數(shù)字濾波器設(shè)計一、實驗?zāi)康?、熟悉iir數(shù)字濾波器的設(shè)計原理與方法。2、掌握數(shù)字濾波器的計算機軟件實現(xiàn)方法。3、通過觀察對實際心電圖信號的濾波作用，學(xué)習(xí)數(shù)字濾波器在實際中的應(yīng)用。二、實驗儀器及材料計算機，matlab軟件三、實驗內(nèi)容及要求1. 設(shè)計巴特沃斯低通數(shù)字濾波器對人體心電信號進行濾波(1) 人體心電圖

6、信號在測量過程中會受到工業(yè)高頻干擾，所以必須經(jīng)過低通濾波處理, 才能作為判斷心臟功能的有用信息。以下為一個實際心電圖信號采樣序列x(n),其中存在 高頻干擾，采樣周期t=1秒。在實驗中，以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對序列x(n)用fft做頻譜分析，生成x(n)的頻譜

7、圖。(2) 用沖激響應(yīng)不變法設(shè)計一個巴特沃斯低通iir數(shù)字濾波器h(z)o設(shè)計指標(biāo)參數(shù)為：在通帶內(nèi)頻率低于0.2兀時，最大衰減小于2db；在肌帶內(nèi)0.371,7i頻率區(qū)間上，最小衰減大于25db。寫出數(shù)字濾波器h的表達式，畫岀濾波器的幅頻響應(yīng)曲線|h(ejo) |。(3) 用所設(shè)計的濾波器對實際心電圖信號采樣序列x(n)進行濾波處理，編寫程序，求 濾波厲的序列y(n),并分別畫出濾波前后的心電圖信號波形圖和頻譜圖。2. 用help查看內(nèi)部函數(shù)cheblord.m及chebyl.m, 了解調(diào)用格式。編程設(shè)計教材習(xí)題62,求模擬濾波器ha(s)的衣達式，畫出該模擬濾波器頻響。3模擬濾波器的數(shù)字化用

8、內(nèi)部函數(shù)impinvar及bilinear實現(xiàn)教材習(xí)題6-5,求數(shù)字濾波器h的衣達式。四、思考題1. 對比實驗內(nèi)容1屮濾波前后心電圖信號x(n), y(n),說明濾波器的作用。2. 雙線性變換與沖激不變法比較有何優(yōu)缺點？五、實驗報告要求1. 列岀本實驗編寫的所有文件及各項實驗結(jié)果曲線，加注必耍的說明；2. 回答思考題；3. 總結(jié)實驗體會及實驗屮存在的問題。六、本實驗所用的部分matlab函數(shù) l=length(x):求序列x長度。 y=fft(x,l):將序列x(n)做l點快速傅立葉變換，結(jié)果賦給序列y(n)。 n,wn = buttord(wp,ws,rprs,0)：計算模擬 butterw

9、orth 濾波器的最小階次 n 和 截止頻率為wn。其中wp,ws為模擬角頻率，單位為rad/s, rp,rs的單位為db。 b,a = butter(n,wn,9j：設(shè)計模擬截止頻率為 wn (rad/s)的口 階 butterworth 低通 濾波器，返冋值為模擬濾波器的系數(shù)。 n,wn = buttord(wp,ws,rp,rs):計算數(shù)字butterworth濾波器的最小階次n和截 止頻率為wn。其中wp,ws為數(shù)字角頻率/pi,例如數(shù)字角頻率為0.2pi,則賦值 時為0.2。rp,rs的單-位為db。 b,a = butter(n,wn)：設(shè)計數(shù)字截止頻率為wn的n階butterwo

10、rth低通濾波器， 返冋值為模擬濾波器的系數(shù)。 y=filter(b,a,x):將序列x(n)通過濾波器濾波后生成序列y(n),濾波器的分母多項 式系數(shù)構(gòu)成a向量，分子多項式系數(shù)構(gòu)成b向量。 bz,az = impinvar(b,a, fs)：沖激響應(yīng)不變法，返回值為數(shù)字濾波器的系數(shù)bz、az。 其中fs為抽樣頻率。 bz,az = bilinear(b,a,fs):雙線性變換，返冋值為數(shù)字濾波器的系數(shù)bz、az。其中 fs為抽樣頻率。 h w=freqz(b,a)：由數(shù)字濾波器濾波器分母多項式系數(shù)構(gòu)成的a向量和分子多項 式系數(shù)構(gòu)成的b向量求系統(tǒng)頻響。 h w=freqs(b,a)：由模擬濾波

11、器濾波器分母多項式系數(shù)構(gòu)成的a向量和分子多項 式系數(shù)構(gòu)成的b向量求系統(tǒng)頻響。實驗三fir數(shù)字濾波器設(shè)計一、實驗?zāi)康? 熟悉線性相位fir數(shù)字濾波器特性。2. 掌握用窗函數(shù)法設(shè)計fir數(shù)字濾波器的原理和方法，理解各種窗函數(shù)對濾波性能 的影響。3. 掌握用頻率抽樣法設(shè)計fir數(shù)字濾波器的原理和方法，理解過渡帶上抽樣點數(shù)對 濾波性能的影響。二、實驗儀器及材料計算機，matlab軟件三、實驗內(nèi)容及方法1. 學(xué)習(xí)firl例程，用窗函數(shù)法設(shè)計fir數(shù)字濾波器1) 設(shè)計線性相位數(shù)字低通濾波器，要求通帶截止頻率fp=.5kh乙，阻帶截止頻率 fst = 2.5khz，通帶最大衰減ap = db ,阻帶最小衰減

12、as =40db 0設(shè)采樣頻率為 fs=okhzf要求濾波器階數(shù)盡量低。選取合適的窗函數(shù)編程設(shè)計，求濾波器的單位脈沖響應(yīng)h(n),并畫出幅頻特性和和頻 特性曲線。2) 設(shè)計線性相位數(shù)字帶阻濾波器，要求阻帶截止頻率costh = 0.65 , ®=0.35/r, 通帶截止頻率co 山=0.83龍，co pi = 0.17兀通帶最大衰減勺=db ,阻帶最小衰減 as =50db o選取合適的窗函數(shù)編程設(shè)計，求濾波器的單位脈沖響應(yīng)h(n),并畫出幅頻特性和相頻 特性曲線。2. 學(xué)習(xí)fk2例程，用頻率抽樣法設(shè)計fir數(shù)字濾波器，指標(biāo)同實驗內(nèi)容1。四、思考題1 .總結(jié)窗函數(shù)法屮窗兩數(shù)主耍參數(shù)對

13、濾波性能的影響。2.頻率抽樣法設(shè)計fir濾波器的優(yōu)缺點是什么？五、實驗報告要求1. 列岀本實驗編寫的所有文件及各項實驗結(jié)果曲線，加注必耍的說明；2. 對實驗結(jié)果作理論計算，解釋實驗結(jié)果；3. 冋答思考題；4. 總結(jié)實驗體會及實驗中存在的問題。六、參考firl implements the classical method of windowed linear-phase fir digital filter design.b = firl(n,wn)b = firl(n,wn；ftype,)b = fir 1 (n,wn, window)b = firl (nvn/ftype*,window)firl例程% design a 48th-order fir bandpass filter with passband 0.35 w 3 w 0.65:b = firl(48,0.35 0.65);freqz(b,l,512)fir2例程% design a 30th-order lowpass filter and overplot the desired f requency response w