脉冲响应不变法设计巴特沃斯低通滤波器.docx

1、脉冲响应不变法设计巴特沃斯低通滤波器燕山大学课程设计说明书题目：脉冲响应不变法设计数字低通滤波 器学院（系）：电气工程学院年级专业：09级检测1班学 号：0901*学生姓名： 郝瑾指导教师： 王娜教师职称： 讲师电气工程学院课程设计任务书课程名称： 数字信号处理课程设讣 基层教学单位：仪器科学与工程系 指导教师：王娜学号0901*学生姓名郝瑾（专业）班级09检测1班设计题目7、脉冲响应不变法设计数字低通滤波器设 计 技 术参 数给定技术指标为：几=100处，人=300股，勺=3B, a,=20dB, 采样频率Fx =1000 Hz o设 计 要 求设计Butterworth低通滤波器，用脉冲响

2、应不变法转换成数字滤波器。参 考 资 料数字信号处理方面资料MATLAB方而资料周次前半周后半周应 完 成 内 容收集消化资料、学习MATLAB软件， 进行相关参数计算编写仿真程序、调试指导教 师签字王娜 嚟囂说w： i.此表一式四份.系、抬导教师、学生各一份.报送院教务科一份。2、学生那份任务书要求装订到课程设计报告前面。电气工程学院教务科第一章摘要滤波器是自动控制、信号处理和通信领域的重要组成部分，广泛地应用于各种系统 中。MATLAB语言是一种简单、髙效的髙级语言，是一种内容丰富、功能强大的分析工具, 其应用范围几乎覆盖了所有的科学和工程计算领域。MATLAB中提供了丰富的用于模拟滤 波

3、器设计的函数，通过编程可以很容易实现低通、髙通、带通、带阻滤波器，并能画岀滤 波器的幅频特性曲线，大大简化了模拟滤波器设计。本文将通过利用MATLAB滤波器设计 函数直接实现巴特沃斯低通滤波器的设汁，找到应用MATLAB来设汁低通滤波器的方法。 介绍了用脉冲响应不变法设计巴特沃斯低通滤波器的基本理论和设计思想，给岀了基于 MATLAB设计巴特沃斯低通滤波器的具体步骤，文中还给出了利用MATLAB产生的滤波器 的幅频及相频特性曲线。详细介绍了在基于MATLAB设计低通滤波器过程中常用到的工具 和命令。关键字：巴特沃斯低通滤波器脉冲响应不变法第二章引言滤波器设讣在电子工程、应用数学和汁算机科学领域

4、都是非常重要的内，。而低通滤 波器在信号处理中的作用等同于其它领域如金融领域中移动平均数（moving average） 所起的作用，低通滤波器有很多种，苴中，最通用的就是巴特沃斯滤波器和切比雪夫滤 波器。巴特沃斯滤波器的特点是内的曲线最大限度平坦，没有起伏，而在则逐渐下降为零。 在振幅的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而 逐步减少，趋向负无穷大。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的 无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在 阻频带振幅衰减速度越快。苴他滤波器髙阶的振幅对角频率图和地结束的振幅对角频率 有不同

5、的形状。第三章基本原理3.1滤波器的基本理论3.1.1滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。经典滤波器是假泄输入信号班“）中的 有用成分和希望取出的成分各自占有不同的频带。这样，当M）通过一个线性系统（即 滤波器）后可讲欲去除的成分有效的去除。现代滤波器理论研究的主要内容是从含有噪 声的数据记录（又称时间序列）中估计岀信号的某些特征或信号本身。经典滤波器分为低通.髙通.带通、带阻滤波器。每一种又有模拟滤波器(AF)和数字滤波器(DF)o对数字滤波器，又有IIR滤波器和FIR滤波器匚IIR DF的转移函数是:.wr-()FIR DF的转移函数是：FIR滤波器可以对给泄的频率特性直接

6、进行设汁，而IIR滤波器目前最通用的方法 是利用已经很成熟的模拟滤波器的设计方法进行设讣。3.1.2 滤波器的技术要求低通滤波器：卩：通带截止频率 阻带下限截止频率勺：通带允许的最大衰减 久：阻带允许的最小衰减(勺，久的单位 dB)p：通带上限角频率阻带下限角频率3.2模拟滤波器的设汁3.2. 1巴特沃思(Butterworth)滤波器的介绍：因为我们设计的滤波器的冲击响应一般都为实数，所以有G(s)Gs) = G(QG(5)|.、.应=|G(jQJ|2这样，如果我们能由Q, 卩，求出Ig(7O)I ,那么就容易得到所需要 的G(s)。不同类型的lG(7Q)l的表达式，代表了几种不同类型的滤波

7、器。巴特沃思 (Butterworth)滤波器为：lGn)|2=TTFWC为待定常数，N为待泄的滤波器阶次。3. 2. 2巴特沃思模拟低通滤波器的设计由于每一个滤波器的频率范用将直接取决于设讣者的应用目的，因此必然是千差万 别。为了使设计规范化，我们需要将滤波器的频率参数作归一化处理。设所给的实际频 率为G (或/)，归一化后的频率为兄，对低通模拟滤波器，令显然，Ap=人又令归一化复数变量为P , P = j入,显然p =以=冋 =slGp所以巴特沃思模拟低通滤波器的设讣可按以下三个步骤来进行。(1)将实际频率归一化求C和NC,N可按下式求解：若令勺=3dB,则C = L这样巴特沃思滤波器的设

8、计就只剩一个参数N,这时（3）确定G）因为p = j入,根据上而公式有由1 + （-1严严=0解得.2k + N J2Nlh =exp(丿一2龙)这样可得G*(p)= 5一 八)5一 PzQ求得G（“）后，用代替变量厂 即得实际需要得GG）。3. 3数字滤波器的设计按一左规则将给出的数字滤波器的技术指标转换魏模拟低通滤波器的技术指标。根据转换后的技术指标设计模拟低通滤波器G（$）：再按一定的规则将G）转换成H。若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的 技术指标，然后按上述步骤2）设计出低通GG）,再将G（转换为所需的H。第四章设计过程4.1设计低通数字滤波器4.

9、1.1数字低通滤波藩设计结构图用脉冲响应不变法转换原理如下图:Ha (s) 拉氏逆变换 ha (t) 等间隔采样 ha(nT)二h(n) Z变换 H (z) 图4. 1结构框图4.1.2 设计步骤设计要求：脉冲响应不变法设计数字低通滤波器，给泄技术指标为：几00., A =300 Hz 勺=3dB av = 20dB 采样频率 F$=1000Hz首先应该得到角频率然后再按如下步骤求解：(1)将数字滤波器的技术要求转换为模拟滤波器的技术要求。(2)设计低通滤波器GG)由几= /p依次求出心，人，再求岀N,可得G(P)G(s) = G(p)亠然后由 F转换成GC)由G(s)求H4.2滤波器的MAT

10、LAB实现4.1.1在MATLAB下，设计巴特沃斯IIR滤波器可使用butter函数。Butter函数可设计低通、高通、带通和带阻的数字和模拟IIR滤波器，其特性为 使通带内的幅度响应最大限度地平坦，但同时损失截止频率处的下降斜度。在期望通带 平滑的情况下，可使用butter函数0 but ter函数的用法为：b, a：二butter (n, Wn, /ftype/) 其中n代表滤波器阶数，Wn代表滤波器的截止频率，这两个参数可使用buttord函数来 确定。buttord函数可在给左滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n, 同时给出对应的截止频率Wn。butt ord函数的用法为

11、：n, Wn二buttord (Wp, Ws, Rp, Rs) 其中Wp和Ws分别是通带和阻带的拐角频率(截I匕频率)，英取值范用为0至1之间。当 其值为1时代表采样频率的一半。Rp和Rs分别是通带和阻带区的波纹系数。不同类型(髙通、低通、带通和带阻)滤波器对应的Wp和恥值遵循以下规则：a.高通滤波器：WpWs为一元矢量且WpWs；b.低通滤波器：Wp和Ws为一元矢量且WpWs：c.带通滤波器:Wp 和 Ws 为二元矢量且 WpWs,如 Wp=0. 1,0. 8,Ws=0. 2, 0. 74.1.2 Mat lab的信号处理工具箱提供了频谱分析函数：freqz函数。freqz函数功能：离散时间

12、系统的频率响应。格式：h, w=freqz (b, a, n)h, f二freqz (b, a, n, Fs)h二freqz (b, a, w)h二freqz (b, a, f, Fs)freqz (b, a, n)说明：freqz用于计算数字滤波器H(Z)的频率响应函数H(ejs)。h,弱二freqz (b, a, n)可得到数字滤波器的n点复频响应值，这n个点均匀地分布在 0, “上,并将这n个频点的频率记录在w中，相应的频响值记录在h中。要求n为大于 零的整数,最好为2的整数次幕，以便采用FFT汁算，提髙速度。缺省时n二512。h, f =freqz(b, a, n, Fs)用于对H(e

13、je)在0, Fs/2上等间隔采样n点，采样 点频率及相应频响值分别记录在f和h中。由用户指左FS (以HZ为单位)值h二freqz (b, a, w)用于对H(ej)在0, 2 n 上进行采样，采样频率点由矢量w指泄。 h二freqz (b, a, f, Fs)用于对H(ejw)在0,FS上采样，釆样频率点由矢量f指定。 freqz (b, a, n)用于在当前图形窗口中绘制幅频和相频特性曲线。4.1. 3采用脉冲响应不变法将模拟滤波器转换为数字滤波器可使用impinvar函数。BZ,AZ二IMPINVAR(B,扎Fs)生成一个数字滤波器。BZ和AZ分别为脉冲值,结果与模拟滤波器的值相当。B

