第六章 多采样率信号处理
在实际应用中,各系统之间的采样率往往是
不同的,例如,在音频范围内,广播系统的
采样率为 32kHz,CD唱机的采样率为 44.1kHz
,而数字录音带( DAT)的采样率为 48kHz。
这就需要设计一类数字滤波器,用于将被处
理信号的采样率转换成与相应系统所要求一
致的采样率。多采样率技术还有许多其他的
应用,过采样技术是其中的一个典型应用实
例。
6.1、采样率降低 —— 整数M倍抽取
?
?? M
T
T
M
f
MTTf
s
s ?
??
?
???
0 5 10 15 20 25 30 35 40-1
-0.5
0
0.5
1 x(n)
0 2 4 6 8 10 12 14 16 18 20-1
-0.5
0
0.5
1 y(m)
? ?
??
?
?
?
?
??
?
其他,
,~
MeH j
??
?
? ? ? ? ? ???
???
??
k
knxkhnw ? ? ? ?Mmwmy ?
x ( n ) < = = > ? ?fTjeX ?2
以± k sf ( k = 0,1,2,… ) 无限地重复
? ? ? ? ? ???
???
??
k
kMmxkhmy ? ? ? ? ? ??
?
???
??
k
kxkMmhmy
? ? ? ?
?
?
?
?
?????
??
其他,
,,,MMnnw
nw
? ? ? ? ? ?
? ? ??????
?
?
?
?
?
?
????
?
?
?
?
?
???
ne
M
nw
kMnnwnw
M
l
Mnlj
k
,
1
1
0
2 ?
?
? ? ? ? ? ?MmwMmwmy ???
? ? ? ? ? ??? ?
???
?
?
???
? ???
m
m
m
m zMmwzmyzY
? ? ? ?
? ?
? ?
? ??
? ?
? ?
?
?
?
?
?
?
?
???
?
?
???
?
?
?
?
???
?
??
?
?
?
?
?
?
?
?
?
?
?
?
?
?
??
1
0
12
1
0
2
1
0
2
1
1
1
M
l
MMlj
M
l m
MmMlmj
m
Mm
M
l
Mlmj
m
Mm
zeW
M
zemw
M
ze
M
mw
zmwzY
?
?
?
? ? ? ? ? ?zXzHzW ?
? ? ? ? ? ??
??
??
????????
M
l
MMljMMlj zeXzeH
M
zY ??
? ? ? ?? ? ? ?? ??
??
??
??????? ??
M
l
MljMljj eXeH
M
eY ?????
? ? ? ? ? ? ? ?? ? ? ?? ?? ?????? ????????? MjMjMjMjj eXeHeXeH
M
eY ???????
6.2, 采样率提高 —— 整数L倍内插
0 2 4 6 8 10 12 14 16 18 20-1
-0.5
0
0.5
1
0 10 20 30 40 50 60-1
-0.5
0
0.5
1
插零的结果 L=4 插零并且滤波的结果,L=4
LT
T ???
ss Lff ?
?
? ?
?
?
?
?
?
?
?????????
?
?
?
?
?
?
其他,
,,,,?LLm
L
m
x
mw
? ? ? ???
???
??
n
nznwzW
? ? ? ??
?
???
??
m
mLzmxzW
? ? ? ?LzXzW ?
? ? ? ?Ljj eXeW ?? ?? ?
? ?
??
?
?
?
???
?????
??
其他,
,~
L
G
eH j
?
??
? ? ? ?
??
?
?
?
??
?
?
?
其他,0
,
L
eGX
eY
j
j
?
???
? ? ? ???
??
?? L
L
Ljj deXGdeYy /
/
,,,,
22
1)0( ?
?
??
?
? ?
???
? ?
L
Gx
L
deXG j )0(
2 ??
???? ? ?
?
如要求 y(0)=x(0),则应有 G=L,即对理想的内插器
要求能恢复抽取前的信号,增益 G 必须等于 L 。
? ? ? ?
? ? ? ? 为整数LkrxrLmh
L
k
xkmhmy
r
k
?????
?
?
?
?
?
?
??
?
?
?
???
?
???
,
? ? ? ? ? ??
?
???
??
k
kwkmhmy
采样率变换一个有理因数
内插器 低通滤波 抽取器
sf
nx
采样率
)}({
sLf
my )}({
sLf
my )}({
sf
M
L
kw )}({
问题,低通滤波的指标如何确定?
6.3 抽取与内插的 FIR结构
一、抽取的 FIR结构
h(0)
h(1)
y(m)
h(N-1)
z-1
z-1
x(n)
M
h(0)
h(1)
y(m)
h(N-1)
x(n)
M
M
M
z-1
z-1
(a) FIR滤波与抽取器级联 (b) 提高运算效率的结构
抽取的 FIR结构
? ? ? ? ? ??
?
?
??
1
0
N
k
kMmxkhmy
抽取器的多相( Polyphase)滤波器结构
k=Mq+i,i=0,1,…M-1,q=0,1,…,N/M-1,N取 M的整数倍
? ? ? ? ? ?
])([)(
1
0
1
0
1
0
iqmMxiMqh
kMmxkhmy
M
i
M
N
q
N
k
????
??
? ?
?
?
?
?
?
?
?
1,,1,0;1,,1,0),()( ?????? MNmMiimMhmh i ??
1,,1,0),()( ???? MiimMxmx i ?
令
如
)}(,),2(),(),0({)(0 MNhMhMhhnh ?? ?
?
? ?
?
?
?
?
?
?
?
??
1
0
1
0
1
0
)(
)()()(
M
i
i
M
i
M
N
k
ii
my
kmxkhmy
y1 (m)
yM-1(m)
y0(m)
xM-1(m)
x0(m)
y(m)x(n)
M
M
M
hM-1 (m)
h1(m)
h0(m)
x1 (m)z-1
z-1
抽取器的多相滤波器结构
二、内插的 FIR结构
h(0)
h(1)
y(m)
h(N-1)
x(n)
L
z-1
z-1
h(0)
h(1)
y(m)
h(N-1)
x(n)
L
L
L
z-1
z-1
(a) 内插器与 FIR滤波级联 (b) 提高运算效率的结构
内插的 FIR结构
6.4,过采样技术
(a) (b)
抗混叠滤波器和 抗镜象后置模拟滤波器 的指
标分配
fs- fcfc0
1
H(f)
fs/2 f fc0
1
H(f)
fs/2 f
数字
低通 L
抽 取 器
模拟
抗混滤
波
A/D
变换器
x(t)
f’s fs
x(n)
过采样 A/D变换器
数字
低通L
内 插 器
后置
模拟滤
波
D/A
变换器
x(t)
fs f’s
x(n)
过采样 D/A变换器
数字
低通 L
抽 取 器
位
一阶 sigma-deltaA/D变换器
D/A
变换器
积分器
1/S
A/D
变换器+
位
'sf
'b b
sf
E(z)e(m)
一阶 Sigma-delta 的等效数学模型
+ ++
x(m)
X(z)
y(m)
Y(z)Z
-1
1
1
1
)( ?
?
?
?
z
zzH
? ? )()()()()( zEzYzXzHzY ???
)(
)(1
1
)(
)(1
)(
)( zE
zH
zX
zH
zH
zY
?
?
?
?
)()()()()( zEzHzXzHzY NSx ??
)(1
)()(
zH
zHzH
x ?? )(1
1)(
zH
zH NS
?
?
1)( ?? zzH
x
11)( ??? zzH
NS
22 )2/s i n (2)( ?? ?j
NS eH
)()1()()( 11 zEzzXzzY ?? ???
)()1()( mmxmy ????
dffH
f
s
s
f
f NS
s
ee
22/
2/'
22 )(1
' ?
?
? ??
'
2
2
' b
b
e
e
?
?
?
?
?
位b
fs
L? 低通数字
位b
fs'
内插器 噪声整形再量化器
环路滤波器
量化器
? ?nWMSB
? ?nWLSB
位'
'
b
fs
AD
b位至 '
? ?ZH
Q
数字输入
)(nw
? ? ? ? ? ?nwnwnw M S BL S B ??
? ?nx
'
sf 位b
? ?ne
? ?nw
? ?ny
? ?zH
'
sf 位'b
? ?ne
? ? ? ? ? ?nenwny ?? ? ? ? ? ? ?nenwny ??
? ? ? ? ? ?zEzWzY ??
? ? ? ? ? ? ? ?zEzHzXzW ??
? ? ? ? ? ?? ? ? ? ? ? ? ? ? ?zEzHzXzEzHzXzY NS????? 1
? ?zHzH NS ?? 1)(
1)( ?? zzH 11)( ??? zzH NS
112)( ?? ?? zzzH 21 )1()( ??? zzH NS
在实际应用中,各系统之间的采样率往往是
不同的,例如,在音频范围内,广播系统的
采样率为 32kHz,CD唱机的采样率为 44.1kHz
,而数字录音带( DAT)的采样率为 48kHz。
这就需要设计一类数字滤波器,用于将被处
理信号的采样率转换成与相应系统所要求一
致的采样率。多采样率技术还有许多其他的
应用,过采样技术是其中的一个典型应用实
例。
6.1、采样率降低 —— 整数M倍抽取
?
?? M
T
T
M
f
MTTf
s
s ?
??
?
???
0 5 10 15 20 25 30 35 40-1
-0.5
0
0.5
1 x(n)
0 2 4 6 8 10 12 14 16 18 20-1
-0.5
0
0.5
1 y(m)
? ?
??
?
?
?
?
??
?
其他,
,~
MeH j
??
?
? ? ? ? ? ???
???
??
k
knxkhnw ? ? ? ?Mmwmy ?
x ( n ) < = = > ? ?fTjeX ?2
以± k sf ( k = 0,1,2,… ) 无限地重复
? ? ? ? ? ???
???
??
k
kMmxkhmy ? ? ? ? ? ??
?
???
??
k
kxkMmhmy
? ? ? ?
?
?
?
?
?????
??
其他,
,,,MMnnw
nw
? ? ? ? ? ?
? ? ??????
?
?
?
?
?
?
????
?
?
?
?
?
???
ne
M
nw
kMnnwnw
M
l
Mnlj
k
,
1
1
0
2 ?
?
? ? ? ? ? ?MmwMmwmy ???
? ? ? ? ? ??? ?
???
?
?
???
? ???
m
m
m
m zMmwzmyzY
? ? ? ?
? ?
? ?
? ??
? ?
? ?
?
?
?
?
?
?
?
???
?
?
???
?
?
?
?
???
?
??
?
?
?
?
?
?
?
?
?
?
?
?
?
?
??
1
0
12
1
0
2
1
0
2
1
1
1
M
l
MMlj
M
l m
MmMlmj
m
Mm
M
l
Mlmj
m
Mm
zeW
M
zemw
M
ze
M
mw
zmwzY
?
?
?
? ? ? ? ? ?zXzHzW ?
? ? ? ? ? ??
??
??
????????
M
l
MMljMMlj zeXzeH
M
zY ??
? ? ? ?? ? ? ?? ??
??
??
??????? ??
M
l
MljMljj eXeH
M
eY ?????
? ? ? ? ? ? ? ?? ? ? ?? ?? ?????? ????????? MjMjMjMjj eXeHeXeH
M
eY ???????
6.2, 采样率提高 —— 整数L倍内插
0 2 4 6 8 10 12 14 16 18 20-1
-0.5
0
0.5
1
0 10 20 30 40 50 60-1
-0.5
0
0.5
1
插零的结果 L=4 插零并且滤波的结果,L=4
LT
T ???
ss Lff ?
?
? ?
?
?
?
?
?
?
?????????
?
?
?
?
?
?
其他,
,,,,?LLm
L
m
x
mw
? ? ? ???
???
??
n
nznwzW
? ? ? ??
?
???
??
m
mLzmxzW
? ? ? ?LzXzW ?
? ? ? ?Ljj eXeW ?? ?? ?
? ?
??
?
?
?
???
?????
??
其他,
,~
L
G
eH j
?
??
? ? ? ?
??
?
?
?
??
?
?
?
其他,0
,
L
eGX
eY
j
j
?
???
? ? ? ???
??
?? L
L
Ljj deXGdeYy /
/
,,,,
22
1)0( ?
?
??
?
? ?
???
? ?
L
Gx
L
deXG j )0(
2 ??
???? ? ?
?
如要求 y(0)=x(0),则应有 G=L,即对理想的内插器
要求能恢复抽取前的信号,增益 G 必须等于 L 。
? ? ? ?
? ? ? ? 为整数LkrxrLmh
L
k
xkmhmy
r
k
?????
?
?
?
?
?
?
??
?
?
?
???
?
???
,
? ? ? ? ? ??
?
???
??
k
kwkmhmy
采样率变换一个有理因数
内插器 低通滤波 抽取器
sf
nx
采样率
)}({
sLf
my )}({
sLf
my )}({
sf
M
L
kw )}({
问题,低通滤波的指标如何确定?
6.3 抽取与内插的 FIR结构
一、抽取的 FIR结构
h(0)
h(1)
y(m)
h(N-1)
z-1
z-1
x(n)
M
h(0)
h(1)
y(m)
h(N-1)
x(n)
M
M
M
z-1
z-1
(a) FIR滤波与抽取器级联 (b) 提高运算效率的结构
抽取的 FIR结构
? ? ? ? ? ??
?
?
??
1
0
N
k
kMmxkhmy
抽取器的多相( Polyphase)滤波器结构
k=Mq+i,i=0,1,…M-1,q=0,1,…,N/M-1,N取 M的整数倍
? ? ? ? ? ?
])([)(
1
0
1
0
1
0
iqmMxiMqh
kMmxkhmy
M
i
M
N
q
N
k
????
??
? ?
?
?
?
?
?
?
?
1,,1,0;1,,1,0),()( ?????? MNmMiimMhmh i ??
1,,1,0),()( ???? MiimMxmx i ?
令
如
)}(,),2(),(),0({)(0 MNhMhMhhnh ?? ?
?
? ?
?
?
?
?
?
?
?
??
1
0
1
0
1
0
)(
)()()(
M
i
i
M
i
M
N
k
ii
my
kmxkhmy
y1 (m)
yM-1(m)
y0(m)
xM-1(m)
x0(m)
y(m)x(n)
M
M
M
hM-1 (m)
h1(m)
h0(m)
x1 (m)z-1
z-1
抽取器的多相滤波器结构
二、内插的 FIR结构
h(0)
h(1)
y(m)
h(N-1)
x(n)
L
z-1
z-1
h(0)
h(1)
y(m)
h(N-1)
x(n)
L
L
L
z-1
z-1
(a) 内插器与 FIR滤波级联 (b) 提高运算效率的结构
内插的 FIR结构
6.4,过采样技术
(a) (b)
抗混叠滤波器和 抗镜象后置模拟滤波器 的指
标分配
fs- fcfc0
1
H(f)
fs/2 f fc0
1
H(f)
fs/2 f
数字
低通 L
抽 取 器
模拟
抗混滤
波
A/D
变换器
x(t)
f’s fs
x(n)
过采样 A/D变换器
数字
低通L
内 插 器
后置
模拟滤
波
D/A
变换器
x(t)
fs f’s
x(n)
过采样 D/A变换器
数字
低通 L
抽 取 器
位
一阶 sigma-deltaA/D变换器
D/A
变换器
积分器
1/S
A/D
变换器+
位
'sf
'b b
sf
E(z)e(m)
一阶 Sigma-delta 的等效数学模型
+ ++
x(m)
X(z)
y(m)
Y(z)Z
-1
1
1
1
)( ?
?
?
?
z
zzH
? ? )()()()()( zEzYzXzHzY ???
)(
)(1
1
)(
)(1
)(
)( zE
zH
zX
zH
zH
zY
?
?
?
?
)()()()()( zEzHzXzHzY NSx ??
)(1
)()(
zH
zHzH
x ?? )(1
1)(
zH
zH NS
?
?
1)( ?? zzH
x
11)( ??? zzH
NS
22 )2/s i n (2)( ?? ?j
NS eH
)()1()()( 11 zEzzXzzY ?? ???
)()1()( mmxmy ????
dffH
f
s
s
f
f NS
s
ee
22/
2/'
22 )(1
' ?
?
? ??
'
2
2
' b
b
e
e
?
?
?
?
?
位b
fs
L? 低通数字
位b
fs'
内插器 噪声整形再量化器
环路滤波器
量化器
? ?nWMSB
? ?nWLSB
位'
'
b
fs
AD
b位至 '
? ?ZH
Q
数字输入
)(nw
? ? ? ? ? ?nwnwnw M S BL S B ??
? ?nx
'
sf 位b
? ?ne
? ?nw
? ?ny
? ?zH
'
sf 位'b
? ?ne
? ? ? ? ? ?nenwny ?? ? ? ? ? ? ?nenwny ??
? ? ? ? ? ?zEzWzY ??
? ? ? ? ? ? ? ?zEzHzXzW ??
? ? ? ? ? ?? ? ? ? ? ? ? ? ? ?zEzHzXzEzHzXzY NS????? 1
? ?zHzH NS ?? 1)(
1)( ?? zzH 11)( ??? zzH NS
112)( ?? ?? zzzH 21 )1()( ??? zzH NS