跳转到内容

2.3 线性时不变系统

杨鸿文 主编(北京邮电大学通信原理课程组)北京邮电大学出版社《通信原理(第5版)》大学本科 · 教材与微课双语助读 (Alt+Y)

从发端到收端,信号需要经历各种各样的变换,有些由模拟电路完成,有些由数字信号处理算法实现。此外,信号所经过的传输媒介也会对信号产生变换作用。无论这些变换的具体细节如何,都可以建模为一个黑箱,它将输入 x(t) 映射为输出 y(t)。

graph LR
    X["x(t)"] --> Block["变换系统"]
    Block --> Y["y(t)"]

把一个信号变成另一个信号的过程构成一个信号系统。信号系统有线性、非线性、时变、时不变之分,本章主要考虑线性时不变系统,也叫线性滤波器或滤波器。

线性时不变(Linear Time-Invariant, LTI)系统同时具有线性与时不变两个特性:

  • 线性(Linearity):多个输入信号线性组合后的输出是各自输出的线性组合;
  • 时不变(Time-Invariance):输入延迟则输出也等量延迟。

线性时不变系统不会产生输入所没有的频率成分,非线性系统或时变系统可能产生新的频率分量。

例 2.3.1

信号 x(t) 经过某线性调制系统后的输出为 y(t) = x(t)cos(2π f₀ t)。分析其系统性质。

解

由于:

[x₁(t) + x₂(t)]cos(2π f₀ t) = [x₁(t)cos(2π f₀ t)] + [x₂(t)cos(2π f₀ t)]

该系统为线性系统。但它不满足时不变性:x(t) 的输出是 y(t) = x(t)cos(2π f₀ t),x(t - t₀) 的输出 x(t - t₀)cos(2π f₀ t) 不等于 y(t - t₀) = x(t - t₀)cos[2π f₀ (t - t₀)]。

若 f₀ = 10 Hz,输入是频率为 1 Hz 的单频信号 cos(2π t),输出为:

y(t) = cos(2π t)cos(20π t) = 1/2cos(18π t) + 1/2cos(22π t)

输出包含两个单频信号,频率分别为 9 Hz 和 11 Hz,产生了输入所没有的新频率。

例 2.3.2

将信号 x(t) 通过一个平方器,输出为 y(t) = x²(t)。分析其系统性质。

解

输入变成 x(t - t₀) 后,输出变成 x²(t - t₀) = y(t - t₀),因此平方器是时不变系统。但将两个信号 x₁(t), x₂(t) 叠加后送入平方器,输出为 [x₁(t) + x₂(t)]² ≠ x₁²(t) + x₂²(t),说明平方器是非线性系统。

若输入是频率为 1 Hz 的单频信号 cos(2π t),输出为:

y(t) = cos²(2π t) = 1/2 + 1/2cos(4π t)

输出包含一个直流信号和一个频率为 2 Hz 的倍频信号。

例 2.3.3

某系统的输入输出关系为 y(t) = x(t) + x(t - T)。分析其系统特性。

解

输入分别为 x₁(t), x₂(t) 时,输出为 y₁(t) = x₁(t) + x₁(t - T) 和 y₂(t) = x₂(t) + x₂(t - T)。输入为 x(t) = x₁(t) + x₂(t) 时,输出为:

y(t) = x(t) + x(t - T) = x₁(t) + x₂(t) + x₁(t - T) + x₂(t - T) = y₁(t) + y₂(t)

当输入为 w(t) = x(t - t₀) 时,输出为:

w(t) + w(t - T) = x(t - t₀) + x(t - t₀ - T) = y(t - t₀)

因此,该系统是线性时不变系统。

若 T = 1/4,输入信号是频率为 1 Hz 的单频信号 cos(2π t),输出为:

y(t) = cos(2π t) + cos(2π t - (π)/(2)) = √(2)cos(2π t - (π)/(4))

输出仍然是频率为 1 Hz 的单频信号,只是幅度和相位发生了变化。

