路漫漫其修远兮,吾将上下而求索。 ——战国 · 屈原 离散时间傅里叶变换是将离散时间信号从时域映射到频域的重要手段,且映射结果有 明确的物理意义,它是数字信号处理技术的基础,在信号分析、系统分析和系统设计等方 面扮演着重要的角色。由于离散时间信号的自变量定义在整数域上,使离散时间傅里叶变 换不同于连续时间傅里叶变换。本章主要讨论离散时间傅里叶变换的基本概念和主要性质。 首先,给出离散时间傅里叶变换及其反变换;然后,论述离散时间傅里叶变换的存在条件和 特殊序列的傅里叶变换;其次,论述离散时间傅里叶变换的运算性质;最后,论述利用离散 时间傅里叶变换分析线性时不变系统的基本方法。 3.1 DTFT 的定义和条件 离散时间傅里叶变换(Discrete Time Fourier Transform,DTFT)是频谱分析、系统分 析、系统设计等的重要工具,它以完备正交函数集 {ejωn| n ∈ Z} 为基础对离散时间信号(序 列)进行正交展开,而序列的时域离散性使 DTFT 有别于连续时间傅里叶变换(Continuous Time Fourier Transform,CTFT)。 3.1.1 DTFT 的定义 序列 x[n] 的离散时间傅里叶变换定义为 X(ejω) = ∞X n=.∞ x[n]e.jωn (3-1) 通常称 X(ejω) 为频谱密度(简称频谱),它是数字频率 ω 的连续函数。X(ejω) 的离散时间 傅里叶反变换(Inverse Discrete Time Fourier Transform,IDTFT)可以表示为 x[n] = 1 2π Zπ .π X(ejω)ejωndω (3-2) 为了证明式(3-2)是式(3-1)的反变换,将式(3-1)两边乘以 ejωm,同时在 [.π, π] 内求积 分,可以得到 1 2π Zπ .π X(ejω)ejωmdω = 1 2π Zπ .π ∞X n=.∞ x[n]e.jωn!ejωmdω (3-3) 如果式(3-1)右边括号中的求和项一致收敛于 X(ejω),即满足一致收敛条件 lim M→∞ X(ejω) . XM(ejω) = 0, XM(ejω) = MX n=.M x[n]e.jωn (3-4) 则根据级数理论,可以交换式(3-3)右边的积分与求和的次序,即 1 2π Zπ .π X(ejω)ejωmdω = ∞X n=.∞ x[n] 1 2π Zπ .π ejω(m.n)dω (3-5) 根据复指数序列的运算性质,式(3-5)右边括号中的积分项为 1 2π Zπ .π ejω(m.n)dω = δ[m . n] =8<: 1, n = m 0, n .= m (3-6) 将式(3-6)代入式(3-5),可以得到 1 2π Zπ .π X(ejω)ejωmdω = ∞X n=.∞ x[n]δ[m . n] = x[m] (3-7) 再将式(3-7)中的变量 m 替换成变量 n,即得到式(3-2)给出的表达式。 因此,式(3-1)和式(3-2)构成了时域和频域之间可逆变换。称式(3-2) 为合成公式,表明 x[n] 是由频率位于 [.π, π] 内的无限多个复指数序列(频率分量)经过叠加运算(积分)得 到。称式(3-1)为分析公式,表明通过分析 x[n] 中包含的频率分量,可以确定合成 x[n] 时, 每个频率分量的相对大小。也就是说,x[n] 由无限小的复指数分量 1 2π ejωndω 叠加而成,而 每个复指数分量的相对大小取决于 X(ejω),且 ω 位于 [.π, π]。 序列 x[n] 的离散时间傅里叶变换 X(ejω) 既是关于 ω 的连续函数,又是以 2π 为周期 的周期函数,即满足条件 X(ejω) = X(ej(ω+2π))。简要证明如下: X(ej(ω+2π)) = ∞X n=.∞ x[n]e.j(ω+2π)n = ∞X n=.∞ x[n]e.jωn · e.j2πn = ∞X n=.∞ x[n]e.jωn · ..ej2πn.1 = X(ejω) · 1 = X(ejω) (3-8) 由于 X(ejω) 是以 2π 为周期的连续函数,因此可以将式(3-1)看作 X(ejω) 的傅里叶级 数展开式,而 x[n] 是傅里叶级数的展开式系数。虽然任何长度为 2π 的区间都是式(3-2)的 可用积分范围,但是将 ω 限定在 [.π, π] 或 [0, 2π] 内会更加简洁。 通常 X(ejω) 是关于 ω 的复函数,可以将 X(ejω) 表示成直角坐标形式,即 42 X(ejω) = XR(ejω) + jXI(ejω) (3-9) 其中:XR(ejω) ∈ R 是实部,是关于 ω 的偶函数;XI(ejω) ∈ R 是虚部,是关于 ω 的奇函数。 与此同时,还可以将 X(ejω) 表示成极坐标形式,即 X(ejω) = X(ejω) ej∠X(ejω) (3-10) 其中:|X(ejω)|(. 0) 为幅度,是关于 ω 的偶函数,反映序列的幅频特性;∠X(ejω) ∈ R 为 相位,是关于 ω 的奇函数,反映序列的相频特性。在实际应用中,称 X(ejω) 为傅里叶谱, 称 |X(ejω)| 为幅度谱,称 |∠X(ejω)| 为相位谱。 特别地,XR(ejω)、XI(ejω)、|X(ejω)| 和 ∠X(ejω) 都是以 2π 为周期的连续函数。在 ∠X(ejω) 上增加 2π 的整数倍,不会影响用式(3-10)表示 X(ejω)。在实际应用中,只能计算 出 ∠X(ejω) 位于 [.π, π] 或 [0, 2π] 内的相位主值,而无法消除 ∠X(ejω) 关于 2π 的模糊性。 3.1.2 一致收敛条件 虽然离散时间傅里叶变换的定义式(3-1)给出了计算频谱的方法,但是并非任意序列 x[n] 都存在频谱,讨论频谱的存在问题即讨论式(3-1)右边级数的收敛问题。如果该级数 收敛,则对于所有的 ω 值,必须满足条件 |X(ejω)| < ∞。根据式(3-1)可以得到 |X(ejω)| . ∞X n=.∞ |x[n]||e.jωn| = ∞X n=.∞ |x[n]| = S (3-11) 式(3-11)表明:如果 S < ∞,则 |X(ejω)| < ∞。也就是说,如果 x[n] 是绝对可求和的,则 X(ejω) 一定存在,即 x[n] 的绝对可求和是 X(ejω) 存在的充分条件。因此,当 S < ∞ 时, 式(3-1)右边级数一致收敛于 X(ejω),即对于所有的 ω 值满足条件 lim M→∞ X(ejω) . XM(ejω) = 0, XM(ejω) = MX n=.M x[n]e.jωn (3-12) 例 3.1 单位脉冲序列的离散时间傅里叶变换:计算单位脉冲序列 x[n] = δ[n] 的离散 时间傅里叶变换 X(ejω),绘制幅度谱 |X(ejω)| 和相位谱 ∠X(ejω)。 解 单位脉冲序列定义为 δ[n] =8<: 1, n = 0 0, n .= 0 (3-13) 将 x[n] = δ[n] 代入式(3-1)可以得到 X(ejω) = ∞X n=.∞ δ[n]e.jωn = δ[0]e.jω0 = 1 (3-14) 根据式(3-14)得到的幅度谱 |X(ejω)| 和相位谱 ∠X(ejω) 如图 3-1所示,在所有频率上 |X(ejω)| = 1 且 ∠X(ejω) = 0。 43 图 3-1 单位脉冲序列 δ[n] 的幅度谱和相位谱 例 3.2 矩形序列的离散时间傅里叶变换:计算长度为 N 的矩形序列 x[n] = RN[n] 的频率响应 X(ejω);给出当 N = 5 时的频率特性,包括 XR(ejω)、XI(ejω)、|X(ejω)| 和 ∠X(ejω)。 解 长度为 N 的矩形序列定义为 x[n] = RN[n] =8<: 1, 0 . n . N . 1 0, 其他 (3-15) 将 x[n] = RN[n] 代入式(3-1)可以得到 X(ejω) = ∞X n=.∞ RN[n]e.jωn = N.1 Xn=0 1 · e.jωn = N.1 Xn=0 (e.jω)n (3-16) 式(3-16)是等比数列的 N 项求和形式,求和结果可以表示为 X(ejω) = 1 . e.jNω 1 . e.jω = e.j(N.1)ω/2 · sin(Nω/2) sin(ω/2) (3-17) 根据式(3-17)可以得到 X(ejω) 的实部和虚部,即 XR(ejω) = cosN . 1 2 ωsin(Nω/2) sin(ω/2) (3-18a) XI(ejω) = .sinN . 1 2 ωsin(Nω/2) sin(ω/2) (3-18b) 以及幅度谱和相位谱,即 |X(ejω)| = sin(Nω/2) sin(ω/2) (3-19a) ∠X(ejω) = . N . 1 2 ω + argsin(Nω/2) sin(ω/2) (3-19b) 44 长度 N = 5 的矩形序列 R5[n] 如图 3-2(a)所示,实函数 XA(ejω) = sin(Nω/2)/ sin(ω/2) 如图 3-2(b)所示。根据式(3-18)得到的实部 XR(ejω) 如图 3-2(c)所示,它关于 ω = 0 呈 偶对称;虚部 XI(ejω) 如图 3-2(d)所示,它关于 ω = 0 呈奇对称。根据式(3-19)得到的幅 频谱 |X(ejω)| 如图 3-2(e)所示,它关于 ω = 0 呈偶对称;相位谱 ∠X(ejω) 如图 3-2(f) 所示,它关于 ω = 0 呈奇对称,且存在着关于 2π 的模糊性。 图 3-2 矩形序列 RN[n] 及其频率特性(N = 5) 例 3.3 指数序列的离散时间傅里叶变换:计算序列 x[n] = αnu[n] 的离散时间傅里 叶变换,绘制当 α = 0.8 和 α = .0.8 时的频谱 X(ejω)。 解 将 x[n] = αnu[n] 代入式(3-1)可以得到 X(ejω) = ∞X n=.∞ αnu[n]e.jωn = ∞Xn=0 αne.jωn = ∞Xn=0 (αe.jω)n (3-20) 式(3-20)是公比为 αe.jω 的等比数列。如果 |αe.jω| < 1,即 |α| < 1,则 x[n] = αnu[n] 是 45 绝对可求和的序列,即存在频谱 X(ejω)。利用等比序列的求和公式,可以得到 X(ejω) = 1 1 . αe.jω , |α| < 1 (3-21) 将欧拉公式 ejθ = cos θ + j sin θ 用于式(3-21),可以得到 X(ejω) 的实部和虚部,即 XR(ejω) = 1 . α cos ω 1 + α2 . 2α cos ω (3-22a) XI(ejω) = .α sin ω 1 + α2 . 2α cos ω (3-22b) 以及幅度和相位,即 |X(ejω)| = 1 [1 + α2 . 2α cos ω]1/2 (3-23a) ∠X(ejω) = arctan .α sin ω 1 . α cos ω (3-23b) 当 α = 0.80 和 α = .0.80 时,x[n] = αnu[n] 的波形分别如图 3-3(a)和图 3-3(b) 所示。实部 XR(ejω) 和虚部 XI(ejω) 分别如图 3-3(c)和图 3-3(d)所示,它们分别关于 ω = 0 呈偶对称和奇对称。幅度 |X(ejω)| 和相位 ∠X(ejω) 分别如图 3-3(e)和图 3-3(f) 所示,它们分别关于 ω = 0 呈偶对称和奇对称。由于幅度-相位的表示方法有明确的物理意 义,因此在工程实践中的应用更加广泛。 图 3-3 序列 x[n] = αnu[n] 及其频率特性(实线代表 α = 0.80,虚线代表 α = .0.80) 46 图 3-3 (续) 3.1.3 均方收敛条件 根据 3.1.2节可知:序列 x[n] 的绝对可求和是离散时间傅里叶变换 X(ejω) 存在的充分 条件,如果 X(ejω) 的傅里叶级数的展开式系数 x[n] = 1 2π Zπ .π X(ejω)ejωndω (3-24) 存在,则 X(ejω) 可以用傅里叶级数表示,但是式(3-1)右边级数可能不是一致收敛的。如果 x[n] 不满足绝对可求和条件,而满足相对松弛的平方可求和条件 ∞X n=.∞ |x[n]|2 < ∞ (3-25) 即 x[n] 在时域的总能量是有限值,则式(3-1)右边的级数按照均方收敛于 X(ejω),即 lim M→∞Zπ .π X(ejω) . XM(ejω) 2 dω = 0, XM(ejω) = MX n=.M x[n]e.jωn (3-26) 特别地,式(3-25)所示序列 x[n] 的平方可求和(总能量是有限值)条件,也是 X(ejω) 存在的充分条件。此外,根据 Cauchy-Schwarz 不等式可以得到 ∞X n=.∞ |x[n]|2 . ∞X n=.∞ |x[n]|!2 (3-27) 根据式(3-27)可知:如果 x[n] 是绝对可求和的,则它一定是平方可求和的。也就是说,如果 级数是一致收敛的,则一定是均方收敛的,反之不一定成立。 例 3.4 理想低通滤波器的均方收敛特性:讨论截止频率是 ωc 的理想低通数字滤波 器 Hlp(ejω) 的均方收敛特性。 解 理想低通滤波器的频率响应表示为 Hlp(ejω) =8<: 1, |ω| < ωc 0, ωc . |ω| . π (3-28) 47 其中:ωc 是截止频率。将式(3-28)代入式(3-2),可以得到理想单位脉冲响应 hlp[n] = 1 2π Zπ .π Hlp(ejω)ejωndω = 1 2π Zωc .ωc ejωndω = 1 2π · 1 jn ..ejωcn . e.jωcn= sin(ωcn) πn , .∞ < n < ∞ (3-29) 根据式(3-29)可知:当 n < 0 时,hlp[n] 不为零,即 hlp[n] 是非因果系统。特别地,由 于 sin(ωcn) 在区间 [.1, 1] 取值,因此当 n → ∞ 时,hlp[n] 以 1/n 形式趋近于零,即 hlp[n] 不是绝对可求和的。因此,hlp[n] 的离散时间傅里叶变换 ∞X n=.∞ hlp[n]e.jωn = ∞X n=.∞ sin(ωcn) πn e.jωn (3-30) 无法一致收敛于 Hlp(ejω),这是 Hlp(ejω) 在 ω = ωc 位置不连续的必然结果。 但是,hlp[n] 是平方可求和的,即总能量是有限的à,即 ∞X n=.∞ |hlp[n]|2 = ∞X n=.∞sin(ωcn) πn 2 = ωc π < ∞ (3-31) 式(3-31)使得式(3-30)右边的级数在均方意义下收敛于 Hlp(ejω),即 lim M→∞Zπ .π H(ejω) . HM(ejω) 2 dω = 0, HM(ejω) = MX n=.M sin(ωcn) πn e.jωn (3-32) 当截止频率 ωc = 0.5π 时,理想低通滤波器 Hlp(ejω) 及其单位脉冲响应 hlp[n] 分别如 图 3-4(a)和图 3-4(b)所示。当 M 取不同值时,HM(ejω) 的特性曲线如图 3-4(c) ~ 图 3-4(f)所示。随着 M 值的增大,HM(ejω) 在 ω = ωc 位置的振荡加剧,但是波纹幅度没 有增加。当 M → ∞ 时,最大振荡幅度(波纹)并不趋近于零,但是振荡行为逐渐收敛于 ω = ±ωc 位置。虽然 HM(ejω) 不能一致收敛于 Hlp(ejω),但是 hlp[n] 是平方可求和的,使 得 HM(ejω) 在均方意义下收敛于 Hlp(ejω)。也就是说,当 M → ∞ 时 .Hlp(ejω) 和 HM(ejω) 仅在 ω = ±ωc 位置存在差异,即 HM(ejω) 在均方误差为零意义上收敛于 Hlp(ejω)。 图 3-4 理想低通滤波器 Hlp(ejω) 的均方收敛特性(ωc = 0.5π) à 利用 Parseval 定理可以证明式(3-31),将在 3.3节论述 Parseval 定理。 48 图 3-4 (续) 3.2 周期序列的 DTFT 根据 3.1节可知:序列 x[n] 的一致收敛和均方收敛都是离散时间傅里叶变换 X(ejω) 存 在的充分条件。如果 x[n] 是周期序列,即对于所有的 n 值,满足条件 x[n] = x[n+N] (N ∈ Z+),则当 n → ∞ 时,x[n] 既不是绝对可求和的,又不是平方可求和的,使得 X(ejω) 既不 是一致收敛的,也不是均方收敛的。为了描述周期序列的频率特性,在频域中引入了冲激 函数 δ(ω)à,使周期序列也存在着离散时间傅里叶变换(DTFT)。本节仅讨论 4 种特殊周 期序列的 DTFT,对于任意周期序列的 DTFT,将在离散傅里叶级数的相关章节中讨论。 3.2.1 复指数和正弦序列 复指数序列 ejω0n 和正弦序列 cos(ω0n) 或 sin(ω0n) 是常见的基础序列,且可以利用欧 拉公式进行相互转换。将复指数和正弦序列纳入离散时间傅里叶变换的分析体系,对深入 了解数字频谱概念和拓展 DTFT 的分析范围非常重要。 1. 复指数序列 频率为 ω0 的复指数序列定义为 x[n] = ejω0n = cos(ω0n) + j sin(ω0n) (3-33) 当 ω0 = 0.25π 时,复指数序列 ejω0n 的实部 xc[n] = cos(ω0n) 和虚部 xs[n] = sin(ω0n) 分别 如图 3-5(a)和图 3-5(b)所示。 à 关于冲激函数的定义和性质,可以参见《信号与系统》教材。 49 图 3-5 复指数序列和正弦序列的频谱结构(ω0 = 0.25π) 复指数序列 x[n] = ejω0n 的离散时间傅里叶变换表示为 X(ejω) = ∞X k=.∞ 2πδ(ω . ω0 + 2πk) (3-34) 根据式(3-34)可知,X(ejω) 是以 ω0 为中心,以 2π 的整数倍为间隔的一组冲激函数,且每 个冲激函数的积分面积是 2π。当 ω0 = 0.25π 时,复指数序列的频谱如图 3-5(c)所示。可 以看出,X(ejω) 在频域上呈现以 2π 为周期的“线谱”结构。 下面证明式(3-33)和式(3-34)可以构成可逆关系。将式(3-34)代入离散时间傅里叶反变 换的定义式(3-2),可以得到 x[n] = 1 2π Zπ .π ∞X k=.∞ 2πδ(ω . ω0 + 2πk)!ejωndω (3-35) 由于式(3-35)限定了积分区间 [.π, π],故只考虑当 k = 0 时的情况,即 x[n] = 1 2π Zπ .π 2πδ(ω . ω0)ejωndω = ejω0n (3-36) 因此,式(3-33)和式(3-34)构成了离散时间傅里叶变换与反变换关系。 2. 正弦序列 频率为 ω0 的正弦序列 cos(ω0n) 和 sin(ω0n) 是复指数序列 ejω0n 的实部和虚部,根据 欧拉公式,用复指数序列将它们表示为 xc[n] = cos(ω0n) = ejω0n + e.jω0n 2 (3-37a) 50 xs[n] = sin(ω0n) = ejω0n . e.jω0n 2j (3-37b) 因此,根据式(3-33)和式(3-34)所示的可逆变换关系,可以得到 Xc(ejω) = π ∞X k=.∞ [δ(ω . ω0 + 2πk) + δ(ω + ω0 + 2πk)] (3-38a) Xs(ejω) = .jπ ∞X k=.∞ [δ(ω . ω0 + 2πk) . δ(ω + ω0 + 2πk)] (3-38b) 根据式(3-38)可知,Xc(ejω) 和 Xs(ejω) 也是以 ±ω0 为中心的,以 2π 的整数倍为间隔的一 组冲激函数。当 ω0 = 0.25π 时,余弦序列 xc[n] = cos(ω0n) 的频谱 Xc(ejω) 如图 3-5(d) 所示,可以看出,Xc(ejω) 在频域上呈现以 2π 为周期的“线谱”结构。 特别地,对于序列 ejω0n、cos(ω0n) 和 sin(ω0n) 而言,只有当 2π/ω0 为有理数时,它们 才是周期序列(周期 N = (2π/ω0)k,k ∈ Z)。由于上述内容没有限定 2π/ω0 是否为有理 数,因此离散时间傅里叶变换结果的适用范围更为广泛。 3.2.2 单位脉冲的衍生序列 任意序列 x[n] 可以表示为单位脉冲序列 δ[n] 的“平移-加权-求和”形式,即 x[n] = ∞X k=.∞ x[k]δ[n . k], .∞ < n < ∞ (3-39) 因此,利用式(3-39)可以构造出特殊的周期序列(如单位值序列等),它们在数字信号的时 域和频域分析中占有独特的地位。 1. 单位值序列 单位值序列是幅值为 1 的序列,如图 3-6(a)所示,可以表示为 x[n] = ∞X k=.∞ δ[n . k] (3-40) 或者表示成更简单的形式 x[n] = 1, .∞ < n < ∞ (3-41) 由于 x[n] = ejω0n|ω0=0 = ej0 = 1,因此单位值序列是频率 ω0 = 0 的复指数序列。将 ω0 = 0 代入式 (3-1) 和式(3-34),可以得到单位值序列 x[n] = 1 (n ∈ Z) 的离散时间傅里 叶变换 X(ejω) = ∞X k=.∞ 1 · e.jωn = ∞X k=.∞ 2πδ(ω + 2πk) (3-42) 根据式(3-42)可知,X(ejω) 是以 ω = 0 为中心,以 2π 为间隔的一组冲激函数,且每个冲激 函数的积分面积为 2π,如图 3-6(b)所示,它同样呈现以 2π 为周期的“线谱”特性。 51 图 3-6 单位值序列及其频谱结构 2. 单位脉冲周期序列 周期为 N 的单位脉冲周期序列,可以表示为 x[n] = ∞X k=.∞ δ[n . kN] (3-43) 当 N = 5 时,单位脉冲周期序列如图 3-7(a)所示。 图 3-7 单位脉冲周期序列及其频谱结构 将式(3-43)代入式(3-1),可以得到 X(ejω) = ∞X n=.∞ ∞X k=.∞ δ[n . kN]!e.jωn = ∞X k=.∞ ∞X n=.∞ δ[n . kN]e.jωn = ∞X k=.∞ e.jωNk (3-44) 如果将式(3-44)的 ωN 看作整体,则式(3-44)表示单位值序列的离散时间傅里叶变换。因此, 根据式(3-42)和式(3-44)可以得到 X(ejω) = ∞X k=.∞ 2πδ(ωN + 2πk) = ∞X k=.∞ 2πδ (N (ω + 2πk/N)) (3-45) 再利用冲激函数的性质 δ(at) = δ(t)/|a|,根据式(3-45)可以得到 52 X(ejω) = ∞X k=.∞ 2π N δ ω + 2πk N (3-46) 由此可知,X(ejω) 是以 ω = 0 为中心,以 ω = 2π/N 为间隔的一组冲激函数,且每个 冲激函数的积分面积为 2π/N。当 N = 5 时,X(ejω) 的频谱结构如图 3-7(b)所示。比 较图 3-6和图 3-7可以看出,在时域和频域上,单位值序列和单位脉冲周期序列存在着对偶 关系。 利用复指数序列、正弦序列、单位值序列以及单位脉冲周期序列的离散时间傅里叶变 换,可以推导出满足 x[n] = x[n + N](.∞ < n < ∞) 关系的任意周期序列的离散时间傅里 叶变换,并可以用离散傅里叶级数(DFS)表示。这些内容超出本章范围,在此不再赘述。 典型序列的离散时间傅里叶变换如表 3-1所示。 表 3-1 典型序列的离散时间傅里叶变换 序列 x[n] 离散时间傅里叶变换 X(ejω) δ[n] 1 δ[n . n0] e.jωn0 1 (.∞ < n < ∞) ∞P k=.∞ 2πδ(ω + 2πk) u[n] 1 1 . e.jω + ∞P k=.∞ πδ(ω + 2πk) anu[n] (|a| < 1) 1 1 . ae.jω (n + 1)anu[n] (|a| < 1) 1 (1 . ae.jω)2 sin ωcn πn X(ejω) =8<: 1, |ω| < ωc 0, ωc < |ω| . π rn sin ω0(n + 1) sin ω0 u[n] (|r| < 1) 1 1 . 2r cos ω0e.jω + r2e.j2ω x[n] =8<: 1, 0 . n . M 0, 其他 sin[ω(M + 1)/2] sin(ω/2) e.jωM/2 ejω0n ∞P k=.∞ 2πδ(ω . ω0 + 2πk) cos(ω0n + .) ∞P k=.∞πej.δ(ω . ω0 + 2πk) + πe.j.δ(ω + ω0 + 2πk) sin(ω0n + .) .j ∞P k=.∞πej.δ(ω . ω0 + 2πk) . πe.j.δ(ω + ω0 + 2πk) 3.3 DTFT 的主要性质 离散时间傅里叶变换(DTFT)和连续时间傅里叶变换(CTFT)有很多类似的性质, 掌握它们可使离散时间信号分析和处理更加快捷。为了论述方便,用符号 F 代表 DTFT 和 53 IDTFT 的可逆变换关系,序列 x[n] 与其频谱 X(ejω) 的对应关系记作 x[n] F ←→ X(ejω)。 3.3.1 基本性质 1. 线性性质 如果 x[n] F ←→ X(ejω),y[n] F ←→ Y (ejω),α 和 β 是常系数,则 αx[n] + βy[n] F ←→ αX(ejω) + βY (ejω) (3-47) 式(3-47)表明,x[n] 和 y[n] 在时域上的线性组合对应着 X(ejω) 和 Y (ejω) 在频域上的线性 组合。简要证明过程如下: ∞X n=.∞ (αx[n] + βy[n])e.jωn = ∞X n=.∞ (αx[n])e.jωn + ∞X n=.∞ (βy[n])e.jωn = α ∞X n=.∞ x[n]e.jωn + β ∞X n=.∞ y[n]e.jωn = αX(ejω) + βY (ejω) 2. 时间平移 如果 x[n] F ←→ X(ejω),则 x[n . m] F ←→ e.jωmX(ejω) (3-48) 式(3-48)表明,x[n] 在时域上平移 m 个单位,对应着 X(ejω) 在频域上产生 .ωm 的相位移 位(简称相移),且相移大小与频率呈正比。简要证明过程如下: ∞X n=.∞ x[n . m]e.jωn = ∞X k=.∞ x[k]e.jω(m+k) = e.jωm ∞X k=.∞ x[k]e.jωk = e.jωmX(ejω) 3. 时间翻转 如果 x[n] F ←→ X(ejω),则 x[.n] F ←→ X(e.jω) (3-49) 式(3-49)表明,x[n] 在时域上的翻转,对应着 X(ejω) 在频域上的翻转。简要证明过程如下: ∞X n=.∞ x[.n]e.jωn = .∞ Xk=∞ x[k]ejωk = .∞ Xk=∞ x[k]e.j(.ω)k = X(e.jω) 如果 x[n] 是实序列,则 x[.n] F ←→ X.(ejω),即存在 X.(ejω) = X(e.jω)。由于 |X(ejω)| = |X.(ejω)|,因此实序列的时间翻转不改变幅度特性,而只改变相位特性——相位值取负号或 相位平移了 π 弧度(ejπ = .1)。 54 4. 序列共轭 如果 x[n] F ←→ X(ejω),则 x.[n] F ←→ X.(e.jω) (3-50) 式(3-50)表明,x[n] 在时域上取共轭,对应着 X(ejω) 在频域上翻转之后取共轭。简要证明 过程如下: ∞X n=.∞ x.[n]e.jωn = ∞X n=.∞ x[n]ejωn!. = ∞X n=.∞ x[n]e.j(.ω)n!. = X.(e.jω) 5. 频率微分 如果 x[n] F ←→ X(ejω),则 nx[n] F ←→ jdX(ejω) dω (3-51) 式(3-51)表明,x[n] 在时域上对下标值 n 的线性加权,对应着 X(ejω) 在频域上取一阶导数 并乘以 j 因子。简要证明过程如下: jdX(ejω) dω = j d dω ∞X n=.∞ x[n]e.jωn!= j ∞X n=.∞ x[n] d(e.jωn) dω = j ∞X n=.∞ x[n](.jn · e.jωn) = ∞X n=.∞ nx[n] · e.jωn 6. 指数乘法 如果 x[n] F ←→ X(ejω),则 anx[n] F ←→ X ejω a (3-52) 式(3-52)表明,x[n] 在时域上与指数序列 an 相乘,对应着 X(ejω) 在频域上用 ejω/a 代替 ejω。简要证明过程如下: ∞X n=.∞ (anx[n])e.jωn = .∞ Xn=∞ x[n]ejω a .n = X ejω a 7. 频率平移 如果 x[n] F ←→ X(ejω),则 ejω0nx[n] F ←→ X ..ej(ω.ω0) (3-53) 式(3-53)表明,x[n] 与复指数序列 ejω0n 在时域上相乘,对应着 X(ejω) 在频域上产生 ω0 的 55 频率平移(简称频移)。简要证明过程如下: ∞X n=.∞ (ejω0nx[n])e.jωn = .∞ Xn=∞ x[n]e.j(ω.ω0)n = X(ej(ω.ω0)) 8. 时域卷积 如果 x[n] F ←→ X(ejω),y[n] F ←→ Y (ejω),则 x[n] . y[n] F ←→ X(ejω)Y (ejω) (3-54) 式(3-54)表明,即 x[n] 和 y[n] 在时域上的卷积,对应着 X(ejω) 和 Y (ejω) 在频域上的乘积。 简要证明过程如下: x[n] . y[n] = ∞X k=.∞ x[k]y[n . k] = ∞X k=.∞ 1 2π Zπ .π X(ejω)ejωkdωy[n . k] = 1 2π Zπ .π X(ejω) ∞X k=.∞ y[n . k]ejωk!dω = 1 2π Zπ .π X(ejω) Y (ejω)ejωndω = 1 2π Zπ .π X(ejω)Y (ejω)ejωndω 9. 频域卷积 如果 x[n] F ←→ X(ejω),y[n] F ←→ Y (ejω),则 x[n]y[n] F ←→ 1 2π X(ejω) . Y (ejω) (3-55) 式(3-55)表明,x[n] 和 y[n] 在时域上的乘积,对应着 X(ejω) 和 Y (ejω) 在频域上的卷积。简 要证明过程如下: ∞X n=.∞ x[n]y[n]e.jωn = ∞X n=.∞ 1 2π Zπ .π X(ejθ)ejθndθy[n]e.jωn = 1 2π Zπ .π X(ejθ) ∞X n=.∞ y[n]e.j(ω.θ)n!dθ = 1 2π Zπ .π X(ejθ)Y (ej(ω.θ))dθ = 1 2π X(ejω) . Y (ejω) 10. Parseval 定理 如果 x[n] F ←→ X(ejω),y[n] F ←→ Y (ejω),则 ∞X n=.∞ x[n]y.[n] = 1 2π Zπ .π X(ejω)Y .(ejω)dω (3-56) 56 式(3-56)表明,x[n] 和 y[n] 在时域上的内积,对应着 X(ejω) 和 Y (ejω) 在频域上的内积。特 别地,当 x[n] = y[n] 时,由于 |x[n]|2 = x[n]x.[n] 且 |X(ejω)|2 = X(ejω)X.(ejω),因此可 以得到 ∞X n=.∞ |x[n]|2 = 1 2π Zπ .π X(ejω) 2 dω (3-57) 式(3-57)表明,x[n] 在时域上的总能量等于 X(ejω) 在频域上的总能量,因此又称 Parserval 定理为能量守恒定理。通常称 |X(ejω)|2/(2π) 为能量谱密度,代表了在频域上的能量分布。 由于频谱函数具有周期为 2π 的特性,因此在任意长度为 2π 区间内都可以求出频域的总能 量。式(3-57)既是普适的能量守恒定律在数字信号处理领域的体现,又是在时域和频域能够 实现等效分析的理论依据。式(3-56)和式(3-57)的简要证明过程如下: ∞X n=.∞ x[n]y.[n] = ∞X n=.∞ x[n] 1 2π Zπ .π Y (ejω)ejωndω. = 1 2π Zπ .π Y .(ejω) ∞X n=.∞ x[n]e.jωn!dω = 1 2π Zπ .π Y .(ejω)X(ejω)dω = 1 2π Zπ .π X(ejω)Y .(ejω)dω 因此,式(3-56)得证;再令 x[n] = y[n],式(3-57)得证。 上述离散时间傅里叶变换的主要性质如表 3-2所示。 表 3-2 离散时间傅里叶变换的基本性质 序 列 离散时间傅里叶变换 x[n . n0] e.jωn0X(ejω) x[.n] X(e.jω) (如果 x[n] ∈ R,则 X.(ejω) = X(e.jω)) anx[n] X(ejω/a) ejω0nx[n] X ej(ω.ω0) nx[n] j dX(ejω) dω ax[n] + by[n] aX(ejω) + bY (ejω) x[n] . y[n] X(ejω)Y (ejω) x[n]y[n] 1 2π ∫π .π X(ejθ)Y (ej(ω.θ))dθ Parseval 定理: ∞X n=.∞ x[n]y.[n] = 1 2π ∫π .π X(ejω)Y .(ejω)dω ∞X n=.∞ |x[n]|2 = 1 2π ∫π .π X(ejω) 2 dω 57 3.3.2 对称性质 在离散时间傅里叶变换过程中,一般假定序列 x[n] 是复序列(将实序列看作复序列的 特殊情况)。将序列 x[n] 表示成实部-虚部形式,即 x[n] = xr[n] + jxi[n] (3-58) 其中 xr[n] = 1 2 (x[n] + x.[n]) (3-59a) xi[n] = 1 2j (x[n] . x.[n]) (3-59b) 通常称满足 xe[n] = x.e [.n] 关系的复序列为共轭对称序列,称满足 xo[n] = .x.o [.n] 关系的复序列为共轭反对称序列。任何复序列 x[n] 都可以表示为共轭对称序列 xe[n] 和共 轭反对称序列 xo[n] 的求和形式,即 x[n] = xe[n] + xo[n] (3-60) 其中:xe[n] 是偶函数,xo[n] 是奇函数,即 xe[n] = 1 2 (x[n] + x.[.n]) = x.e [.n] (3-61a) xo[n] = 1 2 (x[n] . x.[.n]) = .x.o [.n] (3-61b) 同理,序列 x[n] 的离散时间傅里叶变换 X(ejω) 也可以表示成实部-虚部形式,即 X(ejω) = Xr(ejω) + jXi(ejω) (3-62) 其中:Xr(ejω) 是偶函数,Xi(ejω) 是奇函数,即 Xr(ejω) = 1 2 X(ejω) + X.(ejω)= X.r (ejω) (3-63a) Xi(ejω) = 1 2j X(ejω) . X.(ejω)= .X.i (ejω) (3-63b) 还可以将 X(ejω) 表示成共轭对称函数 Xe(ejω) 和共轭反对称函数 Xo(ejω) 的求和形式, 即 X(ejω) = Xe(ejω) + Xo(ejω) (3-64) 其中:Xe(ejω) 是偶函数,Xo(ejω) 是奇函数,即 Xe(ejω) = 1 2 X(ejω) + X.(e.jω)= X.e (e.jω) (3-65a) 58 Xo(ejω) = 1 2 X(ejω) . X.(e.jω)= .Xo(e.jω) (3-65b) 如果将 x[n] 和 X(ejω) 的可逆变换关系表示为 X(ejω) = DTFT(x[n]),则利用式(3-1)可 以得到:X.(e.jω) = DTFT(x.[n]),X(e.jω) = DTFT(x[.n]),X.(ejω) = DTFT(x.[.n])。 比较式(3-59)与式(3-65),式(3-61)与式(3-63),可以得到时域、频域对应关系如下: Xe(ejω) = DTFT(xr[n]), Xo(ejω) = DTFT(jxi[n]) (3-66a) Xr(ejω) = DTFT(xe[n]), jXi(ejω) = DTFT(xo[n]) (3-66b) 式(3-66)给出了离散时间傅里叶变换的时域和频域对偶关系,以此为基础可以推导出更多、 更具体的对应关系,如因果序列的对称关系等à,在此不再赘述。 特别地,实序列 x[n] 的离散时间傅里叶变换 X(ejω) 存在如下关系: Re[X(ejω)] = Re[X(e.jω)], Im[X(ejω)] = .Im[X(e.jω)] (3-67a) |X(ejω)| = |X(e.jω)|, arg[X(ejω)] = .arg[X(e.jω)] (3-67b) 式(3-67)表明,当 x[n] ∈ R 时,X(ejω) 的实部是偶函数,虚部是奇函数;X(ejω) 的幅度 是偶函数,相位是奇函数。在实际工程应用中,经常使用式(3-67)检验 DTFT 结果的正 确性。 离散时间傅里叶变换的对称性质如表 3-3所示,为简化求解过程提供了便捷工具。 表 3-3 离散时间傅里叶变换的对称性质 序号 序列 离散时间傅里叶变换 1 x.[n] X.(e.jω) 2 x.[.n] X.(ejω) 3 Re(x[n]) Xe(ejω) 4 jIm(x[n]) Xo(ejω) 5 xe[n] XR(ejω) = Re ..X(ejω) 6 xo[n] jXI(ejω) = jIm ..X(ejω) 7 x[n] ∈ R X(ejω) = X.(e.jω) 8 x[n] ∈ R XR(ejω) = XR(e.jω) 9 x[n] ∈ R XI(ejω) = .XI(e.jω) 10 x[n] ∈ R |X(ejω)| = |X(e.jω)| 11 x[n] ∈ R ∠X(ejω) = .∠X(e.jω) à 参见程佩青编著的《数字信号处理教程》(第五版),清华大学出版社。 59 3.4 LTI 系统的频率响应 如果离散时间线性时不变系统的输入序列为 x[n],单位脉冲响应为 h[n],则输出序列 是 y[n] = x[n] . h[n]。根据卷积的交换性质(x[n] . h[n]=h[n] . x[n])可知,x[n] 和 h[n] 在 本质上是相同的,都是离散时间信号,即序列。因此,关于序列 x[n] 的离散时间傅里叶变 换 X(ejω) 的收敛性质、运算性质、对称性质等,同样适用于线性时不变(LTI)系统的单 位脉冲响应 h[n],也就是说,离散时间傅里叶变换为 LTI 系统分析提供了理论基础。 3.4.1 频率响应概念 如果 LTI 系统的单位脉冲响应为 h[n],输入序列为 x[n],则输出序列为 y[n] = h[n] . x[n] = ∞X k=.∞ h[k]x[n . k] (3-68) 假定 x[n] 是复指数序列 ejωn,根据式(3-68)可以得到 y[n] = h[n] . ejωn = ∞X k=.∞ h[k]ejω(n.k) = ejωn ∞X k=.∞ h[k]e.jωk = H(ejω)ejωn (3-69) 其中: H(ejω) = ∞X k=.∞ h[k]e.jωk = ∞X n=.∞ h[n]e.jωn (3-70) 式(3-69)和式(3-70)表明:当 LTI 系统的输入是复指数序列 ejωn 时,输出是被 H(ejω) 加权 且频率相同的复指数序列,其中 H(ejω) 是 h[n] 的离散时间傅里叶变换。通常称 ejωn 为 LTI 系统的特征函数,称 H(ejω) 为 LTI 系统的频率响应。利用离散时间傅里叶反变换公 式(3-2),可以依据 H(ejω) 获得单位脉冲响应 h[n] = 1 2π Zπ .π H(ejω)ejωndω (3-71) 通常,LTI 系统的频率响应 H(ejω) 是与单位脉冲响应有关且与输入序列无关的复值函 数,可以将 H(ejω) 表示成实部-虚部形式,即 H(ejω) = HR(ejω) + jHI(ejω) (3-72) 其中:HR(ejω) ∈ R 为实部响应,是偶对称函数;HI(ejω) ∈ R 为虚部响应,是奇对称函数。 还可以将 H(ejω) 表示成幅度-相位形式,即 H(ejω) = |H(ejω)|ej∠H(ejω) (3-73) 60