经典滤波器的设计和使用.doc

1、*实践教学*兰州理工大学计算机与通信学院2015 年春季学期信号处理课程设计报告题 目： 经典滤波器的设计和使用 班 级： 通信工程 12 级（ 1 ）班 姓名： 设计质量（30 分）： 学号： 122501xx 说明书质量（10 分）： 同组成员： 指导教师： 兰 州 理 工 大 学 课 程 设 计 报 告目录摘要 .1一、基本原理 .21.1 IIR 数字滤波器设计原理 .21.2 FIR 数字滤波器设计原理 .3二、滤波器详细设计 .52.1 信号源 .52.2 低通滤波器 .62.3 带通滤波器 .92.3 高通滤波器 .11总结 .13参考文献 .14附录 .15兰 州 理 工 大

2、学 课 程 设 计 报 告1摘要滤波器广泛应用于电子电路中，具有平滑信号、选频、去噪和使交流信号变为直流信号的功能。可分为低通滤波器、高通滤波器、带通滤波器和带阻滤波器，也可分为模拟滤波器和数字滤波器，数字滤波器分为IIR数字滤波器和FIR数字滤波器。设计滤波器之前先根据实际需要确定滤波器类型和指标，通过MATLAB分别设计低通、带通和高通三种滤波器，首先通过三种滤波器对时域信号滤波，然后对信号频域变换后滤波，最后对信号进行时域搬移后滤波。通过MATLAB得出的频谱图对信号的三种形式通过三种滤波器后的结果进行分析，通过分析比较各滤波器的优点和缺点。关键词：MATLAB IIR滤波器 FIR滤波