线性时不变系统不会产生输入所没有的频率分量。当滤波器的输入是频率为 f 的单频信号 e^(j2π f t) 时,输出还是这个单频信号,只是幅度和相位可能发生变化,即输出是输入 e^(j2π f t) 乘上一个表示幅度与相位变化的复系数 Ae^(jφ)。输入频率不同时,幅度和相位的变化可能不同,即复系数 Ae^(jφ) 与频率有关,可记其为 A(f)e^(jφ(f)) = H(f)。此时,滤波器输出为 H(f)e^(j2π f t)。称 H(f) 为滤波器的传递函数(Transfer Function)或频率响应,其幅频特性是 A(f),相频特性是 φ(f)。

按照频谱分析的观点,滤波器输入信号 x(t) 由不同频率的单频信号 e^(j2π f t) 组成,各个频率成分的含量体现为频谱密度 X(f)。滤波器对各频率分量逐一调整幅度与相位,在频率 f 处幅度乘以 A(f),相位增加 φ(f),使滤波器输出信号的频谱密度相应变成:

Y(f) = H(f)X(f)

上式右边是频域乘积,在时域对应卷积。若 H(f) 的傅氏反变换为 h(t),则滤波器输出信号为:

y(t) = x(t) * h(t) = ∫_(-∞)^(∞) x(t - τ)h(τ) dτ

若滤波器输入是单位冲激 x(t) = δ(t),代入上式后,输出是 h(t),即 h(t) 是滤波器对单位冲激的响应,简称冲激响应(Impulse Response)。根据线性时不变特性,若滤波器输入的冲激改变强度与位置,变成 a · δ(t - t₀),输出将变成 a · h(t - t₀)。

滤波器的输入 x(t)、冲激响应 h(t)、输出 y(t) 都可以是复信号。按复信号考虑,若设计 h(t) = g^*(-t),则滤波器输出为:

y(t) = ∫_(-∞)^(∞) x(t - τ)h(τ) dτ = ∫_(-∞)^(∞) x(t - τ)g^*(-τ) dτ = Rₓg(t)

输出是 x(t) 与 g(t) 的互相关函数。说明通过适当设计冲激响应,可以用滤波器来测量两个信号的互相关函数。这是匹配滤波器的核心数学渊源。

理想滤波器是一个频域窗口,它能够理想滤除窗口外的频率成分,完整保留窗口内的分量。图 2.3.2 示出了理想低通滤波器和理想带通滤波器的传递函数。本书提到的低通滤波器(Low-Pass Filter, LPF)和带通滤波器(Band-Pass Filter, BPF)主要指理想滤波器。

图 2.3.2 理想低通滤波器与理想带通滤波器

图 2.3.2 理想低通滤波器与理想带通滤波器(原书影印参考)

根据传递函数可以求出冲激响应:

  • 理想低通滤波器(带宽为 W、增益为 1): H_(LPF)(f) = operatornamerect((f)/(2W)) iff h_(LPF)(t) = 2W · operatornamesinc(2Wt)
  • 理想带通滤波器(带宽为 B、增益为 1,中心频率为 f_c): H_(BPF)(f) = operatornamerect((f - f_c)/(B)) + operatornamerect((f + f_c)/(B)) 其冲激响应为: h_(BPF)(t) &= B · operatornamesinc(Bt)e^(j2π f_c t) + B · operatornamesinc(Bt)e⁻j2π f_c t &= 2B · operatornamesinc(Bt)cos(2π f_c t)

周期信号 y(t) = ∑ₙ=-∞^(∞) x(t - nT) 可以看成是周期冲激序列 δ_T(t) = ∑ₙ=-∞^(∞) δ(t - nT) 通过一个冲激响应为 x(t) 的滤波器的输出。

graph LR
    Delta["δ_T(t) = Σ δ(t - nT)"] --> Filter["系统 x(t)"]
    Filter --> Out["y(t) = Σ x(t - nT)"]

设滤波器的传递函数为 X(f)。周期信号 y(t) 的傅氏变换为:

