1、第6章无限脉冲响应数字滤波器设计说明数字信号处理作业与上机实验(第六章) 班 级: 学 号: 姓 名: 任课老师: 完成时间: 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;a
2、s=40;n,wc=buttord(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=4
3、000*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=4000*2*pi*T;ws=20000*2*pi*T;N,wc=cheb1ord(wp,ws,rp,rs,s);B,A=che
4、by1(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=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,w
5、po,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*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(
6、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=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(相频特
7、性曲线)hp,w=freqz(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)
8、;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);plot(wk/pi,20*log10(abs(Hk(k+1);title(损耗函数);grid onfigure(2)plot(wk/pi,angle(Hk(k+1)/pi);grid ontitle(相频特性曲线);图示如2.1 图2.1
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2