当前位置:文档之家› 设计IIR滤波器的双线性变换法

设计IIR滤波器的双线性变换法


版权所有 违者必究
第五章 第2讲
16
典型例题
② 频率预畸变
c
tg c
2
tg(0.1 ) 0.3249(rad/s)
∴ 数字高通滤波器的系统函数为:
H (z)
HaL (
p)
p
c
1 1
z z
1 1
1
p3
2 p2
2p
1
p
0.32491 1
z z
1 1
(1 z 1 )3
1.859208 3.3358z 1 2.241976z 2 0.57016z 3
模拟域频率变换法
1、归一化模拟低通原型到数字高通滤波器的频率变换
设归一化模拟低通原型滤波器的系统函数为 H aL ( p),
p为模拟域内的拉氏变量 。
在模拟域内从低通到高通变换:即以 p-1代替 p。
HaH ( p) HaL (1 / p)
反后归的一高化通。滤即 波以器的p传输s 函数c 带H(入s)上。式,得到反归一化
均匀的频率点映射到上时变成了非均匀的点,而且随
频率增加越来越密。
双线性变换法除了不能用于线性相位滤波器设计外,仍
然是应用最为广泛的设计IIR数字滤波器的方法。
版权所有 违者必究
第五章 第2讲
7
设计IIR滤波器的双线性变换法
频率预畸变
为了保证各边界频率点为预先指定的频率,在确定模拟
低通滤波器系统函数之前必须按式 c
版权所有 违者必究
第五章 第2讲
17
典型例题
例2:设计一个巴特沃思型数字高通滤波器,3dB数字截
频为c 0.2π弧度,阻带下边频 =s 0.05π弧度,阻带衰
减 As≥48dB,求数字滤波器的系统函数。
解: ① 频率预c 畸t变g 求2c模拟tg低(0通.1的 )截频0.3和2c4阻9(带ra下d/边s) 频 s
H
(s)
H
aH
(
s c
)
H
aL
(
c s
)
双线性变换,得数字高通滤波器的系统函数H(z)。
H (z) H (s) s11zz11
H ( p) aL
p
c
1 1
z z
1 1
版权所有 违者必究
第五章 第2讲
15
设计IIR数字滤波器的频率变换法
变换关系:直接由归一化低通原型变换成数字高通滤波器时:
p
c
1 1
t g(c )
2
进行
频率预畸变;然后将预畸变后的频率代入归一化低通原
型Ha(s) 确定 H(s) Ha (s / c ) ;最后求得数字系统
函数:H (z)
H (s) s
1 1
z z
1 1
Ha (s)
s
1 c
1 1
z z
1 1
频率预畸变示意图
用双线性变换法设计IIR滤波器的过程总结
H (e j ) 频率预畸变 H ( j) H (s) 双线性变换 H (z)
(极点)
u , l :下边频和上边频
sl:下阻带上边频
su :上阻带下边频
版权所有 违者必究
第五章 第2讲
21
设计IIR数字滤波器的频率变换法
满足以上映射关系的变换式推导过程:
s
s(z)
(z
e j0 )( z e j0 ) (z 1)(z 1)
z2
2z cos0
z2 1
1是稳定的
∴ S平面稳定的函数变换到Z平面也是稳定的。
版权所有 违者必究
第五章 第2讲
8
设计IIR滤波器的双线性变换法
返回
版权所有 违者必究
第五章 第2讲
9
设计IIR滤波器的双线性变换法
有关双线性变换公式的说明
有些文献中双线性变换的关系为:
s
2 T
1 1
z 1 z 1
(先前为:s
1 1
z 1 z 1
)
两式存在 2 /T的系数差别
抽样间隔T的选取:
版权所有 违者必究
第五章 第2讲
6
设计IIR滤波器的双线性变换法
双线性变换的频率对应关系
双线性变换法虽然避免了“频率混叠效应”,但出现了
模拟频率与数字频率为一种非线性的关系情形。即:
tg(
)
2
可见:模拟滤波器与数字滤波器的响应在对应的频率关
系上发生了“畸变”,也造成了相位的非线性变化,这
是双线性变换法的主要缺点。具体而言,在上刻度为
满足: s
的任意T值。
T
只要保证将模拟频率在带限之间,不会产生因多值映
射产生频率混迭现象即可。
版权所有 违者必究
第五章 第2讲
10
典型例题
例:利用双线性变换法设计巴特沃斯型数字低通滤波器
设计参数:通带数字截止频率 p 0.25 ,通带内最大 衰减 Ap 0.5dB,阻带数字截止率 s 0.55 ,阻带内
s
c2ctg
s
2
1.341267(rad/s)
② 确定低通滤波器的阶次 N
N 1 lg(100.1As 1) 3.8976 2 lg( s / c )
∴ 取N=4
版权所有 违者必究
第五章 第2讲
18
典型例题
③ 查表得N=4时巴特沃思低通原型滤波器的系统函数
1 HaL ( p) ( p2 0.7654 p 1) ( p2 1.8478p 1)
2
j
面的单位圆
版权所有 违者必究
第五章 第2讲
5
设计IIR滤波器的双线性变换法
1 s
令s
j
, 带入 z
1
得:
s
z 1 j 1 j
(1 )2 2
z
(1 )2 2
当 0 时, z 1
此式表明: S平面的左 半平面轴映 射到了Z平 面的单位圆 内,可保证 系统函数经 映射后稳定 性不变。
版权所有 违者必究
第五章 第2讲
3
设计IIR滤波器的双线性变换法
3、双线性变换法的基本思路:
从频率响应出发,直接使数字滤波器的频率响应 H (e j ) ,
逼近模拟滤波器的频率响应 H ( j) ,进而求得H(z)。
“双线性变换法”设计方法
① 通过正切变换: tg( 1T ) 2
将S平面的jΩ轴压缩到S1平面的jΩ1轴上的
特别地,设 0 < z = r < 1,则有:
s
r
2
2r cos0
r2 1
1
(1
r)2
2r(1 r2 1
cos 0
)
0
此式表明:Z平面单位圆内的极点变到了S平面的左半平面
再令 z e代j入 上式 ,得:
s
e j2
2e j e j2
cos0
1
1
j
cos0 cos sin
j
此式表明:Z平面单位圆变换到了S平面的虚轴
cos(1 2 )
cos 0
cos(1
2
2
)
0
2
0 0.5
版权所有 违者必究
第五章 第2讲
24
典型例题
③ 求模拟低通截止频率和阻带下边频
c
1
cos 0 cos 1 sin 1
cos0.5 cos0.75 sin 0.75
1
s
cos0 cosu sin u
cos(0.5 ) cos(0.925 ) sin( 0.925 )
版权所有 违者必究
第五章 第2讲
20
设计IIR数字滤波器的频率变换法
2、归一化模拟低通原型到数字带通滤波器的频率变换 方法:直接寻求从模拟低通到数字带通之间的映射关系
低通幅度响应
带通幅度响应
0
sb:模拟的阻带下边频
0:中心频率
映 射 关 系
0
0
(零点)
z
e j0
s
0
0, z 1 s j
④ 利用直接变换关系求数字滤波器的系统函数H(z)
H (z)
HaL (
p)
p
c
1 1
z z
1 1
(1 z 1 )4
2.3103 5.4743z 1 5.346z 2 2.4366z 3 0.4324z 4
⑤ 所设计的高通滤波器的幅度特性与相位特性图示
版权所有 违者必究
第五章 第2讲
19
典型例题
0.0662272(1 z1)3 1)(1 0.6762858z1
0.3917468z 2
)
应滤 波 器 的 幅 度 响 应 与 相 位 响
版权所有 违者必究
第五章 第2讲
13
§5 设计IIR数字滤波器的频率变换法
又称为
频率变换的两种基本方法
“模拟域频率变换 法”
法一:从归一化模拟低通原型出发,先在模拟域内经频率
版权所有 违者必究
第五章 第2讲
22
设计IIR数字滤波器的频率变换法
模拟低通与数字带通的频率变换关系:
c
cos0 cosu sin u
,
sb
cos0 cossu sin su
,
cos 0
cos(l u
2
cos(u l
) )
2
当 c 时 1,反归一化处理:
p s z 2 2z cos0 1
求 得 c 0.588148
反归一化得:
H (s) Ha (s / c )
0.203451 (s 0.588148)(s2 0.588148s 0.345918)
相关主题