Skip to content

卷积と线性时不变系统

Info

Linear Time Invariant 线性时不变

Abstract

基本信号应具有以下特点:

  • 本身尽可能简单,并且用它们的线性组合能够表示尽可能广泛的信号
  • LTI 系统对这种信号的响应易于计算且有利于对系统的深入理解

两类有用的基本信号:

  • 复指数信号
  • 单位冲激信号

本章主要考虑以单位冲激信号为基本信号进行信号与系统分析

summary

卷积和

定义

对于任意一个离散时间系统:

离散时间系统

设离散时间系统为 LTI 的,定义单位脉冲响应 h[n] 为以单位冲激信号 δ[n] 为输入时的输出:

δ[n]h[n]

由 LTI 时不变性有:

δ[nk]h[nk]

由 LTI 线性有:

x[n]=k=+x[k]δ[nk]y[n]=k=+x[k]h[nk]

这种求得系统响应的运算关系称为卷积和,即离散时间 LTI 系统的响应可以用卷积和表示为:

y[n]=x[n]h[n]=k=+x[k]h[nk]

⭐更一般的形式:

x[nn0n1]=x[nn0]h[nn1]

⭐只要得到了 LTI 系统对 δ[n] 的 响 应 h[n] ——单位脉冲响应,就可以得到 LTI 系统对任何输入信号 x[n] 的响应。

换言之,一个 LTI 系统可以完全由它的单位脉冲响应 h[n] 来表征。

特殊离散 LTI 系统

延时器

h[n]=δ[nn0]

具有该单位脉冲响应的系统一般情况下存在多种可能性,但是对应的 LTI 系统只有一个:

y[n]=x[nn0]=x[n]δ[nn0]

累加器

h[n]=k=nδ[k]=u[n]

对应的 LTI 系统:

y[n]=x[n]u[n]=k=nx[k]

计算

选择 n 值并固定,求和项视为 k 的函数,步骤如下:

  1. x[n]h[n] 变为 x[k]h[k] ,以 k 为求和变量
  2. h[k] 翻转变为 h[k]
  3. h[k] 平移 n,变为 h[nk]
  4. x[k]h[nk] 相乘;
  5. 对乘积 x[k]h[nk]kZ 取和;
  6. 一般都转化为只对非 0 区间积分。注意非 0 区间上下限!!!
  7. 改变 n 值并重复 3–5 步,直到 n 超过了 x[n] 非零部分定义域

Warning

注意很多题目最后答案会加上一个 u(n) ,是因为只求定义域在 [0,+) 的部分

卷积积分

定义

Abstract

类似积分思想,转化为无限多的离散和的极限

用冲激信号表示连续时间信号:

将输入信号 x(t) 近似表示为移位、放缩后的脉冲之和

脉冲和近似表示输入信号

xΔ(t)=x(kΔ) , kΔ<t<(k+1)Δ

阶梯近似表达过程:

设连续时间系统为 LTI 的,定义(名字待补,与卷积和中的单位脉冲响应对应) hΔ(n) 为以单位冲激信号 δΔ(n) 为输入时的输出:

δΔ(n)hΔ(n)

阶梯近似表达

由 LTI 时不变性有:

δΔ(nkΔ)hΔ(nkΔ)

由 LTI 线性齐次性尺度调整:

x(kΔ)δΔ(nkΔ)x(kΔ)hΔ(nkΔ)

阶梯近似表达2

由 LTI 线性叠加性有:

k=+x(kΔ)δΔ(nkΔ)Δk=+x(kΔ)hΔ(nkΔ)ΔxΔ(t)yΔ(t)

当 Δ → 0 时 :

  • xΔ(t)x(t)
  • yΔ(t)y(t)
  • kΔτ
  • Δdτ
  • δΔ(tkΔ)δ(tτ)
  • hΔ(tkΔ)h(tτ)

从而得到:

x(t)=+x(τ)δ(tτ)dτ=x(t)δ(t)y(t)=+x(τ)h(tτ)dτ=x(t)h(t)
  • 其中第一个式子意义是:用位移的单位冲激信号表示任意的连续时间信号
  • 其中第一个式子意义是:用卷积积分表示连续时间 LTI 系统的响应

阶梯近似图解过程:

阶梯近似图解过程

计算

求卷积积分 x(t)h(t)=+x(τ)h(tτ)dτ 的过程与卷积和计算过程相似。卷积积分的计算步骤如下:

  1. x(t)h(t) 变为 x(τ)h(τ) ,以 τ 为积分变量;
  2. h(τ) 翻转变为 h(τ)
  3. h(τ) 平移 t ,变为 h(tτ)
  4. x(τ)h(tτ) 相乘;
  5. 求区间 (,+) 上乘积 x(τ)h(tτ) 下的面积;
  6. 一般都转化为只对非 0 区间积分。注意非 0 区间上下限!!!
  7. 改变 t 值并重复 3–5 步。

