课 程 设 计 报 告
课程名称: 专业综合课程设计 学生姓名: 段博文 学 号: 201116020227 专业班级: 通信工程11102班 指导教师: 朱明旱 完成时间: 2014年6月15日 报告成绩: 评阅意见: 评阅教师 日期
1
线性相位FIR带通滤波器的设计
1.课程设计目的
1.能够运用本课程中学到的知识,设计基于窗口函数法的FIR数字滤波器。要求掌握数字信号处理的基本方法; 2. FIR滤波器的设计步骤和方法;
3.能够熟练MATLAB进行计算机辅助设计和方针验证设计内容的合理性。
2.课题要求
1.了解数字滤波器的性能要求 2.了解窗函数设计滤波器的基本理论
3.设计原理
3.1 FIR滤波器
数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理,以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展,数字滤波器已可用计算机软件实现,也可用大规模集成数字硬件实时实现。数字滤波器根据其单位冲激响应函数的时域特性分为两种:无限长冲激响应(IIR) 滤波器和有限长冲激响应(FIR)滤波器。
FIR数字滤波器又称有限长单位冲激响应滤波器,它的优点是可以做成具有严格的线性相位,同时又可以具有任意的幅度特性。此外,FIR滤波器的单位抽样响应是有限长的,因而滤波器一定是稳定的。再有,FIR滤波器由于单位冲激响应是有限长的,所以可以用快速傅里叶变换(FFT)算法来实现过滤信号,从而可大大提高运算效率。在滤波器设计中要对理想滤波器抽样响应进行截断. 截
2
断后不可避免的产生了频谱泄漏, 为了尽量减小频谱泄漏, 在设计滤波器时要采用不同的窗函数来满足不同用途的要求.各种窗函数的幅频响应都存在明显的主瓣和旁瓣. 主瓣宽度和旁瓣的幅值衰减特性决定了窗函数的应用. 用于滤波器的窗函数,一般要求窗函数主瓣宽度窄,以获得较好过渡带:旁瓣相对值尽可能小, 以增加通带段的平稳度和增大阻带的衰减.窗函数应满足在0 由线性系统理论可知,在某种适度条件下,输入到线性系统的一个冲击完全可以表征系统。当我们处理有限的离散数据时,线形系统的响应(包括对冲击的响应)也是有限的。若线性系统仅是一个空间滤波器,则通过简单地观察它对冲击的响应,我们就可以完全确定该滤波器。通过这种方式确定的滤波器称为有限冲击响应(FIR)滤波器。FIR滤波器是在数字信号处理(DSP)中经常使用的两种基本的滤波器之一,另一个为IIR滤波器。IIR滤波器是无限冲激响应滤波器。 3.1.1 FIR数字滤波器的线性相位条件 设滤波器单位脉冲相应的长度为N,系统函数为 H(z)??h(n)z?n n?0N?1由此式可见,H(z)是Z?1的(N-1)个零点,原点z=0是(N-1)阶重极点,位于r=1的单位圆内,系统永远稳定。稳定性和线性相位特性是FIR滤波器的突出优点。 FIR滤波器的设计任务是选择有限长度的h(n),使传输函数H(ej?)满足要求。 线性相位条件: 对于长度为N的h(n),传输函数为 H(e)??h(n)e?j?m j?n?0N?1 H(ej?m)?Hg(?)e?j?(?) 3 式中,Hg(?)称为幅度特性,?(?)称为相位特性。线性相位是指相位函数 ?(?)满足如下特性:?(?)???w或?(?)??0??w,?0是起始相位,?为常数,一般称第一种情况为第一类线性相位,称第二种情况为第二类线性相位。 满足第一类线性相位的充要条件是:h(n)为实序列,并且对(N-1)/2偶对称,即h(n)=h(N-n-1) 3.1.2 FIR数字滤波器的基本结构 FIR滤波器的基本结构有一下几种:直接型、级联型、线性相位型、频率采样型。 1.直接型 设FIR滤波器的单位冲击相应h(n)为一个长度为N的序列,则滤波器系统函数为: H(z)??h(n)z?n n?0N?1表示已这一系统输入输出关系的差分方程为 y(n)??h(m)x(n?m) m?0N?1直接由查分方程可得出对应的网络结构如图所示 图3-1 FIR滤波器的直接型结构 直接型结构的优点:简单直观,惩罚运算量较少。 缺点:调整零点较难 2.级联型 4 当需要控制滤波器的传输零点时,可讲H(z)分解为实系数二阶因子的乘积形式: H(z)??(?ok??1kz?1??2kz?2) k?1N/2式中,H(z)为h(n)的z变换,?ok,?1k,?2k为实数。级联型结构如图所示 图3-2 FIR滤波器的级联型结 该结构的优点:调整零点比直接型方便。 缺点:H(z)中的系数比直接型多,因而需要的乘法器多。当H(z)的阶次高时,也不易分解。 3.2窗函数设计法 窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。通常在设计滤波器之前,应该先根据具体的工程应用确定滤波器的技术指标。在大多数实际应用中,数字滤波器常常被用来实现选频操作,所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。 用窗函数法设计FIR滤波器的步骤如下: 1.根据过渡带宽及阻带衰减要求,选择窗函数的类型并估计窗口长度N(或阶数M=N-1)。窗函数类型可根据最小阻带衰减AS独立选择,因为窗口长度N对 5