第9章 多采样率信号处理
多采样率信号处理广泛应用于要求转换采样率,或要求系统工作在多采样率状态的信号处理系统中。如多种媒体——语音、视频、数据的传输,它们的频率很不相同,采样率自然不同,必须实行采样率的转换;又如信号要在两个时钟频率的数字系统中传输时,为了便于信号的处理、编码、传输和存储,要求根据时钟频率对信号的采样率加以转换;再如一种信号处理算法在系统的不同部分采用不同的采样率(如子带编码等),使处理更加有效,等等。
本章首先介绍直接在数字域对离散时间信号进行采样率转换的抽取(Decimation )和内插(Interpolation )方法,然后讨论抽取滤波器与内插滤波器的设计与实现方法。由不同采样率构成的系统称为多采样率系统,大部分多采样率系统使用了滤波器组,以正交镜像滤波器组(QMF )为基础的树状结构滤波器组是一典型的多采样率系统,它与离散小波(Wavelet )变换的关系密切。本章最后由正交镜像滤波器组的概念引入小波变换的基本原理和多分辨率分析的概念,以利开拓思路,为进一步深入学习打下基础。
9.1 离散信号的抽取与内插
9.1.1 抽取与内插的时域描述
离散序列的抽取与内插是多采样率系统中的基本运算,抽取运算将降低信号的采样频率,内插运算将提高信号的采样频率。离散序列x (k ) 的M 倍抽取定义为
x D (k ) =x (M k ), k ∈ (9-1-1)
其中M 为一正整数。抽取运算的框图如图9-1所示。
x (k ) D (k )
图9-1 M倍抽取运算的框图
图9-2画出了M =3时序列抽取的示意图。由图可知,离散序列的抽取表示保留第M 个样本点,而去除两个样本之间的M-1个样本点,设原离散信号x (k ) 的采样周期为T ,经M 倍抽取后的信号x D (k ) 的采样周期为T ' ,满足T ' =MT 。为了强调此概念,在图9-2中,
'
有意将抽取后的序列的间隔画为原序列的3倍。这时新的采样频率f s 为
f s =
'
1T
'
=
1MT
=
f s M
(9-1-2)
式中,f s 为原有的采样频率。
(a )原序列 (b) 3倍抽取后的序列
图9-2 离散序列的抽取
离散序列x (k ) 的L 倍内插定义为
⎧x
(k /L ), k =0, ±L , ±2L ,
x I (k ) =⎨ (9-1-3)
其他⎩0 ,
其中L 为一正整数。内插运算的框图如图9-3所示
x (k ) 图9-3 整数L 倍内插运算框图
I (k )
图9-4画出L =3时序列内插运算的示意图。由图可知,原序列x (k ) 中的所有样本都保留在内插后的序列x I (k
) 中,即内插运算不丢失信息。离散序列的内插运算是在原序列x (k ) 的每两个样本点之间插入L -1个零值样本点。
若对内插后的序列通过一个低通滤波器进行平滑处理,则可将序列中的零值转换为内插值,使得低通滤波后的序列的采样率是原序列的L 倍。如果序列x (k ) 的抽样间隔为T ,则低通滤波后序列的采样间隔为T ' =T /L 。
(a )原序列 (b )3倍内插后的序列
图9-4 离散序列的内插
9.1.2 抽取与内插的变换域描述
以上分析了序列抽取与内插的时域特性,下面从频域和z 域讨论抽取与内插的变换域特
性。由z 变换的定义,M 倍抽取后的序列x D (k ) 的z 变换X D (z ) 为
∞
∞
X D (z ) =
∑
k =-∞
∞
x D (k ) z x (n ) z
-k
=
∑
k =-∞
∞
x (M k ) z
-k
=
∑
n =-∞n =kM
-
n M
=
∑
n =-∞
x (n ) δ M (n ) z
-
n M
(9-1-4)
~
式中δ M (n ) 是周期为M 的单位脉冲序列。由离散傅里叶级数可知,δM (n ) 可以表示为
1
δ M (n ) =
M
M -1
∑e
l =0
j 2πnl /M
(9-1-5)
将(9-1-5)式代入(9-1-4) 式可得
X D (z ) =
1M
∞
M -1
∑∑
n =-∞l =0
x (n ) e
j 2πnl /M
z
-
n M
=
1M
M -1
∑
l =0
X (e
-j
2πM
l
1
z M ) (9-1-6)
将z =e j ω代入式(9-1-6)可得M 倍抽取后序列x D (k ) 的频谱为
X D (e
j ω
) =
1M
M -1
∑
l =0
X (e
j
ω-2πl
M
) (9-1-7)
式(9-1-7)表明,M 倍抽取后序列x D (k ) 的频谱可由下列步骤获得: (1)将X (e j ω) 扩展M 倍得到X (e j ω/M ) ,X (e j ω/M ) 的周期为2πM 。 (2)将X (e j ω/M ) 右移2π的整数倍得到{X (e j (ω-2πl ) /M ); l =0,1, , M -1}。 (3)将(2)中的M 个周期为2πM 的函数相加并乘以因子1/M ,得到周期为2π的M 倍抽取后序列的频谱X D (e
j ω
) 。
图9-3画出了3倍抽取后序列的频谱。由图可知,抽取后序列的频谱没有混叠。一般地,如果低频信号x (n ) 的频谱是带限,即在区间[-π, π]范围内有
X (e
j ω
,|ω|>π/M (9-1-8) ) =0
则M 倍抽取后信号的频谱不会发生混叠。式(9-1-8)称为序列M 倍抽取不混叠的Nyquist 条件。
图9-3 原序列的频谱和3倍抽取后序列的频谱
由z 变换的定义,L 倍内插后的序列x I (n ) 的z 变换X I (z ) 表示为
∞
∞
∞
X I (z ) =
∑
n =-∞
x I n (z )
-n
=
∑
n =-∞
n =m , L ∈ m
x n (L /z ) =
-n
∑
m =-∞
x m z (
-m L
) =X z
L
( ) (9-1-9)
将z =e j ω代入式(9-1-9)可得L 倍内插后的序列x I (n ) 的频谱为
X I (e
j ω
) =X (e
j ωL
) (9-1-10)
式(9-1-9)表明L 倍内插后的序列x I (n ) 的频谱X I (e j ω) 是原序列频谱X (e j ω) 的L 倍压缩。
图9-6 原序列的频谱和4倍内插后序列的频谱
图9-6画出了L =4时,原序列的频谱和内插后序列的频谱。由图可知,内插后序列在区间
[-π
X (e
/4, π/4]内的频谱,是由原信号在区间[-π, π]的频谱压缩4倍得到的。除了相差一个
) 的周期为π/2。内插序列在区间[-π/4, π/4]将在区间[-π, -π/4]和[π/4, π]
尺度因子外,两个频谱的形状保持不变。由于原序列频谱X (e j ω) 的周期为2π,因而
j 4
ω
内重复3次,这些重复的部分称为镜像频谱。
9.2抽取滤波器与内插滤波器
9.2.1 抽取滤波器
通常抽取后的离散序列的频谱将会出现混叠。为了避免混叠,可在信号抽取前利用低通滤波器对信号进行滤波,如图9-7所示。该滤波器称为抽取滤波器(decimation filter)。
x (n ) y (n )
图9-7 M 倍抽取滤波系统
由图9-3可知,抽取滤波器可以是截频为π/M 的理想低通滤波器,其频率响应为:
⎧1⎪⎪j ω
H (e ) =⎨
⎪0⎪⎩
ω|≤|πM
πM
(9-2-1)
如果低通滤波器的单位冲激响应为h (n ) ,则滤波器的输出w (n ) 为
∞
w (n )=
最后的输出y (m ) 为
∑h (k ) x (n -k ) (9-2-2)
k =-∞
y (m ) =w (mM ) (9-2-3)
将式(9-2-2)和式(9-2-3)合起来,y (m)与x (n ) 之间的关系为
∞
∞
y (m )=
∑
k =-∞
h (k x ) M (m -k =)
∑
k =-∞
h M (m -k x k ) ( (9-2-4)
式(9-2-4)表明,在计算M 倍抽取滤波器的输出时,只需计算抽取滤波器每M 个输出中的一个样本。利用
W (z ) =H (z ) X (z ) (9-2-5)
并根据(9-1-6)式,有 Y (z ) =
1M
M -1
1
-j 2πl /M
∑W (e
l =0
z M ) =
1M
M -1
1
-j 2πl /M
1
-j 2πl /M
∑H (e
l =0
z M ) X (e z M ) (9-2-6)
z 在单位圆上取值,即z =e j ω时,可得
Y (e
j ω
) =
1M
M -1
∑
l =0
H (e
j
ω-2πl
M
) X (e
j
ω-2πl
M
) (9-2-7)
式(9-2-7)是对输入信号x (n ) 进行滤波、抽取后的频域表达式。
9.2.2 内插滤波器
信号的内插不会引起频谱的混叠,但会产生镜像频谱,如图9-6所示。为了消除这些镜像频谱,可将内插后的信号通过低通滤波器,如图9-8所示。该低通滤波器称为内插滤波器(interpolation fliter )。由图9-6可知,内插滤波器可以是截止频率为π/L 的理想低通滤波器。该滤波器可以滤除信号w (n ) 频谱中的镜像频谱,仅保留[-π/L , π/L ]范围内的频谱。
内插滤波器的频率响应为
π⎧
G ω ≤ |⎪j ω
H (e ) =⎨L (9-2-8)
⎪0 其他⎩
x (n )
y (n )
图9-8 L 倍内插滤波系统
根据图9-8可知
j ωj ω
Y (e j ω) =H (e ) W (e =)
H (e
j ω
(9-2-9) ) X (e
j ωL
根据式(9-2-8)和式(9-2-9),可近似得出
π⎧j ωL
G X (e ) ω ≤ ||⎪j ω
Y (e ) ≈⎨L (9-2-10)
⎪0 其他⎩
在n =0时刻,有
y (0)=
12πG 2π
⎰π
-
π
Y (e
j ω'
) d ω' =) d ω
=
G 2π
⎰π
-
π/L
/L
X (e
j ω' L
) d ω'
=
⎰πX (e
-
π
j ω
(9-2-11)
G x (0)L
L
因此,如要求y (0)=x (0),则应有G =L ,即对理想的内插器要求能恢复内插前的信号,增益G 必须等于L 。
若设内插滤波器的单位冲激响应为h (n ) ,则
∞
y (n )=
∑
k =-∞
h (n - (9-2-12) k ) w ( k )
由式(9-1-3)和式(9-2-12)式,可以得到
∞
y (n )=
∑
k =-∞
h (n -k ) x (
k L
∞
) =
∑
r =-∞
h (n -rL ) x (r ), k 是L 的整数倍 (9-2-13)
由于内插滤波器的输入信号w (n ) 中每L 个样本中只有一个非0样本,所以内插滤波器的计算量只有常规系统的1/L 。
9.2.3 有理数倍抽样率转换
给定信号x (n ) ,若希望将抽样率转变为L /M 倍,可以通过把M 倍抽取和L 倍内插结合起来得到。一般是先做L 倍的插值,再做M 倍的抽取。这是因为先抽取会使x (n ) 的数据点减少,会产生信息的丢失,并且可能产生频率响应的混叠失真。例如:如果x (n ) 的抽样频率f s =2f h ,此时x (n ) 的基带正好在容许的频带上限之内,即在折叠频率|ω|≤π以内,现在要将x (n ) 的抽样频率转换为3f s /2,此时M =2,L =3,如果先做2倍抽取,则会先丢失掉一些数据,而且信号的数字频带要增加2倍,必然产生混叠失真。为了不产生混叠失真,
必须将防混叠的低通滤波器频带限制在|ω|≤π/2内,这样会丢失很多信息。如果先做3倍插值,使数字频带先缩小3倍,再做2倍抽取,则信号的数字频带变成原信号数字频带的2/3,因而不会产生混叠失真。所以,应先对信号做L 倍的插值,再做M 倍抽取,结构上就是两者的级联。图9-9实现以有理数L /M 来改变采样率的系统,新系统输出信号的采样率为f s ' =Lf s /M 。
图9-9 有理数倍L/M采样率转换框图
图9-9中将内插后的抗镜像滤波器和抽取前的抗混叠滤波器合并为一个数字低通滤波器h (n ) 。由于此滤波器同时用作插值和抽取的运算,因而,它的理想频率响应为 H (e j ω) =⎨
≤|⎧⎪L , ω⎪⎩0, 其他
|{πm i L n π
M /},
(9-2-14)
/
9.2.4 抽取滤波器的FIR 结构和多相结构
从前面的讨论可以知道,抽取器需要一个数字低通滤波器用于滤除可能会引起混叠失真的频谱分量。FIR 滤波器具有绝对稳定,容易实现线性相位特性,特别是连同抽取器一起,采用合理的结构,可以大大提高运算效率。
对于图9-7的M 倍抽取器的直接型FIR 滤波器的实现结构如图9-10(a )所示。这个直接型FIR 滤波器结构概念清楚明白,其实现也很简单。但这是h (n ) 工作在高采样率状态,
x (n ) 的每一个采样点均要和FIR 滤波器的系数相乘。而滤波器输出中,每M 个样值中只抽
取一个作为最终的输出y (m ) ,丢弃了其中M -1个样值,所以该结构效率很低。 为了提高直接型FIR 滤波器结构的运算效率,将图9-10(a )中抽取操作嵌入到FIR 滤波器结构中,如图9-10(b )所示,由级联的延迟器z -1移入各抽头的x (n ) 先做抽取,再和
h (n ), n =0,1, , N -1相乘,由于工作在低采样率状态,系统的运算速率降低了M 倍。FIR
滤波器则可写成
N -1
y (m ) =
∑h (k ) x (Mm
k =0
-k ) (9-2-15)
图9-10 抽取的FIR 结构
(a )FIR 滤波器与抽取器级联 (b )提高运算效率的结构
多相(Polyphase )滤波器结构是按整数因子抽取的另一种高效实现结构。通常取N 是M 的整数倍,若令式(9-2-15)中的k =Mq +i , i =0,1, , M -1, q =0,1, , N /M -1,则
N -1
y (m ) =
∑h (k ) x (M m -k )
k =0
N M -1M
-1
=
i =0
(9-2-16)
h (M q +i ) x [M (m -q ) -i ]
∑∑
q =0
令h i (m ) =h (mM +i ) ,i =0,1, , M -1, m =0,1, , N /M -1为多相滤波器的子滤波器的单位脉冲响应,如h 0(m ) ={h (0),h (M ), h (2M ), , h (N -M )};并且
x i (m ) =x (mM -i ) ,i =0,1, , M -1
则式(9-2-16)可以写成M 个子载波的和的形式
N M -1M
-1
i
i
M -1
i
y (m ) =
N M -1
i
i
∑∑h (k ) x (m -k ) =∑y (m ) (9-2-17)
i =0k =0
i =0
式中,y i (m ) =示。
∑h (k ) x (m -k ) 为子滤波器输出。抽取器的多相滤波器结构如图9-11所
k =0
)
图9-11 抽取器的多相滤波器结构
9.2.5 内插滤波器的FIR 结构和多相结构
根据图9-8得到按整数因子L 内插系统的直接型FIR 滤波器结构如图9-12(a )所示,其中FIR 滤波器采用了转置型结构。该结构中x (n ) 插零后再进行滤波,使得FIR 滤波器要和大量的零值相乘,FIR 滤波器以高采样率运行,该结构效率很低。
类似抽取滤波器的FIR 结构,将图9-12(a )内插器嵌入到FIR 滤波器结构中的N 个乘法器之后,得到如图9-12(b )所示的结构,h (n ) 以低的运算速率与x (n ) 相乘后再插零,运算速率降低了L 倍。
图9-12 内插的FIR 结构
(a) 内插器与FIR 滤波器级联 (b) 提高运算效率的结构
仿照抽取器的多相滤波器结构,取内插因子为L 的低通滤波器h (n ) 的长度N 为L 的整数倍,则可以分解成L 个子滤波器
h i (n ) =h (nL +i ) ,i =0,1, , L -1, n =0,1, , N /L -1 (9-2-18)
由图9-12(b )得到图9-13形式的内插器的多相滤波器结构,图中
N L
-1
i
y i (n ) =
∑h (k ) x (n -k ) ,
k =0
i =0,1, , L -1 (9-2-19)
)
图9-13 内插器的多相滤波器结构
9.3 正交镜像滤波器组
9.3.1 数字滤波器组的一般概念和定义
数字滤波器组在语音分析、数据压缩、信号传输等领域均要用到。M 通道滤波器组的基本结构如图9-14所示。H 0(z ) 、H 1(z ) 、…H M -1(z ) 组成了分析滤波器组,输入信号x (n ) 通过这一组滤波器后,得到的x 0(m ) 、x 1(m ) 、…x M -1(m ) 是x (n ) 的子带信号;而合成过程
ˆ(n ) 表示该综合滤波器的输出信号。一般来说,各路则在滤波器组的综合滤波器中完成,x
信号的采样率与输入或输出信号的采样率是不同的,因此,系统本质上是多采样率的。设计滤波器组的一个重要任务或基本准则是,综合应用分析滤波器组H 0(z ) 、H 1(z ) 、…H M -1(z ) 和综合滤波器组G 0(z ) 、G 1(z ) 、…G M -1(z ) ,抵消或尽可能地抑制混叠失真。
x 0(m x 1(m 9.3.2 正交镜像滤波器组
M -1(m
图9-14 M 通道滤波器组对信号分解和综合示意图
如果M =2,图9.14的M 通道滤波器组就简化为两通道滤波器组,如图9-15(a )所示。利用前述信号抽取与内插理论,可得到两通道滤波器组的输入与输出关系为
ˆ(z ) =1{H (z ) G (z ) +H (z ) G (z ) }X (z ) X 0011
2
+12
(9-3-1)
{H 0(-z ) G 0(z ) +H 1(-z ) G 1(z ) }X (-z )
式中X (-z ) 分量是由抽取过程中的混叠产生的响应,称之为混叠项。式(9-3-1)表明,为了消除混叠项对输出的影响,系统需满足
H 0(-z ) G 0(z ) +H 1(-z ) G 1(z ) =0 (9-3-2)
为满足式(9-3-2)约束条件,存在多种选择方案。一种简单的选择方案为
G 0(z ) =H 1(-z ), G 1(z ) =-H 0(-z ) (9-3-3)
这样即可消除混叠失真对整个系统输出响应的影响。滤波器组输出项中无混叠项的滤波器组称为无混叠滤波器组。
) )
(a )
2
(b )
图 9-15 正交镜像滤波器组
(a )结构原理 (b )镜像对称的幅频响应
对于两通道无混叠滤波器组,若分析滤波器组满足
H 1(z ) =H 0(-z ) (9-3-4)
则称该滤波器组为π/2镜像滤波器组。当h 0(k ) 为实系数时,式(9-3-4)意味着
H 1(e j ω) =H 0(e j (π-ω) ) (9-3-5)
j ωj ω
即它们的幅度响应关于π/2镜像对称,如图9-15(b )所示。如果H 0(e ) 和H 1(e ) 两者j ωj ωj ω
没有重合,即当π/2≤ω≤π时,H 0(e ) ≡0,则H 0(e ) 和H 1(e ) 是正交的,称为正
j ωj ω
交镜像滤波器组(QMFB ,Quadrature Mirror Filter Bank)。在实际应用中,H 0(e ) 和H 1(e )
可以有少量重叠,也称为QMFB 。
令X (e j ω) 、X 0(e j ω) 、X 1(e j ω) 、H 0(e j ω) 和H 1(e j ω) 分别是x (n ) 、x 0(m ) 、x 1(m ) 、
h 0(n ) 和h 1(n ) 的DTFT ,可以得到
X 0(e
j ω
) =
1212
{X (e
j ω/2
) H 0(e
j ω/2
) +X [e
j (ω+2π) /2
]H 0[e
j (ω+2π) /2
]} (9-3-6)
及
X 1(e j ω) =
{X (e
j ω/2
) H 1(e
j ω/2
) +X [e
j (ω+2π) /2
]H 1[e
j (ω+2π) /2
]} (9-3-7)
j ωˆˆ0(m ) 、x ˆ1(m ) 、g 0(n ) 和类似的,令X 0(e ) 、X 1(e j ω) 、G 0(e j ω) 和G 1(e j ω) 分别是x g 1(n ) 的DTFT ,得到综合关系式为
ˆ(e j ω) =X ˆ(e j 2ω) G (e j ω) +X ˆ(e j 2ω) G (e j ω) (9-3-8) X 0011
ˆ0(m ) =x 0(m ) ,x ˆ1(m ) =x 1(m ) ,可以得到滤波 用分析器的输出作为综合器的输入,即x
器组的输入输出频域关系式为
ˆ(e j ω) =1[H (e j ω) G (e j ω) +H (e j ω) G (e j ω)]X (e j ω) +X 0011
2 {H 0[e
2
1
j (ω+π)
(9-3-9)
j (ω+π)
]G 0(e
j ω
) +H 1[e
j (ω+π)
]G 1(e
j ω
)}X [e ]
ˆ(e j ω) 的有用信号变换,第二项则代表不希望出现的这个式中的第一项代表从X (e j ω) 到X
频域混叠分量。
为了去除不希望出现的频域混叠分量,必须满足
H 0[e
j (ω+π)
]G 0(e
j ω
) +H 1[e
j (ω+π)
]G 1(e
j ω
) =0 (9-3-10)
只要仔细选择滤波器,将这个关系式中的前后两项抵消掉,就可以实现这个条件。
通常可先设计一公共低通滤波器h (n ) ,再由它得到所有的分析和综合滤波器,即
h 0(n ) =h (n ) (9-3-11) h 1(n ) =(-1) h (n ) (9-3-12)
等效的,它们的DTFT 满足
j ωj ω
H 0(e ) =H (e ) (9-3-13)
n
H 1(e
j ω
) =H [e
j ω
j (ω+π)
] (9-3-14)
上面的条件表明,滤波器H 0(e (b )所示。
j ω
) 和H 1(e
) 对于频率ω=π/2是镜像对称的,如图9-15
将式(9-3-13)和式(9-3-14)代入式(9-3-10),得
H [e
j (ω+π)
]G 0(e
j ω
) +H (e
j ω
) G 1(e
j ω
) =0 (9-3-15)
下面说明对G 0(e j ω) 和G 1(e j ω) 的要求。因为G 0(e j ω) 必须是一个低通滤波器,令
G 0(e j ω) =2H (e j ω) (9-3-16)
或等效为
g 0(n ) =2h (n ) (9-3-17)
式中2是与内插滤波器有关的增益因子。将式(9-3-16)代入式(9-3-15)中,得
G 1(e j ω) =-2H [e j (ω+π) ] (9-3-18)
可以看出G 1(e j ω) 必须具有高通性质,其冲激响应为
g 1(n ) =-2(-1) n h (n ) (9-3-19)
因此,图9-15中分析和综合滤波器的设计转换为对公共低通滤波器H (e j ω) 的设计。将式(9-3-13)、式(9-3-14)、式(9-3-16)和式(9-3-18)代入式(9-3-9),得到这种正交镜像滤波器组最终的输入和输出关系
ˆ(e j ω) ={H 2(e j ω) -H 2[e j (ω+π) ]}X (e j ω) (9-3-20) X
这里代表混叠的第二项已被消去,这意味着在分析结构中由抽取引起的混叠分量被综合结构中内插引起的镜像分量精确地抵消掉了。由式(9-3-20)可以看到,在低通滤波器H (e j ω) 满足条件
|H 2(e j ω) -H 2[e j (ω+π) ]|=1 (9-3-21)
时,分析和综合滤波器组构成的系统其增益将变成1。一般希望H (e j ω) 逼近理想低通条件
π⎧
1 0≤ω≤⎪⎪2j ω
(9-3-22) H (e ) =⎨
⎪0 π
QMF 公共低通滤波器的设计可以综合考虑式(9-3-21)和式(9-3-22)的约束条件,采用第7章的FIR 滤波器设计方法或借助计算机辅助设计的最优化方法实现。
9.4 离散小波变换
9.4.1 连续小波变换
由3.5节可知,在用DFT 对信号进行频谱分析时,要提高频率分辨率就必须加大信号时间截取的长度,这意味着信号分析只能得到某一时间段上信号所包含的频率成份,截取长度越大,关于某一频率分量的时间定位精确度就越低,即时间分辨率降低。小波变换将时域
信号变换到时-频域中,克服了DFT 的不足,它对于不同的时间和不同的频率具有不同的精确度,可以同时获得信号时域特性和频域特性。
对一个平方可积函数f (t ) 进行小波变换,就是用一个母小波函数ψ(t ) 的伸缩平移函数族{ψa , b (t ) }来对其进行展开,记为
WT {f (t ), ψa , b (t )}=
f (t ), ψa , b (t ) =
⎰
∞-∞
f (t ) ψa , b (t ) dt (9-4-1)
式中ψa , b (t ) =
(
t -b a
) ,b 为时间平移,a 为时间伸缩, , 表示两个函数的内积。母
小波函数ψ(t ) 必须满足一定的要求,实际上存在着许多母小波函数,因而有许多小波变换,并具有不同的特征。显然,式(9-4-1)中参数a 为一折中参数,在时间精度和频率分辨率之间取得折中。
这样一个连续参数a 和b 的小波变换具有相当大的冗余性,相应的反变换也不是唯一的,因此很少有实用性,下面将重点讨论二进小波(Dyadic wavelet)。
9.4.2 多分辨率分析
Mallat 提出了用多分辨率分析(或逼近)概念定义小波,将以前的各种正交小波基构造统一了起来,并给出了分解和重构算法。
令 表示实数域, 表示整数域,f (t ) ∈L 2( ) 。多分辨率分析建立在以下三个基本假设的基础上:
(1)存在着一簇具有不同分辨率的子空间,满足
V -∞⊂ ⊂V -2⊂V -1⊂V 0⊂V 1⊂V 2⊂ ⊂V ∞=L ( ) (9-4-2)
2
每一个子空间具有不同的基向量,它们给定了不同的时间分辨率。当下标i 增加,空间V i 的时间分辨率也增加。
(2)存在着一个尺度函数ϕ(t ) ,以及它的整数平移
ϕk (t ) =ϕ(t -k ), k ∈ (9-4-2)
构成空间V 0的一组正交基{ϕ
k (t ) },满足ϕk , ϕl =δkl , k , l ∈ 。
(3)伸缩规则性
f (t ) ∈V i f (2t ) ∈V i +1 (9-4-3)
利用伸缩规则性和尺度函数产生了V 0空间正交基的特性,可以立即得到所有空间V i 的基的集合
ϕi , k (t ) =2ϕ(2t -k )
i /2i
i , k ∈ (9-4-4)
函数ϕi , k (t ) 代表空间V
i 的一个正交基。因此有ϕi , k , ϕi , l =δkl , i , k , l ∈ 。上式中的2i /2保证模i , k (t ) 总为1,它独立于下标i 与k 。
由于ϕ(t ) ∈V 0,V 0⊂V 1,并且函数21/2ϕ(2t -n ), n ∈ ,在V 1空间ϕ(t ) 可以表示成函数ϕ(2t -n ) 的线性组合
ϕ(t ) =
式中h 0(n ) 为线性组合系数。
∑h
n
(n ) ϕ(2t -n ) (9-4-5)
利用引入的基,一个信号f (t ) ∈V i 可以表示为
f (t ) =
∑α
m
i
(m ) ϕi , m (t ) (9-4-6)
式中αi (m ) =以展开成
f (t ), ϕi , m (t ) 为展开系数。进一步,f (t ) 也是V i +1空间上的分量,从而也可
f (t ) =
∑α
n
i +1
(n ) ϕi +1, n (t ) (9-4-7)
与函数ϕ(2i t ) 比较,函数ϕ(2i +1t ) 沿着时间轴压缩了2倍,因此,在V i +1空间上的信号的时间分辨率是V i 空间上的2倍。由于子空间上的分辨率以2的幂次增加,因而称为二进尺度函数或二进多分辨率。
9.4.3 二进小波
对每一子空间V i ⊂V i +1,可以将V i +1表示成V i 与它的正交补空间的直和,即
V i +1=V i ⊕W i , i ∈ (9-4-8)
和V i 一样,希望找到一个确定的函数ψ(t ) ∈W 0,使得每个i ∈ ,函数系
ψi , k (t ) =2ψ(2t -k )
i /2i
i , k ∈ (9-4-9)
构成空间W i 的正交基,这里ψ(t ) 是前面介绍的母小波函数。因此,对于一个信号f (t ) ∈W i ,可以写成
f (t ) =
∑β
m
i
(m ) ψi , m (t ) (9-4-10)
式中βi (m ) =f (t ), ψi , m (t ) 。
由于W i ⊂V i +1,信号f (t ) ∈W i 也可以在V i +1空间中由尺度函数ϕi +1, k (t ) 展开。对于
i =k =0的母小波函数,也可以在空间V 1中表示成
ψ(t ) =
∑h (n ) ϕ(2t -n ) (9-4-11)
1n
式中h 1(n ) 为线性组合系数。类似地,V i 空间也可分解成直和的形式V i =V i -1⊕W i -1,V i -1也可进一步分解。根据式(9-4-2),信号空间可以表示成
L ( ) =V j ⊕W j ⊕W j +1⊕ ⊕W -1⊕W 0⊕W 1 (9-4-12)
2
式中下标j 是任意的,表示分解的深度。
由式(9-4-6)、式(9-4-10)和式(9-4-12),对于任意平方可积信号f (t ) ∈L ( ) ,都
2
可以用尺度函数与小波函数线性表示为
∞
f (t ) =
∑α
m
j
(m ) ϕj , m (t ) +
∑∑β
i =j
m
i
(m ) ψi , m (t ) (9-4-13)
式中∑αj (m ) ϕj , m (t ) 对应f (t ) 的低频分量,是信号f (t ) 的最粗略近似,由低分辨率的尺
m
∞
度信号ϕj , m (t ) 表示;∑
i =j
∑β
m
i
(m ) ψi , m (t ) 对应于信号f (t ) 的高频分量,由具有更高分辨率
的一系列不同尺度的小波信号ψi , m (t ) 表示,随着尺度j 是的增加,表示了信号f (t ) 中更丰富的细节。也就是说,展开系数αj (m ) 反映了信号f (t ) 中的低频分量的分布情况,而展开系数βi (m ) 反映了信号f (t ) 中的高频分量的分布情况,这些展开系数就是信号f (t ) 的离散小波变换(DWT ,Discrete Wavelet Transform)。
9.4.4 二进小波变换与滤波器组
在信号的小波分析中,实际上并不直接利用小波函数ψ(t ) 和尺度函数ϕ(t ) 计算信号的DWT ,而是利用其对应的小波函数系数h 1(n ) 和尺度函数系数h 0(n ) 。
设f (t ) ∈V i +1空间,有
f (t ) =
∑α
n
i +1
(n ) ϕi +1, n (t ) (9-4-14)
式中系数αi +1(n ) 必须已知。由于V i +1=V i ⊕W i ,这一信号可以唯一表示成两子空间的投影和的形式,这里的投影被展开成各自的基的展开式,即
f (t ) =
∑α
m
i
(m ) ϕi , m (t ) +∑βi (m ) ψi , m (t ) (9-4-15)
m
实际上,投影是由已知系数αi +1(n ) 来计算未知系数αi (m ) 和βi (m ) 。
根据式(9-4-5),V i 空间的基函数ϕi , m (t ) 可以由V i +1空间的基函数ϕi +1, m (t ) 递推得到。将式(9-4-5)代入式(9-4-4),得
ϕi , m (t ) =2ϕ(2t -m ) =2
i /2i i /2
∑h
l
(l ) ϕ(2
i +1
t -2m -l ) (9-4-16)
-1/2
h 0(l ) 代入上式,得 将关系式2m +l =n 和h 0(l ) =2
ϕi , m (t ) =
=
∑h
n
(n -2m ) 2
(i +1) /2
ϕ(2
i +1
t -n )
∑
n
h 0(n -2m ) ϕi +1, n (t )
(9-4-17)
同理,W i 空间中的基函数ψi , m (t ) 可以相应于V i +1空间展开成
ψi , m (t ) =
则展开系数αi (m ) 可写成
∑h (n -2m ) ϕ
1n
i +1, n
(t ) (9-4-18)
αi (m ) =
==
f (t ), ϕi , m (t ) =
f (t ), ∑h 0(n -2m ) ϕi +1, n (t )
n
∑h
n
(n -2m ) f (t ), ϕi +1, n (t ) (n -2m ) αi +1(n )
n =2m
(9-4-19)
∑h
n
=h 0(-n ) *αi +1(n )
式中*表示线性卷积。因此,系数αi (m ) 可通过αi +1(n ) 与序列h 0(-n ) 的卷积,并对其结果2倍抽取得到。类似地,βi (m ) 可以用下式求得。
βi (m ) =h 1(-n ) *αi +1(n )
n =2m
(9-4-20)
这一结果示于图9-16。因此,分解一个信号的系数可以通过单位脉冲响应为h 0(-n ) 和h 1(-n ) 的分析滤波器组求得。将信号f (t ) ∈V i +1投影到其子空间时,它具有两倍于子空间的分辨率。相应的,展开系数αi (m ) 和βi (m ) 是系数αi +1(n ) 的采样率的一半。
αi (n )
i (n )
图 9-16 用分析滤波器计算小波展开系数
联系到伸缩的规则性,V i +1到V i 上的投影相应于低通滤波,V i +1到W i 上的投影相应于高通滤波。因此,h 0(-n ) 是一低通滤波器,而h 1(-n ) 是一互补的高通滤波器。如果重复地将更多的两通道分析滤波器连接到一个两通道滤波器组的低通输出端,将实现式(9-4-13)的信号分解,其结构如图9-17所示。
图 9-17 树状结构的二进小波变换
与上述讨论的将一个信号投影到两个子空间的过程相反,可以将两个子空间V i 和W i 的信号合成为V i +1空间的一个信号,这一过程伴随着信号的时间精度提高了一倍。
这一信号的合成相应于两通道综合滤波器组的运算。在综合运算中,合成信号的未知系数αi +1(n ) ,可由子带信号的已知系数αi (m ) 和βi (m ) 计算得到。为了计算这些系数,需要用基ϕi , m (t ) 和ψi , m (t ) 表示ϕi +1, m (t ) 。因为V 1=V 0⊕W 0,所以存在着以下关系式
ϕ(2t ) =
利用关系式
∑k
(k ) ϕ(t -k ) +∑1(k ) ψ(t -k ) (9-4-21)
k
ϕi +1, n (t ) =2
=
(i +1) /2
ϕ(2
i +1
t -n )
i
i
∑2
k
1/2
0(k ) 2ϕ(2t -
n 2
-k ) +
∑2
k
1/2
1(k ) 2ψ(2t -
i /2i
n 2
-k )
(9-4-21)
令g 0(-2k ) =21/20(k ) ,g 1(-2k ) =21/21(k ) 和n +2k =2m 代入式(9-4-21),得
ϕi +1, n (t ) =
=
∑g
m
(n -2m ) 2
i /2
ϕ(2t -m ) +∑g 1(n -2m ) 2ψ(2t -m )
m
i i /2i
∑g
m
(n -2m ) ϕi , m (t ) +
∑g
m
1
(n -2m ) ψi , m (t )
(9-4-22)
因此,展开系数可写成
αi +1(n ) =
==
f (t ), ϕi +1, n (t ) =
f (t ), ∑g 0(n -2m ) ϕi , m (t ) +
m
∑g
m
1
(n -2m ) ψi , m (t )
∑g
n
(n -2m ) f (t ), ϕi , m (t ) +(n -2m ) αi (m ) +
∑g
n
1
(n -2m ) f (t ), ψi , m (t )
(9-4-23)
∑g
n
∑g
n
1
(n -2m ) βi (m )
式(9-4-23)可以进一步表示为
αi +1(n ) =g 0(n ) *αi () +g 1(n ) *βi () (9-4-24)
2
2
n
n
其中αi () 和βi () 分别为αi (n ) 和βi (n ) 的2倍内插序列。式(9-4-24)可以用图9-18的
2
2
n n
综合滤波器组来实现。由于内插过程,数字滤波器g 0(n ) 和g 1(n ) 的偶数下标才与系数
αi (m ) 和βi (m ) 相乘。式(9-4-24)就是小波重构算法,对应于离散逆小波变换(IDWT: Inverse Discrete Wavelet Transform)。
αi (m βi (m
图 9-18 用综合滤波器组合两个子空间的信号
对于图9-18结构的二进级联,可将不同子空间的信号合成到一起,得到了与图9-17相反的结构,示于图9-19所示。
图 9-19 树状结构的逆二进小波变换
采用二进树结构的分析滤波器组和综合滤波器组对信号进行分解或重构是非常有效的,它甚至不涉及函数ϕ(t ) 和ψ(t ) 的具体形式。
9.5 本章相关的MATLAB 命令及应用
在MA TLAB 信号处理工具箱中,提供了抽取函数decimate ()、内插函数interp ()和重采样函数resample ()等。
1. 抽取函数decimate 格式:y = decimate(x,r)
功能:对离散时间信号向量x 按抽取因子r 抽取,得到信号向量y ,相当于降低了采样率r 倍。向量y 的长度是原信号向量x 长度的1/r倍。
2. 内插函数interp 格式:y = interp(x,r)
功能:对离散时间信号向量x 按内插因子r 内插,得到信号向量y ,相当于提高了采样率r 倍。向量y 的长度是原信号向量x 长度的r 倍。
3. 重采样函数resample 格式:y = resample(x,p,q)
功能:对离散时间信号向量x 按有理数p/q倍的采样率重新采样,得到信号向量y 。p 和q 必须为正整数。向量y 的长度等于ceil(length(x)*p/q)。
4. 离散小波变换dwt
格式:[cA,cD] = dwt(X , ‘wname ’)
功能:对离散时间信号序列X 进行小波分析,cA 是信号离散小波变换对应的近似展开系数,cD 是对应的细节系数,wname 是小波名称。
5. 离散小波逆变换idwt 格式:X = idwt(cA,cD,’wname ’)
功能:实现一维单级离散小波逆变换,是函数dwt 的逆运算。各参数与函数dwt 相同。 6. 分解与重构滤波器组 wfilters
格式:[Lo_D, Hi_D, Lo_R,Hi_R] = wfilters(‘wname ’)
功能:产生用于信号DWT 和IDWT 对应的滤波器组。Lo_D为分解数字低通滤波器的单位脉冲响应,Hi_D为分解数字高通滤波器的单位脉冲响应,Lo_R为重构数字低通滤波器的单位脉冲响应,Hi_R为重构数字高通滤波器的单位脉冲响应。
9.5.1 时域和频域抽样率的改变
[例9-1] 编写程序对时域正弦信号
x (n ) =sin(2π⋅0.12n )
进行3倍内插和3倍抽取操作。 程序代码如下: % 程序9_1
% 对一个正弦信号进行3倍内插和3倍抽取 clf; n = 0:49; m = 0: 50*3 - 1; x1 = sin(2*pi*0.12*m); x = sin(2*pi*0.12*n); y = zeros(1, 3*length(x));
y([1: 3: length(y)]) = x; %3倍内插 y1=x1([1 : 3 : length(x1)]); %3倍抽取 subplot(3,1,1) stem(n,x); title(' 输入序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(3,1,2) stem(n,y(1:length(x))); title(' 内插序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(3,1,3) stem(n,y1); title(' 抽取序列' );
xlabel(' 时间 n' );ylabel(' 幅度' );
程序运行结果如下:
输入序列
幅度
时间 n内插序列
幅度
时间 n幅度
时间 n
图9-20 对正弦信号的时域3倍内插和3倍抽取
[例9-2] 编写程序对一个时域有限长且频域有限带宽的输入序列进行整数倍内插与抽取,分析原输入序列与内插和抽取后序列的频谱变化。
% 程序 9_2
% 分析内插对频谱的影响 freq = [0 0.45 0.5 1]; mag = [0 1 0 0];
x = fir2(99, freq, mag); % 利用fir2函数产生一个有限长序列 % 求取并画入输出谱 [Xz, w] = freqz(x, 1, 512); subplot(2,1,1);
plot(w/pi, abs(Xz)); axis([0 1 0 1]); grid xlabel('\omega/ \pi'); ylabel(' 幅度' ); title(' 输入谱' ); subplot(2,1,2); %产生内插序列
L = input(' 输入内插因子 = '); y = zeros(1, L*length(x)); y([1: L: length(y)]) = x; % 求取并画出输出谱 [Yz, w] = freqz(y, 1, 512);
plot(w/pi, abs(Yz)); axis([0 1 0 1]); grid xlabel('\omega/ \pi'); ylabel(' 幅度' ); title(' 输出谱' );
程序运行结果如下:
输入谱
1
幅度
0.5
00.10.20.30.4
0.50.6ω/ π输出谱
0.70.80.91
1
幅度
0.5
00.10.20.30.4
0.5ω/ π
0.60.70.80.91
图9-21 原序列与3倍内插序列的频谱
9.5.2 抽取器与内插器的设计
[例9-3 用抽取函数decimate 、内插函数interp 和重采样函数resample 编写程序对时域信号
x (n ) =sin(2π⋅0.043n ) +sin(2π⋅0.031n )
分别进行M 倍的抽取、L 倍的内插和L/M倍的分数率抽样率改变操作。
程序代码如下:
% 程序 9_3
% 用抽取函数decimate 、内插函数interp 和重采样函数resample 对信号进行重采样 clf;
M = input(' 抽取因子 = '); L = input(' 内插因子 = '); n = 0:99;
x = sin(2*pi*0.043*n) + sin(2*pi*0.031*n); y = decimate(x,M,'fir' ); y1 = interp(x,L); y2 = resample(x,L,M); subplot(2,2,1); stem(n,x(1:100)); title(' 输入序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,2); m = 0:(100/M)-1; stem(m,y(1:100/M)); title(' 抽取序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,3); m = 0:(100*L)-1; stem(m,y1(1:100*L)); title(' 内插序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,4); m = 0:(100*L/M)-1; stem(m,y2(1:100*L/M)); title(' 分数率抽样序列' ); xlabel(' 时间 n' );ylabel(' 幅度' );
程序运行结果如下:
输入序列
幅度
抽取序列
幅度
时间 n
时间 n
内插序列
幅度
分数率抽样序列
幅度
时间 n
时间 n
图9-22 原序列与2倍抽取、3倍内插、3/2分数倍抽样率得到的序列
9.5.3 滤波器组的设计 对于图9-14中的分析滤波器组,有
∞
∞
k
H k (z ) =
∑h
n =0
(n ) z
-n
=∑h 0(n )(zW M )
n =0
k -n
=H 0(zW M ), k =0, , M -1 (9-5-1)
k
-j 2π/M
式中W M =e 。这些传输函数的频率响应为
H k (e
j ω
) =H 0(e
j (ω-
2πk M
)
), k =0, , M -1 (9-5-2)
即H k (z ) 的频率响应可通过将H 0(z ) 的响应以步长2πk /M 右移得到。式(9-5-1)定义的M 个滤波器H k (z ) 可以作为图9-14中的分析滤波器组中的分析滤波器用,也可在图9-14中的综合滤波器组中作为合成滤波器使用。得到的滤波器称为均匀滤波器组。
[例9-4] 设原型低通滤波器的归一化通带边界为0.2,归一化阻带边界为0.25,试编写程序设计4频带均匀分析或合成滤波器组。
程序代码如下:
% 程序 9_4
% 设计均匀DFT 滤波器组 clf;
b = remez(20, [0 0.2 0.25 1], [1 1 0 0], [10 ]); %用Remez 算法设计原型FIR 低通滤波器 w = 0:2*pi/255:2*pi; n = 0:20; for k = 1:4;
c = exp(2*pi*(k-1)*n*i/4); FB = b.*c;
HB(k,:) = freqz(FB,1,w); end
%画出各子带滤波器的幅频响应¦ subplot(2,2,1) plot(w/pi, abs(HB(1,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 1'); axis([0 2 0 1.1]); subplot(2,2,2) plot(w/pi,abs(HB(2,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 2');axis([0 2 0 1.1]); subplot(2,2,3) plot(w/pi,abs(HB(3,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 3'); axis([0 2 0 1.1]); subplot(2,2,4) plot(w/pi,abs(HB(4,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 4'); axis([0 2 0 1.1]);
程序运行结果如下:
滤波器 No. 1
滤波器 No. 2
度
幅度
幅ω/ π
ω/ π
滤波器 No. 3
滤波器 No. 4
度
幅度
幅0
0.5
11.5
2
00.5
11.52
ω/ π
ω/ π
图9-23 4频带滤波器组
9.5.4 离散小波变换
[例9-5] 试利用Daubechies 系列中的db2小波和函数dwt 对某一维信号进行小波分析。 程序代码如下:
%产生一个基本信号序列,长度=16
s = 2+kron(ones(1,8), [1 -1]) +((1:16).^2)/32+0.2*randn(1,16); % 利用小波db2进行DWT [ca1,cd1] = dwt(s, 'db2'); subplot(211);plot(s); title('原始信号');
subplot(223);stem(ca1,'k.'); title('近似系数');
subplot(224);stem(cd1,'k.'); title('细节系数');
程序运行结果如下:
原始信号
024
近似系数
681012细节系数
1416
0510
0510
图9-24 信号DWT 波形
[例9-6] 试利用Daubechies 系列中的db2小波和函数idwt 对某一维信号进行离散小波逆变换,并计算重构误差。
程序代码如下:
%产生一个基本信号序列,长度=16
s = 2+kron(ones(1,8), [1 -1]) +((1:16).^2)/32+0.2*randn(1,16); % 利用小波db2进行DWT [cA1,cD1] = dwt(s, 'db2');
% 根据dwt 所得尺度展开系数ca1,重构信号低频分量 a1 = idwt(cA1, [], ’db2’); subplot(221);plot(a1); title(‘低频分量’);
% 根据dwt 所得尺度展开系数cd1,重构信号高频分量 d1 = idwt([], cD1, ’db2’); subplot(222);plot(d1); title(‘高频分量’);
% 根据信号分量a1和d1,重构信号ss ss = a1 + d1; Err = norm(s-ss); subplot(212); plot([s;ss];);
title('原始信号与重构信号'); xlabel ([‘重构误差=’, num2str(Err)]);
程序运行结果如下:
低频分量
5
10
15
20
高频分量
051015
20
原始信号与重构信号
重构误差=1.2632e-012
图9-24 信号经IDWT 的重构信号波形
小 结
本章对多采样率信号处理的基本概念与理论知识作了简要介绍。首先介绍了改变离散信号采样率的两种基本操作-抽取与内插,分别从时域与变换域对这两种操作进行了分析。然后讨论了抽取滤波器与内插滤波器,以及它们的FIR 结构和高效的多相滤波器结构。第三讨论了正交镜像滤波器组的结构原理和设计要求。最后扼要介绍了离散小波变换的基本概念,并讨论了它与树状结构滤波器组间的联系。本章所讨论的内容,可以为进一步深入学习软件无线电、数字信道化接收机等理论知识打好基础。
习 题
1. 图P9-1所示系统输入为x (n ) ,输出为y (n ) ,内插模块在输入序列每个样本值间插
入2个零,抽取模块定义为y (n ) =w (5n ) ,其中w (n ) 是抽取模块的输入系列。若输入
sin(ω1n )
x (n ) =
πn 试确定下列ω值时的输出y (n ) :
(1)ω1≤
35
35
(2)ω1>π;
π
图P9-1
2. 如果抽取器的抽取因子M =2k ,k 为整数,可采用抽取因子为2的抽取器级联而成,试给出这一抽取器的框图。假定要求0~π/(2M ) 范围内的信号分量在抽取后无混叠失真,试给出抽取滤波器的技术指标。
3. 已知有理数I/D作抽样率转换的两个系数,如图P9-2所示。 (1)写出X ID 1(z ) 、X ID 2(z ) 、X ID 1(e
j ω
) 、X ID 2(e
j ω
) 的表示式;
(2)若
I =D ,试分析这两个系统是否有x ID 1(n ) =x ID 2(n ) ,请说明理由; (3)若
I ≠D ,请说明在什么条件下x ID 1(n ) =x ID 2(n ) ,并说明理由。
x (n ) x (n )
图P9-2
ID 1(n ) ID 2(n )
4. 设计一采样率转换系统,对采样频率为20kHz 的信号降低采样率以减少数据量,要求仅保留4.1-4.9kHz 的频率分量,在此频带内的频谱失真不大于1dB ,频谱混叠不大于1%。
(1)确定满足要求的最低采样频率和相应的采样率转换因子。
(2)确定采样率转换系统中的带通滤波器的技术指标,并设计相应的带通滤波器。
-1
5. 对于图P9-3所示的分析和综合滤波器组,已知H 0(z ) =(1+z ) /2,求图中的
H 1(z ) 、G 0(z ) 和G 1(z ) ,使得输出Y (z ) 可以完美地重建X (z ) 。
X (z Y (z )
图P9-3
第9章 多采样率信号处理
多采样率信号处理广泛应用于要求转换采样率,或要求系统工作在多采样率状态的信号处理系统中。如多种媒体——语音、视频、数据的传输,它们的频率很不相同,采样率自然不同,必须实行采样率的转换;又如信号要在两个时钟频率的数字系统中传输时,为了便于信号的处理、编码、传输和存储,要求根据时钟频率对信号的采样率加以转换;再如一种信号处理算法在系统的不同部分采用不同的采样率(如子带编码等),使处理更加有效,等等。
本章首先介绍直接在数字域对离散时间信号进行采样率转换的抽取(Decimation )和内插(Interpolation )方法,然后讨论抽取滤波器与内插滤波器的设计与实现方法。由不同采样率构成的系统称为多采样率系统,大部分多采样率系统使用了滤波器组,以正交镜像滤波器组(QMF )为基础的树状结构滤波器组是一典型的多采样率系统,它与离散小波(Wavelet )变换的关系密切。本章最后由正交镜像滤波器组的概念引入小波变换的基本原理和多分辨率分析的概念,以利开拓思路,为进一步深入学习打下基础。
9.1 离散信号的抽取与内插
9.1.1 抽取与内插的时域描述
离散序列的抽取与内插是多采样率系统中的基本运算,抽取运算将降低信号的采样频率,内插运算将提高信号的采样频率。离散序列x (k ) 的M 倍抽取定义为
x D (k ) =x (M k ), k ∈ (9-1-1)
其中M 为一正整数。抽取运算的框图如图9-1所示。
x (k ) D (k )
图9-1 M倍抽取运算的框图
图9-2画出了M =3时序列抽取的示意图。由图可知,离散序列的抽取表示保留第M 个样本点,而去除两个样本之间的M-1个样本点,设原离散信号x (k ) 的采样周期为T ,经M 倍抽取后的信号x D (k ) 的采样周期为T ' ,满足T ' =MT 。为了强调此概念,在图9-2中,
'
有意将抽取后的序列的间隔画为原序列的3倍。这时新的采样频率f s 为
f s =
'
1T
'
=
1MT
=
f s M
(9-1-2)
式中,f s 为原有的采样频率。
(a )原序列 (b) 3倍抽取后的序列
图9-2 离散序列的抽取
离散序列x (k ) 的L 倍内插定义为
⎧x
(k /L ), k =0, ±L , ±2L ,
x I (k ) =⎨ (9-1-3)
其他⎩0 ,
其中L 为一正整数。内插运算的框图如图9-3所示
x (k ) 图9-3 整数L 倍内插运算框图
I (k )
图9-4画出L =3时序列内插运算的示意图。由图可知,原序列x (k ) 中的所有样本都保留在内插后的序列x I (k
) 中,即内插运算不丢失信息。离散序列的内插运算是在原序列x (k ) 的每两个样本点之间插入L -1个零值样本点。
若对内插后的序列通过一个低通滤波器进行平滑处理,则可将序列中的零值转换为内插值,使得低通滤波后的序列的采样率是原序列的L 倍。如果序列x (k ) 的抽样间隔为T ,则低通滤波后序列的采样间隔为T ' =T /L 。
(a )原序列 (b )3倍内插后的序列
图9-4 离散序列的内插
9.1.2 抽取与内插的变换域描述
以上分析了序列抽取与内插的时域特性,下面从频域和z 域讨论抽取与内插的变换域特
性。由z 变换的定义,M 倍抽取后的序列x D (k ) 的z 变换X D (z ) 为
∞
∞
X D (z ) =
∑
k =-∞
∞
x D (k ) z x (n ) z
-k
=
∑
k =-∞
∞
x (M k ) z
-k
=
∑
n =-∞n =kM
-
n M
=
∑
n =-∞
x (n ) δ M (n ) z
-
n M
(9-1-4)
~
式中δ M (n ) 是周期为M 的单位脉冲序列。由离散傅里叶级数可知,δM (n ) 可以表示为
1
δ M (n ) =
M
M -1
∑e
l =0
j 2πnl /M
(9-1-5)
将(9-1-5)式代入(9-1-4) 式可得
X D (z ) =
1M
∞
M -1
∑∑
n =-∞l =0
x (n ) e
j 2πnl /M
z
-
n M
=
1M
M -1
∑
l =0
X (e
-j
2πM
l
1
z M ) (9-1-6)
将z =e j ω代入式(9-1-6)可得M 倍抽取后序列x D (k ) 的频谱为
X D (e
j ω
) =
1M
M -1
∑
l =0
X (e
j
ω-2πl
M
) (9-1-7)
式(9-1-7)表明,M 倍抽取后序列x D (k ) 的频谱可由下列步骤获得: (1)将X (e j ω) 扩展M 倍得到X (e j ω/M ) ,X (e j ω/M ) 的周期为2πM 。 (2)将X (e j ω/M ) 右移2π的整数倍得到{X (e j (ω-2πl ) /M ); l =0,1, , M -1}。 (3)将(2)中的M 个周期为2πM 的函数相加并乘以因子1/M ,得到周期为2π的M 倍抽取后序列的频谱X D (e
j ω
) 。
图9-3画出了3倍抽取后序列的频谱。由图可知,抽取后序列的频谱没有混叠。一般地,如果低频信号x (n ) 的频谱是带限,即在区间[-π, π]范围内有
X (e
j ω
,|ω|>π/M (9-1-8) ) =0
则M 倍抽取后信号的频谱不会发生混叠。式(9-1-8)称为序列M 倍抽取不混叠的Nyquist 条件。
图9-3 原序列的频谱和3倍抽取后序列的频谱
由z 变换的定义,L 倍内插后的序列x I (n ) 的z 变换X I (z ) 表示为
∞
∞
∞
X I (z ) =
∑
n =-∞
x I n (z )
-n
=
∑
n =-∞
n =m , L ∈ m
x n (L /z ) =
-n
∑
m =-∞
x m z (
-m L
) =X z
L
( ) (9-1-9)
将z =e j ω代入式(9-1-9)可得L 倍内插后的序列x I (n ) 的频谱为
X I (e
j ω
) =X (e
j ωL
) (9-1-10)
式(9-1-9)表明L 倍内插后的序列x I (n ) 的频谱X I (e j ω) 是原序列频谱X (e j ω) 的L 倍压缩。
图9-6 原序列的频谱和4倍内插后序列的频谱
图9-6画出了L =4时,原序列的频谱和内插后序列的频谱。由图可知,内插后序列在区间
[-π
X (e
/4, π/4]内的频谱,是由原信号在区间[-π, π]的频谱压缩4倍得到的。除了相差一个
) 的周期为π/2。内插序列在区间[-π/4, π/4]将在区间[-π, -π/4]和[π/4, π]
尺度因子外,两个频谱的形状保持不变。由于原序列频谱X (e j ω) 的周期为2π,因而
j 4
ω
内重复3次,这些重复的部分称为镜像频谱。
9.2抽取滤波器与内插滤波器
9.2.1 抽取滤波器
通常抽取后的离散序列的频谱将会出现混叠。为了避免混叠,可在信号抽取前利用低通滤波器对信号进行滤波,如图9-7所示。该滤波器称为抽取滤波器(decimation filter)。
x (n ) y (n )
图9-7 M 倍抽取滤波系统
由图9-3可知,抽取滤波器可以是截频为π/M 的理想低通滤波器,其频率响应为:
⎧1⎪⎪j ω
H (e ) =⎨
⎪0⎪⎩
ω|≤|πM
πM
(9-2-1)
如果低通滤波器的单位冲激响应为h (n ) ,则滤波器的输出w (n ) 为
∞
w (n )=
最后的输出y (m ) 为
∑h (k ) x (n -k ) (9-2-2)
k =-∞
y (m ) =w (mM ) (9-2-3)
将式(9-2-2)和式(9-2-3)合起来,y (m)与x (n ) 之间的关系为
∞
∞
y (m )=
∑
k =-∞
h (k x ) M (m -k =)
∑
k =-∞
h M (m -k x k ) ( (9-2-4)
式(9-2-4)表明,在计算M 倍抽取滤波器的输出时,只需计算抽取滤波器每M 个输出中的一个样本。利用
W (z ) =H (z ) X (z ) (9-2-5)
并根据(9-1-6)式,有 Y (z ) =
1M
M -1
1
-j 2πl /M
∑W (e
l =0
z M ) =
1M
M -1
1
-j 2πl /M
1
-j 2πl /M
∑H (e
l =0
z M ) X (e z M ) (9-2-6)
z 在单位圆上取值,即z =e j ω时,可得
Y (e
j ω
) =
1M
M -1
∑
l =0
H (e
j
ω-2πl
M
) X (e
j
ω-2πl
M
) (9-2-7)
式(9-2-7)是对输入信号x (n ) 进行滤波、抽取后的频域表达式。
9.2.2 内插滤波器
信号的内插不会引起频谱的混叠,但会产生镜像频谱,如图9-6所示。为了消除这些镜像频谱,可将内插后的信号通过低通滤波器,如图9-8所示。该低通滤波器称为内插滤波器(interpolation fliter )。由图9-6可知,内插滤波器可以是截止频率为π/L 的理想低通滤波器。该滤波器可以滤除信号w (n ) 频谱中的镜像频谱,仅保留[-π/L , π/L ]范围内的频谱。
内插滤波器的频率响应为
π⎧
G ω ≤ |⎪j ω
H (e ) =⎨L (9-2-8)
⎪0 其他⎩
x (n )
y (n )
图9-8 L 倍内插滤波系统
根据图9-8可知
j ωj ω
Y (e j ω) =H (e ) W (e =)
H (e
j ω
(9-2-9) ) X (e
j ωL
根据式(9-2-8)和式(9-2-9),可近似得出
π⎧j ωL
G X (e ) ω ≤ ||⎪j ω
Y (e ) ≈⎨L (9-2-10)
⎪0 其他⎩
在n =0时刻,有
y (0)=
12πG 2π
⎰π
-
π
Y (e
j ω'
) d ω' =) d ω
=
G 2π
⎰π
-
π/L
/L
X (e
j ω' L
) d ω'
=
⎰πX (e
-
π
j ω
(9-2-11)
G x (0)L
L
因此,如要求y (0)=x (0),则应有G =L ,即对理想的内插器要求能恢复内插前的信号,增益G 必须等于L 。
若设内插滤波器的单位冲激响应为h (n ) ,则
∞
y (n )=
∑
k =-∞
h (n - (9-2-12) k ) w ( k )
由式(9-1-3)和式(9-2-12)式,可以得到
∞
y (n )=
∑
k =-∞
h (n -k ) x (
k L
∞
) =
∑
r =-∞
h (n -rL ) x (r ), k 是L 的整数倍 (9-2-13)
由于内插滤波器的输入信号w (n ) 中每L 个样本中只有一个非0样本,所以内插滤波器的计算量只有常规系统的1/L 。
9.2.3 有理数倍抽样率转换
给定信号x (n ) ,若希望将抽样率转变为L /M 倍,可以通过把M 倍抽取和L 倍内插结合起来得到。一般是先做L 倍的插值,再做M 倍的抽取。这是因为先抽取会使x (n ) 的数据点减少,会产生信息的丢失,并且可能产生频率响应的混叠失真。例如:如果x (n ) 的抽样频率f s =2f h ,此时x (n ) 的基带正好在容许的频带上限之内,即在折叠频率|ω|≤π以内,现在要将x (n ) 的抽样频率转换为3f s /2,此时M =2,L =3,如果先做2倍抽取,则会先丢失掉一些数据,而且信号的数字频带要增加2倍,必然产生混叠失真。为了不产生混叠失真,
必须将防混叠的低通滤波器频带限制在|ω|≤π/2内,这样会丢失很多信息。如果先做3倍插值,使数字频带先缩小3倍,再做2倍抽取,则信号的数字频带变成原信号数字频带的2/3,因而不会产生混叠失真。所以,应先对信号做L 倍的插值,再做M 倍抽取,结构上就是两者的级联。图9-9实现以有理数L /M 来改变采样率的系统,新系统输出信号的采样率为f s ' =Lf s /M 。
图9-9 有理数倍L/M采样率转换框图
图9-9中将内插后的抗镜像滤波器和抽取前的抗混叠滤波器合并为一个数字低通滤波器h (n ) 。由于此滤波器同时用作插值和抽取的运算,因而,它的理想频率响应为 H (e j ω) =⎨
≤|⎧⎪L , ω⎪⎩0, 其他
|{πm i L n π
M /},
(9-2-14)
/
9.2.4 抽取滤波器的FIR 结构和多相结构
从前面的讨论可以知道,抽取器需要一个数字低通滤波器用于滤除可能会引起混叠失真的频谱分量。FIR 滤波器具有绝对稳定,容易实现线性相位特性,特别是连同抽取器一起,采用合理的结构,可以大大提高运算效率。
对于图9-7的M 倍抽取器的直接型FIR 滤波器的实现结构如图9-10(a )所示。这个直接型FIR 滤波器结构概念清楚明白,其实现也很简单。但这是h (n ) 工作在高采样率状态,
x (n ) 的每一个采样点均要和FIR 滤波器的系数相乘。而滤波器输出中,每M 个样值中只抽
取一个作为最终的输出y (m ) ,丢弃了其中M -1个样值,所以该结构效率很低。 为了提高直接型FIR 滤波器结构的运算效率,将图9-10(a )中抽取操作嵌入到FIR 滤波器结构中,如图9-10(b )所示,由级联的延迟器z -1移入各抽头的x (n ) 先做抽取,再和
h (n ), n =0,1, , N -1相乘,由于工作在低采样率状态,系统的运算速率降低了M 倍。FIR
滤波器则可写成
N -1
y (m ) =
∑h (k ) x (Mm
k =0
-k ) (9-2-15)
图9-10 抽取的FIR 结构
(a )FIR 滤波器与抽取器级联 (b )提高运算效率的结构
多相(Polyphase )滤波器结构是按整数因子抽取的另一种高效实现结构。通常取N 是M 的整数倍,若令式(9-2-15)中的k =Mq +i , i =0,1, , M -1, q =0,1, , N /M -1,则
N -1
y (m ) =
∑h (k ) x (M m -k )
k =0
N M -1M
-1
=
i =0
(9-2-16)
h (M q +i ) x [M (m -q ) -i ]
∑∑
q =0
令h i (m ) =h (mM +i ) ,i =0,1, , M -1, m =0,1, , N /M -1为多相滤波器的子滤波器的单位脉冲响应,如h 0(m ) ={h (0),h (M ), h (2M ), , h (N -M )};并且
x i (m ) =x (mM -i ) ,i =0,1, , M -1
则式(9-2-16)可以写成M 个子载波的和的形式
N M -1M
-1
i
i
M -1
i
y (m ) =
N M -1
i
i
∑∑h (k ) x (m -k ) =∑y (m ) (9-2-17)
i =0k =0
i =0
式中,y i (m ) =示。
∑h (k ) x (m -k ) 为子滤波器输出。抽取器的多相滤波器结构如图9-11所
k =0
)
图9-11 抽取器的多相滤波器结构
9.2.5 内插滤波器的FIR 结构和多相结构
根据图9-8得到按整数因子L 内插系统的直接型FIR 滤波器结构如图9-12(a )所示,其中FIR 滤波器采用了转置型结构。该结构中x (n ) 插零后再进行滤波,使得FIR 滤波器要和大量的零值相乘,FIR 滤波器以高采样率运行,该结构效率很低。
类似抽取滤波器的FIR 结构,将图9-12(a )内插器嵌入到FIR 滤波器结构中的N 个乘法器之后,得到如图9-12(b )所示的结构,h (n ) 以低的运算速率与x (n ) 相乘后再插零,运算速率降低了L 倍。
图9-12 内插的FIR 结构
(a) 内插器与FIR 滤波器级联 (b) 提高运算效率的结构
仿照抽取器的多相滤波器结构,取内插因子为L 的低通滤波器h (n ) 的长度N 为L 的整数倍,则可以分解成L 个子滤波器
h i (n ) =h (nL +i ) ,i =0,1, , L -1, n =0,1, , N /L -1 (9-2-18)
由图9-12(b )得到图9-13形式的内插器的多相滤波器结构,图中
N L
-1
i
y i (n ) =
∑h (k ) x (n -k ) ,
k =0
i =0,1, , L -1 (9-2-19)
)
图9-13 内插器的多相滤波器结构
9.3 正交镜像滤波器组
9.3.1 数字滤波器组的一般概念和定义
数字滤波器组在语音分析、数据压缩、信号传输等领域均要用到。M 通道滤波器组的基本结构如图9-14所示。H 0(z ) 、H 1(z ) 、…H M -1(z ) 组成了分析滤波器组,输入信号x (n ) 通过这一组滤波器后,得到的x 0(m ) 、x 1(m ) 、…x M -1(m ) 是x (n ) 的子带信号;而合成过程
ˆ(n ) 表示该综合滤波器的输出信号。一般来说,各路则在滤波器组的综合滤波器中完成,x
信号的采样率与输入或输出信号的采样率是不同的,因此,系统本质上是多采样率的。设计滤波器组的一个重要任务或基本准则是,综合应用分析滤波器组H 0(z ) 、H 1(z ) 、…H M -1(z ) 和综合滤波器组G 0(z ) 、G 1(z ) 、…G M -1(z ) ,抵消或尽可能地抑制混叠失真。
x 0(m x 1(m 9.3.2 正交镜像滤波器组
M -1(m
图9-14 M 通道滤波器组对信号分解和综合示意图
如果M =2,图9.14的M 通道滤波器组就简化为两通道滤波器组,如图9-15(a )所示。利用前述信号抽取与内插理论,可得到两通道滤波器组的输入与输出关系为
ˆ(z ) =1{H (z ) G (z ) +H (z ) G (z ) }X (z ) X 0011
2
+12
(9-3-1)
{H 0(-z ) G 0(z ) +H 1(-z ) G 1(z ) }X (-z )
式中X (-z ) 分量是由抽取过程中的混叠产生的响应,称之为混叠项。式(9-3-1)表明,为了消除混叠项对输出的影响,系统需满足
H 0(-z ) G 0(z ) +H 1(-z ) G 1(z ) =0 (9-3-2)
为满足式(9-3-2)约束条件,存在多种选择方案。一种简单的选择方案为
G 0(z ) =H 1(-z ), G 1(z ) =-H 0(-z ) (9-3-3)
这样即可消除混叠失真对整个系统输出响应的影响。滤波器组输出项中无混叠项的滤波器组称为无混叠滤波器组。
) )
(a )
2
(b )
图 9-15 正交镜像滤波器组
(a )结构原理 (b )镜像对称的幅频响应
对于两通道无混叠滤波器组,若分析滤波器组满足
H 1(z ) =H 0(-z ) (9-3-4)
则称该滤波器组为π/2镜像滤波器组。当h 0(k ) 为实系数时,式(9-3-4)意味着
H 1(e j ω) =H 0(e j (π-ω) ) (9-3-5)
j ωj ω
即它们的幅度响应关于π/2镜像对称,如图9-15(b )所示。如果H 0(e ) 和H 1(e ) 两者j ωj ωj ω
没有重合,即当π/2≤ω≤π时,H 0(e ) ≡0,则H 0(e ) 和H 1(e ) 是正交的,称为正
j ωj ω
交镜像滤波器组(QMFB ,Quadrature Mirror Filter Bank)。在实际应用中,H 0(e ) 和H 1(e )
可以有少量重叠,也称为QMFB 。
令X (e j ω) 、X 0(e j ω) 、X 1(e j ω) 、H 0(e j ω) 和H 1(e j ω) 分别是x (n ) 、x 0(m ) 、x 1(m ) 、
h 0(n ) 和h 1(n ) 的DTFT ,可以得到
X 0(e
j ω
) =
1212
{X (e
j ω/2
) H 0(e
j ω/2
) +X [e
j (ω+2π) /2
]H 0[e
j (ω+2π) /2
]} (9-3-6)
及
X 1(e j ω) =
{X (e
j ω/2
) H 1(e
j ω/2
) +X [e
j (ω+2π) /2
]H 1[e
j (ω+2π) /2
]} (9-3-7)
j ωˆˆ0(m ) 、x ˆ1(m ) 、g 0(n ) 和类似的,令X 0(e ) 、X 1(e j ω) 、G 0(e j ω) 和G 1(e j ω) 分别是x g 1(n ) 的DTFT ,得到综合关系式为
ˆ(e j ω) =X ˆ(e j 2ω) G (e j ω) +X ˆ(e j 2ω) G (e j ω) (9-3-8) X 0011
ˆ0(m ) =x 0(m ) ,x ˆ1(m ) =x 1(m ) ,可以得到滤波 用分析器的输出作为综合器的输入,即x
器组的输入输出频域关系式为
ˆ(e j ω) =1[H (e j ω) G (e j ω) +H (e j ω) G (e j ω)]X (e j ω) +X 0011
2 {H 0[e
2
1
j (ω+π)
(9-3-9)
j (ω+π)
]G 0(e
j ω
) +H 1[e
j (ω+π)
]G 1(e
j ω
)}X [e ]
ˆ(e j ω) 的有用信号变换,第二项则代表不希望出现的这个式中的第一项代表从X (e j ω) 到X
频域混叠分量。
为了去除不希望出现的频域混叠分量,必须满足
H 0[e
j (ω+π)
]G 0(e
j ω
) +H 1[e
j (ω+π)
]G 1(e
j ω
) =0 (9-3-10)
只要仔细选择滤波器,将这个关系式中的前后两项抵消掉,就可以实现这个条件。
通常可先设计一公共低通滤波器h (n ) ,再由它得到所有的分析和综合滤波器,即
h 0(n ) =h (n ) (9-3-11) h 1(n ) =(-1) h (n ) (9-3-12)
等效的,它们的DTFT 满足
j ωj ω
H 0(e ) =H (e ) (9-3-13)
n
H 1(e
j ω
) =H [e
j ω
j (ω+π)
] (9-3-14)
上面的条件表明,滤波器H 0(e (b )所示。
j ω
) 和H 1(e
) 对于频率ω=π/2是镜像对称的,如图9-15
将式(9-3-13)和式(9-3-14)代入式(9-3-10),得
H [e
j (ω+π)
]G 0(e
j ω
) +H (e
j ω
) G 1(e
j ω
) =0 (9-3-15)
下面说明对G 0(e j ω) 和G 1(e j ω) 的要求。因为G 0(e j ω) 必须是一个低通滤波器,令
G 0(e j ω) =2H (e j ω) (9-3-16)
或等效为
g 0(n ) =2h (n ) (9-3-17)
式中2是与内插滤波器有关的增益因子。将式(9-3-16)代入式(9-3-15)中,得
G 1(e j ω) =-2H [e j (ω+π) ] (9-3-18)
可以看出G 1(e j ω) 必须具有高通性质,其冲激响应为
g 1(n ) =-2(-1) n h (n ) (9-3-19)
因此,图9-15中分析和综合滤波器的设计转换为对公共低通滤波器H (e j ω) 的设计。将式(9-3-13)、式(9-3-14)、式(9-3-16)和式(9-3-18)代入式(9-3-9),得到这种正交镜像滤波器组最终的输入和输出关系
ˆ(e j ω) ={H 2(e j ω) -H 2[e j (ω+π) ]}X (e j ω) (9-3-20) X
这里代表混叠的第二项已被消去,这意味着在分析结构中由抽取引起的混叠分量被综合结构中内插引起的镜像分量精确地抵消掉了。由式(9-3-20)可以看到,在低通滤波器H (e j ω) 满足条件
|H 2(e j ω) -H 2[e j (ω+π) ]|=1 (9-3-21)
时,分析和综合滤波器组构成的系统其增益将变成1。一般希望H (e j ω) 逼近理想低通条件
π⎧
1 0≤ω≤⎪⎪2j ω
(9-3-22) H (e ) =⎨
⎪0 π
QMF 公共低通滤波器的设计可以综合考虑式(9-3-21)和式(9-3-22)的约束条件,采用第7章的FIR 滤波器设计方法或借助计算机辅助设计的最优化方法实现。
9.4 离散小波变换
9.4.1 连续小波变换
由3.5节可知,在用DFT 对信号进行频谱分析时,要提高频率分辨率就必须加大信号时间截取的长度,这意味着信号分析只能得到某一时间段上信号所包含的频率成份,截取长度越大,关于某一频率分量的时间定位精确度就越低,即时间分辨率降低。小波变换将时域
信号变换到时-频域中,克服了DFT 的不足,它对于不同的时间和不同的频率具有不同的精确度,可以同时获得信号时域特性和频域特性。
对一个平方可积函数f (t ) 进行小波变换,就是用一个母小波函数ψ(t ) 的伸缩平移函数族{ψa , b (t ) }来对其进行展开,记为
WT {f (t ), ψa , b (t )}=
f (t ), ψa , b (t ) =
⎰
∞-∞
f (t ) ψa , b (t ) dt (9-4-1)
式中ψa , b (t ) =
(
t -b a
) ,b 为时间平移,a 为时间伸缩, , 表示两个函数的内积。母
小波函数ψ(t ) 必须满足一定的要求,实际上存在着许多母小波函数,因而有许多小波变换,并具有不同的特征。显然,式(9-4-1)中参数a 为一折中参数,在时间精度和频率分辨率之间取得折中。
这样一个连续参数a 和b 的小波变换具有相当大的冗余性,相应的反变换也不是唯一的,因此很少有实用性,下面将重点讨论二进小波(Dyadic wavelet)。
9.4.2 多分辨率分析
Mallat 提出了用多分辨率分析(或逼近)概念定义小波,将以前的各种正交小波基构造统一了起来,并给出了分解和重构算法。
令 表示实数域, 表示整数域,f (t ) ∈L 2( ) 。多分辨率分析建立在以下三个基本假设的基础上:
(1)存在着一簇具有不同分辨率的子空间,满足
V -∞⊂ ⊂V -2⊂V -1⊂V 0⊂V 1⊂V 2⊂ ⊂V ∞=L ( ) (9-4-2)
2
每一个子空间具有不同的基向量,它们给定了不同的时间分辨率。当下标i 增加,空间V i 的时间分辨率也增加。
(2)存在着一个尺度函数ϕ(t ) ,以及它的整数平移
ϕk (t ) =ϕ(t -k ), k ∈ (9-4-2)
构成空间V 0的一组正交基{ϕ
k (t ) },满足ϕk , ϕl =δkl , k , l ∈ 。
(3)伸缩规则性
f (t ) ∈V i f (2t ) ∈V i +1 (9-4-3)
利用伸缩规则性和尺度函数产生了V 0空间正交基的特性,可以立即得到所有空间V i 的基的集合
ϕi , k (t ) =2ϕ(2t -k )
i /2i
i , k ∈ (9-4-4)
函数ϕi , k (t ) 代表空间V
i 的一个正交基。因此有ϕi , k , ϕi , l =δkl , i , k , l ∈ 。上式中的2i /2保证模i , k (t ) 总为1,它独立于下标i 与k 。
由于ϕ(t ) ∈V 0,V 0⊂V 1,并且函数21/2ϕ(2t -n ), n ∈ ,在V 1空间ϕ(t ) 可以表示成函数ϕ(2t -n ) 的线性组合
ϕ(t ) =
式中h 0(n ) 为线性组合系数。
∑h
n
(n ) ϕ(2t -n ) (9-4-5)
利用引入的基,一个信号f (t ) ∈V i 可以表示为
f (t ) =
∑α
m
i
(m ) ϕi , m (t ) (9-4-6)
式中αi (m ) =以展开成
f (t ), ϕi , m (t ) 为展开系数。进一步,f (t ) 也是V i +1空间上的分量,从而也可
f (t ) =
∑α
n
i +1
(n ) ϕi +1, n (t ) (9-4-7)
与函数ϕ(2i t ) 比较,函数ϕ(2i +1t ) 沿着时间轴压缩了2倍,因此,在V i +1空间上的信号的时间分辨率是V i 空间上的2倍。由于子空间上的分辨率以2的幂次增加,因而称为二进尺度函数或二进多分辨率。
9.4.3 二进小波
对每一子空间V i ⊂V i +1,可以将V i +1表示成V i 与它的正交补空间的直和,即
V i +1=V i ⊕W i , i ∈ (9-4-8)
和V i 一样,希望找到一个确定的函数ψ(t ) ∈W 0,使得每个i ∈ ,函数系
ψi , k (t ) =2ψ(2t -k )
i /2i
i , k ∈ (9-4-9)
构成空间W i 的正交基,这里ψ(t ) 是前面介绍的母小波函数。因此,对于一个信号f (t ) ∈W i ,可以写成
f (t ) =
∑β
m
i
(m ) ψi , m (t ) (9-4-10)
式中βi (m ) =f (t ), ψi , m (t ) 。
由于W i ⊂V i +1,信号f (t ) ∈W i 也可以在V i +1空间中由尺度函数ϕi +1, k (t ) 展开。对于
i =k =0的母小波函数,也可以在空间V 1中表示成
ψ(t ) =
∑h (n ) ϕ(2t -n ) (9-4-11)
1n
式中h 1(n ) 为线性组合系数。类似地,V i 空间也可分解成直和的形式V i =V i -1⊕W i -1,V i -1也可进一步分解。根据式(9-4-2),信号空间可以表示成
L ( ) =V j ⊕W j ⊕W j +1⊕ ⊕W -1⊕W 0⊕W 1 (9-4-12)
2
式中下标j 是任意的,表示分解的深度。
由式(9-4-6)、式(9-4-10)和式(9-4-12),对于任意平方可积信号f (t ) ∈L ( ) ,都
2
可以用尺度函数与小波函数线性表示为
∞
f (t ) =
∑α
m
j
(m ) ϕj , m (t ) +
∑∑β
i =j
m
i
(m ) ψi , m (t ) (9-4-13)
式中∑αj (m ) ϕj , m (t ) 对应f (t ) 的低频分量,是信号f (t ) 的最粗略近似,由低分辨率的尺
m
∞
度信号ϕj , m (t ) 表示;∑
i =j
∑β
m
i
(m ) ψi , m (t ) 对应于信号f (t ) 的高频分量,由具有更高分辨率
的一系列不同尺度的小波信号ψi , m (t ) 表示,随着尺度j 是的增加,表示了信号f (t ) 中更丰富的细节。也就是说,展开系数αj (m ) 反映了信号f (t ) 中的低频分量的分布情况,而展开系数βi (m ) 反映了信号f (t ) 中的高频分量的分布情况,这些展开系数就是信号f (t ) 的离散小波变换(DWT ,Discrete Wavelet Transform)。
9.4.4 二进小波变换与滤波器组
在信号的小波分析中,实际上并不直接利用小波函数ψ(t ) 和尺度函数ϕ(t ) 计算信号的DWT ,而是利用其对应的小波函数系数h 1(n ) 和尺度函数系数h 0(n ) 。
设f (t ) ∈V i +1空间,有
f (t ) =
∑α
n
i +1
(n ) ϕi +1, n (t ) (9-4-14)
式中系数αi +1(n ) 必须已知。由于V i +1=V i ⊕W i ,这一信号可以唯一表示成两子空间的投影和的形式,这里的投影被展开成各自的基的展开式,即
f (t ) =
∑α
m
i
(m ) ϕi , m (t ) +∑βi (m ) ψi , m (t ) (9-4-15)
m
实际上,投影是由已知系数αi +1(n ) 来计算未知系数αi (m ) 和βi (m ) 。
根据式(9-4-5),V i 空间的基函数ϕi , m (t ) 可以由V i +1空间的基函数ϕi +1, m (t ) 递推得到。将式(9-4-5)代入式(9-4-4),得
ϕi , m (t ) =2ϕ(2t -m ) =2
i /2i i /2
∑h
l
(l ) ϕ(2
i +1
t -2m -l ) (9-4-16)
-1/2
h 0(l ) 代入上式,得 将关系式2m +l =n 和h 0(l ) =2
ϕi , m (t ) =
=
∑h
n
(n -2m ) 2
(i +1) /2
ϕ(2
i +1
t -n )
∑
n
h 0(n -2m ) ϕi +1, n (t )
(9-4-17)
同理,W i 空间中的基函数ψi , m (t ) 可以相应于V i +1空间展开成
ψi , m (t ) =
则展开系数αi (m ) 可写成
∑h (n -2m ) ϕ
1n
i +1, n
(t ) (9-4-18)
αi (m ) =
==
f (t ), ϕi , m (t ) =
f (t ), ∑h 0(n -2m ) ϕi +1, n (t )
n
∑h
n
(n -2m ) f (t ), ϕi +1, n (t ) (n -2m ) αi +1(n )
n =2m
(9-4-19)
∑h
n
=h 0(-n ) *αi +1(n )
式中*表示线性卷积。因此,系数αi (m ) 可通过αi +1(n ) 与序列h 0(-n ) 的卷积,并对其结果2倍抽取得到。类似地,βi (m ) 可以用下式求得。
βi (m ) =h 1(-n ) *αi +1(n )
n =2m
(9-4-20)
这一结果示于图9-16。因此,分解一个信号的系数可以通过单位脉冲响应为h 0(-n ) 和h 1(-n ) 的分析滤波器组求得。将信号f (t ) ∈V i +1投影到其子空间时,它具有两倍于子空间的分辨率。相应的,展开系数αi (m ) 和βi (m ) 是系数αi +1(n ) 的采样率的一半。
αi (n )
i (n )
图 9-16 用分析滤波器计算小波展开系数
联系到伸缩的规则性,V i +1到V i 上的投影相应于低通滤波,V i +1到W i 上的投影相应于高通滤波。因此,h 0(-n ) 是一低通滤波器,而h 1(-n ) 是一互补的高通滤波器。如果重复地将更多的两通道分析滤波器连接到一个两通道滤波器组的低通输出端,将实现式(9-4-13)的信号分解,其结构如图9-17所示。
图 9-17 树状结构的二进小波变换
与上述讨论的将一个信号投影到两个子空间的过程相反,可以将两个子空间V i 和W i 的信号合成为V i +1空间的一个信号,这一过程伴随着信号的时间精度提高了一倍。
这一信号的合成相应于两通道综合滤波器组的运算。在综合运算中,合成信号的未知系数αi +1(n ) ,可由子带信号的已知系数αi (m ) 和βi (m ) 计算得到。为了计算这些系数,需要用基ϕi , m (t ) 和ψi , m (t ) 表示ϕi +1, m (t ) 。因为V 1=V 0⊕W 0,所以存在着以下关系式
ϕ(2t ) =
利用关系式
∑k
(k ) ϕ(t -k ) +∑1(k ) ψ(t -k ) (9-4-21)
k
ϕi +1, n (t ) =2
=
(i +1) /2
ϕ(2
i +1
t -n )
i
i
∑2
k
1/2
0(k ) 2ϕ(2t -
n 2
-k ) +
∑2
k
1/2
1(k ) 2ψ(2t -
i /2i
n 2
-k )
(9-4-21)
令g 0(-2k ) =21/20(k ) ,g 1(-2k ) =21/21(k ) 和n +2k =2m 代入式(9-4-21),得
ϕi +1, n (t ) =
=
∑g
m
(n -2m ) 2
i /2
ϕ(2t -m ) +∑g 1(n -2m ) 2ψ(2t -m )
m
i i /2i
∑g
m
(n -2m ) ϕi , m (t ) +
∑g
m
1
(n -2m ) ψi , m (t )
(9-4-22)
因此,展开系数可写成
αi +1(n ) =
==
f (t ), ϕi +1, n (t ) =
f (t ), ∑g 0(n -2m ) ϕi , m (t ) +
m
∑g
m
1
(n -2m ) ψi , m (t )
∑g
n
(n -2m ) f (t ), ϕi , m (t ) +(n -2m ) αi (m ) +
∑g
n
1
(n -2m ) f (t ), ψi , m (t )
(9-4-23)
∑g
n
∑g
n
1
(n -2m ) βi (m )
式(9-4-23)可以进一步表示为
αi +1(n ) =g 0(n ) *αi () +g 1(n ) *βi () (9-4-24)
2
2
n
n
其中αi () 和βi () 分别为αi (n ) 和βi (n ) 的2倍内插序列。式(9-4-24)可以用图9-18的
2
2
n n
综合滤波器组来实现。由于内插过程,数字滤波器g 0(n ) 和g 1(n ) 的偶数下标才与系数
αi (m ) 和βi (m ) 相乘。式(9-4-24)就是小波重构算法,对应于离散逆小波变换(IDWT: Inverse Discrete Wavelet Transform)。
αi (m βi (m
图 9-18 用综合滤波器组合两个子空间的信号
对于图9-18结构的二进级联,可将不同子空间的信号合成到一起,得到了与图9-17相反的结构,示于图9-19所示。
图 9-19 树状结构的逆二进小波变换
采用二进树结构的分析滤波器组和综合滤波器组对信号进行分解或重构是非常有效的,它甚至不涉及函数ϕ(t ) 和ψ(t ) 的具体形式。
9.5 本章相关的MATLAB 命令及应用
在MA TLAB 信号处理工具箱中,提供了抽取函数decimate ()、内插函数interp ()和重采样函数resample ()等。
1. 抽取函数decimate 格式:y = decimate(x,r)
功能:对离散时间信号向量x 按抽取因子r 抽取,得到信号向量y ,相当于降低了采样率r 倍。向量y 的长度是原信号向量x 长度的1/r倍。
2. 内插函数interp 格式:y = interp(x,r)
功能:对离散时间信号向量x 按内插因子r 内插,得到信号向量y ,相当于提高了采样率r 倍。向量y 的长度是原信号向量x 长度的r 倍。
3. 重采样函数resample 格式:y = resample(x,p,q)
功能:对离散时间信号向量x 按有理数p/q倍的采样率重新采样,得到信号向量y 。p 和q 必须为正整数。向量y 的长度等于ceil(length(x)*p/q)。
4. 离散小波变换dwt
格式:[cA,cD] = dwt(X , ‘wname ’)
功能:对离散时间信号序列X 进行小波分析,cA 是信号离散小波变换对应的近似展开系数,cD 是对应的细节系数,wname 是小波名称。
5. 离散小波逆变换idwt 格式:X = idwt(cA,cD,’wname ’)
功能:实现一维单级离散小波逆变换,是函数dwt 的逆运算。各参数与函数dwt 相同。 6. 分解与重构滤波器组 wfilters
格式:[Lo_D, Hi_D, Lo_R,Hi_R] = wfilters(‘wname ’)
功能:产生用于信号DWT 和IDWT 对应的滤波器组。Lo_D为分解数字低通滤波器的单位脉冲响应,Hi_D为分解数字高通滤波器的单位脉冲响应,Lo_R为重构数字低通滤波器的单位脉冲响应,Hi_R为重构数字高通滤波器的单位脉冲响应。
9.5.1 时域和频域抽样率的改变
[例9-1] 编写程序对时域正弦信号
x (n ) =sin(2π⋅0.12n )
进行3倍内插和3倍抽取操作。 程序代码如下: % 程序9_1
% 对一个正弦信号进行3倍内插和3倍抽取 clf; n = 0:49; m = 0: 50*3 - 1; x1 = sin(2*pi*0.12*m); x = sin(2*pi*0.12*n); y = zeros(1, 3*length(x));
y([1: 3: length(y)]) = x; %3倍内插 y1=x1([1 : 3 : length(x1)]); %3倍抽取 subplot(3,1,1) stem(n,x); title(' 输入序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(3,1,2) stem(n,y(1:length(x))); title(' 内插序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(3,1,3) stem(n,y1); title(' 抽取序列' );
xlabel(' 时间 n' );ylabel(' 幅度' );
程序运行结果如下:
输入序列
幅度
时间 n内插序列
幅度
时间 n幅度
时间 n
图9-20 对正弦信号的时域3倍内插和3倍抽取
[例9-2] 编写程序对一个时域有限长且频域有限带宽的输入序列进行整数倍内插与抽取,分析原输入序列与内插和抽取后序列的频谱变化。
% 程序 9_2
% 分析内插对频谱的影响 freq = [0 0.45 0.5 1]; mag = [0 1 0 0];
x = fir2(99, freq, mag); % 利用fir2函数产生一个有限长序列 % 求取并画入输出谱 [Xz, w] = freqz(x, 1, 512); subplot(2,1,1);
plot(w/pi, abs(Xz)); axis([0 1 0 1]); grid xlabel('\omega/ \pi'); ylabel(' 幅度' ); title(' 输入谱' ); subplot(2,1,2); %产生内插序列
L = input(' 输入内插因子 = '); y = zeros(1, L*length(x)); y([1: L: length(y)]) = x; % 求取并画出输出谱 [Yz, w] = freqz(y, 1, 512);
plot(w/pi, abs(Yz)); axis([0 1 0 1]); grid xlabel('\omega/ \pi'); ylabel(' 幅度' ); title(' 输出谱' );
程序运行结果如下:
输入谱
1
幅度
0.5
00.10.20.30.4
0.50.6ω/ π输出谱
0.70.80.91
1
幅度
0.5
00.10.20.30.4
0.5ω/ π
0.60.70.80.91
图9-21 原序列与3倍内插序列的频谱
9.5.2 抽取器与内插器的设计
[例9-3 用抽取函数decimate 、内插函数interp 和重采样函数resample 编写程序对时域信号
x (n ) =sin(2π⋅0.043n ) +sin(2π⋅0.031n )
分别进行M 倍的抽取、L 倍的内插和L/M倍的分数率抽样率改变操作。
程序代码如下:
% 程序 9_3
% 用抽取函数decimate 、内插函数interp 和重采样函数resample 对信号进行重采样 clf;
M = input(' 抽取因子 = '); L = input(' 内插因子 = '); n = 0:99;
x = sin(2*pi*0.043*n) + sin(2*pi*0.031*n); y = decimate(x,M,'fir' ); y1 = interp(x,L); y2 = resample(x,L,M); subplot(2,2,1); stem(n,x(1:100)); title(' 输入序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,2); m = 0:(100/M)-1; stem(m,y(1:100/M)); title(' 抽取序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,3); m = 0:(100*L)-1; stem(m,y1(1:100*L)); title(' 内插序列' );
xlabel(' 时间 n' );ylabel(' 幅度' ); subplot(2,2,4); m = 0:(100*L/M)-1; stem(m,y2(1:100*L/M)); title(' 分数率抽样序列' ); xlabel(' 时间 n' );ylabel(' 幅度' );
程序运行结果如下:
输入序列
幅度
抽取序列
幅度
时间 n
时间 n
内插序列
幅度
分数率抽样序列
幅度
时间 n
时间 n
图9-22 原序列与2倍抽取、3倍内插、3/2分数倍抽样率得到的序列
9.5.3 滤波器组的设计 对于图9-14中的分析滤波器组,有
∞
∞
k
H k (z ) =
∑h
n =0
(n ) z
-n
=∑h 0(n )(zW M )
n =0
k -n
=H 0(zW M ), k =0, , M -1 (9-5-1)
k
-j 2π/M
式中W M =e 。这些传输函数的频率响应为
H k (e
j ω
) =H 0(e
j (ω-
2πk M
)
), k =0, , M -1 (9-5-2)
即H k (z ) 的频率响应可通过将H 0(z ) 的响应以步长2πk /M 右移得到。式(9-5-1)定义的M 个滤波器H k (z ) 可以作为图9-14中的分析滤波器组中的分析滤波器用,也可在图9-14中的综合滤波器组中作为合成滤波器使用。得到的滤波器称为均匀滤波器组。
[例9-4] 设原型低通滤波器的归一化通带边界为0.2,归一化阻带边界为0.25,试编写程序设计4频带均匀分析或合成滤波器组。
程序代码如下:
% 程序 9_4
% 设计均匀DFT 滤波器组 clf;
b = remez(20, [0 0.2 0.25 1], [1 1 0 0], [10 ]); %用Remez 算法设计原型FIR 低通滤波器 w = 0:2*pi/255:2*pi; n = 0:20; for k = 1:4;
c = exp(2*pi*(k-1)*n*i/4); FB = b.*c;
HB(k,:) = freqz(FB,1,w); end
%画出各子带滤波器的幅频响应¦ subplot(2,2,1) plot(w/pi, abs(HB(1,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 1'); axis([0 2 0 1.1]); subplot(2,2,2) plot(w/pi,abs(HB(2,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 2');axis([0 2 0 1.1]); subplot(2,2,3) plot(w/pi,abs(HB(3,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 3'); axis([0 2 0 1.1]); subplot(2,2,4) plot(w/pi,abs(HB(4,:))); xlabel('\omega/ \pi');ylabel(' 幅度' ); title(' 滤波器No. 4'); axis([0 2 0 1.1]);
程序运行结果如下:
滤波器 No. 1
滤波器 No. 2
度
幅度
幅ω/ π
ω/ π
滤波器 No. 3
滤波器 No. 4
度
幅度
幅0
0.5
11.5
2
00.5
11.52
ω/ π
ω/ π
图9-23 4频带滤波器组
9.5.4 离散小波变换
[例9-5] 试利用Daubechies 系列中的db2小波和函数dwt 对某一维信号进行小波分析。 程序代码如下:
%产生一个基本信号序列,长度=16
s = 2+kron(ones(1,8), [1 -1]) +((1:16).^2)/32+0.2*randn(1,16); % 利用小波db2进行DWT [ca1,cd1] = dwt(s, 'db2'); subplot(211);plot(s); title('原始信号');
subplot(223);stem(ca1,'k.'); title('近似系数');
subplot(224);stem(cd1,'k.'); title('细节系数');
程序运行结果如下:
原始信号
024
近似系数
681012细节系数
1416
0510
0510
图9-24 信号DWT 波形
[例9-6] 试利用Daubechies 系列中的db2小波和函数idwt 对某一维信号进行离散小波逆变换,并计算重构误差。
程序代码如下:
%产生一个基本信号序列,长度=16
s = 2+kron(ones(1,8), [1 -1]) +((1:16).^2)/32+0.2*randn(1,16); % 利用小波db2进行DWT [cA1,cD1] = dwt(s, 'db2');
% 根据dwt 所得尺度展开系数ca1,重构信号低频分量 a1 = idwt(cA1, [], ’db2’); subplot(221);plot(a1); title(‘低频分量’);
% 根据dwt 所得尺度展开系数cd1,重构信号高频分量 d1 = idwt([], cD1, ’db2’); subplot(222);plot(d1); title(‘高频分量’);
% 根据信号分量a1和d1,重构信号ss ss = a1 + d1; Err = norm(s-ss); subplot(212); plot([s;ss];);
title('原始信号与重构信号'); xlabel ([‘重构误差=’, num2str(Err)]);
程序运行结果如下:
低频分量
5
10
15
20
高频分量
051015
20
原始信号与重构信号
重构误差=1.2632e-012
图9-24 信号经IDWT 的重构信号波形
小 结
本章对多采样率信号处理的基本概念与理论知识作了简要介绍。首先介绍了改变离散信号采样率的两种基本操作-抽取与内插,分别从时域与变换域对这两种操作进行了分析。然后讨论了抽取滤波器与内插滤波器,以及它们的FIR 结构和高效的多相滤波器结构。第三讨论了正交镜像滤波器组的结构原理和设计要求。最后扼要介绍了离散小波变换的基本概念,并讨论了它与树状结构滤波器组间的联系。本章所讨论的内容,可以为进一步深入学习软件无线电、数字信道化接收机等理论知识打好基础。
习 题
1. 图P9-1所示系统输入为x (n ) ,输出为y (n ) ,内插模块在输入序列每个样本值间插
入2个零,抽取模块定义为y (n ) =w (5n ) ,其中w (n ) 是抽取模块的输入系列。若输入
sin(ω1n )
x (n ) =
πn 试确定下列ω值时的输出y (n ) :
(1)ω1≤
35
35
(2)ω1>π;
π
图P9-1
2. 如果抽取器的抽取因子M =2k ,k 为整数,可采用抽取因子为2的抽取器级联而成,试给出这一抽取器的框图。假定要求0~π/(2M ) 范围内的信号分量在抽取后无混叠失真,试给出抽取滤波器的技术指标。
3. 已知有理数I/D作抽样率转换的两个系数,如图P9-2所示。 (1)写出X ID 1(z ) 、X ID 2(z ) 、X ID 1(e
j ω
) 、X ID 2(e
j ω
) 的表示式;
(2)若
I =D ,试分析这两个系统是否有x ID 1(n ) =x ID 2(n ) ,请说明理由; (3)若
I ≠D ,请说明在什么条件下x ID 1(n ) =x ID 2(n ) ,并说明理由。
x (n ) x (n )
图P9-2
ID 1(n ) ID 2(n )
4. 设计一采样率转换系统,对采样频率为20kHz 的信号降低采样率以减少数据量,要求仅保留4.1-4.9kHz 的频率分量,在此频带内的频谱失真不大于1dB ,频谱混叠不大于1%。
(1)确定满足要求的最低采样频率和相应的采样率转换因子。
(2)确定采样率转换系统中的带通滤波器的技术指标,并设计相应的带通滤波器。
-1
5. 对于图P9-3所示的分析和综合滤波器组,已知H 0(z ) =(1+z ) /2,求图中的
H 1(z ) 、G 0(z ) 和G 1(z ) ,使得输出Y (z ) 可以完美地重建X (z ) 。
X (z Y (z )
图P9-3