FIR滤波器的设计及性能分析.docxVIP

汕头大学工学院 三级项目报告 课程名称： 数字信号处理 课程设计题目：FIR数字滤波器均衡器设计实践及性能分析 指导教师： 姜永权 系另IJ: 电子工程系专业：通信工程 学号：姓名：黄兰凤 合作者 完成时间：2010年11月19日至11月30_H 项目意义与目标 F1R数字滤波器是一类广泛应用的数字滤波器。如自适应横向数字滤波器 就采用FIR数字滤波器，在GSM手机中就用到。在实际DSP系统中能够止确运 用F1R数字滤波器，理解线性相位的含义和应满足的条件。通过恰当举例，通 过亲自编程，给出设计结果，绘制频率响应曲线，掌握FIR数字滤波器设计算 法。FIR数字滤波器有具有精确的线性相位，FIR滤波器总是稳定的。设计方法 通常是线性，可以用硬件快速实现。 二,项目内容 根据频率选择性数字滤波器的性能指标，计算线性相位FIR数字滤波器的 理想单位冲激响应。选择Kaiser窗函数的阶数及形状控制参数，给出频率选择 性线性相位FIR数字滤波器的设计结果，画出FIR数字滤波器的频率响应曲 线。 根据任意频率响应(均衡器)性能要求，计算线性和位FIR数字均衡器的单 位冲激响应。选择Kaiser窗函数的阶数及形状控制参数，给出均衡器的设计结 果，画出均衡器的频率响应曲线。FIR数字滤波器、均衡器的性能分析。 三，项目报告正文 FIR数字滤波器 FIR数字滤波器的设计方法有：窗函数设计法，频率采样法，最优化设计 法。 这里介绍窗函数设计法。设计思想：用设计好的数字滤波器的频率响应 班严)2話7頃“)4珈去逼近理想数字滤波器的频率响应比(ejw) 设计步骤如下： 给定所要求的频率响应函数出(ejw)0 对该频率响应函数做傅里叶变换，求理想滤波器的单位冲激响应 hd (n)o 由过渡带宽和阻带最小衰减等指标要求，选取窗函数类型及窗的 长度。 求的所设计FIR滤波器的单位冲激响应h (n) = hd (n) ? w(n) o 计算频率响应函数II (e),并检验其是否满足性能要求。 设： 那么理想数字滤波器的单位冲激响应为: h, (n)二圭丄：论(*)严(“先为仏 =- 选择Kaiser窗对h (n进行加窗后: h (n) =h (n) 3 (n) n二0,1, 2,…,N 所以FIR数字滤波器的频率响应：班护比霊汀叫⑴片-宀 F1R数字滤波器的性能指标: 窗谱性能指标 加窗后滤波器性能指标 旁瓣峰值 (dB) 主瓣宽度 (X2 ji/N) 过渡带宽 (X2 jt/N) 阻带最小衰 减(dB) 矩形窗 -13 2 0.9 -21 巴特利特窗 -25 4 2. 1 -25 汉宁窗 -31 4 3. 1 -44 海明窗 -41 4 3. 3 -53 彳lj莱克曼窗 -57 6 5.5 -74 凯泽窗 -57 5 5 -80 FIR数字滤波器的理想单位冲激响应的程序如下: N=201; wc二pi/4; n=0:N-l; r=(N-l)/2; hdn=sin(wc*(n? r))/pi./( n?r); figure(l); stem(n,hdn); title(理想单位冲激响应); xlabel(rn); ylabelChdn); 图形： 用Kaiser窗来逼近理想窗。其函数形式如下: IQ IQ (fVl-ti-WCN-l)]2) 式中，I。Cx)为零阶贝塞尔函数；卩为一个可以自由选择的参数，它可以调 节主瓣与旁瓣宽度。当参数P越大，其频谱的旁瓣越小，主瓣宽度也随之增 加，因此0值可以考虑主瓣与旁瓣的影响进行选择。本实验会选择下表中的 A值来实现FIR数字滤波器。 0 过渡带宽As 阻带最小衰减(dB) 2.120 3.00 n/N -30 6.764 8A4VN -70 10.056 -100 程序如下： clc elf N=201; wc=pi/4; n=O:N-l; r=(N-l)/2; hdn=(wc/pi)*sinc(wc*(n-r)); %理想数字滤波器单位冲激响应hdn plot(n,hdn); titled理想单位冲激响应J; xlabel(n); ylabel(hdn); hold on; beta 1=2.120; beta2=6.764; beta3= 10.056; wn 1 =kaiser(length(n),beta 1); wn2=kaiser(length(n),beta2); wn3=kaiser(length(n),beta3); figure subplot(311 ),plot(abs(wn 1)); title(*betal=2.120 的 Kaiser Window Wave*); %Kaiser 窗波形图 xlabel(n); ylabelC/absCwnl)1); subplot