1、數字信號處理作業與上機實驗（第六章） 班 級： 學 號： 姓 名： 任課老師： 完成時間： 2014.11.18 信息與通信工程學院 20132014學年第 2 學期第6章 無限脈沖響應數字濾波器設計1、教材p195：13,14,15,16,17,18,19（1） 代碼如下：clear;close allt=1;fs=4000;wpi=0.45*pi,wpl=0.25*pi;wsu=0.55*pi,wsl=0.15*pi;wpz=0.25,0.45;wsz=0.15,0.55;wp=2/t*tan(wpz/2);ws=2/t*tan(wsz/2);rp=3;as=40;n,wc=buttord

2、(wp,ws,rp,as,s);b,a=butter(n,wc,s);bz,az=bilinear(b,a,fs);nd,wdc=buttord(wpz,wsz,rp,as);bd,adz=butter(nd,wdc);hk=freqz(bd,adz)figure(1)plot(angle(hk)title(相頻特性曲線)hp,w=freqz(bd,adz,4000);figure(2)plot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.1 圖1.1（2） 代碼如下：fs=80000;T=1/fs;rp=0.5;rs=45;wp=4000*2*pi*T;ws=20

3、000*2*pi*T;N,wc=buttord(wp,ws,rp,rs,s);B,A=butter(N,wc,s);hk=freqz(B,A);figure(1)plot(angle(hk)title(相頻特性曲線)hp,w=freqz(B,A);figure(2)plot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.2 圖1.2（3） 代碼如下：fs=80000;T=1/fs;rp=0.5;rs=45;wp=4000*2*pi*T;ws=20000*2*pi*T;N,wc=cheb1ord(wp,ws,rp,rs,s);B,A=cheby1(N,rp,wc,s);

4、fk=0:120000:120000;wk=2*pi*fkhk=freqs(B,A,wk);figure(1)plot(angle(hk)title(相頻特性曲線)hp,w=freqz(B,A,4000);Bz,Az=impinvar(B,A);figure(2)freqz(Bz,Az);圖示如1.3 圖1.3（4） 代碼如下：fs=2500000;T=1/fs;rp=1;rs=40;wp=325000*2*pi*T;ws=225000*2*pi*T;rp=1;rs=40;N,wpo=ellipord(wp,ws,rp,rs);B,A=ellip(N,rp,rs,wpo,s)hk=freqz(

5、B,A)figure(1)plot(angle(hk)title(相頻特性曲線)hp,w=freqz(B,A);figure(2)plot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.4 圖1.4（5） 代碼如下：fs=5000000;T=1/fs;rp=1;rs=50;wp=560000*2*pi*T,780000*2*pi*T;ws=375000*2*pi*T,1000000*2*pi*T;N,wpo=ellipord(wp,ws,rp,rs,s);B,A=ellip(N,rp,rs,wpo,s)hk=freqz(B,A)figure(1)plot(angle(

6、hk)title(相頻特性曲線)hp,w=freqz(B,A);figure(2)plot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.5 圖1.5（6） 代碼如下：fs=5000;T=1/fs;rp=1;rs=40;wp=500000*2*pi*T,2125000*2*pi*T;ws=1050000*2*pi*T,1400000*2*pi*T;N,wpo=ellipord(wp,ws,rp,rs,s);B,A=ellip(N,rp,rs,wpo,s);hk=freqz(B,A)figure(1)plot(angle(hk)title(相頻特性曲線)hp,w=fre

7、qz(B,A);figure(2)plot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.6 圖1.6 （7） 代碼如下：fs=80000;T=1/fs;rp=0.5;rs=45;wp=4000*2*pi*T;ws=20000*2*pi*T;N,wc=buttord(wp,ws,rp,rs,s);B,A=butter(N,wc,s);hk=freqz(B,A)figure(1)plot(angle(hk)title(相頻特性曲線)Bz,Az=impinvar(B,A);figure(2)freqz(Bz,Az);hp,w=freqz(Bz,Az);figure(3)p

8、lot(w/pi,20*log10(hp);title(損耗函數曲線)圖示如1.7 圖1.72、某信號為：，其中設計最低階數IIR數字濾波器，按下圖所示對進行數字濾波處理，實現：1）IIR數字濾波器代碼如下：fs=3800;T=1/fsfp=130;fs=600;rs=50;wp=2*pi*fp*T;ws=2*pi*fs*T;Bt=ws-wp; alph=0.5842*(rs-21)0.4+0.07886*(rs-21); N=ceil(rs-8)/2.285/Bt); wc=(wp+ws)/2/pi; hn=fir1(N,wc,kaiser(N+1,alph); M=1024;Hk=fft(hn,M);k=0:M/2-1;wk=(2*pi/M)*k;fig