信号系统避坑指南：LTI连续系统初始值跃变到底怎么算？（附经典例题解析）

LTI连续系统初始值跃变问题深度解析：从理论到实战的完整指南

在信号与系统课程中，LTI连续系统的时域分析是一个让许多学生感到困惑的难点，尤其是当系统响应在初始时刻出现跃变时。这种困惑往往源于对微分方程右端含有冲激函数δ(t)时，系统响应及其各阶导数在0+时刻如何变化的机制理解不够透彻。本文将通过一个典型二阶系统案例，带你一步步拆解这个"黑箱"过程。

1. 理解初始值问题的本质

初始值问题之所以成为学习难点，是因为它涉及三个关键概念的区分与联系：

• 初始状态（y(j)(0-)）：系统在激励接入前的"历史记忆"，完全由系统自身特性决定
• 初始值（y(j)(0+)）：激励接入后系统响应的起始条件
• 跃变现象：当微分方程右端含有δ(t)时，响应导数在0-到0+时刻的突变

经典误区警示：

1. 混淆0-和0+时刻的值，直接套用初始状态作为求解条件
2. 忽视微分方程右端是否含有δ(t)，错误判断跃变发生的位置
3. 在零输入和零状态响应中采用相同的初始值处理方法

关键规律：当且仅当微分方程右端含有δ(t)时，系统响应的高阶导数才可能出现跃变，且跃变遵循"最高阶导数含δ(t)→次高阶跃变→其余连续"的连锁反应规则。

2. 典型二阶系统案例解析

考虑如下RLC电路系统微分方程：

y''(t) + 3y'(t) + 2y(t) = δ'(t) + 2δ(t) + ε(t)

已知初始状态：y(0-)=1, y'(0-)=2

2.1 判断跃变发生的条件

首先分析方程右端是否含有冲激函数及其导数：

跃变判定流程：

1. 确定方程右端最高阶冲激项：δ'(t)（二阶导级别）
2. 根据跃变传递规则：
   • y''(t) 将包含δ'(t)
   • y'(t) 将包含δ(t)并有跃变
   • y(t) 将连续无跃变

2.2 计算0+时刻初始值

通过积分法求解跃变量：

1. y(t)连续性：

   y(0+) - y(0-) = \int_{0-}^{0+} y'(τ)dτ = 0 \quad ⇒ \quad y(0+) = y(0-) = 1

2. y'(t)跃变计算： 对原方程在[0-,0+]积分：

   \int_{0-}^{0+} [y'' + 3y' + 2y] dt = \int_{0-}^{0+} [δ'(t) + 2δ(t) + ε(t)] dt

   逐项分析：

   • ∫y''dt = y'(0+) - y'(0-)
   • ∫3y'dt = 0 （有限量积分）
   • ∫2ydt = 0 （有限量积分）
   • ∫δ'(t)dt = 0 （δ'积分为0）
   • ∫2δ(t)dt = 2
   • ∫ε(t)dt = 0

   得到方程：

   y'(0+) - y'(0-) = 2 ⇒ y'(0+) = 2 + 2 = 4

3. y''(t)中的冲激： 虽然不需要计算y''(0+)，但需知道y''(t)包含δ'(t)项，这会影响后续全响应的求解。

2.3 零输入与零状态响应的初始值差异

零输入响应：

• 仅由初始状态引起
• 满足齐次方程，无跃变：

  y_zi(0+) = y(0-) = 1 y_zi'(0+) = y'(0-) = 2

零状态响应：

• 仅由激励引起
• 初始状态为零，但可能有跃变：

  y_zs(0+) = 0 （连续） y_zs'(0+) = 2 （来自跃变计算）

验证：

y(0+) = y_zi(0+) + y_zs(0+) = 1 + 0 = 1 ✔ y'(0+) = y_zi'(0+) + y_zs'(0+) = 2 + 2 = 4 ✔

3. 系统响应求解全流程

3.1 零输入响应求解

方程为：y'' + 3y' + 2y = 0

1. 特征方程：

   r² + 3r + 2 = 0 ⇒ r = -1, -2

2. 通解形式：

   y_zi(t) = C1e^{-t} + C2e^{-2t}

3. 代入初始值：

   y_zi(0+) = C1 + C2 = 1 y_zi'(0+) = -C1 -2C2 = 2

   解得：

   C1 = 4, C2 = -3

   因此：

   y_zi(t) = 4e^{-t} - 3e^{-2t}, t ≥ 0

3.2 零状态响应求解

方程为：y'' + 3y' + 2y = δ'(t) + 2δ(t) + ε(t)

1. 齐次解同前：yh(t) = A1e^{-t} + A2e^{-2t}
2. 特解形式：
   • 对于t>0，方程变为：y'' + 3y' + 2y = 1
   • 设特解为常数B，代入得：2B = 1 ⇒ B = 0.5
3. 全解形式：

   y_zs(t) = A1e^{-t} + A2e^{-2t} + 0.5, t ≥ 0+

4. 代入初始值：

   y_zs(0+) = A1 + A2 + 0.5 = 0 y_zs'(0+) = -A1 -2A2 = 2

   解得：

   A1 = -3, A2 = 0.5

   因此：

   y_zs(t) = -3e^{-t} + 0.5e^{-2t} + 0.5, t ≥ 0+

3.3 全响应合成

y(t) = y_zi(t) + y_zs(t) = e^{-t} - 2.5e^{-2t} + 0.5, t ≥ 0

响应成分分析：

4. 实战检查清单与常见错误

4.1 初始值问题解决五步法

1. 判冲激：检查方程右端是否含有δ(t)或其导数
2. 定跃变：根据冲激最高阶数确定哪些导数会跃变
3. 算跳量：通过积分法计算各阶导数的跃变量
4. 分响应：区分零输入和零状态响应的初始条件
5. 验结果：验证各部分初始值之和是否符合全响应

4.2 典型错误案例

错误1：忽视δ(t)导致跃变判断失误

方程：y'' + 2y' = δ(t) 错误做法：直接认为y'(0+) = y'(0-) 正确解法：y'应有跃变，y连续

错误2：零状态响应初始值处理不当

错误：y_zs(0-) = 0 ⇒ 直接设y_zs(0+) = 0 正确：需考虑可能的跃变，如本例中y_zs'(0+)=2

错误3：冲激平衡分析不完整

方程：y'' + y = δ'(t) 错误：仅考虑y''含δ'(t)，忽略y'需含δ(t) 正确：y''含δ'(t) ⇒ y'含δ(t)并有跃变 ⇒ y连续

4.3 MATLAB验证代码

% 系统定义 num = [1 2 0]; % δ'(t)+2δ(t)对应的系数 den = [1 3 2]; % 微分方程左端系数 sys = tf(num, den); % 计算阶跃响应（对应ε(t)部分） t = 0:0.01:10; yzs_step = step(sys, t); % 计算冲激响应（对应δ(t)部分） yzs_impulse = impulse(sys, t); % 手工计算的理论曲线 y_theory = exp(-t) - 2.5*exp(-2*t) + 0.5; % 绘图比较 plot(t, yzs_step + yzs_impulse, 'b', t, y_theory, 'r--') legend('MATLAB仿真', '理论计算') title('全响应验证') xlabel('时间t') ylabel('响应y(t)')

运行此代码将显示仿真结果与理论计算曲线完美重合，验证了我们推导的正确性。