數(shù)字信號處理實驗 濾波器在語音信號中的應(yīng)用 (2)

1、數(shù)字信號處理實驗姓名：彭金金學(xué)號：071041219同組成員：王春香實驗三IIR數(shù)字濾波器的設(shè)計及濾波實現(xiàn)一、程序如下：wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;ripple=10(-Rp/20);Attn=10(-As/20);T=0.02*pi;Fs=1/T;Omgp=(2/T)*tan(wp/2);Omgs=(2/T)*tan(ws/2);n,Omgc=buttord(Omgp,Omgs,Rp,As,s);bal,aal=butter(n,Omgc,s);bd,ad=bilinear(bal,aal,Fs);sos,g=tf2sos(bd,ad);H,w=freqz(

2、bd,ad);dbH=20*log10(abs(H)+eps)/max(abs(H) );plot(w/pi,dbH);title(濾波器的幅頻響應(yīng));ylabel(dB);xlabel(頻率(pi);axis(0,0.5,-40,5);set(gca,XTickMode,manual,XTick,0,0.2,0.3,0.5);set(gca,YTickMode,manual,YTick,-50,-15,-1,0);gridn=0:55;x=-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

3、-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;figureplot(n,x,k);xlabel(n);ylabel(x);title(心電圖信號線圖);gridfigure stem(n,x,m)xlabel(n);ylabel(x);title(心電圖離散序列)gridw1=linspace(-2*pi,2*pi,500);N=length(x);n=0:N-1;X=x*exp(-j*n*w1);figureplot(w1,abs(X),r)xlabel(w1),ylabel(X)t

4、itle(心電圖幅頻特性)gridh=ifft(H);n=0:length(h)-1;figurestem(n,h,g)axis(0,50,-0.3,0.3)xlabel(n),ylabel(h)title(數(shù)字濾波器時域響應(yīng))gridhx=conv(x,h);n2=0:length(hx)-1;figurestem(n2,hx,g)xlabel(n),ylabel(hx)title(濾波后的心電圖信號)axis(0,60,-50,10)gridn1=0:length(hx)-1;HX=hx*exp(-j*n1*w1);figureplot(w1,abs(HX)xlabel(w1)ylabel

5、(HX)title(濾波后心電圖信號的幅頻響應(yīng))grid二、各信號的波形如下：三、思考題：用雙線性變換發(fā)設(shè)計數(shù)字濾波器過程中，變換公式中T是我取值對設(shè)計結(jié)果有無影響？為什么？解：用雙線性變換發(fā)設(shè)計數(shù)字濾波器過程中，變換公式中T是我取值對設(shè)計結(jié)果沒有影響，因為T只是一個參數(shù)，在設(shè)計濾波器的過程中其值可以被約去，對濾波器的特性沒有任何影響。四、實驗結(jié)論： 實驗中對心電圖信號進(jìn)行濾波有兩種方法，一種是頻域相乘，一種是時域卷積，由濾波器的時域特性可以看出它是一個N值很小的離散信號，采用反傅里葉變換將其轉(zhuǎn)換為時域與心電圖信號卷積更為方便。 對一次濾波后的信號進(jìn)行2次3次濾波得到的信號波形如下：hhx=conv(h,hx);n3=0:length(hhx)-1;figurestem(n3,hhx)axis(0,70,-40,10)xlabel(n3),ylabel(hhx),title(二次濾波后的心電圖信號),gridhhhx=conv(h,hhx);n4=0: