脉冲响应不变法设计巴特沃斯低通滤波器.docx
《脉冲响应不变法设计巴特沃斯低通滤波器.docx》由会员分享，可在线阅读，更多相关《脉冲响应不变法设计巴特沃斯低通滤波器.docx（14页珍藏版）》请在冰豆网上搜索。
脉冲响应不变法设计巴特沃斯低通滤波器
燕山大学
课程设计说明书
题目：
脉冲响应不变法设计数字低通滤波器
学院（系）：
电气工程学院
年级专业：
09级检测1班
学号：
0901********
学生姓名：
郝瑾
指导教师：
王娜
教师职称：
讲师
电气工程学院《课程设计》任务书
课程名称：
数字信号处理课程设讣
基层教学单位：
仪器科学与工程系指导教师：
王娜
学号
0901********
学生姓名
郝瑾
（专业）班级
09检测1班
设计题目
7、脉冲响应不变法设计数字低通滤波器
设计技术
参数
给定技术指标为：
几=100处，人=300股，勺=3〃B,a,=20dB,采样频率Fx=1000Hzo
设计要求
设计Butterworth低通滤波器，用脉冲响应不变法转换成数字滤波器。
参考资料
数字信号处理方面资料
MATLAB方而资料
周次
前半周
后半周
应完成内容
收集消化资料、学习MATLAB软件，进行相关参数计算
编写仿真程序、调试
指导教师签字
王娜嚟囂
说w：
i.此表一式四份.系、抬导教师、学生各一份.报送院教务科一份。
2、学生那份任务书要求装订到课程设计报告前面。
电气工程学院教务科
第一章摘要
滤波器是自动控制、信号处理和通信领域的重要组成部分，广泛地应用于各种系统中。
MATLAB语言是一种简单、髙效的髙级语言，是一种内容丰富、功能强大的分析工具,其应用范围几乎覆盖了所有的科学和工程计算领域。
MATLAB中提供了丰富的用于模拟滤波器设计的函数，通过编程可以很容易实现低通、髙通、带通、带阻滤波器，并能画岀滤波器的幅频特性曲线，大大简化了模拟滤波器设计。
本文将通过利用MATLAB滤波器设计函数直接实现巴特沃斯低通滤波器的设汁，找到应用MATLAB来设汁低通滤波器的方法。
介绍了用脉冲响应不变法设计巴特沃斯低通滤波器的基本理论和设计思想，给岀了基于MATLAB设计巴特沃斯低通滤波器的具体步骤，文中还给出了利用MATLAB产生的滤波器的幅频及相频特性曲线。
详细介绍了在基于MATLAB设计低通滤波器过程中常用到的工具和命令。
关键字：
巴特沃斯低通滤波器脉冲响应不变法
第二章引言
滤波器设讣在电子工程、应用数学和汁算机科学领域都是非常重要的内，。
而低通滤波器在信号处理中的作用等同于其它领域如金融领域中移动平均数（movingaverage）所起的作用，低通滤波器有很多种，苴中，最通用的就是巴特沃斯滤波器和切比雪夫滤波器。
巴特沃斯滤波器的特点是内的曲线最大限度平坦，没有起伏，而在则逐渐下降为零。
在振幅的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负无穷大。
巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。
只不过滤波器阶数越高，在阻频带振幅衰减速度越快。
苴他滤波器髙阶的振幅对角频率图和地结束的振幅对角频率有不同的形状。
第三章基本原理
3.1滤波器的基本理论
3.1.1滤波器的分类
滤波器分两大类：
经典滤波器和现代滤波器。
经典滤波器是假泄输入信号班“）中的有用成分和希望取出的成分各自占有不同的频带。
这样，当M"）通过一个线性系统（即滤波器）后可讲欲去除的成分有效的去除。
现代滤波器理论研究的主要内容是从含有噪声的数据记录（又称时间序列）中估计岀信号的某些特征或信号本身。
经典滤波器分为
低通.髙通.带通、带阻滤波器。
每一种又有模拟滤波器（AF）和数字滤波器（DF）o对
数字滤波器，又有IIR滤波器和FIR滤波器匚
IIRDF的转移函数是:
.w
r-（）
FIRDF的转移函数是：
FIR滤波器可以对给泄的频率特性直接进行设汁，而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设讣。
3.1.2滤波器的技术要求
低通滤波器：
"卩：
通带截止频率阻带下限截止频率
勺：
通带允许的最大衰减久：
阻带允许的最小衰减（勺，久的单位dB）
°p：
通带上限角频率阻带下限角频率
3.2模拟滤波器的设汁
3.2.1巴特沃思（Butterworth）滤波器的介绍：
因为我们设计的滤波器的冲击响应一般都为实数，所以有
G（s）G"s）=G（QG（—5）|.、.应=|G（jQJ|2
这样，如果我们能由Q",°卩，求出Ig（7O）I,那么就容易得到所需要的G（s）。
不同类型的lG（7'Q）l的表达式，代表了几种不同类型的滤波器。
巴特沃思（Butterworth）滤波器为：
lG°n）|2=TTFW
C为待定常数，N为待泄的滤波器阶次。
3.2.2巴特沃思模拟低通滤波器的设计
由于每一个滤波器的频率范用将直接取决于设讣者的应用目的，因此必然是千差万别。
为了使设计规范化，我们需要将滤波器的频率参数作归一化处理。
设所给的实际频率为G（或/），归一化后的频率为兄，对低通模拟滤波器，令
显然，Ap=\人又令归一化复数变量为P,P=j入,显然
p=以=冋％=slGp
所以巴特沃思模拟低通滤波器的设讣可按以下三个步骤来进行。
（1）将实际频率°归一化
⑵求C和N
C,N可按下式求解：
若令勺=3dB,
则C=L这样巴特沃思滤波器的设计就只剩一个参数N,这时
（3）确定G"）
因为p=j入,根据上而公式有
由1+（-1严严=0解得
.2k+N—J
2N
lh=exp（丿一——2龙）
这样可得
G*（p）=
5一八）5一PzQ
求得G（“）后，用代替变量厂即得实际需要得GG）。
3.3数字滤波器的设计
按一左规则将给出的数字滤波器的技术指标转换魏模拟低通滤波器的技术指标。
根据转换后的技术指标设计模拟低通滤波器G（$）：
再按一定的规则将G"）转换成H⑵。
若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标，然后按上述步骤2）设计出低通GG）,再将G（°转换为所需的H⑵。
第四章设计过程
4.1设计低通数字滤波器
4.1.1数字低通滤波藩设计结构图
用脉冲响应不变法转换原理如下图:
Ha（s）拉氏逆变换ha（t）等间隔采样ha（nT）二h（n）Z变换H（z）
►►►
图4.1结构框图
4.1.2设计步骤
设计要求：
脉冲响应不变法设计数字低通滤波器，给泄技术指标为：
几"00^.,A=300Hz勺=3dBav=20dB采样频率F$=1000Hz
首先应该得到角频率然后再按如下步骤求解：
（1）将数字滤波器的技术要求转换为模拟滤波器的技术要求。
（2）设计低通滤波器GG）
由几=°/°p依次求出心，人，再求岀N,可得G（P）
G（s）=G（p）亠
然后由F转换成GC）
⑶由G（s）求H⑵
4.2滤波器的MATLAB实现
4.1.1在MATLAB下，设计巴特沃斯IIR滤波器可使用butter函数。
Butter函数可设计低通、高通、带通和带阻的数字和模拟IIR滤波器，其特性为使通带内的幅度响应最大限度地平坦，但同时损失截止频率处的下降斜度。
在期望通带平滑的情况下，可使用butter函数0butter函数的用法为：
［b,a：
二butter（n,Wn,/ftype/）其中n代表滤波器阶数，Wn代表滤波器的截止频率，这两个参数可使用buttord函数来确定。
buttord函数可在给左滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n,同时给出对应的截止频率Wn。
buttord函数的用法为：
［n,Wn］二buttord（Wp,Ws,Rp,Rs）其中Wp和Ws分别是通带和阻带的拐角频率（截I匕频率），英取值范用为0至1之间。
当其值为1时代表采样频率的一半。
Rp和Rs分别是通带和阻带区的波纹系数。
不同类型（髙通、低通、带通和带阻）滤波器对应的Wp和恥值遵循以下规则：
a.高通滤波器：
Wp^Ws为一元矢量且Wp>Ws；
b.低通滤波器：
Wp和Ws为一元矢量且Wpc.带通滤波器:
Wp和Ws为二元矢量且Wpd.带阻滤波器:
Wp和Ws为二元矢量且Wp>Ws,如Wp=[0.1,0.8],Ws=[0.2,0.7]«
4.1.2Matlab的信号处理工具箱提供了频谱分析函数：
freqz函数。
freqz函数功能：
离散时间系统的频率响应。
格式：
[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（eje）在[0,Fs/2]上等间隔采样n点，采样点频率及相应频响值分别记录在f和h中。
由用户指左FS（以HZ为单位）值°
h二freqz（b,a,w）用于对H（ej）在[0,2n]上进行采样，采样频率点由矢量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和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）;%采样点数
[NzWn]=buttord（WpfWszRpzRs,'s1）;
[b,a]=butter（Nw=linspace（1,400,100）*2*pi;
H=freqs（b,a,w）;%在[0,2n]上进行采样，采样频率点由矢量w指定figure
（1）;
plot（w/（2*pi）,20*logl0（abs（H）））;
title（'巴特沃斯模拟滤波器幅频特性•）;
xlabel（1频率/Hz1）;
ylabel（1幅度/db'）;
[bzzaz]=impinvar（b,a,Fs）;
figure
（2）;
zplane（bzzaz）;
title（'脉冲响应不变法零极点分布图1）;
figure（3）;
Subplot（2,lz1）;
freqz（bz,az,Nn,Fs）;对H（ej<«））在[0,Fs/2]上等间隔采样Nn点gridon;
title（'巴特沃斯数字滤波器幅频特性•）;
xlabel「幅度/Hz'）;ylabel（1频率/db1）;
Subplot（2,lz2）;
title（'巴特沃斯数字滤波器相频特性•）;
xlabel（1幅度/Hz1）;ylabel（1频率/db'）;gridon;
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幅值
gridon;
figure（5）;
y=filter（bzzaz,x）;$对原信号进彳亍滤波…
plot（n/fs,y）;
title（'滤波前输入信号图像，）；
xlabel「时间（s）1）;
ylabel（•幅值■）;
5.2特性曲线仿真图
图5.1巴特沃斯模拟滤波器幅频特性曲线
QFigure2〔og
FileEditViewInsertToolsDesktopWindowHelp
□凰
■Q
脉冲响应不变;去的專扳点團
-10500.51
RealPari
图5.2脉冲响应不变法的零极点图
图5.3椭圆低通滤波器的仿貞•图
图5・4滤波前输入信号图像
图5.4滤波后输岀信号图像
第六章结语
低通滤波器的应用十分广泛，利用Matlab语言，可以很容易地设il•岀巴特沃斯数字低通滤波器，与英他髙级语•言的程序设计相比,MATLAB环境下更方便、快捷，节省大量的编程时间，提高编程效率,且参数的修改也十分方便.还可以进一步进行优化设讣。
巴特沃斯低通滤波器的振幅对角频率单调下降，无论阶数，振幅对角频率曲线都保持同样的形状的滤波器可以用较少的阶数获得很高的选择特性，在衰减特性要求不敏感的场合，适合用巴特沃斯型滤波器。
心得体会
通过一周的设计过程，让我感触最深的就是MA77XB功能的强大性与掌握编程各种函数和语句的重要性。
当然首先要了解所要编程运行的对象的原理。
在课程设计的过程中，我深深的感受到我们所学的东四太少了，需要学习的东西太多了，真的是学无止境。
学习的过程是艰辛的，但是同时也是快乐的，让我们大家朝着自己各自的目标努力奋斗。
最后感谢老师给我们这次学习的机会，让我发现自己的不足，让我知道以后该在实践方而更多地锻炼自己，提髙自己的动手能力。
参考文献
1谢平、王娜、林洪彪主编.信号处理原理及应用.机械工业出版社，2008年
2陈亚勇等编著.matlab信号处理详解.人民邮电岀版社，2001年
3宁彦卿等译.电子滤波器设计.科学出版社.2008年
燕山大学课程设计评审意见表
指导教师评语:
成绩：
指导教师：
年月日
答辩小组评语：
成绩：
评阅人：
年月日
课程设计总成绩：
答辩小组成员签字：
年月日
升级会员 