Y(f) = X(f) · 1/T∑ₘ=-∞^(∞) δ(f - m/T) = 1/T∑ₘ=-∞^(∞) X(m/T)δ(f - m/T)

两边做傅氏反变换可得到周期信号 y(t) 的傅里叶级数展开式为:

∑ₙ=-∞^(∞) x(t - nT) = 1/T∑ₘ=-∞^(∞) X(m/T)e^(j2π m/T t)

通信信号传输的理想情形是信号不发生失真,即波形不变。例如,在电脑上播放音乐,调节音量大小不会造成失真,播放时间略有推迟也不是失真。实信号通过无失真系统,其输入输出关系是:

y(t) = c · x(t - t₀)

其中 c 是实常数。两边做傅氏变换得到:

Y(f) = c · e⁻j2π f t₀ X(f)
无失真传输系统的充要条件

无失真系统的传递函数在拟传输信号的频带范围内必须满足:

H(f) = c · e⁻j2π f t₀
  1. 幅频特性是常数:A(f) = |H(f)| = c(全通,无幅度失真);
  2. 相频特性是过原点的直线:φ(f) = angle H(f) = -2π f t₀(线性相位,无相位失真)。

在通信系统中,如果携带信息的信号经过的某个环节是无失真系统,那么原则上说,无论 c, t₀ 具体是多少,该环节对通信系统的原理没有影响。因此,当不需要特别关注幅度与时延时,对于无失真系统(如理想信道),我们经常假设 c = 1, t₀ = 0。


2.3.2 滤波器输出的能量及功率谱密度

Section titled “2.3.2 滤波器输出的能量及功率谱密度”

滤波器会改变信号中各个频率分量的强度,进而改变信号的能量谱密度或功率谱密度。

设滤波器的输入 x(t) 是能量信号,其能量谱密度是 Eₓ(f) = |X(f)|²。通过传递函数为 H(f) 的滤波器之后,输出信号 y(t) 的傅氏变换为 Y(f) = H(f)X(f),其能量谱密度为:

E_y(f) = |Y(f)|² = |H(f)|² · |X(f)|² = |H(f)|² Eₓ(f)

即输出能量谱密度是输入能量谱密度乘以传递函数的模平方。

对 E_y(f) 做傅氏反变换,频域乘积对应时域卷积,因此滤波器输出的自相关函数为:

R_y(τ) = Rₓ(τ) * h(τ) * h^*(-τ)

同理,设滤波器的输入 x(t) 是功率信号,其功率谱密度为 Pₓ(f)。输入为截短信号 x_T(t) 时,滤波器输出信号的频谱密度为 H(f)X_T(f),能量谱密度为 |H(f)|² · |X_T(f)|²。除以 T 后取极限,得到滤波器输出的功率谱密度为:

P_y(f) = |H(f)|² Pₓ(f)

滤波器输入输出的谱密度与相关函数关系总结于表 2.3.1 与图 2.3.4。

信号类型谱密度关系相关函数关系
能量信号E_y(f) = |H(f)|² Eₓ(f)
Eₓy(f) = H^*(f)Eₓ(f)
E_(yx)(f) = H(f)Eₓ(f)
R_y(τ) = Rₓ(τ) * h(τ) * h^*(-τ)
Rₓy(τ) = Rₓ(τ) * h^*(-τ)
R_(yx)(τ) = Rₓ(τ) * h(τ)
功率信号P_y(f) = |H(f)|² Pₓ(f)
Pₓy(f) = H^*(f)Pₓ(f)
P_(yx)(f) = H(f)Pₓ(f)
R_y(τ) = Rₓ(τ) * h(τ) * h^*(-τ)
Rₓy(τ) = Rₓ(τ) * h^*(-τ)
R_(yx)(τ) = Rₓ(τ) * h(τ)

图 2.3.4 滤波器输入输出的相关函数与谱密度关系

图 2.3.4 滤波器输入输出的相关函数与谱密度关系(原书影印参考)

