数字信号处理-丁玉美 高西全 编著-第7章

(7.2.16)
图7.2.4中各点信号的时域表示式归纳如下：
第７章
多采样率数字信号处理
线性滤波器输出序列为
l x v(l ) = I 0
线性滤波器输出序列为
l = 0, ± I , ±2 I , ±3I , L
(7.2.17)
其他
w(l ) =
k =−∞
∑ h(l − k )v(k ) = ∑ h(l − kI )x(k )
y (m) = v(m) ∗ hI (m) =
k =− ∞
∑
∞
hI (m − k )v(k )
(7.2.9)
因为除了在I的整数倍点v(kI)=x(k)以外， v(k)=0， 所以
y ( m) =
k =− ∞
∑ h (m − kI )x(k )
I
∞
(7.2.10)
第７章
多采样率数字信号处理
变换域输入输出关系式：
V ( z) =
m =− ∞
∑
∞
v ( m) z − m =
m =− ∞
∑
∞
v( Im) z − Im =
m =− ∞
∑
∞
x(m) z − Im = X ( z I )
(7.2.11)
计算单位圆上的V(z)得到v(m)的频谱为
V (e jω ) = V ( z ) z = e ω = X (e jI ω )
7.2.2 整数因子 内插器 整数因子I内插器 内插器
1. 整数因子 内插器原理框图 整数因子I内插器原理框图 内插器原理框图 按整数因子I对x(n)内插的原理方框图如图7.2.2所示。
图7.2.2
第７章
多采样率数字信号处理
2． 整数因子内插器的功能 ． 整数因子内插器的功能 整数因子内插器的功能表现为： 输出端信号采样频率Fy 提高为输入端信号采样频率Fx的I倍， 即Fy=IFx。 3. 知识要点及重要公式 知识要点及重要公式 经过整数倍零值内插使采样率升高至I倍， 不会引起新 I 的频谱混叠失真， 但是会产生I-1个镜像频谱（见图 7.2.3(b)）。 从时域考虑， 零值内插器输出v(m)的两个非零 值之间有I-1个零样值， 不是我们所期望的提高采样率后的 采样序列y(m)。 所以， 必须进行低通滤波， 滤除镜像频谱， 得到y(m)。
图7.2.4
第７章
多采样率数字信号处理
2． 有理数因子 采样率转换系统的功能 ． 有理数因子I/D采样率转换系统的功能 采样率转换系统的功能 有理数因子I/D采样率转换系统首先对输入序列x(n)按 整数因子I内插， 然后再对内插器的输出序列按整数因子 D抽取， 达到按有理数因子I/D的采样率转换。 如果仍用 Fx=1/Tx和Fy=1/Ty分别表示输入序列x(n)和输出序列y(m)的 采样频率， 则Fy=(I/D)Fx。应当注意， 先内插后抽取才能 最大限度地保留输入序列的频谱成分(请读者解释为什么)。
第７章
多采样率数字信号处理
整数因子抽取器的输入x(n)和输出y(m)的关系式如下：
y(m) = v( Dm) =
∑
k =0
∞
hD (k ) x( Dm − k )
(7.2.3)
D− 1 D−1 Y ( z) = H D (e − j2 πk / D z1/ D ) X (e− j2 πk / D z1/ D ) D k =0
第７章
多采样率数字信号处理
与抗混叠滤波器情况类似， 根据具体要求确定镜像 频谱滤波器其他三个指标参数（通带截止频率ωp、 通带最 大衰减αp、 阻带最小衰减αs）。 镜像频谱滤波器的ωp（或 过渡带宽度）和αp取决于内插系统对有用信号频谱的保真 度要求， 而αs取决于内插系统对镜像高频干扰的限制指标。 理想情况下， 镜像频谱滤波器的频率响应函数为
k =−∞
∞
∞
(7.2.18)
第７章
多采样率数字信号处理
整数因子D抽取器输出序列y(m)为
y (m) = w( Dm) =
k =− ∞
∑ h( Dm − kI )x(k )
∞
(7.2.19)
式（7.2.19）就是有理数因子I/D采样率转换系统的输入输出 时域关系。 如果线性滤波器用FIR滤波器实现， 则式 （7.2.19）为有限项之和， 所以可以直接按式（7.2.19）编 程序计算输出序列y(m)。 当然， 也可以采用教材上介绍的 各种高效实现结构以硬件或软硬结合来实现。
第７章
多采样率数字信号处理
下面对三种常用的采样率转换基本系统的重要知识点 及相关公式进行归纳总结， 以便读者复习巩固。 值得注意， 要理解采样率转换原理， 必须熟悉时域 采样概念、 时域采样信号的频谱结构、 时域采样定理。 此外， 时域离散线性时不变系统的时域分析和变换（Z变 换、 傅里叶变换）域分析理论是本章的分析工具。 只有熟 练掌握上述基础知识， 才能掌握本章的知识， 否则， 无 法理解本章内容。
I X (e jI ω ) jω Y (e ) = 0
π 0≤ ω < I (7.2.15) π ≤ ω ≤ π I
第７章
多采样率数字信号处理
7.2.3 有理数因子 采样率转换系统 有理数因子I/D采样率转换系统 采样率转换系统
1． 有理数因子 采样率转换系统原理框图 ． 有理数因子I/D采样率转换系统原理框图 采样率转换系统原理框图 有理数因子I/D采样率转换的原理框图如图7.2.4所示。
第７章
多采样率数字信号处理
Fx fs = = Hz 2 2D
相应的数字阻带截止频率为
Fy
(7.2.1)
2πf s π ωs = = Fx D
(7.2.2)
应当注意， 由于抗混叠滤波器工作于输入信号采样频 率Fx， 所以， 式（7.2.2）中用Fx换算得到相应的数字截 止频率， 绝对不能用Fy换算。
第７章
多采样率数字信号处理
抗混叠滤波器的通带截止频率（或过渡带宽度）取决于 抽取系统对信号频谱的失真度要求。 设计时根据具体要求确 定抗混叠滤波器的其他三个指标参数（通带截止频率ωp、 通 带最大衰减αp、 阻带最小衰减αs）。 例如， 要求抽取过程中频带［0, fp］上幅频失真小于1% 0, 1% （显然fp<fs）， 由于抽取引起的频谱混叠失真不超过0.1%。 这时， 抗混叠滤波器的通带截止频率为ωp=2πfp/Fx， 通带最 大衰减为αp=－20 lg(1－1%)=0.0873 dB， 阻带最小衰减为αs= －20 lg(0.1%)=60 dB。
第７章
多采样率数字信号处理
7.2 学习要点及重要公式
本章要求学生熟悉采样率转换的基本概念和种类， 了 解采样率转换的应用价值和适用场合， 掌握三种常用的采 样率转换基本系统（整数因子D抽取、 整数因子I插值和有 理数因子I/D采样率转换）的基本原理、 构成原理方框图及 其各种高效实现方法（FIR直接实现、 多相滤波器实现)， 每种实现方法的特点。 这些专业基础知识对进一步学习、 设计、 开发工作在多采样率状态的各种复杂系统是极其重 要的。 采样率转换的基本概念和种类， 以及应用价值和 适用场合在教材中已有较详细的介绍， 这里不再重复。
第７章
多采样率数字信号处理
镜像频谱滤波器hI(n)的指标也有两项： 阻带截止频率和 数字域阻带截止频率。 根据采样理论， 采样频率提高到Fy时， 采样信号序列y(m)的频谱以Fy为周期。 由于输入端信号x(n)的 采样频率为Fx， 所以x(n)的频带宽度不会超过Fx/2（对应的数 字频率为π）。 因此， 整数因子I内插器输出信号y(n)的频谱 一定是带宽为Fx/2、 重复周期为Fy的周期谱。 我们将零值内 插器输出的信号v(m)的频谱中除了y(m)频谱以外的其他频谱称 为镜像频谱。 镜像频谱滤波器的作用就是让y(m)频谱尽量无 失真通过， 滤除v(m)中的镜像频谱。 当I=3时， 图7.2.2中各 点信号的频谱示意图如图7.2.3所示。
第７章
多采样率数字信号处理
由于在数字域改变采样频率完全用软件实现， 所以， 采样频率可以任意选择， 又避免了对各种不同的采样频 率设计各种不同的抗混叠模拟滤波器， 从而使系统复杂 度大大降低。 采样率转换在现代通信、 信号处理和图像处理等领域 都有广泛的应用， 应用实例举不胜举， 本科大学生能理 解的几种实例在教材8.1节有介绍， 这里不再重复。
第７章
多采样率数字信号处理
图7.2.3
第７章
多采样率数字信号处理
综上所述可知， 镜像频谱滤波器的阻带截止频率为
Fx fs = Hz 2
相应的数字域阻带截止频率为
(7.2.6)
2πf s 2πFx / 2 π ωs = = = Fy IFx I
(7.2.7)
应当注意， 由于镜像频谱滤波器工作于输出信号采样 频率Fy， 所以， 式（7.2.7）中用Fy换算得到相应的数字截 止频率， 绝对不能用Fx换算。
j
(7.2.12)
I
Y ( z ) = H I ( z )V ( z ) = H I ( z ) X ( z )
(7.2.13) (7.2.14)
Y (e jω ) = H I (e jω )V (e jω ) = H I (e jω ) X (e jI ω )
第７章
多采样率数字信号处理
理想情况下， HI(ejω)由式（7.2.8）确定， 所以，
I jω H I (e y ) = 0 π 0≤ ω y < I π ≤ ωy ≤ π I
(7.2.8)
第７章
多采样率数字信号处理
通带内幅度常数取I是为了确保在m=0, ±I, ±2I, ±3I, …时， 输出序列y(m)=x(m/I)。 整数因子I内插系统的时域输入输出关系式如下：