椭圆IIR高通数字滤波器

南华大学电气工程学院

数字信号处理课程设计

设计题目：椭圆高通HR数字滤波器

专业：10级电子02班

学生姓名：号

起迄日期:2013年12月期日一2014年1月11日

指导教师：

一、总体设计思路

(1)根据rand函数产生的0-1之间随机值计算数字高通滤波器的技术

指标。

(2)将数字高通滤波器的边界频率转换成模拟高通滤波器的边界频率。

(3)将模拟高通滤波器的的技术指标转换成归一化的椭圆模拟低通滤波

器的技术指标。

(4)设计归一化椭圆低通滤波器原型Go(s)。

(5)通过频率转换将GLP(S)转换成过渡的模拟高通滤波器G"(s)。

(6)采用双线性变换法将过渡的G/(s)转换成数字高通滤波器

HHP(Z)。

二、手工计算完成椭圆HR数字高通滤波器的初始设计

1.数字高通滤波器的设计指标

(1)通带截止频率o)pc=321Trad,

(2)阻带截止频率o)sc=31nrad,

(3)通带最大衰减内)=ldB,

(4)阻带最小衰减%=60dB

其中，曲，M为两个用MATLAB的方〃</函数随机产生的0-1之间的值，且

有如<to2。如图LL1所示侬1=0.2785,a)2=0.5489

(1)通带截止频率：

o)pc=6)2nrad=0.5469iTrad=l.7181

(2)阻带截止频率：

u)sc=rad=0.278511rad=0.8749

(3)通带最大衰减：与=ldB

(4)阻带最小衰减：as=60dB

2.将数字高通滤波器设计指标转换成模拟高通滤波器设计指标

91

频率转换关系为：n=^tani(n,为了计算方便，令T=2s计算得：

JL4

(1)Cph=1.1594rad/s

(2)%h=0.4677rad/s

(3)%=ldB

(4)CLS=60dB

3.转换成归一化椭圆滤波器设计指标

(1)归一化边界频率

椭圆滤波器的归一化低通原型一般选择关于通带边界频率Cph归一化

的低通系数。

在虚轴上低通到高通滤波器的映射关系为如下频率变换公式

4pCph

A=--------

综合上述，得:

归一化通带边界频率4P=1

归一化阻带边界频率As=—=2.4789

(2)确定椭圆滤波器最小阶数

定义频率的选择性因数k为:

Ap

fc==2.4789

令:

i

_1l-(l-k2)4_

q。—-------1"

l+(l-k2)4

5913

q=q0+2q0+15q0+15q0=

b2=10°laP=106

£2=io°las=io01

椭圆滤波器的最小阶次N可表示为:

经计算，最小阶数为N=4。

4.设计归一化椭圆低通滤波器原型G'p(s)

椭圆滤波器的幅频响应函数公式为：

1

|G(jC)/=----------------

1+5E/I)

其中，£是与通带衰减有关的参数，Cs为截止频率，EN(生)为雅可比椭

%

圆函数。这里:

6.模拟高通滤波器转换为数字高通滤波器

从Z平面上的一点映射到S平面一点双线性变换为

1-Z-1

1+Z-1;

因此,数字传输函数H“p(Z)和原模拟函数％p(s)之间的关系为

HHP(Z)—G“p(s)|2/i-z~1\

