DSP第二次作业.DOCVIP

DSP第二次作业 （题52）考虑习题51中的低通巴特沃兹数字滤波器设计，各截止频率和衰减要求为： Wp=0.2*pi,Rp=1dB;Ws=0.4*pi,Rs=15dB,T=1s. a. 利用本章讨论的双线性变化技术原理以及bilinear函数，画出对数幅度响应比较模拟原型和数字滤波器的脉冲响应。 b. 利用butter函数进行设计，并与a中的设计进行比较。 a. Wp=0.2*pi;Rp=1;Ws=0.4*pi;Rs=15;T=1; Fs=1/T; omegap=(2/T)*tan(Wp/2);%通带截止频率预畸 omegas=(2/T)*tan(Ws/2);%阻带截止频率预畸 [N,Wn]=buttord(omegap,omegas,Rp,Rs,s); [B,A]=butter(N,Wn,s); [b,a]=bilinear(B,A,Fs); [h,w]=freqz(b,a,256); h1=20*log10(abs(h));%增益(dB) plot(w/pi,h1); grid xlabel(Digital Frequency in pi units) ylabel(Gain in dB) axis([0 1 -80 10]) freqz(b,a,256,Fs) [ha,x,t]=impulse(b,a); plot(t,ha);title(单位脉冲响应ha(t)); b. Wp=0.2*pi;Rp=1;Ws=0.4*pi;Rs=15; %设计滤波器 [N,Wn]=buttord(Wp,Ws,Rp,Rs,s) [z,p,k]=buttap(N); [B,A]=butter(N,Wn,s) f1=linspace(0,Wp,5); f2=linspace(Wp,Ws,15); f3=linspace(Ws,2*pi*10,30); h1=20*log10(abs(freqs(B,A,f1))); h2=20*log10(abs(freqs(B,A,f2))); h3=20*log10(abs(freqs(B,A,f3))); plot([f1 f2 f3]/(2*pi),[h1,h2,h3]); grid xlabel(Frequency in Hz) ylabel(Gain in dB) 其中 B = 0 0 0 0 0.4506 A = 1.0000 2.1410 2.2919 1.4372 0.4506 因此要满足所要求的指标需要4 阶巴特沃斯滤波器，其系统函数为： H(s)=0.4506/(s4+2.1410 s3+2.2919s2+1.4372s+0.4506) （题57）用汉宁窗设计技术设计一个带阻滤波器，技术指标为： 低阻带边缘：0.4pi 高阻带边缘：0.6pi As=40dB 低通带边缘：0.2pi 高通带边缘：0.8pi Rp=1dB 画出设计的滤波器的脉冲响应和幅度响应（dB值）。 解：根据窗函数最小阻带衰减特性，汉宁窗可达到44dB的最小阻带衰减，它提供的过渡带宽为6.2π/N。 wp1=0.2*pi;wr1=0.4*pi;wp2=0.8*pi;wr2=0.6*pi; tr_width=min((wr1-wp1),(wp2-wr2)); N=ceil(6.2*pi/tr_width); n=0:1:N-1 wc1=(wr1+wp1)/2;wc2=(wr2+wp2)/2; hd=ideal_bs(wc1,wc2,N); w_ham=(hanning(N)); h=hd.*w_ham; [db,mag,pha,w]=freqz_m(h,[1]); delta_w=2*pi/1000; Ap=-(min(db(1:1:wp1/delta_w+1))); Ar=-round(max(db(wr1/delta_w+1:1:wr2/delta_w+1))); stem(n,hd);title(理想单位脉冲响应hd(n)); stem(n,w_ham);title(汉宁窗w(n)); stem(n,h);title(实际单位脉冲响应h(n)); plot(w/pi,db);title(幅度响应(db)); 其中freqz_m(b,a)函数如下 function[db,mag,pha,w]=freqz_m(b,a);%计算绝对和相对的幅度响应 [H,w]=freqz(b,a,1000,whole); H=(H(1:1:501)); w=(w(1:1:501)); mag=abs(H);