数字信号论文数字信号处理新技术.doc

分 数格式内容总分数字信号处理多元考核 小论文数字信号处理新技术 专 业： 电子信息工程 班 级： 牛森 姓 名： 1301 学 号： 20130506118 目 录摘要.11.引言.12.原理分析.1 2.1切比雪夫1低通滤波器的原理.2 2.2冲激响应不变法. 2 2.3双线性变换法. 33.目标分析.34.实例分析.4结束语.5致谢.5参考文献.6摘 要IIR数字滤波器在设计上可以借助成熟的模拟滤波器的 成果，如巴特沃斯，切比雪夫和椭圆滤波器等。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的设计方法设计一个满足数字滤波器指标的模拟滤波器，然后通过一定的变换，将模拟滤波器转换成数字滤波器。本文介绍了IIR数字滤波器设计，主要分析了其在上述要求下的工程方法，冲击响应不变法和双线性变换法。通过理论与实践的分析对比，表明了两种方法的优弊端。关键词 IIR数字滤波器 Matlab冲击响应不变法 双线性变换法1.引言 IIR数字滤波器的通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了相当成熟的技术和方法，有完整的设计公式，还有比较完整的图表可以查询，因此设计数字滤弄起可以充分利用这些丰富的资源来进行。对于IIR数字滤波器的设计具体步骤:(1)按照一定的规则给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。(2)根据转换后的技术指标设计模拟低通滤波器G（s）（G（s）是低通滤波器的传递函数）。（3）再按照一定的规则将G（s）转换成H（z）（H（z）是数字滤波器的传递函数）。下面从原理分析在到实践来分析设计用冲击响应不变法和双线性变换法来设计IIR数字低通滤波器的过程。2.原理分析 我们书本上重点介绍了巴特沃斯低通滤波器的设计，但巴特沃斯滤波器的频率特性无论在通带还是阻带都随频率而单调变化，因此，如果在通带边缘满足指标，则在通带内肯定会有富裕量，也就是超过指标的要求，因而并不经济。所以，更有效果的办法是将精度要求均匀的分布在通带内，或均匀分布在阻带内，或同时均匀在通带与阻带内。这样在同样通带，阻带性能要求下，就可以设计出阶数较低的滤波器。2.1切比雪夫1低通滤波器的原理设计IIR滤波器的基础是设计模拟滤波器的原型，其中切比雪夫1就是一种。切比雪夫1低通滤波器的原理：幅度平方函数(2.1.1)式中，为小于1的正数，表示通带内1幅度波动的程度，愈大，波动幅度也愈大。称为通带截止频率。令，称为对的归一化频率。CN(x)称为N阶切比雪夫多项式，定义为 (2.1.2)切比雪夫的特性：(1)切比雪夫多项式的过零点在|x|1的范围内；(2)当|x|1时，|CN(x)|1,在|x|1时，CN(x)是双曲线函数，随x单调上升。Matlab信号工具箱提供了几个直接设计IIR数字滤波器的函数，直接调用这些函数就可以很方便地对滤波器进行设计。2.2冲激响应不变法 顾名思义，所谓冲激响应不变法，就是使离散时间系统的冲激响应等于连续时间系统冲激响应的抽样值，即(2.1.3)这样，如果连续时间系统的冲激响应为 则根据式子(2.1.3)离散时间冲激响应就等于(2.1.4)由式子(2.1.4)可得连续时间系统的系统函数可见，在(2.1.4)中，pi是系统函数的极点，或者系统微分方程的特征根。这里，我们假设pi都是一阶极点，而Ai是与每个极点相对应的系数，其值为(2.1.5)而从式子(2.1.5)可求得离散时间系统的系统函数(2.1.6)根据s平面和z平面之间的映射关系s（2.1.7）平面上的j轴将映射成z平面上的单位园;而左半s平面将映射为单位园内部分；右半s平面将映射为单位园外部分,而且，从连续时间到离散时间频率是一个多值映射。根据这种映射关系可知，左半s平面上的极点将映射到z平面上的单位园内，由此而可得到一个基本结论，即冲激响应不变法可以将一个稳定的连续时间系统转换为一个稳定的离散时间系统。此外，由于可以将离散时间系统的冲激响应看作是连续时间系统的抽样，于是，根据抽样定理和连续时间频率与离散时间频率之间的关系，可以求得离散时间系统的频率响应和连续时间系统的频率响应满足下述关系 （2.1.8此式表明，如果连续时间系统的频率特性不是一个限带频谱，则转换为离散时间系统后会出现频谱混叠，从s平面到z平面的映射来看，多值映射是造成频谱混叠的原因。 2.3 双线性变换法双线性变化法是在冲激响应不变法的基础上采用不同的变换法，分为两步：第一步：将整个s平面压缩到s1平面的一条横带里；第二步：通过标准变换关系将此横带变换到整个z平面上去。为了将s平面的j轴压缩到s1平面j轴上的段上，可以通过以下的正切变换实现：（2.3.1）经过这样的频率变换，当1有时，将这一解析关系延拓整个s平面，平面的映射关系再将s1平面通过标准变换关系映射到z平面，即单值映射关系,令（2.3.2）最后得s平面与z平面的关系称为双线性变换3.目标分析 冲激响应不变法是一种简单的转换方法，它可以将一个稳定的连续时间系统转换为一个稳定的离散时间系统，而且可以使离散时间系统的冲激响应形状和连续时间系统相同，但是，由于连续时间系统的频率响应往往不是带限频谱，因此，这种方法存在着频谱混叠的现象。这种混叠将造成频谱高端的严重失真，从而使得这种方法只适用于低通滤波器或者频谱限带的高通或带通。与冲激响应不变法相比，双线性变换的主要优点：靠频谱的严重非线性关系得到S平面与Z平面的单值一一对应关系，整个j轴单值对应于单位圆一周，这个关系就是上式所表示的，其中和为非线性关系。如图图中看到，在零频率附近，-接近于线性关系，进一步增加时,增长变得缓慢，（终止于折叠频率处），所以双线性变换法不会出现于高频部分超过折叠频率而混淆到低频部分去的现象。双线性变换法的缺点：与的非线性关系，导致数字滤波器的幅度响应相对与模拟滤波器的幅频响应有畸变，（使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变）。4.实例分析下面通过实例来分析两种方法设计IIR低通滤波器。实例分析设计低通数字滤波器，要求在通带内频率低于0.2rad时，允许幅度误差在1dB以内，在频率0.3rad到rad之间的阻带衰减大于15dB。用脉冲响应不变法设计数字滤波器。T1。要求利用切比雪夫滤波器原型来设计IIR低通滤波器。需要对信号做频谱分析，为了简单对比，对信号进行频谱分析时，这里统一采用无纲量的归一化频率单位，即模拟频率对采样频率归一化；模拟角频率对采样角频率归一化；数字频率对2归一化。方法步骤先将技术指标转换为模拟低通滤波器的技术指标，再求出阶数N，查表，得出模拟低通滤波器G(s)，按一定规则G(s)转换成H(z)，并作出频谱分析的图像，利用MATLAB实现上述分析过程的程序如下：（1）冲激响应不变法：wp=0.2*pi;wr=0.3*pi;Ap=1;Ar=15;T=1;Omegap=wp/T;Omegar=wr/T;cs,ds=afd_chb1(Omegap,Omegar,Ap,Ar)%切比雪夫1型C,B,A=sdir2cas(cs,ds);%模拟滤波器级联型db,mag,pha,Omega=freqs_m(cs,ds,pi);%模拟滤波器响应subplot(224);plot(Omega/pi,mag);title(模拟滤波器幅度响应|Ha(jOmega|);b,a=imp_invr(cs,ds,T);%脉冲响应不变法设计C,B,A=dir2par(b,a)%数字低通滤波器并联db,mag,pha,w=freqz_m(b,a);%数字滤波器响应subplot(2，2，1);plot(w/pi,mag);title(数字滤波器幅度响应|H(ejomega|);subplot(2，2，2);plot(w/pi,db);title(数字滤波器幅度响应(dB);subplot(2，2，3);plot(w/pi,pha/pi);title(数字滤波器相位响应);（2）双线性变换法：wp=0.2*pi;wr=0.3*pi;Ap=1;Ar=15; T=1;Omegap=(2/T)*tan(wp/2);Omegar=(2/T)*tan(wr/2);cs,ds=afd_chb1(Omegap,Omegar,Ap,Ar)%切比雪夫原型C,B,A=sdir2cas(cs,ds);%模拟滤波器级联型db,mag,pha,Omega=freqs_m(cs,ds,pi);%模拟滤波器响应subplot(224);plot(Omega/pi,mag);title(模拟滤波器幅度响应|Ha(jOmega|);b,a=bilinear(cs,ds,T);%双线性变换法设计C,B,A=dir2cas(b,a);%数字低通滤波器级联型db,mag,pha,w=freqz_m(b,a);%数字滤波器响应subplot(2，2，1);plot(w/pi,mag);title(数字滤波器幅度响应|H(ejomega|);subplot(2，2，2);plot(w/pi,db);title(数字滤波器幅度响应(dB);subplot(2，2，3);plot(w/pi,pha/pi);title(数字滤波器相位响应);delta_w=2*pi/1000; Ap=-(min(db(1:1:wp/delta_w+1)Ar=-round(max(db(wr/delta_w+1:1:501)程序设计清晰，容易理解，脉冲响应不变法设计的切比雪夫1型数字滤波器，双线性变换法设计的切比雪夫1型数字滤波器，运行结果如图1,2。从结果图像1中可以看到脉冲响应不变法设计会造成频响的混叠，不宜用来设计高通，带阻滤波器，使用于基本上市限带的滤波器，如低通和带通滤波器，当强调以控制时间响应为目的来设计滤波器，采用这种方法比较适合。从结果图像2可以看到双线性变设计克服了频响的混叠，但频率变换关系不是线性的。结束语本次论文的设计是用脉冲响应不变法和双线性不变法设计IIR滤波器，认真阅读数字信号处理课本的相关内容，但在设计的过程中遇到一些问题，通过查阅资料， 度，也锻炼了我的动手能力，提高了自己的独立思考问题，解决问题的能力，所以，整个过程我觉得学到不少的东西。致谢 本论文是在王刚教授的悉心指导下完成的，无论是论文的选题，框架的设计，还是修改定稿，王老师都倾注了大量的心血。导师深厚的理论素养，丰富的实践知识，开拓了我的视野，使我的科研能力有了很大的提高。导师的平易近人，求真务实，积极进取的人生态度也深深影响着我。在生活上，王老师更是给了我无微不至的关怀。借此机会，想王老师表示衷心的感谢。 此外，我还得到了物理系各位领导及老师的帮助，在生活和学习上也得到了各位同学的关心和帮助，在此对大家表示感谢。 特别感谢我的父母和亲朋好友，他们的深切关怀和大力支持是我能够完成学业的