S=T[

_0.0514-0.1654z-1+0.2325z-2-0.1654z-3+0.0514z-4

⑶=1+1.1174Z-1+1.3475z-2+0.7501z-3+0.2673z-4

三、检验

1.检验上述变量计算结果

为了检验上述变量是否计算正确，我编写了从第一步到第三步所有变量计

算的M文件脚本。M文件和计算结果如图1.1和图1.2所示：

图L1M文件脚本

图2.2运行结果

结论：经检验，手工计算与计算机计算得出的结果完全相符。

2.检验最终得出的数字高通滤波器传递函数

根据数字信号处理教材上关于椭圆滤波器的一些介绍，我编写

了直接设计椭圆数字滤波器的M文件，源程序和运行结果如图2.3和图2.4所

图2.3

图2.4

结论；通过matlab函数直接设计椭圆数字高通滤波器检验，结果与上述手

工计算相符。

四、滤波器的不同结构对性能指标的影响

1、利用直接型结构构建数字滤波器

在MATLAB中口J以利用FDATOOL工具箱构建不同类型的数字滤波器。在此

为了使比照效果明显，我们不妨先将将上述初步设计的椭圆数字HR带阻滤波

器的设计参数的字长（即转移函数中分子、分母各项前的系数）进行保存小数

点后3位的进一步的缩减。缩减后的参数如下：

den=[0.051,-0.165,0.235,-0.165,0.051]

num=[1.000,1,117,1,347,0.750,0.267]

图4.1.1filtercoefficients工具工作界面

直接型的结构流图如下图：

y(n)

点击窗口下方

的ImportFilter按钮，构建直接1型结构的椭圆高通滤波器，结果如下图:

图4.1.3Direct-FormI型结构的漉波器幅频响应图

读图可以得Direct-FormT结构的滤波器技术指标（ws,wp单位为

mad/sample;as,av,单位为dB）如表1所不：

表Direct-FormI结构滤波器对性能指标的影响

性能指标初始设计指标Direct-FormI△

3s0.27850.28002930.0015293

3P0.54690.5319824-0.0149176

&6059.69624-0.30376

%11.0206320.020632

分析：由图和表4.1.1可以看出，ws上升0.0015293mad/sample,wp下

降了0.0149176^rad/sample,as下降0.30376db,ap上升0.020632db。阻带的幅

频响应曲线更加平滑，Direct-FormI造成性能指标的误差很大，不能忽略。

2、利用级联结构构建数字滤波器

级联型的结构流图如图4.2.1所示:

x(n)y(n)

agj%

。2匹

图4.2.级联型的结构流图

利用matlab信号工具箱提供的tf2sos函数，可以将直接型结构转换到级

联型结构，具体操作如图4.2.2和图4.2.3所示：

图4.2.2M文件脚本

图4.2.3运行结果

在importassecond-ordersections选项前打勾,在SOS栏和Gain栏

填入上面生成的S和Go

filtercoefficients工具工作界面

选择Edit下拉菜单中点击ConverttoSecond-orderSections选项，将

构建好的Direct-FormI结构的椭圆数字IIR带阻滤波器转换为级联滤波器，

结果如图4.2.5所示。

图级联型结构的滤波器幅频响应图

读图4.2.5可以得级联结构的滤波器技术指标（fs,fp单位为wad/sample；

单位为dB）如表2所示：

表级联结构滤波器对性能指标的影响

性能指标初始设计指标级联结构△

3s0.27850.28088380.00238838

O)p0.54690.5319824-0.0149176

&6060.005740.00574

%10.9999995-0.0000005

分析：由图4.2.5和表可以看出，ws上升0.00238838m-ad/sample,wp下

降了0.0149176/arad/sample,as上升0.00574B,ap下降0.0000005dB。与上面相

比误差减小，级联结构造成性能指标误差较Direct-FomiI更小，误差可以忽

略。

3、两种结构滤波器对指标影响比拟与原因分析

比拟表4.1.1和表4.2.1发现：在参数字长仅保存了小数点后3位的情况

下，两种结构的滤波器较初始设计在性能指标方面均有误差。但是直接型误差

比级联型更大，受有限参数字长影响更大，主要表现在：直接型ws,wp,a和

%与设计要求相应的性能间的差的绝对值普遍大于级联型。此外，直接型和级

联型的幅频响应曲线的通带的波动均不稳定，但是级联型的稳定性要好于直接

型。

造成这一现象的原因是：直接型滤波器的系数不是直接决定单个零极点，

因而不能很好的进行滤波器性能的控制；此外直接型滤波器的极点对参数的变

化过于敏感，从而使得系统的频率响应对参数的变化也特别敏感，也就是对参

数的有限字长运算过于灵敏，容易出现不稳定或产生较大误差。

相比之下，级联型滤波器每个二阶节系数单独控制一对零、极点，有利于

控制频率响应；此外级联结构中后面的网络输出不会再流到前面，运算误差的

积累相对直接型就小。

三、参数字长对性能指标的影响

