数字butterworth滤波器的设计与实现

学生姓名: 专业班级:学生姓名: 专业班级:工作单位:工作单位:指导教师:题目：数字Butterworth滤波器的设计初始条件：1.Matlab软件;2. 课程设计辅导资料：“2. 课程设计辅导资料：“Matlab理与实现”、“Matlab及在电子信息课程中的应用”等;3. 先修课程：信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1.课程设计时间：1周;2.课程设计内容：数字Butterworth 滤波器的设计，具体包括：基本 Butterworth滤波器的设计,数字高通、带通滤波器的设计，冲激响应不变法和双线性变换法的应用等;3.本课程设计统一技术要求： 研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总1.课程设计时间：1周;2.课程设计内容：数字Butterworth 滤波器的设计，具体包括：基本 Butterworth滤波器的设计,数字高通、带通滤波器的设计，冲激响应不变法和双线性变换法的应用等;3.本课程设计统一技术要求： 研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结，按要求进行实验演示和答辩等;4. 课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括:目录;与设计题目相关的理论分析、归纳和总结;与设计内容相关的原理分析、建模、推导、可行性分析;程序设计框图、程序代码（含注释） 、程序运行结果和图表、实验结果分析和总结;⑤ 课程设计的心得体会（至少 ⑤ 课程设计的心得体会（至少 500字）;时间安排:参考文献（近5年的、不少于5篇）;其它必要内容等。附——具体设计内容:1.低通巴特沃斯模拟滤波器设计。 设计一个低通巴特沃斯模拟滤波器： 指标如下：通带截止频率：fp=3400HZ,通带最大衰减： 时间安排:参考文献（近5年的、不少于5篇）;其它必要内容等。附——具体设计内容:1.低通巴特沃斯模拟滤波器设计。 设计一个低通巴特沃斯模拟滤波器： 指标如下：通带截止频率：fp=3400HZ,通带最大衰减： Rp=3dB 阻带截至频率：AS=40dB2.模拟低通转换为数字低通滤波器已知一模拟滤波器的系数函数为fs=4000HZ，阻带最小衰减:Ha（s）忌分别用冲激响应不变法和双线性变换法将Ha（S）转换成数字滤波器系统函数H（Z），并图示Ha（s）和H（Z）的幅度相应曲线。分别取采样频率Fs=1000Hz和Fs=500Hz,分析冲激响应法中存在的频率混叠失真和双线性变换法存在的非线性频率失真等。3.设计一个10阶的Butterworth 滤波器，通带为100〜250Hz,采样频率为1000Hz，绘出滤波器的单位脉冲响应。指导教师签名:系主任（或责任教师）签名:22设计总体思路低通巴特沃斯模拟滤波器设计目录2.1MATLAB提供的设计函数2.1.1巴特沃斯滤波器阶数选择函数2.1.2完整的巴特沃斯滤波器滤波器设计函数2.2设计步骤2.3MATLAB程序及运行结果2.4结果分析模拟低通转换为数字低通滤波器设计3.1冲激响应不变法3.2双线性变换法3.3MATLAB提供的设计函数3.4MATLAB程序及运行结果3.5结果分析设计10阶巴特沃斯滤波器4.1设计原理分析4.2matlab设计程序及运行结果4.3结果分析心得体会参考文献1..0 1..1.1..1.1..2.1..3 1..4..1..5 1设计总体思路IIR数字滤波器的设计，主要采用间接法，即：首先设计出低通模拟滤波器H(S);进行频率变换，将其转换为高通、带通、带阻滤波器；再用脉冲响应不变法或双线性变换法从模拟滤波器转换为数字滤波器。对单极点的N阶H(S)用部分展开式:NAH(S)K1(S-Sk)冲激响应不变法取H(S)的单阶极点Sk的指数函数"作为H(Z)的极点H⑵N—A—K1(1-e'Fsz1)21Z1双线性变换法是用S-厂门代换H(S)中的S得到H(Z)，双线性变换法可完全消除频率混叠失真但存在非线性频率失真，而冲激响应不变法存在混叠失真。在不同的设计阶段MATLAB的信号处理工具箱都给出了相应的滤波器设计函数，这些函数代表了不同类型的逼近函数的滤波器，常用的有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器类型。归一化模拟模拟域模拟高通带通冲激响应不变法f数字高通带通低通原型 归一化模拟模拟域模拟高通带通冲激响应不变法f数字高通带通低通原型 频率变换带阻滤波器j双线性变换法带阻滤波器*图1.1设计整体思路2低通巴特沃斯模拟滤波器设计2.1MATLAB提供的设计函数2.1.1巴特沃斯滤波器阶数选择函数求最小阶数N的函数buttord调用格式1:[N,Wn]=BUTTORD(Wp,Ws,Rp,Rs,'s')功能：求出巴特沃斯模拟滤波器的最小阶数N和频率参数Wn说明：Wp,Ws,为通带截止频率和阻带截止频率如式所示;WsWs2¥2 f^ass,其中fs为时域采样频率Rp,Rs为通带最大衰减Rp,Rs为通带最大衰减Apass(dB)和阻带最小衰减Astop(dB)如式Ap20lgH(ej0)H(ej2fpass)passAp20lgH(ej0)H(ej2fpass)pass(dB)As20lgH(ej0)H(ej2fstop)(dB)20lgH(ej2fpass)20lgH(ej2fstop)'s''s'对应模拟滤波器。(2)(2)模拟低通滤波器设计函数buttaP调用格式：[Z,P,K]=BUTTAP(N)功能：按给定的阶数N设计出巴特沃斯模拟低通滤波器说明：Z,P,K为返回的N阶模拟滤波器的零、极点和增益系数。如:

