深度解析单位冲激响应求解全过程,覆盖连续/离散系统、一阶/二阶系统、稳定性与因果性判断,结合历年考研真题与高频考点,助你精准突破核心难点。
立即掌握求解方法冲击响应是系统对单位冲激信号δ(t)或δ(n)的响应,是表征线性时不变系统(LTI系统)特性的核心函数。掌握其求解方法是信号与系统考研真题得分的关键。
适用于微分方程或差分方程描述的系统,核心在于:设定初始静止条件(零初始状态),求解齐次方程并匹配特解。
系统微分方程:
RC·dy(t)/dt + y(t) = x(t)
令x(t)=δ(t),初始静止:y(0⁻)=0。对方程两边从0⁻到0⁺积分:
RC·[y(0⁺)−y(0⁻)] + ∫₀⁻⁰⁺ y(t)dt = ∫₀⁻⁰⁺ δ(t)dt = 1因y(t)在t=0处有界,第二项→0,故y(0⁺)=1/(RC)。
对t>0,方程变为齐次:RC·y'(t)+y(t)=0,通解:y(t)=A·e⁻ᵗ/ᴿᶜ。
代入初始条件:A = y(0⁺) = 1/(RC),得:
h(t) = (1/RC)·e⁻ᵗ/ᴿᶜ·u(t)其中u(t)为单位阶跃函数,确保因果性。
系统方程:y(n) − 0.6y(n−1) = x(n)
令x(n)=δ(n),零初始状态:y(−1)=0。
计算前几项:
归纳得:h(n) = (0.6)ⁿ·u(n)
通过变换将微分方程转为代数方程,求系统函数H(s)后,再进行拉普拉斯逆变换。
系统微分方程:
d²y(t)/dt² + 5·dy(t)/dt + 6y(t) = 2x(t)
取拉氏变换:
s²Y(s) + 5sY(s) + 6Y(s) = 2X(s)
得系统函数:
H(s) = Y(s)/X(s) = 2/(s² + 5s + 6) = 2/[(s+2)(s+3)]
部分分式分解:
H(s) = A/(s+2) + B/(s+3)
解得:A=2,B=−2,故:
H(s) = 2/(s+2) − 2/(s+3)
查表得冲激响应:
h(t) = [2e⁻²ᵗ − 2e⁻³ᵗ]·u(t)
本例极点为s=−2、s=−3(左半平面实极点),故h(t)为衰减响应,系统稳定。
与拉普拉斯法原理一致,将差分方程转为Z域代数方程。
系统差分方程:
y(n) − 0.5y(n−1) + 0.06y(n−2) = x(n)
取Z变换:
Y(z) − 0.5z⁻¹Y(z) + 0.06z⁻²Y(z) = X(z)
系统函数:
H(z) = Y(z)/X(z) = 1 / (1 − 0.5z⁻¹ + 0.06z⁻²)
整理为有理分式:
H(z) = z² / (z² − 0.5z + 0.06) = z² / [(z−0.2)(z−0.3)]
为便于分解,写为H(z)/z:
H(z)/z = z / [(z−0.2)(z−0.3)] = A/(z−0.2) + B/(z−0.3)
解得:A=2,B=−1,故:
H(z) = 2z/(z−0.2) − z/(z−0.3)
查Z变换表:z/(z−a) ↔ aⁿu(n),得:
h(n) = [2·(0.2)ⁿ − (0.3)ⁿ]·u(n)
极点z=0.2、z=0.3均在单位圆内(|z|<1),故系统稳定,h(n)随n增长而衰减至零。
适用于一阶、二阶系统,尤其适合编程实现或手算前几项验证解析解。
系统方程:y(n) = 0.8y(n−1) + x(n)
设x(n)=δ(n),初始静止:y(−1)=0。
| n | x(n) | y(n)计算 | h(n) |
|---|---|---|---|
| −1 | − | y(−1)=0 | − |
| 0 | δ(0)=1 | 0.8×0 + 1 = 1 | 1 |
| 1 | 0 | 0.8×1 + 0 = 0.8 | 0.8 |
| 2 | 0 | 0.8×0.8 = 0.64 | 0.64 |
| 3 | 0 | 0.8×0.64 = 0.512 | 0.512 |
观察规律:h(n) = (0.8)ⁿ·u(n),与解析法一致。
基于近十年全国重点高校考研真题(如清华、浙大、电子科大、西电等),归纳为四大高频题型,覆盖95%以上考题。
常以RC/RL电路、一阶机械系统为背景,考察时域直接法或拉氏变换法。
真题示例:电子科大2021年
已知系统微分方程:2y'(t)+3y(t)=δ(t),求h(t)。
解析:齐次解y_h(t)=Ae⁻³ᵗ/²,积分得y(0⁺)=1/2,故A=1/2,
h(t) = (1/2)e⁻³ᵗ/²·u(t)
高频考点!涉及欠阻尼、临界阻尼、过阻尼三种形态,极点位置决定响应形式。
真题示例:浙大2022年
系统函数H(s)=ωₙ²/(s²+2ζωₙs+ωₙ²),当ζ=0.5时,求h(t)。
解析:共轭复极点s=−ζωₙ±jωₙ√(1−ζ²)=−0.5ωₙ±j(√3/2)ωₙ,
h(t) = ωₙe⁻⁰·⁵ωₙᵗ·sin(√3/2 ωₙ t)·u(t)
差分方程求解,极点在Z平面位置决定稳定性。
真题示例:西电2020年
差分方程:y(n)−y(n−1)+0.25y(n−2)=x(n),求h(n)。
解析:特征根s=0.5(二重根),
h(n) = (A + Bn)(0.5)ⁿ·u(n),代入初始条件得A=1, B=0.5,
h(n) = (1 + 0.5n)(0.5)ⁿ·u(n)
通过h(t)或h(n)判断:因果性、稳定性、线性性、时不变性。
真题示例:清华2023年
已知h(t)=eᵗ·u(−t),判断系统因果性与稳定性。
解析:
• 因果性:当t>0时h(t)=0,但t<0时非零 → 非因果
• 稳定性:∫₋∞⁰ eᵗ dt = 1 < ∞ → 稳定
冲击响应不仅是数学表达式,更是系统物理本质的体现。深入理解其特性,是解决复杂问题的关键。
对线性时不变系统(LTI),其冲击响应具有以下性质:
反例:若系统含时变元件(如可变电阻),则h(t,τ)依赖于延迟量τ,不再是单一函数,系统非时不变。
核心在于:δ(t)的积分突变特性。对微分方程两边从0⁻到0⁺积分:
∫₀⁻⁰⁺ [aₙy⁽ⁿ⁾(t) + ... + a₀y(t)] dt = ∫₀⁻⁰⁺ [bₘx⁽ᵐ⁾(t) + ... + b₀x(t)] dt由于y(t)及其导数在t=0处最多有有限跳变(除最高阶导数外),左边仅aₙ[y⁽ⁿ⁻¹⁾(0⁺)−y⁽ⁿ⁻¹⁾(0⁻)]非零;右边若m
关键结论:对δ(t)输入,y(t)的(n−1)阶导数发生跳变,而y(t)本身连续。
应用示例:二阶系统y''+3y'+2y=δ(t),可得:
结合初始静止条件y(0⁻)=0, y'(0⁻)=0,得y(0⁺)=0, y'(0⁺)=1,用于求解t>0的齐次解。
在拉普拉斯域,H(s)是h(t)的变换;在Z域,H(z)是h(n)的变换。二者完全等价,但各有优势:
重要关系:
• 稳定性:连续系统H(s)所有极点在左半平面;离散系统H(z)所有极点在单位圆内;
• 因果性:h(t)=0 (t<0) ⇔ H(s)为真有理分式;h(n)=0 (n<0) ⇔ H(z)为真有理分式。
者通过积分/微分关联:
应用价值:当直接求h(t)困难时(如含高阶导数),可先求阶跃响应s(t),再求导。
真题技巧:若系统方程含x'(t),直接求h(t)会涉及δ'(t),此时先求s(t)更简便。
掌握求解只是起点,灵活应用才是高分关键。冲击响应是连接理论与工程的桥梁。
充要条件:∫|h(t)|dt<∞(BIBO稳定)
快速判断法: