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;

3、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('相頻特性曲線')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

4、=4000*2*pi*T;ws=20000*2*pi*T;N,wc=cheb1ord(wp,ws,rp,rs,'s');B,A=cheby1(N,rp,wc,'s');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

5、=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(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*

6、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(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=

7、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=freqz(B,A);figure(2)plot(w/pi,20*log10(hp);title('損耗函數曲線')圖示如1.6 圖1.6 （7） 代碼如下：fs=800

8、00;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)plot(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;figure(1)