N(ZZj)j1H(Z)=k (ZPi)i1其中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<W<W2调用格式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阶零、极点矢量.幅频响应FREQS调用格式：[H,W]=FREQS(B,A,M)功能：此函数可以求出系统频率响应的数值解，并列出此系统的幅频及相频响应曲线。2.2设计步骤(1)确定滤波器设计指标通带内截止频率wp、阻带内截止频率WS、通带内最大衰减rp减rp、阻带内最大衰减rs。(2)运用函数buttord计算巴特沃斯滤波器的阶数N和归一化3db截止频率。运用函数butter求得低通滤波器的系统函数的分子分母多项式形式。

4)作图显示滤波器的幅频特性和相频特性。2.3MATLAB程序及运行结果％巴特沃兹滤波器的幅频响应图subplot(1,2,1)；subplot(1,2,1)；％分两个窗口，幅频图在第一个窗口wp=1000;ws=2000;rp=3;rs=40;％设置指标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]=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，和对应的频率点[h,w]=freqs(B,A,1024);％计算1024点模拟滤波器频率响应h，和对应的频率点plot(w,20*log10(abs(h)/abs(h(1))))％画频率响应幅度图grid;xIabelC频率Hz');ylabel('幅度(dB)');%给x轴和y轴加标注title('巴特沃斯幅频响应title('巴特沃斯幅频响应')%给图形加标题axis([0,3000,-40,3]);line([0,2000],[-3,-3]);%line([0,2000],[-3,-3]);%为方便观察加上线条line([1000,1000],[-40,3]);逼Figure1- □ X文件迥潟辑反1 插入①IMCE 窗拜葡凹 爼Q出日◎R\ ® 則□目B□图2.1低通巴特沃斯模拟滤波器2.4结果分析巴特沃斯低通滤波器的特点是具有通带内最大平坦的幅度特性， 而且随着频率的升高而单调地下降。它的平方幅度响应为|H(j)|21(―)2nc其中，n为整数，称为滤波器的阶数,n值越大通带和阻带的近似性越好，过渡带也越陡。c为低通滤波器的截止频率。该滤波器具有一些特殊的性质①对所有的n，都有当=0时，IH(j)|2=1；1②对所有的n，都有当c时，|H(j)|2=-，即在c处有②对所有的n，都有当2

