利用MATLAB實(shí)現(xiàn)數(shù)字低通濾波器的設(shè)計.doc

西 南 石 油 大 學(xué) 實(shí) 驗(yàn) 報 告 課 程數(shù)字信號處理實(shí)驗(yàn)項(xiàng)目利用MATLAB實(shí)現(xiàn)數(shù)字低通濾波器的設(shè)計成 績專業(yè)年級通信工程09學(xué) 號0907050228指導(dǎo)教師姓 名王延禮同組人姓名實(shí)驗(yàn)日期一 實(shí)驗(yàn)?zāi)康模?學(xué)習(xí)用Matlab直接設(shè)計模擬濾波器和數(shù)字濾波器。2學(xué)習(xí)用沖激響應(yīng)不變法和雙線性變換法的Matlab的實(shí)現(xiàn)。 二 實(shí)驗(yàn)內(nèi)容：設(shè)計滿足下列指標(biāo)的數(shù)字低通濾波器： Wp=0.2*pi, Rp=1db Ws=0.5*pi Rs=20db Fs=1khz1. 利用B、C1型設(shè)計出模擬低通濾波器，采用沖激響應(yīng)不變法、雙線性發(fā)轉(zhuǎn)換成數(shù)字低通濾波器。 2. 直接設(shè)計出B、C1型數(shù)字低通濾波器。三 實(shí)驗(yàn)步驟：程序1Wp=2*pi*0.1*1000;Ws=2*pi*0.25*1000;Rp=1;Rs=20;N,Wn=buttord(Wp,Ws,Rp,Rs,s);z,p,k=buttap(N);B,A=butter(N,Wn,s);freq1=linspace(0,Wp,5);freq2=linspace(Wp,Ws,15);freq3=linspace(Ws,10*pi*2,25);h1=20*log10(abs(freqs(B,A,freq1);h2=20*log10(abs(freqs(B,A,freq2);h3=20*log10(abs(freqs(B,A,freq3);plot(freq1 freq2 freq3/(2*pi),h1,h2,h3);grid;Xlabel(Frequency in Hz);Ylabel(gain in DB);圖一程序2wp=0.2*pi;ws=0.5*pi;rp=1;rs=20;fs=1000;omegap=wp*fs;omegas=ws*fs;N,Wn=buttord(omegap,omegas,rp,rs,s);B A=butter(N,Wn,s);b,a=impinvar(B,A,fs);h,w=freqz(b,a,256);h=20*log10(abs(h);plot(w/pi,h);圖二程序3wp=0.2*pi;ws=0.5*pi;rp=1;rs=20;fs=1000;omegap=2*fs*tan(wp/2);omegas=2*fs*tan(ws/2);N,Wn=cheb1ord(omegap,omegas,rp,rs,s);B A=cheby1(N,rp,Wn,s);b,a=bilinear(B,A,fs);h,w=freqz(b,a,256);h=20*log10(abs(h);plot(w/pi,h);圖三程序4wp=0.2*pi;ws=0.5*pi;rp=1;rs=20;N,Wn=buttord(wp/pi,ws/pi,rp,rs);B A=butter(N,Wn);h,w=freqz(B,A,256);h=20*log10(abs(h);plot(w/pi,h);圖四程序5 Wp=0.2*pi;Ws=0.5*pi;Rp=1;Rs=20;T=0.001;Fs=1000;omegap=(2/T)*tan(Wp/2);omegas=(2/T)*tan(Ws/2);N,Wn=cheb1ord(omegap,omegas,Rp,Rs,s);B,A=cheby1(N,Rp,Wn,s);b,a=bilinear(B,A,Fs);h,w=freqz(b,a,256);h1=20*log10(abs(h);plot(w/pi,h1);grid;xlabel(Digital Frequency in pi units);ylabel(Gain in DB);axis(0 1 -50 10);圖五程序6 Wp=0.2;Ws=0.5;Rp=1;Rs=20;disp(I)N,Wn=cheb1ord(Wp,Ws,Rp,Rs)B,A=cheby1(N,Rp,Wn);disp();fprintf(%.4en,B);disp();fprintf(%.4en,A);w=linspace(0,0.8*pi,50);h1=20*log10(abs(freqz(B,A,w);plot(w/pi,h1);grid;xlabel(Normalized frequency);ylabel(Gain in DB );axis(0 0.8 -50 1);圖六 四、實(shí)驗(yàn)小