3、器 频域变换 时域搬移兰 州 理 工 大 学 课 程 设 计 报 告2一、基本原理数字滤波器是一个离散系统，其系统函数一般可表示为 的有理多项式形式，即1zNaiiMjjzbzH1)(0当a i；i=1，2，N都为0时，上式描述的系统称为有限脉冲响应数字滤波器，简称FIR(Finite-Impulse Response)数字滤波器。当系数a i；i=1，2，N中至少有一个是非0时，上式描述的系统称为无限脉冲响应数字滤波器，简称IIR(Infinite-Impulse Response)数字滤波器。对于IIR数字滤波器，一般满足 MN ，这时将系统称为N阶IIR数字滤波器。对于FIR数字滤波器，

4、系统函数中 的有理多项式的最高次幂M就是其阶数。1z在设计数字滤波器时，首先根据具体的滤波成分，确定待设计数字滤波器的技术指标，再求解出满足设计指标的离散系统的系统函数 H(z)。数字滤波器的技术指标通常用通带截频 ，p通带最大衰减A p，阻带截频 ，阻带最小衰减A S表示。s1.1 IIR 数字滤波器设计原理IIR滤波器的设计就是根据给定的数字滤波器技术指标，确定上式中滤波器的阶数N和系数a i，b j 。在满足技术指标的条件下，滤波器的阶数应尽可能低。应为滤波器阶数越低，实现滤波器的成本就越低。由于模拟滤波器的设计技术已经非常成熟，且可得闭合形式的解，因此在设计IIR滤波器时，一般是通过模

5、拟滤波器来设计滤波器。设计方法是先将数字滤波器技术指标转换为对应的模拟滤波器技术指标，然后设计满足技术指标的模拟滤波器H(s)，再将模拟滤波器H(s)转换为对应的数字滤波器H(z)。因此，在IIR数字滤波器设计中，模拟滤波器设计是基础，模拟滤波器到数字滤波器的转换是核心。而模拟滤波器的设计都是通过设计模拟低通滤波器来实现，比较常见的模拟低通滤波器有巴特沃思和切比雪夫等。将模拟滤波器变换为数字滤波器的主要方法有脉冲响应不变法和双线性变换法。兰 州 理 工 大 学 课 程 设 计 报 告3图1-1 IIR数字滤波器的设计过程1.2 FIR 数字滤波器设计原理表1-1 四种线性相位FIR滤波器的特性

6、类型 I II III IV阶数M 偶数 奇数 偶数 奇数Hk的对称性 偶对称 偶对称奇对称 奇对称0.5M 0.5M 0.5M 0.5M0 0 0.50.5A( )关于 =0的对称性偶对称 偶对称奇对称 奇对称A( )关于 = 的对称性偶对称 奇对称奇对称 偶对称A( )的周期 24 24A(0) 任意 任意 0 0A( ) 任意 0 0 任意可适用的滤波器类型LP，HP，BP，BS等LP，BP微分器，Hilbert(希尔伯特)变换器微分器，Hilbert(希尔伯特)变换器，HP窗函数法设计FIR滤波器的基本思想是在时域逼近理想滤波器的单位脉冲相应。首先根据待逼近的理想滤波器的频率响应 ，由

7、IDTFT求出理想滤波器的单位脉冲响应 ，再将无限)(jdeH khd长的 加窗截断得到有限长序列 。为了获得线性相位FIR滤波器，在窗函数法设计FIRkhd khd滤波器的过程中，需要将线性相位因子 加入理想滤波器的频率响应 ，具体步)5.0(Mj )(jdeH骤如下：（1） 根据所需设计的滤波器，确定线性相位滤波器的类型（I型，II型，III型，IV型） 。（2） 确定理想滤波器的幅度函数 。)(dA数字滤波器 H(z)模拟滤波器 H(s)待设计数字滤波器指标模拟滤波器指标频率转换设计模拟滤波器S 到 z 域转换兰 州 理 工 大 学 课 程 设 计 报 告4（3） 确定理想滤波器的相位

8、。)(dM5.0对I型和II型线性相位FIR滤波器 ，对III型和IV型线性相位FIR滤波器 。0 2/（4） 计算 。khd deAkhjkjddd)(21（5） 截断 得khd khkd其中wk是长度N=M+1的矩形窗。兰 州 理 工 大 学 课 程 设 计 报 告5二、滤波器详细设计2.1 信号源产生一个连续信号，包含低频、中频、高频分量，对其进行采样，进行频谱分析分别设计三种高通、低通、带通滤波器对信号进行滤波处理，观察滤波后信号的频谱。合成信号频谱程序如下:t=(1:100)/2000;x=sin(2*pi*t*20)+sin(2*pi*t*200)+sin(2*pi*t*500);

9、plot(t,x);n=1:100;t=n/2000;X=fft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*20)+sin(2*pi*t*200)+sin(2*pi*t*500);plot(w,abs(X(1:256);xlabel(Hz);ylabel(频率响应幅度);title(合成信号频谱图);图 2-1 合成信号波形图兰 州 理 工 大 学 课 程 设 计 报 告6图 2-2 合成信号频谱图2.2 低通滤波器根据数字滤波器的设计原理，首先将数字域的指标转化为模拟域的指标设计模拟低通滤波器，然后应用双线性变换法将模拟滤波器转化为数字滤波器将指标转换成归

10、一化模拟低通滤波器的指标，通过归一化的模拟低通滤波器阶数 N 和 3dB 截止频率的计算，将模拟域频率变换成模拟低滤波器 H(s)，并用双线性变换法将 H(s)转换成数字低通滤波器 H(z)，由此得到低通数字滤波器。IIR 低通滤波器的设计程序为:Ft=8000;Fp=1000;Fs=1200;As=100;Ap=1;wp=2*pi* Fp/Ft;ws=2*pi*Fs/Ft;兰 州 理 工 大 学 课 程 设 计 报 告7fp=2*Fp*tan(wp/2);fs=2*Fs*tan(ws/2);n11,wn11=buttord(wp,ws,1,50,s);b11,a11=butter(n11,wn11,s);num11,den11=bilinear(b11,a11,0.5);h,w=freqz(num11,den11);plot(w*8000*0.5/pi,abs(h);title(低通滤波器);图 2-3 IIR 低通滤波器读入信号对信号傅里叶变换进行频谱分析，利用 filter 函数滤波，对得到的信号傅里叶变换进行频谱分析，与滤波前进行比较。图 2-4 滤波前信号波形兰 州 理 工 大 学 课 程 设 计 报 告8图 2-5 滤波后信号波形图 2-6 滤波后信号的频谱图 2-7 滤波前信号的频谱