特殊连续 LTI 系统

积分器

单位冲激信号为输入下系统的响应为:

h(t)=u(t)

对任意激励信号 x(t) ,系统响应为:

y(t)=x(t)h(t)=+x(τ)h(tτ)dτ=tx(τ)dτ

卷积和与卷积积分的性质

交换律

变量替换:

x[n]h[n]=k=+x[k]h[nk]=l=+x[nl]h[l]=h[n]x[n]
x(n)h(n)=+x(τ)h(tτ)dτ=+x(tτ)h(τ)dτ=h(t)x(t)

交换律

s[n]=u[n]h[n]=h[n]u[n]=k=nh[k]

利用 u[n] 是累加器的单位脉冲响应以及交换律可以得到:

该系统的单位脉冲响应是 h[n]u[n] 单位阶跃信号作为输入,则其单位阶跃响应 s[n] 等价于 以 h[n] 作为输入时,累加器的响应

分配律

x[n](h1[n]+h2[n])=x[n]h1[n]+x[n]h2[n]x(t)(h1(t)+h2(t))=x(t)h1(t)+x(t)h2(t)

分配律

结合律

(x[n]h1[n])h2[n]=x[n](h1[n]h2[n])(x(t)h1(t))h2(t)=x(t)(h1(t)h2(t))

结合律

级联的 LTI 系统其总的单位冲激(脉冲)响应等于各子系统单位冲激(脉冲)响应的卷积

级联可交换

卷积运算满足结合律和交换律,因此,系统级联的先后次序可以调换。

x[n]h1[n]h2[n]=x[n]h2[n]h1[n]x(t)h1(t)h2(t)=x(t)h2(t)h1(t)

级联可交换

级联次序可交换的前提条件:

  • 系统是 LTI 的
  • 所有涉及到的卷积运算必须收敛

微分、积分及时移特性

卷积积分满足微分、积分及时移特性。若 x(t)h(t)=y(t),则:

微分特性

x(t)h(t)=x(t)h(t)=y(t)
x(m)(t)h(n)(t)=y(m+n)(t)

Warning

t 求导而不是对 τ 求导

积分特性

(tx(τ)dτ)h(t)=+(τ1x(τ)dτ)h(tτ1)dτ1=x(t)(th(τ)dτ)=ty(τ)dτ

时移特性

x(tt0)h(tt1)=x(tt1)h(tt0)=y(tt0t1)

差分、求和及时移特性

卷积和满足差分、求和及时移特性。若 x[n] * h[n] = y[n],则:

差分特性

x[n]x[n1]h[n]=x[n](h[n]h[n1])=y[n]y[n1]

求和特性

(k=nx[k])h[n]=x[n](k=nh[k])=k=ny[k]

时移特性

x[tt0]h[tt1]=x[tt1]h[tt0]=y[tt0t1]

LTI 系统的性质

记忆性

对于离散时间 LTI 系统:

y[n]=k=+x[k]h[nk]

如果系统是无记忆的,则在任何时刻 ny[n] 都只能和 n 时刻的输入有关,等式右边的和式中只能有 k=n 时的一项非零,即必须有 h[nk]=0kn,也即:

h[n]=0,   n0

所以,离散时间无记忆系统的单位脉冲响应为:

h[n]=kδ[n]

即为离散时间系统无记忆充要条件

类似地,连续时间系统无记忆的充要条件为:

h(t)=kδ(t)

k=1 时系统是恒等系统。

可逆性

如果 LTI 系统是可逆的,则其逆系统也是 LTI 系统。

