数字butterworth滤波器的设计与实现.docx
- 文档编号:11105595
- 上传时间:2023-02-25
- 格式:DOCX
- 页数:16
- 大小:164.29KB
数字butterworth滤波器的设计与实现.docx
《数字butterworth滤波器的设计与实现.docx》由会员分享,可在线阅读,更多相关《数字butterworth滤波器的设计与实现.docx(16页珍藏版)》请在冰豆网上搜索。
数字butterworth滤波器的设计与实现
课程设计任务书
学生姓名:
专业班级:
指导教师:
工作单位:
题目:
数字Butterworth滤波器的设计
初始条件:
1.Matlab软件;
2.课程设计辅导资料:
“Matlab语言基础及使用入门”、“信号与系统”、“数字信号处理原理与实现”、“Matlab及在电子信息课程中的应用”等;
3.先修课程:
信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。
要求完成的主要任务:
(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)
1.课程设计时间:
1周;
2.课程设计内容:
数字Butterworth滤波器的设计,具体包括:
基本Butterworth滤波器的设计,数字高通、带通滤波器的设计,冲激响应不变法和双线性变换法的应用等;
3.本课程设计统一技术要求:
研读辅导资料对应章节,对选定的设计题目进行理论分析,针对具体设计部分的原理分析、建模、必要的推导和可行性分析,画出程序设计框图,编写程序代码(含注释),上机调试运行程序,记录实验结果(含计算结果和图表),并对实验结果进行分析和总结,按要求进行实验演示和答辩等;
4.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写,具体包括:
1目录;
2与设计题目相关的理论分析、归纳和总结;
3与设计内容相关的原理分析、建模、推导、可行性分析;
4程序设计框图、程序代码(含注释)、程序运行结果和图表、实验结果分析和总结;
5课程设计的心得体会(至少500字);
6参考文献(近5年的、不少于5篇);
7其它必要内容等。
时间安排:
1周
附——具体设计内容:
1.低通巴特沃斯模拟滤波器设计。
设计一个低通巴特沃斯模拟滤波器:
指标如下:
通带截止频率:
=3400HZ,通带最大衰减:
=3dB阻带截至频率:
=4000HZ,阻带最小衰减:
=40dB
2.模拟低通转换为数字低通滤波器已知一模拟滤波器的系数函数为
分别用冲激响应不变法和双线性变换法将
转换成数字滤波器系统函数
,并图示
和
的幅度相应曲线。
分别取采样频率
=1000Hz和
=500Hz,分析冲激响应法中存在的频率混叠失真和双线性变换法存在的非线性频率失真等。
3.设计一个10阶的Butterworth滤波器,通带为100~250Hz,采样频率为1000Hz,绘出滤波器的单位脉冲响应。
指导教师签名:
年月日
系主任(或责任教师)签名:
年月日
1设计总体思路
IIR数字滤波器的设计,主要采用间接法,即:
首先设计出低通模拟滤波器H(S);进行频率变换,将其转换为高通、带通、带阻滤波器;再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤波器。
对单极点的N阶H(S)用部分展开式:
冲激响应不变法取H(S)的单阶极点Sk的指数函数
作为H(Z)的极点
双线性变换法是用
代换H(S)中的S得到H(Z),双线性变换法可完全消除频率混叠失真但存在非线性频率失真,而冲激响应不变法存在混叠失真。
在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的滤波器设计函数,这些函数代表了不同类型的逼近函数的滤波器,常用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。
模拟域冲激响应不变法
频率变换双线性变换法
图1.1设计整体思路
2低通巴特沃斯模拟滤波器设计
2.1MATLAB提供的设计函数
2.1.1巴特沃斯滤波器阶数选择函数
(1)求最小阶数N的函数buttord
调用格式1:
[N,Wn]=BUTTORD(Wp,Ws,Rp,Rs,'s')
功能:
求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn
说明:
Wp,Ws,为通带截止频率和阻带截止频率如式所示;
,
Rp,Rs为通带最大衰减Apass(dB)和阻带最小衰减Astop(dB)如式
's'对应模拟滤波器。
(2)模拟低通滤波器设计函数buttap
调用格式:
[Z,P,K]=BUTTAP(N)
功能:
按给定的阶数N设计出巴特沃斯模拟低通滤波器
说明:
Z,P,K为返回的N阶模拟滤波器的零、极点和增益系数。
如:
其中k为增益常数,zj和pi为系统函数H(Z)的N个零点和N个极点。
2.1.2完整的巴特沃斯滤波器滤波器设计函数
(1)合为一步的数字滤波器设计函数butter
调用格式1:
[B,A]=BUTTER(N,Wn)
功能:
设计N阶低通数字滤波器,系数矢量B,A按Z的降幂排列,当Wn为二元矢量时,Wn=[W1W2],返回一个2N阶带通滤波器,通带范围:
W1 调用格式2: [B,A]=BUTTER(N,Wn,'high')设计高通数字滤波器. 调用格式3: [B,A]=BUTTER(N,Wn,'stop')设计带阻滤波器Wn=[W1,W2]. 调用格式4: [Z,P,K]=BUTTER(N,Wn)返回低通数字滤波器的N阶零、极点矢量. (2)幅频响应FREQS 调用格式: [H,W]=FREQS(B,A,M) 功能: 此函数可以求出系统频率响应的数值解,并列出此系统的幅频及相频响应曲线。 2.2设计步骤 (1)确定滤波器设计指标通带内截止频率wp、阻带内截止频率ws、通带内最大衰减rp、阻带内最大衰减rs。 (2)运用函数buttord计算巴特沃斯滤波器的阶数N和归一化3db截止频率。 (3)运用函数butter求得低通滤波器的系统函数的分子分母多项式形式。 (4)作图显示滤波器的幅频特性和相频特性。 2.3MATLAB程序及运行结果 %巴特沃兹滤波器的幅频响应图 subplot(1,2,1);%分两个窗口,幅频图在第一个窗口 wp=1000;ws=2000;rp=3;rs=40;%设置指标 [N,wn]=buttord(wp,ws,rp,rs,'s')%计算巴特沃斯低通滤波器的阶数和3dB截止频率 [B,A]=butter(N,wn,'s');%代入N和Wn设计巴特沃斯模拟低通滤波器 [Z,P,K]=buttap(N);%计算滤波器的零、极点 [h,w]=freqs(B,A,1024);%计算1024点模拟滤波器频率响应h,和对应的频率点w plot(w,20*log10(abs(h)/abs(h (1))))%画频率响应幅度图 grid; xlabel('频率Hz');ylabel('幅度(dB)');%给x轴和y轴加标注 title('巴特沃斯幅频响应')%给图形加标题 axis([0,3000,-40,3]); line([0,2000],[-3,-3]);%为方便观察加上线条 line([1000,1000],[-40,3]); 图2.1低通巴特沃斯模拟滤波器 2.4结果分析 巴特沃斯低通滤波器的特点是具有通带内最大平坦的幅度特性,而且随着频率的升高而单调地下降。 它的平方幅度响应为 其中,n为整数,称为滤波器的阶数,n值越大通带和阻带的近似性越好,过渡带也越陡。 为低通滤波器的截止频率。 该滤波器具有一些特殊的性质: ①对所有的n,都有当 =0时, =1; ②对所有的n,都有当 = 时, = ,即在 处有3dB的衰减; ③ 是的单调递减函数,即不会出现幅度响应的起伏; ④当n 时,巴特沃斯滤波器趋向于理想的低通滤波器; ⑤在 =0处平方幅度响应的各级导数均存在且等于0,因此 在该点上取得最大值,且具有最大平坦特性。 3模拟低通转换为数字低通滤波器设计 3.1冲激响应不变法 冲激响应不变法是使数字滤波器的单位冲激响应序列 模仿模拟滤波器的单位冲激响应 ,将模拟滤波器的单位冲激响应加以等间隔抽样,使 正好等于 的抽样值,即满足: 式中: T为抽样周期。 冲激不变法把稳定的 转换为稳定的 。 由此方法可得到一阶系统的最基本的转换关系为: 3.2双线性变换法 双线性变换法式使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法,为了克服多值映射这一缺点,先将s平面压缩成s1平面上一个宽度为2π/T的水平带状区域,然后通过 将这个带状区域映射到整个z平面,实现s平面到z平面的单值映射。 图3.1S平面到Z平面映射关系 其S平面到Z平面的单值映射关系如下: 3.3MATLAB提供的设计函数 模拟数字变换函数-双线性变换函数bilinear或脉冲响应不变法函数impinvar 调用格式: [NUMd,DENd]=BILINEAR(NUM,DEN,Fs) ✧功能: 把模拟滤波器系数为NUM,DEN变为近似等价的数字滤波器系数为NUMd,DENd 调用格式: [BZ,AZ]=IMPINVAR(B,A,Fs) ✧功能: 设计出数字滤波器其单位抽样响应h(n)为模拟滤波器的单位冲激响应h(t)的采样值,采样频率为Fs ✧说明: B,A为模拟滤波器的系统函数H(S)的分子和分母多项式的系数,如下所示: B=[b1,b2,…,bm,bm+1],A=[a1,a2,…,an,an+1] 3.4MATLAB程序及运行结果 b=1;a=[1,1000]; w=[0: 1000]*2*pi;%模拟频率为2f,其中f取0~1000Hz [h,w]=freqs(b,a,w);%计算模拟滤波器的频率响应 subplot(2,2,1) plot(w/2/pi,abs(h)/abs(h (1)));grid;%画模拟滤波器幅频特性 title('模拟频率响应'); xlabel('f(Hz)');ylabel('幅度'); fs0=[1000,500]; form=1: 2; Fs=fs0(m); [bz,az]=impinvar(b,a,Fs);%冲激响应不变法设计数字滤波器 [bzl,azl]=bilinear(b,a,Fs);%双线性变换法设计数字滤波器 wz=[0: pi/512: pi]; hz1=freqz(bz,az,wz); hz2=freqz(bzl,azl,wz); subplot(2,2,2);plot(wz/pi,abs(hz1)/hz1 (1));grid;%画出冲激响应不变法滤波器的幅频图, holdon; axis([0,1,0,1])%数字频率wz归一化为0~1 title('冲激响应不变法数字频率响应') subplot(2,2,3);plot(wz/pi,abs(hz2)/hz2 (1));grid;%画出双线性变换法滤波器的幅频图, holdon; axis([0,1,0,1])%数字频率wz归一化为0~1 title('双线性变换法数字频率响应'); end; 图3.1模拟滤波器转换为数字滤波器 3.5结果分析 冲激响应不变法使得数字滤波器的h(n)能完全模仿模拟滤波器的h(t),即时域逼近良好;而且模拟频率和数字频率之间呈线性关系,因此频率之间不存在失真。 但是由于映射 简单的代数映射,根据奈奎斯特抽样定理,只有当模拟滤波器的频率响应是严格限带的,且带限于折叠频率 以内时,才能使数字滤波器的频率响应在折叠频率以内,重现模拟滤波器的频率响应而不产生混叠失真。 但是,任何一个实际的模拟滤波器响应都不是严格限带的,变换后都会产生周期延拓分量的频谱交叠,即产生频率响应的混叠失真,因此模拟滤波器的频率响应在折叠频率以上衰减越大、越快,变换后频率响应混叠失真就越小。 冲激响应不变法的主要优点是h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好;线性相位模拟滤波器转变为线性相位数字滤波器。 缺点为对时域的采样会造成频域的“混叠效应”,故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大,故不能用来设计高通和带阻滤波器。 只适用于限带的低通、带通滤波器。 双线性变换法最大的优点是避免了频率响应的混叠失真。 如仿真波形图可知,在零频率附近, 的频率变换关系接近线性关系外,当 增加到高频段时,变换关系就是非线性的了,也就是说, 与 之间存在严重的非线性关系。 这就是双线性变换法存在的非线性频率失真。 它的主要优点是避免了频率响应的混迭失真现象,在特定数字滤波器和特定模拟滤波器处,频率响应是严格相等的,它可以较准确地控制截止频率的位置。 此外,由于它是一种简单的代数关系,故设计十分方便。 但缺点在于,在高频段, 与 之间严重非线性,即线性相位模拟滤波器变为非线性相位数字滤波器;其次这种非线性关系要求模拟滤波器的幅频响应为分段常数型,不然会产生畸变,对于分段常数型AF滤波器,经双线性变换后,仍得到幅频特性为分段常数的DF.但在各个分段边缘的临界频率点产生畸变,这种频率的畸变,可通过频率预畸变加以校正。 4设计10阶巴特沃斯滤波器 4.1设计原理分析 由题目可知,要求设计一个10阶的带通的巴特沃斯滤波器,理论上所用的方法可以是先设计出模拟低通滤波器,再转变成数字带通滤波器。 然而,这里我们可以直接利用butter函数直接设计IIR数字带通滤波器。 用于设计高通、带通或带阻滤波器,并由参数ftype确定滤波器的形式。 当ftype (2)=band(默认值)时,为带通滤波器,此时Wn=[wlwh],则返回的[b,a]所构成的滤波器是阶次为2N的带通滤波器,带通范围为: 。 4.2matlab设计程序及运行结果 n=10; Wn=[100250]/1000; [b,a]=butter(n,Wn); [y,t]=impz(b,a,101); stem(t,y) 图4.1十阶巴特沃斯数字滤波器 4.3结果分析 在阶数已知的情况下可以直接利用butter函数来设计巴特沃斯滤波器,滤波器类型则由其中的ftype参数控制。 5心得体会 本次数字信号处理的课程设计主要是关于MATLAB的运用,利用MATLAB来进行巴特沃斯滤波器的设计。 MATLAB虽然早在信号与系统课程中就已接触到,但对其运用的并不是十分娴熟,对其函数库的理解和使用并不十分熟练。 虽然MATLAB函数库提供课很多关于巴特沃斯滤波器设计的函数,但要想完成本次设计,仅有这些是不够的。 首先要对设计的原理和理论有充分的把握,期间反反复复看了书上关于IIR滤波器设计的内容,对IIR滤波器设计的方法(冲激响应不变法、双线性变换法、频率变换法)有了充分的认识。 其次通过网上查阅MATLAB中巴特沃斯滤波器设计的相关函数(butter、buttap、buttord、freqs等)及其使用方法,最后通过编写完整的程序完成任务要求。 本次课程设计内容并不算多,但时间的规划和处理也很重要,三个设计内容上是递进的关系,一步一步的达到用模拟低通滤波器来设计数字滤波器的目的。 这样的安排由易到难,使我对设计的过程映像更加的深刻。 6参考文献 [1]刘泉,阙大顺.数字信号处理原理与实现.电子工业出版社.2008 [2]丁伟.精通MATLABR2014a.清华大学出版社.2015 [3]阙大顺.数字信号处理学习辅导与考研指导.武汉理工大学出版社.2007 [4]赵健.数字信号处理.清华大学出版社.2011 [5]李正周.MATLAB数字信号处理与应用.2008
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字 butterworth 滤波器 设计 实现