在实际的数字滤波器的设计中，由于计算机或DSP芯片等的字长和存储空

间有限，所以也只能对设计参数取有限的字长进行设计。然而，如果字长太短，

那么设计的滤波器误差就会太大，造成滤波效果不佳。卜面就将以上述性能指

标为依据，初始设计的椭圆数字HR带阻滤波器为例，研究不同参数字长对性

能指标的影响。并为适宜参数字长确实定探索规律。

运用FDATOOL工具左下侧上数第三个的Setquantizationparameters

按钮，在filterarithmetic下拉菜单下选择FixedPoint选项，进入工作

界面。

通过改变coefficientwordlength的值便可以改变参与构建滤波器的参

数字长。

1.参数字长取2位对性能指标的影响

将coefficientwordlength的值改为2,点击下方的Apply按钮,此时

设计的滤波器幅频响应由线和性能指标如下图：

图参数字长取2位时的滤波器幅频响应曲线图

图中的虚线为供参考的理想字长下生成的滤波器的幅频响应曲线，图中实

线为参数字长取为2位时的滤波器幅频响应曲线。从图中可以看出：字长为2

位时，滤波器的各项性能指标离设计指标偏差很大，滤波器失真明显，滤波效

果很差，远远不能满足设计指标的要求。

2.参数字长取4位对性能指标的影响

将coefficientwordlength的值改为4,点击下方的Apply按钮，此时

设计的滤波器幅频响应曲线和性能指标如图8和表3所示；

图参数字长取4位时的滤波器幅频响应曲线图

表3参数字长取4位时的实际性能指标一览表

性能指标初始设计指标级联结构△

3s0.27850.30163570.0231357

0)p0.54690.62048340.073584

&60**

%1**

由图8和表3可以看出，当参数字长取为4位时，幅频曲线失真度较2位

时有明显改良，但仍很明显：通、阻带截止频率分别有较大的偏移，使得截止

频率坡度非常平缓，滤波效果很差。通带和阻带波动平缓，最小和最大衰减频

率不明显，通带甚至出现了正增益。离设计指标差距仍然很大。

3.参数字长取6位对性能指标的影响

将coefficientwordlength的值改为6,点击下方的Apply按钮,此时

设计的滤波器幅频响应曲线和性能指标如图9和表4所示：

图4.3.1参数字长取6位时的滤波器幅频响应曲线图

表4.3.1参数字长取6位时的实际性能指标一览表

性能指标初始设计指标级联结构△

3s0.27850.27941890.0009189

3P0.54690.53918450.0077155

%6060.557650.55765

11.9214430.921443

由图和表4.3.1可以看出，当参数字长取为6位时，幅频曲线失真进一步

减小，已经初具带阻滤波器的形制。性能指标也明显接近设计要求：通、阻带

截止频率与设计要求间的误差已经缩小至±0.02的范围内，截止频率的坡度也

已经十分陡峭；但是通带最大衰减频率和阻带最小衰减频率与设计指间标误差

还是很大，离设计要求距离还是较远。

4.参数字长取8位对性能指标的影响

将coefficientwordlength的值改为8,点击下方的Apply按钮,此时

设计的滤波器幅频响应由线和性能指标如图10和表5所示:

图4.4.1参数字长取8位时的滤波器幅频响应曲线图

表参数字长取8位时的实际性能指标一览表

性能指标初始设计指标级联结构△

0.27850.28027340.0017734

0.54690.5440674-0.0028326

%6060.016070.016607

%11.3030120.303012

由图和表可以看出，当参数字长取为8位时，幅频曲线失真进一步减小,

但仍可以看出。性能指标与设计要求的差值继续减小，也已不明显；通带最大

衰减频率和阻带最小衰减频率与设计指间标误差还是有一定的差距，距离设计

要求仍可以进一步接近。

5.参数字长取10位对性能指标的影响

将coefficientwordlength的值改为10,点击下方的Apply按钮,此

时设计的滤波器幅频响应曲线和性能指标如图11和表6所示：

图参数字长取10位时的滤波器幅频响应曲线图

表参数保存小数点后10位时的性能指标一览表

性能指标初始设计指标级联结构