ARTICLE

单位冲激响应

单位冲激响应 单位冲激响应(Unit Impulse Response,简称冲激响应)是线性时不变系统 (LTI系统)在受到单位冲激函数(Dirac Delta 函数或 Kronecker Delta 序列)激励时的输出响应。它构成了连续时间和离散时间信号与系统理论中最核心的概念之一,因为一个 LTI 系统的全部动态特性——包括稳定性、因果性、频率响应和传递

浏览 0 更新 2025-07-11

单位冲激响应

单位冲激响应(Unit Impulse Response,简称冲激响应)是线性时不变系统 (LTI系统)在受到单位冲激函数(Dirac Delta 函数或 Kronecker Delta 序列)激励时的输出响应。它构成了连续时间和离散时间信号与系统理论中最核心的概念之一,因为一个 LTI 系统的全部动态特性——包括稳定性、因果性、频率响应和传递函数——都可由其单位冲激响应唯一刻画。从控制系统到通信工程,从数字信号处理到经济时间序列分析,单位冲激响应无处不在。

连续时间系统中的定义

在连续时间域,单位冲激响应 h(t)h(t) 定义为:当系统输入为单位冲激函数 δ(t)\delta(t)、初始状态为零时,系统的输出。数学上写作:

δ(t)LTI系统h(t)\delta(t) \longrightarrow \text{LTI系统} \longrightarrow h(t)

其中 Dirac Delta 函数 δ(t)\delta(t) 满足 δ(t)dt=1\int_{-\infty}^{\infty} \delta(t)\,dt = 1,且在 t0t \neq 0 处为零,严格而言是广义函数或分布意义下的函数。

单位冲激响应的根本重要性源自 LTI 系统的叠加性质(线性性)和时不变性。对于任意输入信号 x(t)x(t),可将其表示为无穷多个加权、时移冲激函数的叠加(即冲激分解):

x(t)=x(τ)δ(tτ)dτx(t) = \int_{-\infty}^{\infty} x(\tau)\delta(t-\tau)\,d\tau

根据线性时不变性,系统对 x(t)x(t) 的响应 y(t)y(t)卷积积分(Convolution Integral)给出:

y(t)=x(τ)h(tτ)dτ=(xh)(t)y(t) = \int_{-\infty}^{\infty} x(\tau)h(t-\tau)\,d\tau = (x * h)(t)

这一简单而优美的卷积关系意味着:一旦知道了系统的单位冲激响应 h(t)h(t),就能计算任意输入下的输出。这正是系统分析中"时域完全描述"的含义。

离散时间系统中的定义

在离散时间域中,单位冲激响应 h[n]h[n] 是系统对单位脉冲序列 δ[n]\delta[n] 的零状态响应。离散单位脉冲定义为:

\delta[n] = \begin{cases}

1, \& n = 0 \\ 0, \& n \neq 0

\end{cases}

离散 LTI 系统的输入-输出关系由卷积和(Convolution Sum)描述:

y[n]=k=x[k]h[nk]=(xh)[n]y[n] = \sum_{k=-\infty}^{\infty} x[k]\,h[n-k] = (x * h)[n]

离散冲激响应在数字滤波器设计、ARMA模型(自回归滑动平均模型)和有限脉冲响应 (FIR)无限脉冲响应 (IIR)滤波器的分类中扮演着关键角色:若 h[n]h[n] 在有限区间外为零,则为 FIR 滤波器,具有固有稳定性;若 h[n]h[n] 无限延伸,则为 IIR 滤波器,通常需要更少的系数但可能面临稳定性问题。

与传递函数和频率响应的关系

对单位冲激响应 h(t)h(t) 进行拉普拉斯变换,即得系统的传递函数

H(s)=L{h(t)}=0h(t)estdtH(s) = \mathcal{L}\{h(t)\} = \int_{0}^{\infty} h(t) e^{-st}\,dt

类似地,对 h(t)h(t) 进行傅里叶变换,即得系统的频率响应

H(jω)=F{h(t)}=h(t)ejωtdtH(j\omega) = \mathcal{F}\{h(t)\} = \int_{-\infty}^{\infty} h(t) e^{-j\omega t}\,dt

在离散系统中,Z 变换h[n]h[n] 映射为系统的传递函数 H(z)=Z{h[n]}H(z) = \mathcal{Z}\{h[n]\}。因此,时域的冲激响应与频域的系统函数互为傅里叶/拉普拉斯变换对——这一对偶关系构成了经典系统分析的理论支柱。事实上,拉普拉斯变换将时域的卷积运算简化为频域的代数乘法:Y(s)=H(s)X(s)Y(s) = H(s)X(s),这极大地简化了系统分析和综合的数学操作。

