数字信号低通滤波器课程设计.doc

数字信号处理课程设计 题目:有限冲激响应数字滤波器设计 班级: 学号: 姓名: 一、实验目的:加深对数字滤波器的常用指标和设计过程的理解。 二、实验原理:低通滤波器的常用指 标: P P P for H Ω≤Ω+≤Ω≤-,1(1δδ π δ≤Ω≤Ω≤ΩS S for H ,( 数字滤波器技术指标: 通带边缘频率P ω,阻带边缘频率S ω ,通带最大衰减1δ,最小阻带衰减2δ 通带峰值起伏1(log 201101αδ--=, 最小阻带衰减(log 20210 2αδ-=。 设计步骤: ? 给定理想的频率响应函数及技术指标 ? 根据阻带衰减选择窗函数 ? 根据过渡带宽度确定N 值 ? 求所设计的FIR 滤波器的单位抽样响应 ? 作出图形 三、设计内容: 设计FIR 低通滤波器,通带边界频率πω 2.0=p ,阻带边界频率 πω3.0=s ,阻带衰减2δ不小于50dB 。选择一个合适的窗函数,确定脉冲响应,并给出所设计的滤波器的频率响应图。(要求设计完后,再选择其他的窗函数,并比较结果,分析原因 一.低通滤波器 1.求N ,(d h n S Ω-P Ω-P S ΩPassband Stopband Transition band Fig 1 Typical magnitude specification for a digital LPF 根据公式推导计算:. 0.30.20.16.66.666 s p s p w w w w N N w w πππππ?=-=-=?== =- 12 N τ-= 1 (sin[(],( d s h n w n n n ττ πτ=-≠- 海明窗 2.设计程序如下 N=66 wn=hamming(66; nn=[0:1:65]; alfa=(66-1/2; hd=sin(0.25*pi*(nn-alfa+eps./(pi*(nn-alfa+eps; h=hd.*wn; [h1,w1]=freqz(h,1; plot(w1/pi,20*log10(abs(h1; axis([0,1,-100,10]; grid; xlabel(归一化频率/pi; ylabel(幅度/dB; 3. 图形如下 0.1 0.2 0.3 0.40.50.60.7 0.8 0.9 1 -100 -90-80-70-60-50-40-30 -20-100 10归一化频率/pi 幅度/d B 布拉克曼窗 2 5.50.1w N ππ??= =,解得110N =,所以154.5 2 N τ -==。 布拉克曼窗 2424([0.420.5cos(0.08cos( ]([0.420.5cos( 0.08cos( ]( 1 1 109 109 N N n n n n w n R n R n N N ππππ=-+=-+-- 所以: sin[0.25(54.5] 24(((([0.420.5cos( 0.08cos( ]( 54.5 109 109 d N n n n h n h n w n w n R n n ππ-=?= ?=-+- 程序: wn=blackman(110; nn=[0:1:109]; alfa=(110-1/2; hd=sin(0.25*pi*(nn-alfa+eps./(pi*(nn-alfa+eps; h=hd.*wn; [h1,w1]=freqz(h,1; plot(w1/pi,20*log10(abs(h1; axis([0,1,-100,10]; grid; xlabel(归一化频率/pi; ylabel(幅度/dB; 0.1 0.2 0.3 0.40.50.60.7 0.8 0.9 1 -100 -90-80-70-60-50-40-30-20-100 10归一化频率/pi 幅度/d B 二.高通滤波器 海明窗: 22([0.540.46cos(]([0.540.46cos( ](1 65 N N n n w n R n R n N ππ=-=-- 其脉冲响应为: 11sin[(]sin[(]sin[(32.5]sin[0.25(32.5] 2 2 (132.5 ( 2c d N N n w n n n h n N n n πππ--- -----= = --- 因为对于高通滤波器中N 为奇数,即N 取67, 所以: sin[(32.5]sin[0.25(32.5] 2((([0.540.46cos( ]( 32.5 65 d N n n n h n h n w n R n n ππ---=?= ?-- 程序： wn=hamming(67; nn=[0:1:66]; alfa=(67-1/2; hd=(sin(pi*(nn-alfa+eps-sin(0.25*pi*(nn-alfa+eps./(