实信号 x(t) 通过一个冲激响应为 (1)/(π t) 的滤波器,其输出 hatx(t) 称为 x(t) 的希尔伯特变换(Hilbert Transform):

hatx(t) = mathscrH[x(t)] = x(t) * (1)/(π t) = (1)/(π)∫_(-∞)^(∞) (x(τ))/(t - τ) dτ

希尔伯特变换的传递函数为:

H(f) = -j · operatornamesgn(f)

图 2.3.5 示出了希尔伯特变换的冲激响应与传递函数示意图。

图 2.3.5 希尔伯特变换的冲激响应与传递函数

图 2.3.5 希尔伯特变换的冲激响应与传递函数(原书影印参考)
希尔伯特变换的六大核心性质
  1. 全频带 90° 移相器:传递函数在 f > 0 时为 -j,在 f < 0 时为 +j。即正频率后移 90°,负频率前移 90°。例如 cos(20π t) 经希尔伯特变换后是 sin(20π t);
  2. 级联与反变换:一次希尔伯特变换移相 90°,两次希尔伯特变换移相 180°,也就是信号反相:hathatx(t) = -x(t)。四次变换移相 360° 回到原信号。因此希尔伯特反变换相当于连续做三次希尔伯特变换(移相 270° 或 -90^circ);
  3. 奇偶性反转:偶函数的希尔伯特变换是奇函数,奇函数的希尔伯特变换是偶函数;实偶函数变换为实奇函数,实奇函数变换为实偶函数;
  4. 能量与功率守恒:传递函数的模平方 |-j · operatornamesgn(f)|² = 1,故 hatx(t) 与 x(t) 有完全相同的能量谱密度或功率谱密度,能量与功率相同,自相关函数相同;
  5. 互相关函数:信号的希尔伯特变换与原信号的互相关函数是原信号自相关函数的希尔伯特变换:R_(hatx)x(τ) = hatRₓ(τ);
  6. 正交性:信号与其自身的希尔伯特变换在时域严格相互正交: ∫_(-∞)^(∞) hatx(t)x(t) dt = R_(hatx)x(0) = hatRₓ(0) = 0

实信号的频谱正负频率对称。用如图 2.3.6(a) 所示的滤波器去掉负频率,得到一个只有正频率的信号,这种信号称为解析信号(Analytic Signal)。

图 2.3.6 实信号通过滤波器形成解析信号

图 2.3.6 实信号通过滤波器形成解析信号(原书影印参考)

设有实信号 x(t),其傅氏变换是 X(f)。x(t) 通过滤波器的传递函数为:

H(f) = 1 + operatornamesgn(f) = 2, & f > 0 0, & f < 0

输出信号 z(t) 的傅氏变换为:

Z(f) = X(f)H(f) = 2X(f), & f > 0 0, & f < 0

由于 H(f) = 1 + operatornamesgn(f) = 1 + j[-j · operatornamesgn(f)],于是:

Z(f) = X(f)H(f) = X(f) + j[-j · operatornamesgn(f)X(f)]

两边做傅氏反变换得到:

解析信号的定义公式
z(t) = x(t) + j · hatx(t)

解析信号的虚部是实部的希尔伯特变换。解析信号的傅氏变换是原信号正频率部分的 2 倍。能量谱密度与功率谱密度是正频率部分的 4 倍。

解析信号 z(t) 的共轭为 z^*(t) = x(t) - jhatx(t),对应的传递函数为 1 - operatornamesgn(f),只有负频率。z(t) 与 z^*(t) 没有共同的频率分量,它们的互能量或互功率为零,互能量谱密度或互功率谱密度均为零。

例 2.3.4

设 z(t) = m(t)e^(j2π f_c t) = m(t)cos(2π f_c t) + jm(t)sin(2π f_c t),其中 m(t) 是带宽为 W 的实基带信号,f_c > W。

解