因果性与稳定性判据

单位冲激响应直接蕴含了两个关键的系统特性:

因果性(Causality):一个 LTI 系统是因果的当且仅当其冲激响应在负时间上为零,即 h(t)=0, t<0h(t) = 0,\ \forall t < 0(连续)或 h[n]=0, n<0h[n] = 0,\ \forall n < 0(离散)。这意味着系统的输出只依赖于当前和过去的输入,而不依赖于未来的输入——这是所有物理可实现系统的必要条件。

BIBO 稳定性(有界输入有界输出稳定性):一个 LTI 系统是 BIBO 稳定的当且仅当其冲激响应绝对可积(连续)或绝对可和(离散):

h(t)dt<n=h[n]<\int_{-\infty}^{\infty} |h(t)|\,dt < \infty \quad\text{或}\quad \sum_{n=-\infty}^{\infty} |h[n]| < \infty

这是因为若冲激响应绝对可积,则对于任意有界输入 x(t)M|x(t)| \leq M,输出满足 y(t)Mh(τ)dτ<|y(t)| \leq M \int |h(\tau)|d\tau < \infty,从而输出也有界。这一简洁的充要条件使得工程师仅通过检查冲激响应即可判断系统的稳定性。

经济时间序列中的脉冲响应分析

单位冲激响应的思想已深入经济学和计量经济学。在向量自回归 (VAR)模型中,研究者常计算脉冲响应函数(Impulse Response Function, IRF),即给某个结构冲击(如货币政策的一个标准差正冲击)施加一个单位脉冲后,各内生变量在后续各期的动态响应路径。这本质上正是将单位冲激响应的概念移植到随机经济系统中。通过Cholesky 分解符号约束等方法识别结构冲击后,脉冲响应函数直观展示了经济变量之间的短期和长期传导机制——例如,一个紧缩性货币政策冲击如何通过利率渠道传导至投资和产出,以及这种效应的衰减速度。这与控制系统工程中通过单位冲激响应分析系统动态特性的思路如出一辙。

小结

单位冲激响应是 LTI 系统理论中最精炼而完备的时域描述工具。它直接对应系统的传递函数和频率响应,是卷积运算的基础,也是因果性和 BIBO 稳定性的直接判据。无论是在连续时间系统(电路、机械振动)还是离散时间系统(数字滤波器、经济时间序列)中,理解和计算单位冲激响应都是分析动态系统行为的起点。

常见系统的单位冲激响应示例

不同物理系统呈现出形态各异的单位冲激响应。一阶 RC 低通电路的冲激响应为指数衰减函数 h(t)=(1/RC)et/RCu(t)h(t) = (1/RC)e^{-t/RC}u(t),其中衰减时间常数 RCRC 决定了系统记忆的长度。二阶 RLC 电路的冲激响应根据阻尼比 ζ\zeta 的不同可表现为欠阻尼(振荡衰减)、临界阻尼(最快无振荡衰减)和过阻尼(缓慢指数衰减)三种形态。在力学系统中,质量-弹簧-阻尼器系统的冲激响应同样由二阶微分方程控制,其冲激响应直接反映了系统的固有频率和阻尼特性。在数字信号处理中,移动平均(MA)滤波器的冲激响应是有限长的系数序列,而自回归(AR)滤波器的冲激响应则无限延伸并呈指数衰减——这正是 FIR 和 IIR 滤波器在时域上的本质区别。

格林函数与偏微分方程

单位冲激响应的思想在偏微分方程理论中以格林函数(Green's Function)的形式出现。对于线性偏微分算子 L\mathcal{L},格林函数 G(r,r)G(\mathbf{r}, \mathbf{r}') 定义为方程 LG(r,r)=δ(rr)\mathcal{L}G(\mathbf{r}, \mathbf{r}') = \delta(\mathbf{r} - \mathbf{r}') 在齐次边界条件下的解。一旦求得格林函数,任意源项 f(r)f(\mathbf{r}) 对应的解可表示为叠加积分 u(r)=G(r,r)f(r)dru(\mathbf{r}) = \int G(\mathbf{r}, \mathbf{r}') f(\mathbf{r}')\,d\mathbf{r}'。这一框架在电磁学、热传导、量子力学和流体力学中均有广泛应用——格林函数本质上就是多维空间中的单位冲激响应。