切比雪夫I型数字带通滤波器
课 程 设 计
课程名称: 数字信号处理 题目编号: 0801 题目名称:切比雪夫I型带通IIR数字滤波器 专业名称: 电子信息工程 班 级: 1203班 学 号: [1**********] 学生姓名: 段超 任课教师: 陈 忠 泽
2015年08月30日
目 录
1. 数字滤波器的设计任务及要求(编号0801) ................... 2 2. 数字滤波器的设计及仿真 ................................... 2
2.1数字滤波器(编号0801)的设计 ................................... 2 2.2数字滤波器(编号0801)的性能分析 ............................... 4
3. 数字滤波器的实现结构对其性能影响的分析 ................... 6
3.1数字滤波器的实现结构一(0801)及其幅频响应...................... 7 3.2数字滤波器的实现结构二(0801)及其幅频响应...................... 9 3.3 数字滤波器的实现结构对其性能影响的小结 ........................ 12
4. 数字滤波器的参数字长对其性能影响的分析 .................. 12
4.1数字滤波器的实现结构一(0801)参数字长及幅频响应特性变化...... 14 4.2数字滤波器的实现结构二(0801)参数字长及幅频响应特性变化....... 17 4.3 数字滤波器的参数字长对其性能影响的小结 ........................ 20
5. 结论及体会 .............................................. 20
5.1 滤波器设计、分析结论 .......................................... 21 5.2 我的体会 ...................................................... 21 5.3 展望 .......................................................... 21
1. 数字滤波器的设计任务及要求
(1)切比雪夫I型带通IIR数字滤波器各项指标如下: 阻带下截止频率s1=32
ln(d)
rad
;
(id)rad;
20
(id)rad;
20
通带下截止频率p1=ln(d)log10
32
通带上截止频率p2=1-ln(d)log10
32
ln(d)1-rad
32 阻带上截止频率s2=; 通带最大衰减p=1dB; 阻带最小衰减s=60dB.
其中的id为抽到的题目的四位数编号,我的题目编号是0801,所以取 id为801.
(2)滤波器的初始设计通过手工计算完成;
(3)在计算机辅助计算基础上分析滤波器结构对其性能指标的影响(至少两种 合适的滤波器结构进行分析);
(4)在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响;
(5)以上各项要有理论分析和推导、原程序以及表示计算结果的图表;
(6)课程设计结束时提交说明书。
2. 数字滤波器的设计及仿真 2.1数字滤波器的设计
(1)根据任务,确定性能指标:在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:(手工计算完成)
通带下截止频率p1=0.354π 通带上截止频率p2=0.646π 阻带下截止频率s1=0.208π 阻带上截止频率s2=0.792π 阻带最小衰减 s=60dB 通带最大衰减 p=1dB 抽样频率Fs=8kHz
(2)有所给的技术参数可直接得到各模拟频率。
p1=p1*Fs,p2=p2*Fs, s1=s1*Fs, s2=s2*Fs p1=2832, p2=5168, s1=1664, s2=6336
(3) 根据教程表7.8的变换关系,设归一化(p)低通滤波器的阻带衰减用
_
表示,则有 Bp=p2-p1=2336, p0=
p1p2=3825.6732
st1p0
st1Bp
2
2
st2=
_
st2p0
st2Bp
22
=1.723 st1=
=-3.053
取 st=min(| st1|,| st2st2|)=1.723,可满足阻带衰减要求。
(4)求波纹系数
=
_
0.1p
1=0.50885
(5)根据st及给定的p,s,利用教程(7.5.65)式,可得切比雪夫I型低通滤波器的阶次为
0.2s1arcch(1965.2147)8.2765arcch =arcch(1.723)=1.1398=7.26114
arcchst/p
取N=8.
(6)查教程表7.5,可得N=8,p=1dB时的归一化原型切比雪夫I型低通滤波器的系统函数Han(s)
Han(s)=
s80.9198s72.4230s61.6552s51.8369s
2
4
0.8468s30.4478s20.1073s10.01723
(7)按教表7.8的相应变换关系,求出模拟带通滤波器系统函数为:
(8)按冲激响应不变法的修正公式(7.78)式,求所需数字滤波器的系统
Hbp(s)Han(s)
s
s2p02sBp
4144
8.010s7.172210 1.01061s1.89551065
函数
H(z)
0.0001z30.0003z40.0015z50.0015z60.0013z7
10.5316z15.7911z22.6014z315.7331z45.9831z526.0350z68.2727z728.5847z8
0.0038z80.0024z90.0005z100.0012z110.0004z12
7.3787z921.2727z104.2324z1110.4690z121.4457z133.1159z140.2275z150.4301z16
2.2 数字滤波器的性能分析
MATLAB 程序如下:
>> Rs=1;As=60;Fs=8000;T=1/Fs;
wp1=0.354*pi; wp2=0.646*pi; ws1=0.208*pi; ws2=0.792*pi; >>
Omegap1=wp1*Fs;Omegap2=wp2*Fs;Omegap=[Omegap1,Omegap2];Omegas1=ws1*Fs;Omegas2=ws2*Fs;Omegas=[Omegas1,Omegas2]; >> bw=Omegap2-Omegap1;w0=sqrt(Omegap1*Omegap2); >> [N,OmegaC]=cheb1ord(Omegap1,Omegas1,Rp,As,'s') N = 8 OmegaC=
8.8970e+03 >>
[z0,p0,k0]=cheb1ap(N,Rp);ba1=k0*real(poly(z0));aa1=real(poly(p0)); >> [ba,aa]=lp2bp(ba1,aa1,w0,bw); [bd,ad]=impinvar(ba,aa,Fs); %(C,B,A)=dir2par(bd,ad)
[H,w]=freqz(bd,ad);dbH=20*log10((abs(H)+eps)/max(abs(H))); >> subplot(2,2,1),plot(w/pi,abs(H)); >> subplot(2,2,2),plot(w/pi,angle(H)/pi); >> subplot(2,2,3),plot(w/pi,dbH); >> subplot(2,2,4),zplane(ba,ad); 相应的图形
图1频率特性曲线和零极点图
3. 数字滤波器的实现结构对其性能影响的分析
在理想状态下,对于同一个传递函数几乎对应着无数种等效结构,然而这些 结构却并不一定都能实现。在无限参数字长的情况下,所有能实现传递函数的结构之间,其表现完全相同。
然而在实际中,由于参数字长有限的限制,各实现结构的表现不同。下面我们就对比直接型和级联型两种结构在本设计中对性能指标的影响。
在MATLAB中可以利用FDATOOL工具箱建构不同类型的数字滤波器。再次为了使对比效果明显,将上述初步设计的切比雪夫带通IRR数字滤波器的设计参数字长(即转移函数中分子分母各项前的系数)进行保留小数点后3位的进一步缩减。缩减参数如下:
Nun=[0 0 0 0 0 0.002 -0.002 -0.001 0.004 -0.002 -0.001 0.001 0 0 0 0] Den=[1 -0.532 5.791 -2.601 15.733 -5.983 26.035 -8.273 28.585 -7.379
21.273 -4.232 10.469 -1.446 3.116 -0.228 0.430] 将上述参数输入FDATOOL中的filter coefficients 工具中
3.1数字滤波器的实现结构一及其幅频响应
直接型结构滤波器的转移函数为:
H(z)
0.002z50.002z60.001z7
10.532z15.791z22.601z315.733z45.983z526.035z68.273z728.585z8
0.004z80.002z90.001z100.001z11
7.379z921.273z104.232z1110.469z121.446z133.116z140.228z150.430z16
差分方程:
y(n)0.002x(n5)0.002x(n6)0.001x(n7)0.004x(n8)0.002x(n9)0.001x(n10)
0.001x(n11)0.532y(n1)5.791y(n2)2.601y(n3)15.733y(n4)5.983y(n5)
26.035y(n6)8.273y(n7)28.585y(n8)7.379y(n9)21.273y(n10)4.232y(n11)10.469y(n12)1.446y(n13)3.116y(n14)0.228y(n15)0.430y(n16)
直接一型的结构流图如下:
y(n)
图2 Direct-Form I滤波器结构流程图
选择filter structure 选项框中的Direct-Form I选项,点击窗口下方的Import Filter 按钮,构建直接1型结构的椭圆带通IRR数字滤波器,结果如图所示。
图 3 Direct-Form I型结构滤波器幅频响应图
表1 Direct-Form I 结构滤波器性能指标影响(p1,p2,s1,s2单
位为rad/sample;s,p,单位为dB)
由图3和表1可以看出:
(1)滤波器幅频曲线在通带和阻带内波动幅度不均匀; (2) 阻带最小衰减与通带最大衰减都存在误差; (3) wp和ws分别较初始设计变化了。
3.2数字滤波器的实现结构一及其幅频响应 级联型结构滤波器的转移函数:
Z2Z2Z11.1283Z2
H(
z)0.002()()()
10.0328Z10.8103Z210.5222Z10.8638Z210.6715Z10.8759Z210.6304Z111()()()12121210.4108Z0.8809Z10.3369Z0.8924Z10.8774Z0.9608Z11.4550Z10.6328Z211.3048Z11.1108Z2()10.7747Z10.9629Z210.8642Z10.9643Z2
图
4 级联型的结构流图
图
5级联型结构滤波器幅频响应图
表2 级联型结构滤波器性能指标影响(p1,p2,s1,s2单位为
rad/sample;s,p,单位为dB)
从上图和表2中可以看出,采用级联结构的滤波器 (1)滤波器幅频曲线在通带和阻带内波动较均匀;
(2)出现了阻带最小衰减上升了0.45133dB,通带最大衰减下降了0.1324799dB的
误差;
(3)wp1和ws2分别较初始设计减小-0.0026818和-0.0289624rad/sample ,
(4)wp2和ws1分别较初始设计增大了0.003292和0.0897295rad/sample
3.3 数字滤波器的实现结构对其性能影响的小结
比较表1和表2发现:在参数字长仅保留了小数点后三位的情况下,两种
结构的滤波器较初始设计在性能指标方面均有误差。但是直接型误差比级联型更大,受有限参数字长影响更大,其各项性能指标与设计要求间的差的绝对值普遍大于级联型。此外,直接型和级联型的幅频响应曲线的通带波动均不稳定,但是级联型的稳定性要好于直接型。
造成这一现象的原因是:直接型滤波器的系数不是直接决定单个零极点,因而不能很好的进行滤波器性能的控制;此外直接性滤波器对参数的变化过于敏感,从而使得系统的响应对参数的变化也特别敏感,对参数的有限字长运算过于灵敏,容易出现不稳定或产生较大误差。
而级联型滤波器每个二阶节系数单独控制一对零极点,有利于控制频率响应;此外级联结构中后面的网络输出不会再流到前面,运算误差的积累相对直接型就
4.数字滤波器参数字长对其性能指标的影响的分析
在实际的数字滤波器的设计中,由于计算机或DSP芯片字长和存储空间有限,
所以也只能对设计参数取有限字长进行设计。然而,如果字长太短,则设计的滤波器误差就会太大,造成的滤波效果不佳。下面就将以上述性能指标为依据,初始设计的切比雪夫I型带通IRR数字滤波器为例,研究不同参数字长对性能指标的影响。并为合适的字长确定探索规律。
将计算获得的滤波器系数输入filter coefficients 工具中,并点击Import Filter 按钮,生成数字滤波器。运用FDATOOL工具中Set
quantizationparameters 按钮,在filter arithmetic下拉菜单下选择Fixed
Point选项,进入图 所示面。 通过改变 coefficient word length 的值便可以改变参与构建滤波器的参数字长
图6 Set quantization parameters gong工作界面
4.1 直接型结构滤波器参数字长及其幅频响应变化 1. 直接型结构滤波器(2
位)参数字长及其扶贫响应变化
图 7参数字长取2位时的滤波器幅频响应曲线(直接)
由图9可以看到当参数字长取2位时,直接型滤波器严重失真,误差相当大,远远不能满足设计要求。
2. 直接型结构滤波器(15位)参数字长及其幅频响应变化
图 8 参数字长取15位时的滤波器幅频响应曲线(直接)
表3参数字长取15位时的直接型结构滤波器实际性能指标一览表(p1,p2,s1,s2单位为rad/sample;s,p,单位为dB)
由图10和表3可以看出当参数字长取15位时,幅频曲线失真明显减小,已具有带通滤波器形制。性能指标也接近设计要求, 但通带最大衰减和阻带最小衰减频率与指标误差还较大。
3. 直接型结构滤波器(20位)参数字长及其幅频响应变化
图 9 参数字长取22位时的滤波器幅频响应曲线(直接)
表4:参数字长取20位时的直接型结构滤波器实际性能指标一览表(p1,p2,s1,s2单位为rad/sample;s,p,单位为dB)
由表4和图9可知,当参数字长取22时,失真程度已经非常小,通带最大衰减和阻带最小衰减和设计要求几乎无差距,基本上满足设计要求。
因直接型结构滤波器根据参数字长的变化跨度很大,未能详细的对更多字长 位数进行分析,下面对级联型结构滤波器进行详细分析。
4.2级联型滤波器参数字长及幅频响应变化
1. 级联型滤波器(2位)参数字长及幅频响应变化
图 10参数字长取2位时的滤波器幅频响应曲线图(级联型)
图10中的虚线为提供参考的理想字长生成的滤波器的幅频响应曲线,实线为参数字长为2位时的滤波器幅频响应曲线。由图可知:字长为2位时,滤波器的各项性能指标偏差很大,滤波器失真明显,滤波效果很差,远远不能满足设计指标的要求。
图 11参数字长取5位时的滤波器幅频响应曲线图(级联型)
表5:参数字长取5位时的级联型结构滤波器实际性能指标一览表 (p1,p2,s1,s2单位为rad/sample;s,p,单位为dB)
由图11和表5可看出,幅频曲线失真变得更小,程度变小,通带最大衰减和阻带最小衰减仍不可忽略。
图12 参数字长取12位时的滤波器幅频响应曲线图(级联型)
表6:参数字长取12位时的级联型结构滤波器实际性能指标一览表 (p1,p2,s1,s2单位为rad/sample;s,p,单位为dB)
由图12和表5可以看出,当参数字长取为12时,幅频曲线失真进一步减小,但程度已经很小,失真几乎可以忽略。截止频率已与设计要求相差无几,改善的程度也变得很小;通带最大衰减和阻带最小衰减和设计要求几乎无差距。
7.级联型结构滤波器(12位)参数字长及幅频响应变化
图13 参数字长取12位时的滤波器幅频响应曲线图
参数为12位及以上时,幅频曲线失真可看作为零,设计的曲线与要求的曲线几乎完全重合。各项指标几乎完全相同,达到设计要求。
4.3 数字滤波器的参数字长对其性能影响的分析
参数的字长对数字滤波器的性能有影响,字长越短,影响越大,误差也越大,会导致滤波器失真,直接型滤波器字长为22位以上时达到设计要求,级联型滤波器字长为12位以上时达到设计要求,故设计时要保有合适的字长。
5.结论及体会
5.1滤波器设计、分析结论
本次数字滤波器课程设计过程中,发现滤波器的级联型结构比直接型结构误差小,幅频特性曲线的通带的波动稳定性要比直接型好。对于参数字长来说,字长越长,滤波器的实际指标就越能达到设计指标的要求,但不是越长越好
5.2我的体会
采用MATLAB设计滤波器,使原来繁琐复杂的程序设计变成了简单的函数调用,为滤波器的设计和实现开辟了广阔的天地,尤其是MATLAB工具箱使各个领域的研究人员可以直观方便地进行科学研究于工程应用。其中的信号处理工具箱,图像处理工具箱,小波工具箱等更是为数字滤波器研究的蓬勃发展提供了可能。MATLAB信号处理工具箱为滤波器设计和分析提供了非常优秀的辅助设计工具,在设计数字滤波器时,善于应用MATLAB进行辅助设计,能够大大提高设计效率。
5.3展望
通过这次切比雪夫i型设计带通IIR数字滤波器,对冲激响应法这一方法有了较为深入的认知,但其他方法还是没能有进一步的理解与掌握,希望以后在数字滤波器的设计方面,去更多的尝试多种方法,有一个全面的提高。