实验7\8基于MATLAB勺数字滤波器设计实验目的:加深对数字滤波器的常用指标和设计过程的理解。
实验原理:低通滤波器的常用指标:1 一6P 兰G(e^) ≤ 1 + 6P , for 国≤ ωPG(J") ≤ 6s, for 国s ≤ ⑷≤ ∏通带边缘频率:'P ,阻带边缘频率:'s,通带起伏:J P,通带峰值起伏:C(P= —20 IOg io (^-OP )【d B 】阻带起伏.冠SPaSSband StOPbandTran Siti onband Fig 7.1 TyPiCaI magn itude SPeCifiCati On for a digital LPF:S = -20 log ιo(r)[dB 】O数字滤波器有IIR和FlR两种类型,它们的特点和设计方法不同。
在MATLAB^,可以用[b , a]=butter ( N,Wr)等函数辅助设计IIR数字滤波器,也可以用b=fir1(N,Wn, 'type ')等函数辅助设计FIR数字滤波器。
实验内容:利用MATLAB编程设计一个数字带通滤波器,指标要求如下:通带边缘频率:∙∙P1=0.45^,∙∙ P2=0∙65 二,通带峰值起伏:[dB】O阻带边缘频率:'s10.3…,'s2 0.75…,最小阻带衰减:-S 4°[dB] O分别用IIR和FlR两种数字滤波器类型进行设计。
实验要求:给出IIR数字滤波器参数和FIR数字滤波器的冲激响应,绘出它们的幅度和相位频响曲线,讨论它们各自的实现形式和特点。
实验内容:IRR代码:wp=[0.45*pi,0.65*pi];ws=[0.3*pi,0.75*pi];Ap=1;A S=40;[N,Wc]=buttord(wp∕pi,ws∕pi,Ap,As);[b,a]=butter(N,Wc)%[b,a] = butter( n, Wn,'ftype')最小阻带衰减:I -W llrreqz(b-a=SUbP-Of(211=FnagHabS(H)-P-Of(WHLmag)-x ωb e 一(->5(W))y ωb e一(-≡M 一 HamW)口≡e (--R 一 HamW)一) grid On-SUbPOf(212=PhaSeHang-e(H=P -OfWPLPhaSe=x ωb e 一(->5(W))y ωb e -(- W⅛) ≡e (-一R盘a ≡m sB )≡-x ⅛x -幅度IHaOMI≡相拉SS 5 CUS S S 0⅛口 g >酉tt 2 =R ⅛⅛⅛J ≡B ⅛O S F NJ £Q 4S S U7 CJCD S >⅛⅛≡F-RV VWPl HO∙4*prWP2H0∙6*prV V WSl H 0∙3*prws2 H0∙7*prVV =Γl w i d f h Hmin((wprwsu(ws2lwp2))>> tr_width =0.3142>> M = ceil(6.2*pi/tr_width) + 1>> M = 63>> n=[0:1:M-1];>> wc1 = (ws1+wp1)/2; wc2 = (wp2+ws2)/2; >> wc=[wc1/pi,wc2/pi];>> window= hanning(M);>> [h1,w]=freqz(window,1);>> figure(1);>> subplot(2,1,1)>> stem(window);>> axis([0 60 0 1.2]);>> grid;>> xlabel('n');>> title('Hanning 窗函数');>> subplot(2,1,2)>> plot(w/pi,20*log(abs(h1)/abs(h1(1)))); >> axis([0 1 -350 0]);>> grid;>> xlabel('w/\pi');>> ylabel('幅度(dB)');>> title('Hanning 窗函数的频谱');>> hn = fir1(M-1,wc, hanning (M));>> [h2,w]=freqz(hn,1,512);>> figure(2);>> subplot(2,1,1)>> stem(n,hn);>> axis([0 60 -0.25 0.25]);>> grid;>> xlabel('n');>> ylabel('h(n)');>> title('Hanning 窗函数的单位脉冲响应'); >> subplot(2,1,2)>> plot(w/pi,20*log(abs(h2)/abs(h2(1)))); >> grid;>> xlabel('w/\pi');>> ylabel(' 幅度(dB)');>> figure(3);>> phase=angle(h1);>> plot(phase);>> axis([1 pi -1 0]);>> xlabel('w/\pi');>> ylabel(' 线性相位');>> title('Hanning 窗函数相位特性曲线')回FiIe Edit VieW InSert TOOlS DeSktOP WindOW HeIP◎ A 聾紳⑥毘謠▼层□ Ξ∣ ■ 0M FigUre 1 Hanning 窗函数OC S I I * iHanning®函数的频谱-100-200 -300 00.1 0.2 0.3 0.40.5 0.6 0.7 08 0.9 1w/x(8P)置Q FigUre 2 口 回 耳iHanni叩窗函数相位特性曲线O I I I I I I I。