互逆的两个系统的满足(充要条件卷积和等于单位冲激信号

h(t)g(t)=δ(t)h[n]g[n]=δ[n]

从而有:

x(t)h(t)g(t)=x(t)δ(t)=x(t)x[t]h[n]g[n]=x[t]δ[n]=x[t]

可逆性

例如:延时器是可逆的 LTI 系统,h(t) = \delta(t - t0),其逆系统是 g(t) = \delta(t + t0),显然有:

h(t)g(t)=(tt0)δ(t+t0)=δ(t)

又如:累加器是可逆的 LTI 系统,其 h[n] = u[n],其逆系统是 g[n] = \delta[n] - \delta[n - 1],显然也有:

h[n]g[n]=u[n](δ[n]δ[n1])=u[n]u[n1]=δ[n]

因果性

对于离散时间 LTI 系统:

y[n]=k=+x[k]h[nk]

因果性意味着上式等号右边和式中所有 k > n 的项都必须为零,即 h[n - k] = 0,k > n,也即:

h[n]=0,   n<0

同理对连续时间因果系统有:

h(t)=0,   t<0

即为 LTI 系统具有因果性的充要条件

对于线性系统而言,因果性与初始松弛条件等价。所谓系统满足初始松弛条件是指该系统具有如下表现:对于任意时刻,只要该时刻之前系统的输入为 0,则该时刻之前系统的输出也必然为 0。

(参考教材习题 1.44(a))

– 该等价性的前提是系统为线性的。否则不一定成立。例如 y n = 2x n + 3 是因果的,但显然不满足初始松弛条件。

稳定性

稳定性意味着有界输入必然产生有界输出。以离散时间 LTI 系统为例,设 |x[n]| 有界,| x[n - k]| ≤ A。有:

|y[n]|=|k=+h[k]x[nk]|k=+|h[k]||x[nk]|Ak=+|h[k]|

离散 LTI 系统稳定的充分必要条件为:

n=+|h[n]|<+

类似地有连续 LTI 系统稳定的充分必要条件为:

+|h(t)|dt<+

单位阶跃响应

在工程实际中,也常用单位阶跃响应来描述 LTI 系统。单位阶跃响应就是系统对 u(t) 或 u[n] 所产生的响应,记为 s(t)或 s[n]。因此有:

s[n]=u[n]h[n]=k=nh[k]
h[n]=s[n]s[n1]
s(t)=u(t)h(t)=th(τ)dτ
h(t)=ds(t)dt=s(t)

用微分和差分方程描述的因果 LTI 系统

在工程实际中有相当普遍的一类系统,其数学模型可以用线性常系数微分方程或线性常系数差分方程来描述。

在时域内分析这类 LTI 系统,就是要求解线性常系数微分方程或差分方程。

线性常系数微分方程

一般情况 NM , N 阶非其次常系数微分方程的一般形式:

k=0Nakdky(t)dtk=k=0Mbkdky(t)dtk

其中 ak, bk 均为常数

求解该微分方程(求 y( t ) 的显式表达)的步骤:

  1. 求出齐次解通解)yh(t) ,也称为系统的自然响应自由响应
  2. 建立特征方程求出特征根 𝜆k
  3. 根据特征根得到齐次解形式
  4. 根据附加条件确定齐次解形式中的待定常数 Ck
  5. 特解 y𝑝(t),也称系统的受迫响应强迫响应,其形式与输入信号有关
  6. 根据输入信号形式得到特解形式
  7. 将输入信号与特解代入微分方程确定特解形式中的待定常数
  8. 得到完全解 y(t) = yh(t) + y𝑝(t)

Warning

附加条件常给出为初始松弛条件,即 t<0 时若 x(t)=0 则有 y(t)=0 ,或 n<0 时若 x[n]=0 则有 y[n]=0

对于线性系统,初始松弛条件与满足因果性是等价的

解题时,可以直接认为就是有 y(0)=0x[0]=0

线性常系数差分方程

一般情况 NMN 阶非其次常系数差分方程的一般形式

k=0Naky[nk]=k=0Mbkx[nk]

其中 ak, bk 均为常数

还可以将其改写为(正向)递归形式

y[n]=1a0(k=0Mbkx[nk]k=0Naky[nk])

还可将差分方程改写反向递归求解形式

y[nN]=1aN(k=0Mbkx[nk]k=0N1aky[nk])

ak=0,k0 时,差分方程变为一个卷积和的形式即非递归方程

y[n]=k=0Mbka0x[nk]=+h[n] x[nk]

其中:

h[n]={bna0,  0nM0,    others
  • 由于 h[n] 是有限长的,当 ak=0,k0非递归方程描述的系统称为有限脉冲响应( Finite Impulse Response,FIR)系统
  • 而其他情况下递归方程描述的系统结合初始松弛条件具有无限长的单位脉冲响应,因而称为无限脉冲响应(Infinite Impulse Response,IIR)系统

  • 与微分方程类似,差分方程的齐次解通解)与输入信号的形式无关,称为系统的自然响应自由响应

  • 特解具有与输入信号相同的函数形式,称为受迫响应强迫响应

具有非零起始状态的线性常系数微分或差分方程对应增量线性系统,其响应由两部分组成:

  • 零输入响应,与输入信号无关,属于自然响应
  • 零状态响应,既与输入信号有关,也与系统特性有关,它包含了一部分受迫响应,也包含有一部分自然响应

方框图表示

差分方程描述系统

对于正向递归形式的 LTI 系统差分方程,其方程中包括三种基本运算:

  • 移位(延迟)——对应延迟器
  • 乘系数——对应数乘器
  • 相加——对应加法器

令:

w[n]=k=0Mbkx[nk]

正向递归形式可写为:

y[n]=1a0(w[n]k=0Naky[nk])

差分-直接I型

将其级联起来,就成为线性常系数差分方程描述的系统,它具有与差分方程完全相同的运算功能。可以调换级联次序,并将移位单元合并,于是得到直接 II 型