14、和A为采样区间.采样周期为1/Fs。当Fs缺省时，默认采样频率为1H氛第五章程序和仿真图5.1设计程序fp=100;fs=300;Rp=3;Rs=20;Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;Nn=256;n=(0:100-l) ;%采样点数Nz Wn=buttord(Wpf Wsz Rpz Rs, s1);b,a=butter(NWn,1s1);w=linspace(1,400,100) *2*pi;H=freqs (b, a, w) ; %在0, 2 n 上进行采样，采样频率点由矢量w指定 figure (1);plot(w/(2*pi),20*logl0(abs(H)

15、;title (巴特沃斯模拟滤波器幅频特性);xlabel (1 频率/Hz 1);ylabel (1 幅度/db);bzz az = imp invar (b, a, Fs);figure (2);zplane(bzz az);title (脉冲响应不变法零极点分布图1);figure (3);Subplot(2,lz1);f reqz (bz, az,Nn, Fs);对H(ej )在0,Fs/2上等间隔采样Nn点 grid on;title (巴特沃斯数字滤波器幅频特性);xlabel幅度/Hz ) ; ylabel (1 频率/db 1);Subplot(2,lz 2);title (巴

16、特沃斯数字滤波器相频特性);xlabel (1 幅度/Hz 1) ; ylabel (1 频率/db); grid on;figure (4);x=cos (2*pi*50*n/Fs)+cos(2*pi*500*n/Fs); plot(n,x);title (滤波前输入信号图像);xlabel (1 时间(s) 1);ylabel (1 幅值grid on;figure (5);y=filter (bzz az, x) ; $对原信号进彳亍滤波plot(n/fs, y);title (滤波前输入信号图像，)；xlabel时间(s) 1);ylabel (幅值);5.2 特性曲线仿真图图5. 1

17、巴特沃斯模拟滤波器幅频特性曲线Q Figure 2 o gFile Edit View Insert Tools Desktop Window Help凰 Q脉冲响应不变;去的專扳点團-1 05 0 0.5 1Real Pari图5. 2脉冲响应不变法的零极点图图5. 3椭圆低通滤波器的仿貞图图54滤波前输入信号图像图5. 4滤波后输岀信号图像第六章结语低通滤波器的应用十分广泛，利用Matlab语言，可以很容易地设il岀巴特沃斯数字 低通滤波器，与英他髙级语言的程序设计相比,MATLAB环境下更方便、快捷，节省大量的 编程时间，提高编程效率,且参数的修改也十分方便.还可以进一步进行优化设讣。巴

18、特沃斯低通滤波器的振幅对角频率单调下降，无论阶数，振幅对角频率曲线都保 持同样的形状的滤波器可以用较少的阶数获得很高的选择特性，在衰减特性要求不敏感 的场合，适合用巴特沃斯型滤波器。心得体会通过一周的设计过程，让我感触最深的就是MA77XB功能的强大性与掌握编程各 种函数和语句的重要性。当然首先要了解所要编程运行的对象的原理。在课程设计的过 程中，我深深的感受到我们所学的东四太少了，需要学习的东西太多了，真的是学无止 境。学习的过程是艰辛的，但是同时也是快乐的，让我们大家朝着自己各自的目标努力 奋斗。最后感谢老师给我们这次学习的机会，让我发现自己的不足，让我知道以后该在 实践方而更多地锻炼自己，提髙自己的动手能力。参考文献1谢平、王娜、林洪彪主编.信号处理原理及应用.机械工业出版社，2008年2陈亚勇等编著.matlab信号处理详解.人民邮电岀版社，2001年3宁彦卿等译.电子滤波器设计.科学出版社.2008年燕山大学课程设计评审意见表指导教师评语:成绩： 指导教师： 年 月 日答辩小组评语：成绩： 评阅人： 年 月 日课程设计总成绩：答辩小组成员签字：年 月 日