③|H(j)|2③|H(j④当n④当n时，巴特沃斯滤波器趋向于理想的低通滤波器⑤在=0⑤在=0处平方幅度响应的各级导数均存在且等于0，因此|H(j)|2在该点上取得最大值，且具有最大平坦特性。得最大值，且具有最大平坦特性。itj itj 3模拟低通转换为数字低通滤波器设计3.1冲激响应不变法冲激响应不变法是使数字滤波器的单位冲激响应序列 h(n)模仿模拟滤波器的单位冲激响应ha(t)，将模拟滤波器的单位冲激响应加以等间隔抽样， 使h(n)正好等于ha(t)的抽样值，即满足:h(n)hh(n)ha(nT)式中：T为抽样周期。冲激不变法把稳定的H冲激不变法把稳定的Ha(S)转换为稳定的H(z)。由此方法可得到一阶系统的最基本的转换关系为:3.2双线性变换法双线性变换法式使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法，为了克服多值映射这一缺点，先将S平面压缩成s1平面上一个宽度为2n的水平带状区域，然后通过玄状区域，然后通过玄=总将这个带状区域映射到整个z平面，实现S平面到z平面的单值映射。iiri切-\riiri切-\rw=[0:1000]*2*pi;% w=[0:1000]*2*pi;% 模拟频率为2f,其中f取0~1000Hz图3.1S平面到Z平面映射关系其S平面到Z平面的单值映射关系如下:2£二―T2—十日一2——ST3.3MATLAB提供的设计函数模拟数字变换函数—双线性变换函数 bilinear或脉冲响应不变法函数impinvar调用格式：[NUMdQENd]=BILINEAR(NUM,DEN,Fs)功能：把模拟滤波器系数为NUMQEN变为近似等价的数字滤波器系数为NUMd,DENd调用格式：[BZ,AZ]=IMPINVAR(B,A,Fs)功能：设计出数字滤波器其单位抽样响应h(n)为模拟滤波器的单位冲激响应h(t)的采样值，采样频率为Fs说明：B,A为模拟滤波器的系统函数H(S)的分子和分母多项式的系数，如下所示:H(S)鸞H(S)鸞A(S)a1s a2SbmS bm1anS an1B=[b1,b2,…,bm,bm+1],A=[a1,a2, …,an,an+1]3.4MATLAB程序及运行结果b=1;a=[1,1000];

[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~1title（'[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~1title（'冲激响应不变法数字频率响应'）subplot(2,2,3);plot(wz/pi,abs(hz2)/hz2(1));grid;%画出双线性变换法滤波器的幅频图，holdon;axis([0,1,0,1])%数字频率wz归一化为0~1titleC双线性变换法数字频率响应');end;甫Figure1 — □X件曲福lE蛋君迪廷入tu工麹肯助奁awo辛曲(H)□日日&fe鱼严®C摇，a「口®I3口檯加叔車11宾屉00.501300'OSflf(Hz)収线性空檢沃茨檯加叔車11宾屉00.501300'OSflf(Hz)収线性空檢沃茨1^頻率响咸图3.1模拟滤波器转换为数字滤波器3.5结果分析冲激响应不变法使得数字滤波器的h(n)能完全模仿模拟滤波器的h(t)，即时域逼近良好；而且模拟频率和数字频率之间呈线性关系，因此频率之间不存在失真。但是由于映射nV简单的代数映射，根据奈奎斯特抽样定理，只有当模拟滤波器的频率响应是严格限带的，且带限于折叠频率［s/2，s/2］以内时，才能使数字滤波器的频率响应在折叠频率以内，重现模拟滤波器的频率响应而不产生混叠失真。但是，任何一个实际的模拟滤波器响应都不是严格限带的，变换后都会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真，因此模拟滤波器的频率响应在折叠频率以上衰减越大、越快，变换后频率响应混叠失真就越小。冲激响应不变法的主要优点是h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好；线性相位模拟滤波器转变为线性相位数字滤波器。缺点为对时域的采样会造成频域的混叠效应”，故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大，故不能用来设计高通和带阻滤波器。只适用于限带的低通、带通滤波器。双线性变换法最大的优点是避免了频率响应的混叠失真。如仿真波形图可知，在零频率附近，arctan(/2)的频率变换关系接近线性关系外，当增加到高频段时，变换关系就是非线性的了，也就是说，与之间存在严重的非线性关系。这就是双线性变换法存在的非线性频率失真。它的主要优点是避免了频率响应的混迭失真现象，在特定数字滤波器和特定模拟滤波器处，频率响应是严格相等的，它可以较准确地控制截止频率的位置。此外，由于它是与之间严重非线性，种简单的代数关系，故设计十分方便。但缺点在于，在高频段，即线性相位模拟滤波器变为非线性相位数字滤波器；其次这种非线性关系要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变，对于分段常数型AF滤波器，经双线性变换与之间严重非线性，后，仍得到幅频特性为分段常数的 DF.但在各个分段边缘的临界频率点产生畸变，这种频率的畸变，可通过频率预畸变加以校正。设计10阶巴特沃斯滤波器4.1设计原理分析由题目可知，要求设计一个10阶的带通的巴特沃斯滤波器，理论上所用的方法可以是先设计出模拟低通滤波器，再转变成数字带通滤波器。然而，这里我们可以直接利用

butter函数直接设计IIR数字带通滤波器。[b,a]butter(N,Wn,'ftype')用于设计高通、带通或带阻滤波器，并由参数 ftype确定滤波器的形式。当ftype(2)=band(默认值)时，为带通滤波器，此时Wn=[wlwh]，则返回的[b,a]所构成的滤波器是阶次为2N的带通滤波器，带通范围为：l4.2matlab 设计程序及运行结果n=10;Wn=[100250]/1000;[b,a]=butter(n,Wn);[y,t]=i