实验报告基于matlab的数字滤波器设计.pdfVIP

实验 7\8 基于 MATLAB 的数字滤波器设计 实验目的：加深对数字滤波器的常用指标和设计过程的理解。 实验原理：低通滤波器的常用指标： j  1 P  G(ej  )  1 P , f or   P 1  G(e ) P 1  P G(ej  )  S , f orS      s   通带边缘频率：P ，阻带边缘频率：S  S  P P S Passband Stopband  Transition 通 带 起 伏 ： P ，通 带 峰 值 起 伏 ： band Fig 7.1 Typical magnitude p  20 log10 (1  p )[dB ] ，阻带起伏： s ，最 specification for a digital LPF   20 log ( )[dB ] 小阻带衰减： S 10 s 。 数字滤波器有 IIR 和 FIR 两种类型，它们的特点和设计方法不同。 在 MATLAB 中，可以用[b a]=butter （N,Wn）等函数辅助设计 IIR 数字滤波器, 也可以 用 b=fir1(N,Wn,’ftype’) 等函数辅助设计 FIR 数字滤波器。 实验内容： 利用 MATLAB 编程设计一个数字带通滤波器，指标要求如下： 通带边缘频率：P 1  0.45 P 2  0.65 通带峰值起伏：p  1[dB ] 。 阻带边缘频率：S 1  0.3 S 2  0.75 最小阻带衰减： S  40[dB ] 。 分别用 IIR 和 FIR 两种数字滤波器类型进行设计。 实验要求：给出 IIR 数字滤波器参数和 FIR 数字滤波器的冲激响应，绘出它们的幅度和相位 频响曲线，讨论它们各自的实现形式和特点。 实验内容： IRR 代码： wp=[0.45*pi,0.65*pi]; ws=[0.3*pi,0.75*pi]; Ap=1; As=40; [N,Wc]=buttord(wp/pi,ws/pi,Ap,As) ; [b,a]=butter(N,Wc)%[b,a] = butter(n,Wn,ftype) [H,w]=freqz(b,a); subplot(211); mag=abs(H); plot(w/pi,mag); xlabel(角频率 (w)); ylabel(幅度 |Ha(j\w)|); title(IIR|Ha(j\w)|); grid on; subplot(212); phase=angle(H);