【信号与系统】3.1 从电路到方程：LTI系统微分方程的建模与经典解法

1. 从电路到微分方程：LTI系统的建模实战

当你第一次看到RLC电路时，那些弯弯曲曲的符号和连线可能让人望而生畏。但有趣的是，这些看似复杂的物理系统，都可以用数学语言精确描述。我刚开始学信号与系统时，最震撼的时刻就是发现电路图居然能变成微分方程。

以最常见的二阶RLC串联电路为例。假设我们有一个电压源vi(t)连接着电阻R、电感L和电容C。选择电容两端的电压vc(t)作为系统响应，现在我们要建立这个系统的微分方程模型。

关键步骤其实很直观：

1. 根据基尔霍夫电压定律(KVL)，总电压等于各元件电压之和：vi(t) = vR(t) + vL(t) + vC(t)
2. 用元件特性方程替换各项：
   • 电阻电压vR(t) = R·i(t)
   • 电感电压vL(t) = L·di(t)/dt
   • 电容电流i(t) = C·dvc(t)/dt
3. 把这些关系代入KVL方程，经过整理就得到：
   LC·d²vc/dt² + RC·dvc/dt + vc = vi(t)

这个推导过程中最需要注意的就是变量统一。我刚开始经常犯的错误是保留多个变量（比如既有i(t)又有vc(t)），导致方程无法闭合。后来发现，选定一个主变量（这里是vc(t)）后，其他量都要用它的导数或积分表示。

2. 微分方程的标准形式与物理意义

观察我们得到的RLC电路方程，可以写成更通用的形式：
y''(t) + a₁y'(t) + a₀y(t) = b₀x(t)

这个标准形式揭示了LTI系统的核心特征：

• 阶数：最高导数项决定系统阶数（这里是二阶）
• 系数：a₁=1/RC，a₀=1/LC，这些系数直接反映电路参数
• 动态特性：二阶导数项对应惯性（电感），一阶导数项对应阻尼（电阻），零阶项对应弹性（电容）

在实际工程中，我经常用这个标准形式快速判断系统特性。比如：

• 当R很大时，a₁项主导，系统呈现强阻尼
• 当L很小时，a₀项很大，系统响应会更快
• 当C很大时，a₀很小，系统会有更长的振荡周期

3. 经典解法：拆解齐次解与特解

解这类微分方程，工程师们最常用的就是经典法——将解拆分为齐次解和特解。这个方法就像解数学题时的"分而治之"策略。

3.1 齐次解：系统的自由响应

齐次方程（右边=0）的解反映了系统的固有特性。以RLC电路为例：

1. 写出特征方程：s² + (R/L)s + 1/LC = 0
2. 求特征根：s = [-R/L ± √((R/L)² - 4/LC)] / 2

这里会出现三种情况，对应不同的物理现象：

• 过阻尼（实根）：(R/L)² > 4/LC，响应缓慢衰减
• 临界阻尼（重根）：(R/L)² = 4/LC，最快回到稳态
• 欠阻尼（共轭复根）：(R/L)² < 4/LC，振荡衰减

我在实验室调试电路时，经常通过调整电阻值来观察这三种状态。特别是欠阻尼情况下，能看到优美的衰减振荡波形。

3.2 特解：系统的强迫响应

特解对应系统在外部激励下的响应。选择特解形式时有个实用技巧：跟随激励的脚步。常见配对：

• 常数激励 → 常数特解
• e^at → Ae^at（a不是特征根时）
• sin(ωt) → Asin(ωt) + Bcos(ωt)

有个容易踩的坑：当激励形式与齐次解某部分重合时，需要在特解前乘以t。比如激励是e^t，而e^t已经是齐次解，那么特解应该用Ate^t。

4. 完整求解案例：RLC电路响应分析

让我们通过一个具体数值例子，体验完整的建模求解过程。假设：

• R=6Ω，L=1H，C=1/8F
• 输入电压vi(t)=10V（阶跃输入）
• 初始条件vc(0)=1V，vc'(0)=2V/s

步骤1：建立微分方程
代入参数得到：v''c(t) + 6v'c(t) + 8vc(t) = 80

步骤2：求齐次解
特征方程s² + 6s + 8 = 0 → s₁=-2, s₂=-4
齐次解：vc_h(t) = C₁e^(-2t) + C₂e^(-4t)

步骤3：求特解
假设常数特解vc_p(t)=A
代入得8A=80 → A=10

步骤4：确定常数
全解vc(t)=10 + C₁e^(-2t) + C₂e^(-4t)
利用初始条件：
vc(0)=10+C₁+C₂=1
v'c(0)=-2C₁-4C₂=2
解得C₁=-13, C₂=4

最终解：
vc(t) = 10 - 13e^(-2t) + 4e^(-4t) (t≥0)

这个结果清晰地展示了电路响应如何从初始状态过渡到稳态值10V。在实际示波器上观察这个过程，你会看到电压快速上升后略微过冲，然后平稳收敛到10V，这正是-2和-4两个指数项共同作用的结果。

5. 工程应用中的实用技巧

经过多个项目的实践，我总结了一些微分方程建模与求解的实用经验：

建模阶段：

• 选择合适的状态变量（通常选电容电压和电感电流）
• 检查方程量纲是否一致（防止推导错误）
• 无量纲化可以简化计算（特别是时间常数归一化）

求解阶段：

• 特征根实部决定稳定性（全部负实部才稳定）
• 复根的虚部对应振荡频率
• MATLAB的dsolve函数可以验证手算结果

参数设计：

• 临界阻尼电阻R=2√(L/C)常用于需要快速响应又不要振荡的系统
• 品质因数Q=√(L/C)/R反映系统选择性
• 带宽与特征根位置直接相关

在最近的一个电源滤波电路设计中，我正是通过这些方法，仅用两次迭代就确定了最优的LC参数组合，避免了反复试错的时间消耗。