数字切比雪夫滤波器的设计综述.docx
《数字切比雪夫滤波器的设计综述.docx》由会员分享,可在线阅读,更多相关《数字切比雪夫滤波器的设计综述.docx(23页珍藏版)》请在冰豆网上搜索。
数字切比雪夫滤波器的设计综述
电气工程学院《课程设计》任务书
课程名称:
数字信号处理课程设计
基层教学单位:
仪器科学与工程系指导教师:
王娜
学号
学生姓名
(专业)班级
设计题目
21切比雪夫滤波器设计
设
计
技
术
参
数
采样频率100Hz,低频、中频、高频信号频率分别为5Hz、15Hz、30Hz
设
计
要
求
产生一个连续信号,包含低频,中频,高频分量,对其进行采样,进行频谱分析。
设计低通,高通,带通滤波器对信号进行滤波处理,观察滤波后信号的频谱。
分析该类型滤波器与其他类型低通滤波器(如Butterworth)优势及特点。
参
考
资
料
数字信号处理方面资料
MATLAB方面资料
周次
前半周
后半周
应
完
成
内
容
收集消化资料、学习MATLAB软件,进行相关参数计算
编写仿真程序、调试
指导教
师签字
基层教学单位主任签字
说明:
1、此表一式四份,系、指导教师、学生各一份,报送院教务科一份。
2、学生那份任务书要求装订到课程设计报告前面。
电气工程学院教务科
摘要
随着信息和数字时代的到来,数字信号处理已成为当今一门极其重要的学科和技术领域。
在现代通信系统中,由于信号中经常混有各种复杂成分,因此很多信号的处理都是基于滤波器而进行的。
所以,数字滤波器在数字信号处理中起着举足轻重的作用。
数字滤波器是数字信号处理的重要工具之一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。
故本课题使用MATLAB信号处理箱和运用切比雪夫法设计数字滤波器。
数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。
它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。
数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与l两个电平状态)、灵活性强等优点。
第一信号处理设计原理-----------------------------------------4
1.1MATLAB介绍------------------------------------------------------4
1.2数字滤波器的介绍---------------------------------------------------4
1.3滤波器的分类--------------------------------------------------------5
1.4数字滤波器的设计方法-------------------------------------------5
1.5切比雪夫滤波器简介------------------------------------------------5
1.6模拟低通滤波器---------------------------------------------------6
第二用MATLAB设计切比低通雪夫滤波器----------------7
2.1MATLAB相关函数------------------------------------------------7
2.2数字高通滤波器的设计-------------------------------------------7
2.3数字带通滤波器的设计-----------------------------------------8
2.4双线性变换---------------------------------------------------------9
第三切比雪夫数字滤波器的设计过程及程序,仿真-------11
3.1切比雪夫型低通数字滤波器----------------------------------11
3.2切比雪夫型高通数字滤波器----------------------------------14
3.3切比雪夫型带通数字滤波器----------------------------------17第四总结----------------------------------------------------------20
第五设计结果及分析---------------------------------------------21
5.1滤波器与其他类型低通滤波器优势及特点--------------22
5.2心得体------------------------------------------------------------22
第六参考文献----------------------------------------------------23
1信号处理设计原理
1.1MATLAB介绍
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
MATLAB和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件中在数值计算方面首屈一指。
MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。
MATLAB产品族可以用来进行以下各种工作:
数值分析、数值和符号计算、工程与科学绘图、控制系统的设计与仿真、数字图像处理技术、数字信号处理技术、通讯系统设计与仿真、财务与金融工程。
MATLAB的应用范围非常广,包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。
附加的工具箱(单独提供的专用MATLAB函数集)扩展了MATLAB环境,以解决这些应用领域内特定类型的问题。
1.2数字滤波器的介绍
数字滤波器处理模拟信号时,首先须对输入模拟信号进行限带、抽样和模数转换。
数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍,其频率响应具有以抽样频率为间隔的周期重复特性,且以折叠频率即1/2抽样频率点呈镜像对称。
为得到模拟信号,数字滤波器处理的输出数字信号须经数模转换、平滑。
数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。
数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。
1.3滤波器的分类
1.从功能上分;低、带、高、带阻。
2.从实现方法上分:
FIR、IIR
3.从设计方法上来分:
Chebyshev(切比雪夫),Butterworth(巴特沃斯)
4.从处理信号分:
经典滤波器、现代滤波器
1.4数字滤波器的设计方法
它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。
应用最广的是线性、时不变数字滤波器,以及FIR滤波器。
由模拟滤波器设计IIR数字滤波器,必须建立好s平面和z平面的映射关系。
使模拟系统函数
变换成数字滤波器的系统函数
,通常采用冲激相应不变法和双线性变换法。
冲激相应不变法存在频谱混叠现象,双线性变换法消除了这一现象,在IIR数字滤波器的设计中得到了更广泛的应用。
s平面和Z平面的映射关系为
将
和
带入数字频率和
等效的模拟频率之间的映射关系:
,由于二者不是线性关系,所以称为预畸变。
1.5切比雪夫滤波器简介
切比雪夫滤波器特点:
误差值在规定的频段上等波纹变化。
切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹波动的滤波器。
在通带波动的为“I型切比雪夫滤波器”,在阻带波动的为“II型切比雪夫滤波器”。
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
这种滤波器来自切比雪夫多项式,因此得名,用以纪念俄罗斯数学家巴夫尼提•列波维其•切比雪夫切比雪夫滤波器特点:
误差值在规定的频段上等波纹变化。
为了从模拟滤波器出发设计IIR数字滤波器,必须先设计一个满足技术指标的模拟滤波器,亦即要把数字滤波器的指标转换成模拟滤波器的指标,因此必须先设计对应的模拟原型滤波器。
1.6模拟低通滤波器
巴特沃兹滤波器在通带内幅度特性是单调下降的,如果阶次一定,则在靠近截止
处,幅度下降很多,或者说,为了使通带内的衰减足够小,需要的阶次N很高,为了克服这一缺点,采用切比雪夫多项式来逼近所希望的
。
切比雪夫滤波器的
在通带范围内是等幅起伏的,所以在同样的通常内衰减要求下,其阶数较巴特沃兹滤波器要小。
切比雪夫滤波器的振幅平方函数为
(1-4)
式中Ωc为有效通带截止频率,
表示与通带波纹有关的参量,
值越大通带不动愈大。
VN(x)是N阶切比雪夫多项式,定义为
(1-5)
N为偶数,cos2(
)=1,得到min,
(1-6)N为奇数,cos2(
,得到max,
(1-7)
2用MATLAB设计切比低通雪夫滤波器
2.1MATLAB相关函数
(1)
该格式用于计算N阶切比雪夫I型归一化模拟低通滤波器系统的零极点和增益因子。
返回长度为N的列向量Z和P,分别给出N个零点和极点的位置。
Rs是阻带最小衰减。
(2)
该格式用于计算切比雪夫I型数字滤波器的阶数N和阻带截止频率wso。
调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。
(3)
该格式用于计算切比雪夫I型模拟滤波器的阶数N和阻带截止频率wso。
wp和ws是实际模拟角频率。
(4)
该格式用于计算N阶切比雪夫I型数字滤波器系统函数的分子和分母多项式系数向量B和A。
调用参数N和wso分别为切比雪夫I型数字滤波器的阶数和阻带截止频率的归一化值。
(5)
该格式用于计算N阶切比雪夫I型模拟滤波器系统函数的分子和分母多项式系数向量B和A。
调用参数N和wso分别为切比雪夫I型模拟滤波器的阶数和阻带截止频率的归一化值。
2.2数字高通滤波器的设计
对高通滤波器技术指标进行归一化处理,通常对切比雪夫滤波器以及其等波动通带截止频率为归一化因子
将高通的技术指标变换成低通的技术指标
根据转换出的低通滤波器的技术指标,按照前面介绍的低通滤波器设计方法设计满足技术指标的低通滤波器对设计出的归一化低通滤波器的系统函数进行得到归一化高通滤波器的系统函数
将设计的归一化高通滤波器进行反归一化处理,得到实际的高通滤波器
2.3数字带通滤波器的设计
运用双线性变换法基于MATLAB设计一个IIR带通滤波器
(1)确定性能指标
在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:
通带截止频率wp1,wp2;阻带截止频率ws1,ws2;阻带最小衰减As和通带最大衰减Ap;中心频率wp0。
(2)频率预畸变
用Ω=2/T*tan(w/2)对带通数字滤波器H(z)的数字边界频率预畸变,得到带通模拟滤波器H(s)的边界频率主要是通带截止频率Wp1,Wp2;阻带截止频率Ws1,Ws2的转换。
双线性变换法一般T=2s。
通带截止频率Wp1=(2/T)*tan(wp1/2)
Wp2=(2/T)*tan(wp2/2)
阻带截止频率Ws1=(2/T)*tan(ws1/2)
Ws2=(2/T)*tan(ws2/2)
(3)模拟带通性能指标转换成模拟低通性能指标
BW=Wp2-Wp1;%带通滤波器的通带宽度
W0=Wp1*Wp2;
WP=1;%归一化处理
WS=WP*(W0^2-Ws1^2)/(Ws1*BW);
(4)模拟低通滤波器的构造
借助切比雪夫(Chebyshev)滤波器得到模拟低通滤波器的传输函数Ha(s)。
(5)模拟低通滤波器转换成模拟高通滤波器
调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。
(6)模拟带通滤波器转换成数字带通滤波器
利用双线性变换法将模拟带通滤波器Ha(s)转换成数字带通滤波器H(z)。
(7)输入信号检验滤波器性能
输入不同频率的正弦波,观察输出波形,检验滤波器性能。
2.4双线性变换
目的:
将模拟带通滤波器转换成数字高通滤波器
为克服冲激响应法可能产生的频率响应的混叠失真,这是因为从S平面到Z平面是多值的映射关系所造成的。
为了克服这一缺点,可以采用非线性频率压缩方法,将整个频率轴上的频率范围压缩到-π/T~π/T之间,再用z=esT转换到Z平面上。
也就是说,第一步先将整个S平面压缩映射到S1平面的-π/T~π/T一条横带里;第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。
这样就使S平面与Z平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象,映射关系如图1-3
图1-3双线性变换的映射关系
为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上,可以通过以下的正切变换实现(1-13)
式中,T仍是采样间隔。
当Ω1由-π/T经过0变化到π/T时,Ω由-∞经过0变化到+∞,也即映射了整个jΩ轴。
将式(1-9)写成
(1-14)
将此关系解析延拓到整个S平面和S1平面,令jΩ=s,jΩ1=s1,则得
(1-15)
再将S1平面通过以下标准变换关系映射到Z平面z=es1T
从而得到S平面和Z平面的单值映射关系为:
(1-16)(1-17)
式(1-10)与式(1-11)是S平面与Z平面之间的单值映射关系,这种变换都是两个线性函数之比,因此称为双线性变换
式(1-9)与式(1-10)的双线性变换符合映射变换应满足的两点要求。
首先,把z=ejω,可得
(1-18)
即S平面的虚轴映射到Z平面的单位圆。
其次,将s=σ+jΩ代入式(1-12),得
因此有(1-19)
双线性变换法优缺点:
双线性变换法与脉冲响应不变法相比,其主要的优点是避免了频率响应的混叠现象。
这是因为S平面与Z平面是单值的一一对应关系。
S平面整个jΩ轴单值地对应于Z平面单位圆一周,即频率轴是单值变换关系。
这个关系如式(1-12)所示,重写如下:
(1-20)
上式表明,S平面上Ω与Z平面的ω成非线性的正切关系,如图1-4所示。
由图1-4看出,在零频率附近,模拟角频率Ω与数字频率ω之间的变换关系接近于线性关系;但当Ω进一步增加时,ω增长得越来越慢,最后当Ω→∞时,ω终止在折叠频率ω=π处,因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象,从而消除了频率混叠现象。
3三种切比雪夫数字滤波器的设计过程及程序
3.1切比雪夫型低通滤波器
3.1.1设计要求
设计一个切比雪夫型低通滤波器,指标如下:
通带边界频率:
wp=6Hz,通带最大衰减:
阻带截止频率:
,阻带最小衰减:
实验程序:
3.1.2低通滤波器仿真图:
3.1.3低通数字滤波器对连续信号的处理:
实验程序:
仿真图:
3.2切比雪夫型高通数字滤波器
3.2.1设计要求
设计一个高通切比雪夫型型数字滤波器,要求达到的指标是:
wp=30Hz,ws=28Hz,Fs=100Hz,rp=1db,rs=45db.
实验程序:
3.2.2高通滤波器仿真图:
3.2.3高通数字滤波器对连续信号的处理:
实验程序:
仿真图:
3.3带通切比雪夫数字滤波器
3.3.1设计要求
设计一个带通切比雪夫型型数字滤波器,要求达到的指标是:
wp=8Hz,ws=21Hz,Fs=100Hz,rp=1db,rs=45db.
实验程序:
3.3.2带通滤波器仿真图:
3.3.3带通数字滤波器对连续信号的处理:
实验程序:
仿真图:
4总结
这次课程设计过程,看了一遍《数字信号处理》课程中关于数字滤波器的设计的内容,再通过利用参考文献与网络,完成了用Matlab进行数字信号处理课程设计。
通过课程设计的内容,加深了对课堂学习中抽象概念——切比雪夫定理的理解,并且巩固了课堂上所学的理论知识,从而能很好地理解与掌握数字信号处理中的基本概念、基本原理、基本分析方法。
同时在使用MATLAB软件的过程中,掌握了matlab的编程方法和以及实际中解决问题的技巧。
与其他高级语言的程序软件相比,MATLAB环境下可以更方便、快捷地设计出具有严格线性相位的切比雪夫滤波器,节省了编程时间,提高了编程效率,并且相应参数的修改也十分方便,还可以进一步根据图像的变化特点进行优化设计。
相信随着版本的不断提高,MATLAB软件在数字滤波器等技术中必将发挥更大的作用。
同时,用MATLAB计算有关数字滤波器的设计参数,如H(z)、h(n)等,对我而言是在数字滤波器的硬件实现方面提供了一条简单而更加准确的途径和依据。
我还通过实验了解到,切比雪夫滤波器包括两种,即切比雪夫1型滤波器和2型滤波器,切比雪夫滤波器频率响应的特点是:
在带通内是等幅的纹波,在阻带内单调衰减,过度迅速。
更重要的是,我知道了,在实际工程中,需要设计高通、带通、和带阻滤波器时,通常将设计好的低通滤波器,如巴特沃斯低通滤波器或切比雪夫低通滤波器等,在传递函数H(s)中通过频率变化,转换成为其他类型的滤波器。
这些对我来说都是难能可贵的,并且从书本上难以获得的知识。
5设计结果及分析
5.1滤波器与其他类型低通滤波器优势及特点
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
如果需要快速衰减而允许通频带存在少许幅度波动,可用第一类切比雪夫滤波器;如果需要快速衰减而不允许通频带存在幅度波动,可用第二类切比雪夫滤波器。
切比雪夫滤波器的幅度特性具有等纹波特性,I型切比雪夫滤波器:
通带内是等纹波的,在阻带内是单调的;II型切比雪夫滤波器:
通带内是变化的,阻带内等纹波变化的。
切比雪夫滤波器频率响应的特点是:
在带通内是等幅的纹波,在阻带内单调衰减,过度迅速。
在实际工程中,需要设计高通、带通、和带阻滤波器时,通常将设计好的低通滤波器,如巴特沃斯低通滤波器或切比雪夫低通滤波器等,在传递函数H(s)中通过频率变化,转换成为其他类型的滤波器。
而且切比雪夫滤波器的截频衰减陡度比同阶数巴特沃斯特性更陡度比同阶数程时的衰减就超过6NdB.在阶数一定时,波纹越大,截频衰减陡度越陡。
相位响应也是非线性,但较之巴特沃斯为差。
5.2心得体会
在课程设计刚刚开始的时候,觉得很无助,于是请教同学、上网搜资料、去图书馆查找等等,终于功夫不负有心人,在不断的资料搜寻当中我渐渐的明白了我们将要做的设计的基本原理和相关知识。
本次课程设计,我选题为切比雪夫数字滤波器的设计,通过设计,我具体的了解了
(无限脉冲响应)数字滤波器设计的原理和其相关的知识。
首先,我们要了解设计数字滤波器的原理方法,从模拟滤波器设计
数字滤波器在工程上常用的有两种:
脉冲响应不变法、双线性变换法。
其设计过程都是由模拟滤波器的系统函数
去变换出相应的数字滤波器的系统函数
。
然后,我们要知道需要用函数来编写程序通过MATLAB来实现滤波器的设计。
我们知道的现有的逼近函数如巴特沃斯、切比雪夫可供使用。
切比雪夫滤波器是使通带内误差分布均匀的滤波器。
由此函数设计出的模拟低通滤波器通过频率变换
可得到高通、带通、带阻模拟滤波器。
最后通过所掌握的
(无限脉冲响应)数字高通滤波器设计原理和方法将其应用到实际问题当中进行处理。
通过问题分析了离散经典数字滤波器分类和设计方法。
这次课程设计让我收获很多。
6参考文献
[1]谢平,王娜,林洪彬编著信号处理原理及应用——北京:
机械工业出版社,2008.10
[2]郭仕剑,王宝顺,贺志国《MATLAB数字信号处理》.人民邮电出版社,2007.11