双线性变换巴特沃斯IIR数字带通滤波器

双线性变换巴特沃斯IIR数字带通滤波器 ------------------------------------------------------------------------------------------------ 双线性变换巴特沃斯IIR数字带通滤波器 课 程 设 计 报 告 课程名称:专业综合课程设计学生姓名:学 号:专业班级:指导教师: 完成时间: 报告 软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载 成绩: IIR数字带通滤波器 1.课程设计目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法;熟悉用双线性变换法设计IIR 数字带通滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法。掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 2.课题要求 采用双线性变换法设计一数字带通滤波器，抽样频率为 求为:通带范围从250 480HzHzfs?1kHz，性能要Hz到400Hz，在此两频率处衰减不大于3dB,在150和频率处衰减不小于20dB，采用巴特沃思型滤波器。 3.设计原理 3.1 数字滤波器介绍 滤波器，顾名思义，其作用是对输入信号起到滤波作用。数字滤波器(DF，Digital Filter)在数字信号处理中起着重要作用。数字滤波—————————————————————————————————————— ------------------------------------------------------------------------------------------------ 器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。 数字滤波器有低通(LP,Low pass)、高通(HP,High Pass)、带通(BS,Band PASS)、带阻(BS,Band Stop)和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。应用最广的是线性、时不变数字滤波器，以及FIR滤波器。 数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语音信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。 3.2 巴特沃思的原理 巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。 3.3 双线性变换原理 为了克服冲激响应法可能产生的频率响应的混叠失真，这是因为从S平面到,平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T,π/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T,π/T一条横带里;第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，—————————————————————————————————————— ------------------------------------------------------------------------------------------------ 也就消除了频谱混叠现象，映射关系如图1-1 S平面S1平面Z平面 图1-1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现 ????T?tan?1?T ?2?2(1) 式中,T仍是采样间隔。 当Ω1由-π/T经过0变化到π/T时，Ω由-?经过0变化到+?，也即映射了整个jΩ轴。将式写成 j?? 2T ? ee j?1T/2j?1T/2 ?e?e j?1T/2 (2) ?j?1T/2 将此关系解析延拓到整个S平面和S1平面，令jΩ=s，jΩ1=s1，则得 s? 2T?ee —————————————————————————————————————— ------------------------------------------------------------------------------------------------ s1T/2s1T/2 ?e?e ?s1T/2?s1T/2 21?e1?s1T? ?tanh?????sT T2T1?e1?? 2 ?sT (3) 再将S1平面通过以下标准变换关系映射到Z平面 z=es1T 从而得到S平面和Z平面的单值映射关系为: s? 21?zT1?z ?1?1 (4) 1?z? T2T2 ss 2?T2T ?s (5) ?s —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 1? S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换 4.设计思路 4.1设计步骤 (1).数字通带滤波器的技术指标转换为模拟低通滤波器的技术指标。 (2).设计模拟低通滤波器G(p)。 (3).将模拟低通滤波器G(p)转换成数字带通滤波器H(z)。 4.2 设计过程 (1)确定技术指标。 a.数字带通滤波器的技术指标(T ?1fs )。 2*pi*fp1*T?0.5?2*pi*fp3*T?0.8? 通带截止频率:?1? ?3? 阻带截止频率: ?sl ?sh ?2*pi*fsl*T?0.3? ?2*pi*fsh*T?0.96? 通带最大衰减:?p 最小阻带衰减:?p 抽样频率: ?3dB ?20dB —————————————————————————————————————— ------------------------------------------------------------------------------------------------ fs?1kHz b.数字带通滤波器转换成模拟带通滤波器技术指标。 通带截止频率:?1? ??1?tan???2000radsT?2?2 ??3?tan???6155radT2??2 ??sl?tan???1019radT?2?2 ??sh?tan???31789radT2??2 s s ?3? 阻带截止频率:?sl ?s ?sh? 中心频率:?2 2 ??1?3?1231000 通带带宽:?BW c.模拟低通技术指标. ??3??1?4155rads 以?BW为参考频率将Ω归一化，得 ?1? ?1?BW?3?BW —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ? 2000415561554155 ?0.4813 ?3???1.4813 ?sl? ?sl?BW?sh?BW ? 10194155317894155 ?0.2452 ?sh???7.6507 显然 ?BW 2 ??3??1?1 ?2??1?3?0.7129 2 由? ? ???2 ? 2 ,将?转换为低通滤波器G?p?的归一化频率?, 可知?p?1，但?s不唯一，即 ??s? —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?sl??2 ?sl 2 22 ??2.6622 ?s? ?sh??2 ?sh 2 22 ?7.5575 由于求出的两个?s差别较大，为了保证滤波器的衰减性，应取?s 为最小者，固取?s? 2.7。 (2).设计巴特沃思模拟低通滤波器G?p?。 根据?p ?1，?s?2.7， ?p?3dB ，?s? 20dB ，设计G?p?。 ?s2 lg??10?1? ?3?10 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?plg?? 10?1? N? ?1??? ???100.?1??2lg??s? ? 2lg?2.7? ?2.3 故整数N=3。 此时三阶巴特沃思模拟低通滤波器的系统函数为 G?p?? 1 p3 ?2p2 ?2p?1 (3).将 G(p)转换为数字带通滤波器Hd(z)。 由 p ? ?s 2 ??1?3? s —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 2 ??2 2? s??3??1? ? ?s? , BW 而 s ? 21-z -1T1?z -1 则p与z之间转换关系，即 6 ?12 2?12 p? 4?10 ?1?z???2 ?1?z? 2?103 1?z?2 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ? BW 将?22 ，?BW 带入上式，得 1p? 0.4813?1?z ?1?2 ?1.4813?1?z ?? 2 1?z?2 再将所得的结果代入G?p?中即可得数字带通滤波器H?z?。 0.04953?1?z ?64 H?z?? ?3z ??3z ?2? 1?0.1391z ?6 ?0.5704z ?5 ?0.9130z?4 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ ?1.9184z?3 ?1.981z ?2 ?2.1405z ?1 5.实验程序及结果 5.1实验程序 clear;close all t=0.001; fs=1000; wpu=0.8*pi,wpl=0.5*pi; wsu=0.96*pi,wsl=0.3*pi; wpz=[0.5,0.8]; wsz=[0.3,0.96]; wp=2/t*tan(wpz/2);ws=2/t*tan(wsz/2); rp=3;as=20; [n,wc]=buttord(wp,ws,rp,as,'s');%计算带通滤波器阶数N和3dB截止频率wc [b,a]=butter(n,wc,'s'); %计算带通滤波器系统函数分子分母多项式系数向量b,a [bz,az]=bilinear(b,a,fs); [nd,wdc]=buttord(wpz,wsz,rp,as); [bd,adz]=butter(nd,wdc); —————————————————————————————————————— ------------------------------------------------------------------------------------------------ hk=freqz(bd,adz) figure(1) plot(abs(hk)); grid on; xlabel('Hz'); ylabel(' 幅度'); title('幅度函数曲线') figure(2) plot(angle(hk)); grid on; xlabel('\omega/\pi'); ylabel('phi(omega)'); title('相频特性曲线') 5.2 MATLAB结果 \ 6.心得体会 我这次设计的是采用双线性变换法设计数字带通滤波器。经过这段时间的学习设计，能够参考 书 关于书的成语关于读书的排比句社区图书漂流公约怎么写关于读书的小报汉书pdf 籍拓展出来的东西运用到自己的课题上，这是目的的学习。说实在，这东西看这就不易，比如书纸上的东西对我这次设计的促进不够实际，也就是将不详尽。所以只得参考其他的资料，资料的繁多外加冗杂迫使我本来不太清楚的问题，愈加暴露。这就要学会选择，精炼而实际有用的知识。同时老师及同学帮辅也是不可缺少的一部分，因为自身认知的局限导致不可避免的盲区出现，使得设计进度及过程很迟缓。最终，还是完成了。 —————————————————————————————————————— ------------------------------------------------------------------------------------------------ 7.参考资料 [1] 胡广书，数字信号处理:理论、算法与实现.北京:清华大学出版社，1997 [2] 飞科科技产品研发中心，MATLAB7辅助信号处理机应用.北京:电子工业出版社，2005 [3] 高西全、丁玉美、阔永红，数字信号处理:原理、实现.北京，电子工业出版社，2010 [4] 吴镇扬，信号处理原理.北京:高等教育出版社，2012 ——————————————————————————————————————