第九章 多采样率数字信号处理
多采样率系统
第3页/共33页
Enjoy Science
数字信号处理
这么一来，x(n)、w(n)和y(m)三个序列的关系就是
y(m) w(Dm) x(Dm) (D是正整数 )
(9.7)
将这个关系应用到公式(9.5)，得到
Y (z) w(Dm)z m (变量代换Dm n和m n / D)
m
w(n)(z1/ D )n (利用公式(9.6)))
D
e
j
2 D
k
)
D k0
(9.13)
第6页/共33页
Enjoy Science
数字信号处理
只要将z=ejω代入上式，就可以得到抽取的频谱关系
Y (e j )
1
D1
X
(e
j /
eD
j 2 D
k
)
1
D1
j 2k
X (e D )
D k0Βιβλιοθήκη D k0(9.14)
借鉴X(ejω)=X(ω)的关系，还能将抽取的频谱关系(9.14) 变为简单的形式
Y ()
1
D1 2k
X(
)
D k0
D
( y(m)的采样率 f y
fx ) D
(9.15)
该式 说明 ：按 照时 序间 隔 D对x(n)抽 取后得 到 序列 y(m)，它的频谱Y(ω)是D个X(ω)变形后相加的结果。
第7页/共33页
Enjoy Science
数字信号处理
（2）从模拟域的角度观看
X
s2
(
)
CTFT
[ xs2
(t)]
1 Ts2
X a (
j
s2 j)
(9.17)
实验七 多采样率数字信号处理
实验七 多采样率数字信号处理一.实验目的1.掌握信号抽取和插值的基本原理和实现2.掌握信号的有理数倍速率转换 二.实验内容例7.1 对信号进行抽取，使采样率为原来的1/4倍。
MATLAB 程序： t=0:.00025:1;x=sin(2*pi*30*t)+sin(2*pi*60*t); y=decimate(x,4);figure,subplot(2,2,1),stem(x(1:120)); title('原始信号时域图'),xlabel('(a)');subplot(2,2,2),plot(abs(fft(x))),title('原始信号频域图'),xlabel('(b)'); subplot(2,2,3),stem(y(1:30))title('抽样后的信号时域图'),xlabel('(c)'); subplot(2,2,4),plot(abs(fft(y)));title('抽样后的信号频域图'),xlabel('(d)');050100150-2-1012原始信号时域图(a)2000400060000100020003000原始信号频域图(b)0102030-2-1012抽样后的信号时域图(c)050010001500200400600抽样后的信号频域图(d)例 7.2 信号()⎪⎪⎭⎫ ⎝⎛=s f f nn x π3cos ，采样频率16/1=sf f，现将采样率提高为原来的4倍。
MATLAB 程序： n=0:30;x=cos(3*pi*n/16); y=interp(x,4);figure,subplot(2,2,1),stem(x);title('原始信号时域图'),xlabel('(a)');subplot(2,2,2),plot(abs(fft(x))),title('原始信号频域图'),xlabel('(b)'); subplot(2,2,3),stem(y(1:30))title('抽样后的信号时域图'),xlabel('(c)');subplot(2,2,4),plot(abs(fft(y)));title('抽样后的信号频域图'),xlabel('(d)');10203040-1-0.500.51原始信号时域图(a)102030405101520原始信号频域图(b)0102030-1-0.500.51抽样后的信号时域图(c)05010015020406080抽样后的信号频域图(d)例7.3 序列()()()2.12.1cos 5.0sin ++=n n n x ππ,调用resample 函数对想按因子3/7进行采样率转换，并绘出图形。
第九章多采样率数字信号处理
图中， 是第i级整数因子Ii内插系统的 镜F 像iIiF i 1,i 1 ,2 ,L,L
返回
2020/6/20 样率仍满足抽样定理要求时，才不会
回到本节
整数因子抽取特点： (1)已抽样序列x(n)和抽取序列y(n)的频
谱差别在频 率尺度上不同。
(2)抽取的效果使原序列的频谱带宽扩展 。
(3)为避免在抽取过程中发生频率响应的 混叠失真，原序列x(n)的频谱就不能 占满频带(0-π).
2020/6/20
x (n )
h (0 )
y(m )
↓D
z 1
h (1 )
z 1
↓D
M
M
M
z 1
h(M 2)
↓D
z 1
h(M 1)
↓D
2020/6/20
返回
回到本节
x(n) h D ( n )
v (n )
D
X (e jw ) H D ( e j ) V ( e jw )
y (n ) Y ( e jw )
M 1
v(n)x(n)hD (n) hD (k)v(nk) k0
返回
9.3整数因子内插
整数因子I内插的目的将原信号采样频
率提高I倍
Fx
1 Tx
采样频率：
xa(t)x(n)
Fx
Fy IFx
整数因子内插：将x(n)的抽样频率
x(n增) 加 I
v(m )
y(m)xa(mTy)
倍，即为II倍插值结果 h I ( m )
y(m)xa(mTy) TyTxI
2020/6/20
I
hI (m )
整数因子I内插系统的直接型FIR滤波
数字信号处理 第九章多采样率数字信号处理
y(0) 1
2
Y (e jy
)dy
1
2
I
I
CX
(e
jIy
)d
y
C
2 I
X
(e jx
)dx
C I
x(0)
C I
*
时域关系
x(n)
v(m)
y(m) xa (mTy )
I
hI (m)
y(m) v(m) hI (m) hI (m k)v(k) k v(kI ) x(k), v(k) 0, k 0, I , 2I ,
I
y
2
*
如何实现 加滤波器
y(m) xa (mTy ) Ty Tx I
H
I
(e
j y
)
C
,
0,
y I I y
y x I
镜像滤波器
Y
(e
j y
)
CV 0,
(e
j y
) I
CX (e jIy
y
),
y I
C=？
x(n)
v(m)
y(m) xa (mTy )
I
hI (m)
*
*
x(n) X (e jw )
hD (n)
v(n)
H D (e j ) V (e jw )
D
y(n) Y (e jw )
H
D
(e
j
)
1,
0,
D D
V (z) Hd (z)X (z)
Y (e jy )
1
D1
j (y 2 k )
V (e D D )
D k0
1 D1
j (y 2 k )
第9章多采样率信号管理方案计划
第9章 多采样率信号处理多采样率信号处理广泛应用于要求转换采样率，或要求系统工作在多采样率状态的信号处理系统中。
如多种媒体——语音、视频、数据的传输，它们的频率很不相同，采样率自然不同，必须实行采样率的转换；又如信号要在两个时钟频率的数字系统中传输时，为了便于信号的处理、编码、传输和存储，要求根据时钟频率对信号的采样率加以转换；再如一种信号处理算法在系统的不同部分采用不同的采样率（如子带编码等），使处理更加有效，等等。
本章首先介绍直接在数字域对离散时间信号进行采样率转换的抽取（Decimation ）和内插（Interpolation ）方法，然后讨论抽取滤波器与内插滤波器的设计与实现方法。
由不同采样率构成的系统称为多采样率系统，大部分多采样率系统使用了滤波器组，以正交镜像滤波器组（QMF ）为基础的树状结构滤波器组是一典型的多采样率系统，它与离散小波（Wavelet ）变换的关系密切。
本章最后由正交镜像滤波器组的概念引入小波变换的基本原理和多分辨率分析的概念，以利开拓思路，为进一步深入学习打下基础。
9．1 离散信号的抽取与内插9.1.1 抽取与内插的时域描述离散序列的抽取与内插是多采样率系统中的基本运算，抽取运算将降低信号的采样频率，内插运算将提高信号的采样频率。
离散序列()x k 的M 倍抽取定义为()(),D x k x Mk k =∈¢ (9-1-1)其中M 为一正整数。
抽取运算的框图如图9-1所示。
()x k ()D k图9-1 Ｍ倍抽取运算的框图图9-2画出了M =3时序列抽取的示意图。
由图可知，离散序列的抽取表示保留第M 个样本点，而去除两个样本之间的M-1个样本点，设原离散信号()x k 的采样周期为T ，经M 倍抽取后的信号()D x k 的采样周期为'T ，满足'T MT =。
为了强调此概念，在图9-2中，有意将抽取后的序列的间隔画为原序列的3倍。
这时新的采样频率's f 为Mf MT T f ss ===11'' (9-1-2) 式中，s f 为原有的采样频率。
DSP的多采样率数字信号处理及其应用
目录1.背景 12．具体过程 22.1 整数因子抽取 22.2 整数因子内插 22.3 I/D的采样率转换 22.4多采样率数字信号处理的应用 23．实验过程 23.1整数倍抽取实验 23.2整数倍插值实验 23.3用有理因子I/D的采样率转换进行的实验 2 4．实验结果 24.1信号的整数倍抽取 24.2信号的整数倍插值 24.3用有理因子I/D的采样速率转换 25．结论 25.1整数因子抽取 25.2整数因子插值 25.3有理因子I/D的采样速率转换 26.心得体会与总结 21.背景现在实际系统中，经常要求一个数字系统能工作在多采样率状态，例如：在数字电视系统中，图像采集系统一般按4:4:4标准或4:2:2标准采集数字电视信号，再根据不同的电视质量要求将其转换成其它标准的数字电视信号（如4:2:2,4:1:1,2:1:1）进行处理。
在数字电话系统中，传输的信号既有语音信号又有传真信号，甚至有视频信号。
这些信号的频域成分相差甚远。
因此该系统应具有多种采样率，并能根据所传输的信号自动完成采样率转换。
对一个非平稳随机信号（如语音信号）做频谱分析或编码时，对不同的信号段可根据其频域成分的不同而采用不同的采样率，已到达既满足采样定理，又最大限度的减少数据量的目的。
如果以高采样率采集的数据存在冗余，这时就希望在该数字信号的基础上降低采样率。
多采样率数字信号处理是建立在单抽样率信号处理基础上的一类信号处理。
在传输信号时，由于语音﹑图像、视频信号的中心频率相差很大，所以需要以多种抽样频率来对信号采样来满足各种传输类型的需要。
2．具体过程2.1 整数因子抽取信号的抽取是实现频率降低的方法。
在第二章曾经讨论过，当采样频率大于信号最高频率的2倍时，不会产生混叠失真。
显然，当采样频率远高于信号最高频率时，采样后的信号就会有冗余数据。
此时，通过信号的抽取来降低采样频率，同样不会产生混叠失真。
Xd(n)整数因子抽取原理图：设x(n)=x(t)|t=nTs，欲使fs减少D倍，最简单的方法就是从x(n)中每D个点中抽取一个，依次组成一个新的序列xd(n),即xd(n)=x(Dn)因为是舍去部分点，故可引入冲激函数来进行抽样，得到xd(n)与x(n)之间的表达式：xd(n)=x(n) D(n)其中为周期单位脉冲序列，当且仅当n为D的整数倍时， D(n)的值为1,n为其他值时为零。
《数字信号处理》第8章 多采样率数字信号处理
－ Ωc 0 Ωc
Ωsa 1
Ω
(a)
～
Λ (k )
－3 －2 －1 0
1
2
3
4
5
k
(b)
X^(e jΩ T1)
0 Ωsa 1/D
Ωsa 1
Ω
(c)
Y(e jΩ T2)
0 Ωsa 2 ＝Ωsa 1/D
Ω
(d)
图 8.1.8 在Ωc<Ωsa2/2时， 抽取前后信号的时域和频域关系示意图
第8章 多采样率数字信号处理
T2 = DT1
第8章 多采样率数字信号处理
其中n1和n2分别表示x(n1T1)和x(n2T2)序列的序号， 于是有
y(n2T2) = x(n2DT1) 当n1=n2D时，
y(n2，T2) = x(n1，T1) 或
y(n) = x(Dn) D倍抽取就是每隔D-1个点抽取一个。
第8章 多采样率数字信号处理
建立在采样率转换基础上的“多采样率数字信号处理” 已成为数字信号处理学科的主要内容之一，在语音信 号处理、图像处理、通信系统等领域有着广泛的应用。
第8章 多采样率数字信号处理
采样率转换通常分为“抽取”和“内插”。 抽取：是降低采样率以去掉多余数据的过程。
内插：是提高采样率以增加数据的过程。 本章重点讨论抽取和内插的概念以及其基本实现方
法。
第8章 多采样率数字信号处理
8.1 信号的整数倍抽取
1. 信号的整数倍抽取的时域描述 设x(n1, T1)是连续信号xa(t)的采样序列，采样率
F1=1/T1(Hz)，T1称为采样间隔，单位为秒，即 x(n1T1) = xa(n1T1)
若将采样率降低到原来的1/D（D为大于1的整数，称为 抽取因子），采样间隔为T2，采样率F2=1/T2(Hz)，组 成的新序列为y(n2, T2)，则有
多采样率信号处理
多采样率信号处理1．绪论随着数字信号处理的发展，信号的处理、编码、传输和存储等工作量越来越大。
为了节省计算工作量及存储空间，在一个信号处理系统中常常需要不同的采样率及其相互转换，在这种需求下，多速率数字信号处理产生并发展起来。
它的应用带来许多好处，例如：可降低计算复杂度、降低传输速率、减少存储量等。
在信号处理领域，多速率信号处理最早于20世纪70年代提出，由其引出的多速率滤波在数学领域里基于多格算法解决了大量的微分等式。
在多速率数字信号处理发展中，一个突破点是70年代两通道正交镜像滤波器组应用于语音信号的压缩。
在该方法中，信号通过分析滤波器组被分成低通和高通两个子带，每个子带经过2倍抽取和量化后再进行压缩，之后可以通过综合滤波器组近似地重建出原始信号，重建的近似误差一部分源于子带信号的压缩编码，一部分是由分析和综合滤波器组产生的误差，其中最主要的误差是混叠误差，它是由分析滤波器组不是理想带限而引起的。
在很多应用系统中，混叠误差存在一定程度的影响，因此就需要对其进行改进。
多速率系统应用于通信、语音信号处理、谱分析、雷达系统和天线系统，以及在数字音频系统、子带编码技术( 用于声音和图像的压缩) 和模拟语音个人系统(如标准电话通信) 等方面的应用。
另外还应用于多相理论和多速率系统在一些非传统领域，包括：高效率信号压缩的多速率理论；高效窄带滤波器的脉冲响应序列的编码新技术的推导；可调整的多级响应FIR滤波器的设计等。
基于上述研究的发展，从20世纪80年代初开始，多速率数字信号处理技术在工程实践中得到广泛的应用，主要用于通信系统、语音、图像压缩、数字音频系统、统计和自适应信号处理、差分方程的数值解等。
多速率信号处理在基础理论和应用领域的蓬勃发展，也促进了整个数字信号处理界的发展。
2．采样率转换基础理论实现采样率转换的方法有三个：一是若原模拟信号x (t)可以再生，或是己记录下来了的话，那么可重新抽样；二是将x(n)通过D/A变成模拟信号x(t)后，对x (t)经A/D再抽样；三是发展一套算法，对抽样后的数字信号x(n)在“数字域”作采样率转换，以得到新的抽样。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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
实际分解和重建滤波器
理想分解重建滤波器属于非因果系统，物理上无法实现