1
0.5
0 -3
-2
-1
0
1
2
3
Nomalized frequency (/pi)
Digital Signal Processing
序列信号的抽取运算
降采样后频谱扩展了M倍
必须在降采样前增加抗混叠低通滤波防止扩展后频谱的混叠 抽取运算
x D (n)
k
x( k ) h

D
(nM k )
60
80
100
120
2
decimated Signal
1 0 -1 -2
0
5
10
15
20
25
30
Digital Signal Processing
非整数采样率转换
f s' f s L / M
1, [ , ], N max( L, M ) H (e j ) N N 0, 其他
Matalb非整数采样率运算函数：Y= resample (x,r,k,alpha)
Digital Signal Processing
9．2内插和抽取滤波器的设计和实现
内插和抽取滤波的特殊性
hI (n)内插运算中的数字抗镜像低通滤波器
hD (n)抽取数字抗混叠低通滤波器
多采样率系统中，FIR数字滤波器的计算效率优于IIR滤波器
由抗混叠条件得：
G0 ( z ) H ( z) 1 G0 ( z ) C ( z ) H1 ( z ), G1 ( z ) C ( z ) H 0 ( z ), C ( z )有理分式 G1 ( z ) H 0 ( z )
由纯延迟条件得：
C( z) H0 ( z)H1 ( z) H0 ( z)H1 (z) cz k
15
20
25
30
35
40
Digital Signal Processing
降采样序列和原始序列的频谱关系 x' (n)
X DN ( z )
x(n), n kM

0, n kM
n
,k Z
k / M
n
x
1/ M

DN
( n) z
n

n M 1 k 0
Digital Signal Processing
序列信号的降采样
M倍降采样：每隔M-1点抽取一点 xDN (n) x(nM )
1
original sample
0.5 0 -0.5 -1
0
20
40
60
80
100
120
140
160
1
4X downsample
0.5 0 -0.5 -1
0
5
10
Digital Signal Processing
内插数字抗镜像低通滤波器设计考虑
过渡带要求
通带和阻带纹波提出很高要求
上例仅考虑过渡带要求，并采用balckman窗： N 1713
Digital Signal Processing
抽取滤波器的设计
f s f s1 48KHz f p 1.8KHz f stop f s2 f p 4 1.8 2.6 KHz 2 ( f stop f p ) / f s / 3
i
Mi
h( Mi 1) z
i 1

( Mi 1)
... h( Mi M 1) z ( Mi M 1)
i

Mi
z
i
i
h(Mi 1) z

Mi
... z
( M 1)
i
FIR型内插滤波实现
Digital Signal Processing
FIR型抽取滤波实现
Digital Signal Processing
内插和抽取的级联实现
2级内插过程
f s f s2 2.4 KHz f p 0.18KHz f stop f s1 f p 0.4 0.18 0.22 KHz 2 ( f stop f p ) / f s / 30
E j ( z)
i
h(Mi j) z

i
H ( z ) 某一多相分解满足： Ek ( z M ) , k固定
, n 0 h(Mn k ) 0, 其他 n
结论： h(n) 等间隔地存在零值，计算量减少
Digital Signal Processing
双通道滤波器组精确重建的数学分析
Y ( z) 1 H 0 ( z )G0 ( z ) H1 ( z )G1 ( z ) X ( z ) 2 1 H 0 ( z )G0 ( z ) H1 ( z )G1 ( z ) X ( z ) 2
精确重建 Y ( z ) X ( z )条件 抗混叠条件
分析和综合滤波器的多相结构
H k ( z ) z Ekj ( z ), Gk ( z ) z ( M 1 j ) R jk ( z M )
j M j 0 j 0 M 1 M 1
精确重建条件：
E ( z ) R( z ) I
Digital Signal Processing
Digital Signal Processing
内插和抽取数字滤波器的实现结构
级联等效
Digital Signal Processing
滤波器的多相分解
h(k ), k ~

H ( z)
n
h( n) z

n
i
h(Mi) z h(Mi) z
Digital Signal Processing
9．1序列信号的采样率变换
序列信号的升采样
L倍升采样
x(n / L), n kL xUP (n) 0, n kL
k Z
Digital Signal Processing
升采样序列和原始序列的频谱关系
X UP ( z )
j 0
M 1
Digital Signal Processing
II型多相分解结构
R j ( z) EM 1 j ( z)
H ( z ) z ( M 1 j ) R j ( z M )
j 0
M 1
Digital Signal Processing
Nyquist或1/M带滤波器
1, [ , ] j H D (e ) M M 0, Others
Digital Signal Processing
Matalb用于抽取函数运算：Y= decimate (x,r,k,alpha)
2
original signal
1 0 -1 -2
0
20
40
Digital Signal Processing
序列信号的内插运算
升采样后产生L-1个镜像频谱
必须增加抗镜像低通滤波滤去L-1个镜像频谱 内插运算
x I ( n)
k
x(k )h (n kL)
I

L, , j H I (e ) L L 0, Others
h(Mi M 1) z

( Mi M 1)
E j ( z)
i h ( Mi j ) z
H ( z) z j E j ( z M )
j 0
M 1
Digital Signal Processing
I型多相分解结构
H ( z) z j E j ( z M )
Digital Signal Processing
双通道QMF滤波器组
设定两个分解滤波器为：
H1 ( z) H0 ( z)
C ( z) 1
H1 (e j ) H 0 (e j ( ) )
由H1(z)和H0(z)的抗混叠条件得
G1 ( z) H1 ( z) H0 ( z)
双通道CQF滤波器组
设定两个分解滤波器为：
x(nM ) z
1/ M k WM )

n

n
x '(nM ) z


k
x '(k ) z
j 2k M
X '( z
1.5
1 ) M
X (z
X DN (e ) X (e
j
)
original
1
0.5
0 -3
-2
-1
0
1
2
3
1.5
4X downsample
n


xUP (n) z
n

k


xUP (kL) z kL X ( z L )
X UP (e j ) X (e jL )
1.5
original signal
1
0.5
0 -3
-2
-1
0
1
2
3
1.5
4X upsample
1
0.5
0 -3
-2
-1
0
1
2
3
Nomalized frequency (/pi)
双通道滤波器组及其精确重建
理想 分解 和重 建滤 波器
Digital Signal Processing
实际分解和重建滤波器
理想分解重建滤波器属于非因果系统，物理上无法实现