差分-直接II型

微分方程描述系统

微分方程也包括三种基本运算:

  • 微分
  • 乘系数
  • 相加

但微分器不仅在工程实现上有困难,而且对误差及噪声极为敏感,因此,通常使用积分器而不用微分器

N 阶微分方程两边同时积分 N 次,即可得到对应 N 阶积分方程:

k=0Nak y(Nk)(t)=k=0Mbk x(Nk)(t)

对此积分方程完全按照差分方程的办法有其递归方程:

y(t)=1a0(k=0Mbkx(Nk)(t)k=0Naky(Nk)(t))

令:

w(t)=k=0Mbkx(Nk)(t)

可得:

y(t)=1a0(w(t)k=0Naky(Nk)(t))

积分-直接I型

通过交换级联次序,合并积分器可得直接 Ⅱ 型

积分-直接II型

卷积定义奇异函数

单位冲激信号u0(t)

Abstract

δ(t)=ddtu(t) 定义下的单位冲激。这是不严密的,因为 u(t) 在 t=0 处不连续。

进而采用极限的观点,将 δ(t) 视为 δΔ(t)Δ0 时的极限。但这种定义或描述\delta(t)的方法在数学上仍然是不严格的,因为存在无穷多个不同的函数在 Δ0 时都表现为与δΔ(t) 有相同的特性。

极限定义举例1

极限定义举例2

上述现象出现的原因在于\delta(t)是一个理想化的非常规函数,被称为奇异函数,通常采用卷积或积分运算下表现的特性来定义奇异函数

从系统的角度,将 δ(t) 定义为恒等系统单位冲激响应,也就是通过卷积定义 δ(t) ,对任意 x(t)δ(t) 应满足:

x(t)=x(t)δ(t)

::: primary 单位冲激函数相当于卷积运算系统的单位元 根据定义可以得出 δ(t) 的如下性质:

1=1δ(t)=+δ(τ)dτ

卷积定义还有一个等价版本——积分定义

g(t)=g(t)δ(t)=+g(tτ)δ(τ)dτg(0)=+g(t)δ(τ)dτ

偶函数

单位冲激信号 δ(t) 是偶函数:

δ(t)=δ(t)

证明:

+δ(t)f(t)dt=+δ(τ)f(τ)d(τ)=+δ(τ)f(τ)dτ=f(0)

时间尺度变换

δ(at)=1|a|δ(t)

验证信号下方面积:

a>0 时 :

+δ(aτ)dτ=1a+δ(aτ)d(aτ)=1a=1|a|

a<0 时 :

+δ(aτ)dτ=1a+δ(aτ)d(aτ)=1a+δ(t)dt=1a+δ(t)dt=1a=|1a|

单位冲激偶u1(t)

系统角度定义为理想微分器单位冲激响应,即 δ(t) 的微分:

u1(t)=ddtδ(t)

微分器

单位冲激偶

卷积定义,从 LTI 系统分析的角度有:

x(t)u1(t)=ddtx(t)

该信号可视为短脉冲求导的理想化结果。

短脉冲求导的理想化

积分性质

令输入 x(t)=1 有:

+u1(τ)dτ=+x(tτ)u1(τ)dτ=x(t)u1(t)=ddtx(t)=0

积分意义下的定义

+g(τt)u1(τ)dτ=g(t)u1(t)=ddtg(t)

t=0 有:

g(0)=+g(τ)u1(τ)dτ

该积分即可作为 u1(t) 在积分意义下的定义

δ(t)的积分u1(t)

从系统的角度,定义 u1(t)理想积分器单位冲激响应

u1(t)=tδ(τ)dτ=u(t)

运算定义为:

x(t)u1(t)=tx(τ)dτ

单位斜坡函数u2(t)

系统的角度,可以定义为:(卷积)

u2(t)=u1(t)u1(t)=tu(t)dt=tu(t)

运算定义为二重积分:

x(t)u2(t)=x(t)u2(t)=x(t)u1(t)u1(t)=(tu(σ)dσ)u1(t)=t(τu(σ)dσ)dτ=tτu(σ) dσdτ

卷积定义推广

例如 u2(t) 可以定义为:u2(t)=u1(t)u1(t) :

x(t)u2(t)=d2dt2x(t)=x(t)u1(t)u1(t)

以此类推:

u2(t)=u1(t)u1(t)u2(t)=u1(t)u1(t)u1(t)......uk(t)=u1(t)u1(t)...u1(t)k个

单位阶跃函数的高阶积分推广:

u3(t)=u1(t)u1(t)u1(t)=12t2u(t)u4(t)=u1(t)u1(t)u1(t)u1(t)=16t3u(t)......uk(t)=u1(t)u1(t)...u1(t)k个=tk1(k1)!u(t)

Info

实际上, \delta(t)的各次积分已经是常规函数了