若 m(t) 的频谱为 M(f),则 z(t) 的频谱是 M(f - f_c),如图 2.3.7 所示。从图中可以看出,z(t) 的频谱 Z(f) 只有正频率,说明 z(t) 是解析信号,其实部 m(t)cos(2π f_c t) 的希尔伯特变换正是其虚部 m(t)sin(2π f_c t)。

图 2.3.7 解析信号的频谱示意图

图 2.3.7 解析信号的频谱示意图(原书影印参考)

实信号 m(t) 的频谱 M(f) 正负频率对称。令 M_+(f) 和 M_-(f) 分别表示 M(f) 的正频率部分和负频率部分,其傅氏反变换分别为 m_+(t) 和 m_-(t),则 M(f) = M_+(f) + M_-(f), m(t) = m_+(t) + m_-(t)。m(t) 对应的解析信号 m(t) + jhatm(t) 的频谱是 M(f) 正频率部分的 2 倍,即 2M_+(f),因此:

m_+(t) = 1/2[m(t) + jhatm(t)], m_-(t) = m(t) - m_+(t) = 1/2[m(t) - jhatm(t)] = m_+^*(t)

物理可实现的滤波器必须满足因果性:响应不能发生在激励之前。默认假设滤波器初始状态为零。若在零时刻输入单位冲激 δ(t),滤波器不能在零时刻之前就开始响应。因此,物理可实现滤波器的冲激响应 h(t) 必须满足:

h(t) = 0, t < 0

理想低通滤波器的冲激响应 h_(LPF)(t) = 2Woperatornamesinc(2Wt) 和希尔伯特变换的冲激响应 (1)/(π t) 在 t < 0 处均不为零,都是非因果系统。这些系统在实际实现中需要做因果化处理。

图 2.3.8(a)(b) 分别示出了两个系统的冲激响应,其中系统 1 非因果,系统 2 因果,二者是延迟关系:h₂(t) = h₁(t - t₀)。对于相同的输入 x(t),系统 1 的输出是 y(t),系统 2 的输出是 y(t - t₀),系统 2 等价于系统 1 的输出推迟送出。延迟是无失真关系,因此在无失真意义下,系统 2 和系统 1 是等价的。此时,h₂(t) 就是 h₁(t) 的因果化。我们可以按系统 1 设计通信系统,而真正实现的是系统 2。

图 2.3.8 冲激响应因果化

图 2.3.8 冲激响应因果化(原书影印参考)

对于无限持续时间的非因果系统,如理想低通滤波器的冲激响应 operatornamesinc(t/T) 或希尔伯特变换冲激响应 (1)/(π t),实际实现时可采取近似的方法:先将 h(t) 截短为 h_T(t) = operatornamerect(t/T)h(t),然后对持续时间受限的 h_T(t) 进行因果化延迟,如图 2.3.9 所示。

图 2.3.9 截短后因果化

图 2.3.9 截短后因果化(原书影印参考)

理想滤波器的传递函数在截止频率处直接跳到零,是锐降,如图 2.3.10) 所示。能够实现的实际滤波器则是滚降,传递函数逐渐变到零,存在一个过渡带,如图 2.3.10(a) 所示。能够实现的实际滤波器则是滚降,传递函数逐渐变到零,存在一个过渡带,如图 2.3.10(b) 所示。

图 2.3.10 锐降特性的可实现问题

图 2.3.10 锐降特性的可实现问题(原书影印参考)

实际滤波器是对理想滤波器的近似实现。这种近似能否被接受,取决于通过滤波器的信号频谱在截止频率处有多大:

  • 如果滤波器输入频谱如图 2.3.10(c) 所示,在截止频率处的能量很小,那么通过实际滤波器后的输出与通过理想滤波器后的输出差别不大,可以认为此时的理想滤波器是可实现的;
  • 如果滤波器输入频谱如图 2.3.10(d) 所示,在截止频率处有丰富的能量,那么通过实际滤波器后的输出与通过理想滤波器后的输出相比就会有显著失真,此时理想滤波器不可实现。
共 0 题