ARTICLE DETAIL

资讯详情

深耕网站视觉设计与运营推广的一线实战洞察。

连续系统时域分析:从微分方程到冲激响应与卷积积分

连续系统时域分析:从微分方程到冲激响应与卷积积分 1. 这一章为什么是分水岭连续系统时域分析的整体设计思路连续系统的时域分析说白了就是一句话不借助任何变换工具直接盯着时间轴看一个连续系统对输入到底做了什么。信号与系统的整门课本质上就两条腿一条是时域一条是变换域。第1章把信号本身拆开揉碎了讲讲的是原料到了第2章主角换成了系统讲的是原料怎么被加工。很多人学到这里会突然卡住不是因为公式难而是因为前面的思路还没切换过来——第一章你在分析一个函数第二章你要分析一个映射。这一章要建立的核心能力有三块第一能把一个实际的电路、机械结构或者任何物理过程翻译成一个带初始条件的线性常系数微分方程第二能从这个方程里把系统本身的固有属性和外部激励带来的响应分开算清楚第三能用冲激响应 h(t)和卷积积分这套组合拳直接算出任意激励下系统的零状态响应。这三块能力是后面拉普拉斯变换、傅里叶分析的底座底座没打牢后面的部分基本就是背公式。这篇内容适合三类人正在跟着教材啃第2章的学生、准备考研需要把这一章彻底吃透的人、还有工作里偶尔要处理动态系统建模的工程人员。我不会按教材目录照抄一遍而是按我自己踩过坑之后的思路顺序来讲——先讲清楚为什么这么设计再讲具体怎么算最后把常见错法和考场节奏摊开来说。1.1 时域分析到底在解决什么问题先给一个最朴素的问题给你一个系统再给你一个输入你怎么知道输出长什么样时域分析的答案非常直接——把系统的微分方程写出来解它。听起来很笨但正因为笨它把所有细节都摊在明面上每一个跳变、每一个冲激、每一段积分限你都看得见这也是为什么考研大题特别偏爱这一章它没法像变换域那样靠背表蒙过去。时域分析的另一个价值在于它给出了系统的身份证这个概念。一个线性时不变系统它的全部外部特性都可以被一个函数 h(t) 完整描述。你只要拿到了 h(t)任意输入下的零状态响应都是 h(t) 和输入的卷积。这个结论是整个信号与系统课程里最漂亮的结果之一也是后面系统函数 H(s) 与 h(t) 是一对拉氏变换的伏笔。所以这一章不是让你学会解微分方程——解微分方程是大一高数的事。它真正想让你理解的是系统的固有响应长什么样、外部激励怎么改造它、这两者怎么叠加成最终的输出。想清楚这个层次很多题目的解法就自然出来了。1.2 三条解题路线的取舍经典法、零输入零状态法、卷积法同一个微分方程教材上一共给了三条路很多人学完就乱了不知道该用哪条。我按自己的使用习惯说一下取舍逻辑。第一条是经典法也就是齐次解 特解 待定系数。它的优点是通用任何线性常系数方程都能用缺点是待定系数必须用 0 时刻的初始条件而题目通常给的是 0- 时刻的初始状态中间差一个跳变这一步最容易翻车。第二条是零输入响应 零状态响应分解法。这条路的好处是把问题劈成两半一半只看初始储能用齐次解形式加上 0- 条件就能定另一半只看输入用卷积算。两条路各自内部不打架不会出现初始条件到底用 0- 还是 0的纠结。第三条是卷积法专门劈给零状态响应。它的物理意义最清楚也最适合处理分段函数、方波、指数衰减这类输入。我的建议是用第二条做主框架第一条只用来处理强迫响应这一小部分第三条专门负责零状态响应。三条路线不是并列关系是分工关系。真到了考场上一道大题往往需要三条路配合着用而不是死抱着一条。1.3 学习顺序建议与刷题节奏顺序上我强烈建议倒过来学先学冲激响应和卷积再回头学经典法。原因很简单冲激响应 h(t) 的形式本身就是齐次解的形状你先搞懂 h(t) 怎么来的再回头看齐次解、特征根、待定系数会有一种原来你早就在这儿等着我的感觉。反过来先啃经典法的待定系数很容易陷入机械计算而不能理解意义。刷题节奏上这一章我建议按三遍法来第一遍只做方程到 h(t) 的转换不做卷积第二遍专攻卷积积分限的确定拿方波、三角波、指数函数互相卷第三遍再上综合大题把零输入、零状态、全响应、自由响应、强迫响应五个量一次性算全。第三遍做完这一章基本就稳了。大理大学这类学校的考研真题里这一章通常是一道 15 到 20 分的大题分值占比不低值得花这个时间。2. 第一步把物理系统翻译成微分方程时域分析的所有计算都建立在先把方程写对这个前提上。方程写错了后面解得再漂亮也是零分。我见过太多人卷积算得很溜结果 KVL 列错一个符号整道题白做。所以这一节单独拎出来讲建模不是为了凑篇幅是因为它真的值得。2.1 电路系统建模的固定套路电路建模的套路非常固定只有三步第一步对每个独立节点列 KCL对每个独立回路列 KVL第二步把电容、电感的伏安关系代进去——电容是 i C·du/dt电感是 u L·di/dt第三步消去中间变量只留下你关心的那个响应量和一个激励量整理成标准形式。标准形式长这样a_n · d^n y/dt^n ... a_1 · dy/dt a_0 · y b_m · d^m f/dt^m ... b_1 · df/dt b_0 · f左边是响应及其各阶导数右边是激励及其各阶导数系数都是常数对线性时不变系统而言。这里有个经验整理完以后先把最高阶的系数归一化也就是两边同除以 a_n让最高阶导数前面的系数变成 1。这么做有两个好处一是后面求特征根时形式统一二是右边激励的系数会直接变成 h(t) 里的系数省得来回倒。举个典型的 RC 电路电压源 f(t) 通过电阻 R 给电容 C 充电取电容电压 u_C 为响应。由 KVL 有 R·i u_C f再由 i C·du_C/dt 代入得到 RC·du_C/dt u_C f。归一化后是 du_C/dt (1/RC)·u_C (1/RC)·f时间常数 τ RC。这个例子非常小但它的结果会反复出现在后面所有的一阶分析里。2.2 机械系统与其他领域的建模与类比机械系统的建模逻辑和电路一模一样只是换了一套物理量。弹簧的胡克定律 F k·x阻尼器的粘滞摩擦 F c·dx/dt质量的牛顿第二定律 F m·d²x/dt²。把它们按受力分析拼起来得到的就是一个二阶微分方程形式和 RLC 电路完全同构。这种同构不是巧合而是线性系统理论的普适性体现。下面这张表是我自己整理时常用的对应关系理解了它你就明白为什么教材讲电路就够了——本质上是同一套数学电路量机械平移量数学角色电压 u力 F势变量跨变量电流 i速度 v流变量通变量电阻 R阻尼 c耗能元件电容 C弹簧柔度 1/k储能元件势电感 L质量 m储能元件流并联并联同势节点分析串联串联同流回路分析这张表的价值在于你只要会一种建模剩下的靠类比就能平移过去。实际工作中遇到热传导、流体、甚至某些经济系统的一阶惯性过程绝大多数都能塞进这个框架。2.3 建模阶段的三个高频陷阱第一个陷阱是响应量选错。电路题有时候问的是电流而不是电压这时候如果你一开始就按电压列方程最后还得再求一次导多一步就多一个出错点。我的习惯是题目问什么我就让什么当 y中间变量全部用 y 和激励表示最后一起消。第二个陷阱是符号方向。KVL 里到底是 u_C 还是 -u_C取决于你规定的参考方向。很多人列方程时符号乱本质是参考方向没先标好。我的做法是动笔前先在草稿纸上把每个元件的电压电流参考方向画出来画完了再列方程这一步花三十秒能省十分钟的复查时间。第三个陷阱是把非线性元件当线性。二极管、饱和电感、变容二极管这些都不是线性的一旦出现微分方程就不再是常系数线性形式本章的所有方法全部失效。考试里如果真出现了要么是题目给了小信号线性化条件要么就是你理解错了题意。提示方程列完以后先别急着解。拿 t → ∞ 的稳态情况代进去验一遍看看量纲和数量级对不对。这一步能拦住大部分低级错误。3. 经典法求解微分方程齐次解、特解与待定系数经典法是这一章的基本功虽然实际解题中我更愿意用分解法但齐次解的形式、特解的猜法、待定系数的确定这三件事你必须清清楚楚否则冲激响应那一节你根本看不懂它是怎么来的。3.1 齐次解的形式由特征根决定把方程右边置零得到齐次方程。设解的形式为 e^{λt}代进去就能得到特征方程解出特征根 λ。齐次解的形式完全由特征根的类型决定一共三种情况我整理成一张表特征根类型齐次解对应项单实根 λC·e^{λt}r 重实根 λ(C_1 C_2·t ... C_r·t^{r-1})·e^{λt}共轭复根 α ± jβe^{αt}·(C_1·cos βt C_2·sin βt)这里要说清楚一件事齐次解里出现的这些系数还没定。它们要靠初始条件来确定而初始条件到底是 0- 还是 0是这一章最大的争议点我在 3.3 节专门处理。关于重根我必须强调一下很多人看到二重根就只写一项 C·e^{λt}忘了乘以 t。判断方法很简单——特征方程是 n 阶的齐次解里就必须有 n 个独立的待定系数少一个就说明你漏了重根的多项式因子。这个自查方法非常有效我每次做完齐次解都会数一遍系数个数。3.2 特解怎么猜激励与特解形式对照表特解的形式由方程右边的激励函数形状决定这就是所谓的待定系数法。规则的核心是特解的形式和激励的形状一致除非这个形状和齐次解撞车撞车了就乘以 t 升阶直到不撞为止。激励 f(t) 的形式特解 y_p(t) 的设法备注常数 AP常数0 不是特征根时可直接用多项式 t^mP_m t^m ... P_1 t P_00 是 k 重根时整体乘 t^ke^{at}P·e^{at}a 是 k 重特征根时乘 t^kcos ωt 或 sin ωtP_1·cos ωt P_2·sin ωt±jω 是特征根时整体乘 te^{at}·cos ωte^{at}(P_1 cos ωt P_2 sin ωt)同理撞车就乘 t以 f(t) e^{-t} 为例如果 -1 不是特征根特解直接设 P·e^{-t}如果 -1 是单重特征根就要设 P·t·e^{-t}如果是二重根设 P·t²·e^{-t}。这个乘以 t 的幂次的动作就是所谓的共振在微分方程里的体现和物理上讲的外力频率等于固有频率时振幅线性增长是同一件事。注意特解里的系数是在代入原方程之后比较同类项定出来的不能拿初始条件去定。初始条件只管齐次解那部分的系数两个系数的确定渠道完全不同混用必错。3.3 从 0- 到 0初始条件跳变的处理这是整章最容易失分的地方我单独用一小节讲透。0- 时刻的初始状态 y(0-)、y(0-) 是电路或系统在激励接入之前就已经存在的储能状态比如电容上已经有的电压、电感里已经有的电流。这些量是物理事实题目一般直接给。0 时刻的起始条件 y(0)、y(0) 才是经典法待定系数时真正要用的东西。问题在于这两个量有时候不相等——当激励在 t 0 时刻含有冲激或者冲激的导数时系统在 0 点会发生瞬间跳变0 的值会偏离 0- 的值。判断要不要跳变看方程右边激励的最高阶导数次数 m 和左边响应的阶数 n情况0 与 0- 的关系处理方式右侧不含 δ 及其导数m ≤ 0 且激励有界一般相等直接令 0 0-右侧含 δ(t)低阶导数可能跳变用匹配法求跳变量右侧含 δ(t) 或更高阶更多阶导数跳变用匹配法逐阶求右侧含 ε(t)阶跃一般不跳变阶跃本身有界积分后连续一个实用的判断口诀方程两边同时对 t 在 [0-, 0] 上做积分右边积分不出来的项就对应左边的跳变项。这个口诀在处理冲激响应时尤其好用。3.4 冲激函数匹配法实操匹配法的具体操作是假设响应及其各阶导数在 0 附近含有冲激项设出形式代入方程比较 δ 及其各阶导数的系数把跳变量解出来。我拿一道典型题走一遍。设系统的微分方程为y(t) 3y(t) 2y(t) f(t) 2f(t)求冲激响应 h(t)。第一步令 f(t) δ(t)则 h(t) 满足h(t) 3h(t) 2h(t) δ(t) 2δ(t)第二步判断 h(t) 里含不含 δ。因为右边最高阶导数是 1 阶左边阶数是 2 阶m n所以 h(t) 本身不含 δ 项但 h(t) 和 h(t) 在 0 点可能有跳变。第三步设 h(t) 中含 a·δ(t) b·δ(t) 两项。往回推h(t) 中含 a·δ(t) b·ε(t)h(t) 中含 a·ε(t)把这些代回方程左边a·δ(t) b·δ(t) 3a·δ(t) 3b·ε(t) 2a·ε(t)比较系数δ(t) 的系数给出 a 1δ(t) 的系数给出 b 3a 2代入 a 1 得 b -1。第四步把跳变量翻译成起始条件。a·ε(t) 意味着 h(t) 在 0 处跳变 a 1即 h(0) - h(0-) 1而 h(0-) 0所以 h(0) 1。同理 b·ε(t) 意味着 h(0) - h(0-) b -1所以 h(0) -1。第五步h(t) 在 t 0 时满足齐次方程特征根是 -1 和 -2所以h(t) (C_1·e^{-t} C_2·e^{-2t})·ε(t)代入起始条件C_1 C_2 1-C_1 - 2C_2 -1。解得 C_2 0C_1 1。所以h(t) e^{-t}·ε(t)第六步验算。把 h e^{-t}ε、h -e^{-t}ε δ、h e^{-t}ε - δ δ 代回左边得 δ 2δ与右边一致结果正确。我在第一次做这道题时把 h(0) 直接当成了 0算出来 h(t) 0当场蒙了。后来才明白右边那个 f(t) 就是罪魁祸首正是它逼着 h(0) 从 0 跳到了 -1。记住一个反直觉的结论冲激响应的起始条件往往不为零而零输入响应的起始条件就取 0- 的值两者不能混。顺便提一个验证技巧用算子法算一遍。把方程写成算子形式 H(p) (p 2)/(p² 3p 2)因式分解后是 (p 2)/[(p 1)(p 2)] 1/(p 1)对应 h(t) e^{-t}ε(t)。两种方法对上心里就有底了。4. 零输入响应与零状态响应把两块拼图分开算如果只让我教一种解法我会教这一种。零输入零状态分解法的最大好处是它把系统自身和外部激励这两个独立的因果源彻底分开了各自有各自的确定系数的渠道不会互相污染。4.1 零输入响应只跟初始储能有关零输入响应 y_zi(t) 的定义是激励为零仅由 0- 时刻的初始状态引起的响应。它的形式一定是齐次解的形式因为方程右边是零。计算步骤只有三步写出齐次解的通式代入 0- 时刻的初始状态解线性方程组定系数。注意这里用的是0- 而不是 0因为激励为零系统在 0 点不会发生跳变0 等于 0-。举一个例子。系统方程 y 3y 2y f(t)初始状态 y(0-) 1y(0-) 2。特征根是 -1 和 -2零输入响应形式为y_zi(t) (C_1·e^{-t} C_2·e^{-2t})·ε(t)代入初始状态C_1 C_2 1-C_1 - 2C_2 2。解得 C_1 4C_2 -3所以y_zi(t) (4e^{-t} - 3e^{-2t})·ε(t)这里补一个物理直觉零输入响应的每一条指数项都是衰减的对稳定系统而言衰减快慢由特征根的实部决定振荡与否由虚部决定。这跟系统的固有频率、时间常数、阻尼比是同一回事。4.2 零状态响应只跟激励有关零状态响应 y_zs(t) 的定义是初始状态为零仅由激励引起的响应。它的标准算法有两种一种是先用经典法求出全响应再减去零输入响应另一种是直接做卷积 y_zs f(t) * h(t)。我强烈推荐第二种因为 h(t) 一旦求出来后面所有输入都是机械操作。零状态响应里有一个坑值得提虽然初始状态是零但零状态响应在 0 时刻的起始条件一般不为零原因和 3.4 节讲的一样——激励在 0 点的冲激成分会强制产生跳变。所以算零状态响应时千万不要拿0- 为零去定系数那是零输入响应的定法。4.3 两种分解方式的关系教材上还给了一种分解全响应 自由响应 强迫响应。这两种分解容易混淆我用一张表把它们的关系说清楚分解方式分法依据分出来的两块系数由谁定零输入 零状态激励是否为零y_zi、y_zs分别由 0- 状态、激励定自由 强迫是否为齐次解形式自由响应、强迫响应特解自由响应的系数由整体初始条件定关键结论是自由响应 零输入响应 零状态响应中的齐次解部分也就是说零状态响应里也含有齐次解形式的成分它和零输入响应形状一样只是系数不同。很多题目要求分别指出自由响应和强迫响应这时候你先算全响应再把非齐次形式的那部分挑出来当强迫响应剩下的就是自由响应。别硬拆硬拆很容易漏项。提示题目如果只要求求全响应用零输入加零状态最快如果明确问自由响应和强迫响应就先算出全响应再按形式归类。看清题目是哪一种能省不少无谓的计算。5. 冲激响应与阶跃响应拿到系统的身份证到了这一节这一章的重心开始从解方程转向认系统。冲激响应 h(t) 是线性时不变系统最核心的概念没有之一。5.1 冲激响应的定义与求法定义很直接系统初始状态为零时单位冲激函数 δ(t) 作为输入所产生的响应记作 h(t)。它的物理意义是系统对瞬间猛烈一击的自然反应——这一击之后输入就消失了所以剩下的完全是系统自己的固有振动。求法上我总结成两条路径路径一是解方程法。令方程右边只剩 δ(t)如果原方程右边还有 f 的各阶导数那就换成 δ 的各阶导数用 3.4 节的匹配法求出 h(0)、h(0)、h(0) 等起始条件然后套齐次解形式定系数。这条路通用缺点是匹配那一步容易出错。路径二是算子法部分分式展开。把微分方程写成算子形式得到传输算子 H(p) N(p)/D(p)对它做部分分式展开每一项查表反变换回时域。以下是常用对照1/(p - λ) → e^{λt}·ε(t) 1/(p - λ)^k → t^{k-1}/(k-1)! · e^{λt}·ε(t) 1/(p α)² β² → (1/β)·e^{-αt}·sin(βt)·ε(t) (p α)/[(pα)²β²] → e^{-αt}·cos(βt)·ε(t)用 3.4 节那道题验证一下H(p) (p 2)/(p² 3p 2) 1/(p 1)查表得 h(t) e^{-t}ε(t)和解方程法的结果完全一致。实际做题我一般先用算子法快速出结果再用解方程法验算一遍两道防线基本不会错。还有一个结论要记牢h(t) 里一定带 ε(t)对因果系统而言因为 t 0 时系统还没被激发响应必须为零。如果你算出来的 h(t) 里没带 ε(t)八成是漏了。5.2 阶跃响应与冲激响应的互推关系阶跃响应 g(t) 的定义是初始状态为零时单位阶跃 ε(t) 作为输入产生的响应。它和 h(t) 之间有两条互推关系非常好用g(t) ∫_{-∞}^{t} h(τ) dτ h(t) dg(t)/dt这两条关系的直观解释是δ(t) 是 ε(t) 的导数而线性时不变系统对输入求导输出也求导对输入积分输出也积分。据此还能推出一个非常有用的技巧如果求 h(t) 的匹配法太麻烦可以先求 g(t)对阶跃输入的响应通常更容易处理再求导得到 h(t)。阶跃响应本身也有物理意义它描述的是系统对突然接到一个恒定输入的过渡过程工程上讲的上升时间、超调量、调节时间全都是从 g(t) 的曲线上读出来的。这也是为什么控制系统里更喜欢用阶跃响应来做性能指标。5.3 常见系统的 h(t) 清单下面这张表是我整理的常考结果建议直接背下来考场上能省一半时间系统方程特征根h(t)y a·y f-ae^{-at}·ε(t)y a·y f-aδ(t) - a·e^{-at}·ε(t)y 2y y f-1二重t·e^{-t}·ε(t)y 2y 5y f-1 ± 2j(1/2)·e^{-t}·sin(2t)·ε(t)y 3y 2y f-1, -2(e^{-t} - e^{-2t})·ε(t)y 2y y f-1二重(1 - t)·e^{-t}·ε(t)最后一行值得单独说右边出现 f 时h(t) 的结构会发生明显变化可能出现 δ(t) 项。一般规律是如果方程右边的最高阶导数阶数 m 大于等于左边阶数 nh(t) 里就会出现 δ(t) 及其各阶导数。这条规律是判断 h(t) 结构是否合理的快速方法。6. 卷积积分时域方法的心脏如果这一章只能留一个知识点我会留卷积。因为它是把 h(t) 从一个函数变成一把万能钥匙的那一步。6.1 卷积的物理意义与定义卷积的定义式是y(t) f(t) * h(t) ∫_{-∞}^{∞} f(τ)·h(t - τ) dτ它的物理意义可以用一句话概括把输入信号看成无数个在不同时刻出现的、强度不同的冲激的叠加每个冲激在系统里激起一个按 h(t) 形状演化的响应把所有时刻的响应在观察时刻 t 上叠加起来就是输出 y(t)。这个思路就是分解—叠加的时域版本。理解了这个意义很多性质就变成显而易见的了。比如交换律 f * h h * f意味着把输入拆成冲激和把系统拆成冲激两条路等价比如 f(t) * δ(t) f(t)意味着单位冲激是卷积的单位元比如 f(t) * δ(t - t_0) f(t - t_0)意味着延迟的冲激只是把响应整体平移。6.2 图解法的标准四步与积分限确定图解法是考试里最常用的手算方式标准流程四步换元把 f 和 h 的自变量都换成 τ画出两者的波形。翻转把 h(τ) 关于纵轴翻成 h(-τ)。平移把 h(-τ) 沿 τ 轴平移 t得到 h(t - τ)。相乘积分计算两个波形重叠区域内 f(τ)·h(t - τ) 的积分作为 t 时刻的输出。这里面唯一有技术含量的就是第四步的重叠区间判断。我拿一道典型题完整走一遍。已知 f(t) ε(t) - ε(t - 2)宽度为 2 的矩形脉冲h(t) e^{-t}·ε(t)求 y(t) f(t) * h(t)。直接按定义写y(t) ∫_{-∞}^{∞} f(τ)·h(t-τ) dτ ∫_{-∞}^{∞} f(τ)·e^{-(t-τ)}·ε(t-τ) dτ两个约束条件f(τ) 非零要求 0 ≤ τ ≤ 2ε(t - τ) 非零要求 τ ≤ t。两个区间的交集决定了积分上下限按 t 的位置分三段讨论。第一段t 0。此时 τ ≤ t 与 0 ≤ τ ≤ 2 没有交集积分为零所以 y(t) 0。第二段0 ≤ t 2。交集是 0 ≤ τ ≤ t于是y(t) ∫_0^t e^{-(t-τ)} dτ e^{-t}·(e^t - 1) 1 - e^{-t}第三段t ≥ 2。此时 τ ≤ t 的约束在 [0, 2] 内自动满足交集就是整个 0 ≤ τ ≤ 2于是y(t) ∫_0^2 e^{-(t-τ)} dτ e^{-t}·(e^2 - 1) e^{-(t-2)} - e^{-t}在 t 2 处验算一下第二段给出 1 - e^{-2}第三段给出 e^{0} - e^{-2} 1 - e^{-2}两边一致说明分段点接得上结果可信。这道题里我踩过一次坑把第三段的积分上限定成了 t忘了 f(τ) 在 τ 2 之后是零。积分限的确定永远是两个非零区间的交集谁窄听谁的。如果你只盯着 h(t - τ) 的约束忽略了 f 本身的有限支撑结果就会在 t 2 之后继续增长明显与物理直觉不符——输入早就结束了输出不可能一直涨。再补一个更短的例子练手感f(t) e^{-t}ε(t)h(t) e^{-t}ε(t)。两个函数的非零区间在 τ ≥ 0 和 τ ≤ t 上取交集t 0 时为零t 0 时为y(t) ∫_0^t e^{-τ}·e^{-(t-τ)} dτ e^{-t}·∫_0^t dτ t·e^{-t}所以 y(t) t·e^{-t}·ε(t)。这个结果和拉氏变换里 1/(s1)² 反变换的结果一致可以互相印证。6.3 卷积的运算性质与常用卷积表卷积的性质很多但真正高频使用的就那么几条我挑出来列成表性质表达式使用场景交换律f * h h * f哪个好翻转就翻哪个分配律f * (h_1 h_2) fh_1 fh_2输入可分解时拆开算结合律(f * h_1) * h_2 f * (h_1 * h_2)多级系统级联与冲激f(t) * δ(t) f(t)单位元时移f(t - t_1) * h(t - t_2) y(t - t_1 - t_2)延迟输入直接平移结果微分f * h f * h y(t)避开发散积分积分f * ε(t) ∫_{-∞}^{t} f(τ) dτ阶跃卷积就是前置积分尺度变换f(at) * h(at) (1/|a|)·y(at)波形压缩后结果要乘系数尺度变换这条最容易被忽略。很多人卷完发现结果和答案差一个系数十有八九是忘了这个 1/|a|。它背后的道理是自变量压缩会把积分区间也压缩区间长度变成原来的 1/|a|所以积分值要相应补偿。常用的一小张卷积表考前扫一眼就够ε(t) * ε(t) t·ε(t) t·ε(t) * ε(t) (t²/2)·ε(t) e^{-at}ε(t) * ε(t) (1/a)(1 - e^{-at})·ε(t) e^{-at}ε(t) * e^{-bt}ε(t) [1/(b-a)]·(e^{-at} - e^{-bt})·ε(t)a ≠ b e^{-at}ε(t) * e^{-at}ε(t) t·e^{-at}·ε(t)使用时注意一点前四条都是针对因果信号带 ε(t)的如果信号是非因果的比如 e^{-at}ε(-t)必须老老实实从定义式出发算不能套这些表。6.4 零状态响应 激励卷积冲激响应把前面几节串起来这一章的完整闭环就出来了第一步由微分方程求冲激响应 h(t) 第二步零输入响应 y_zi 由 0- 初始状态定 第三步零状态响应 y_zs f(t) * h(t) 第四步全响应 y(t) y_zi(t) y_zs(t)这条链路是这一章的骨架任何一道大题都跑不出这四步。我在做综合题时习惯先把这条链路写在草稿纸最上面每算完一步打个勾最后一步之前严格复查 h(t) 和卷积限基本就不会丢分。顺便说一句h(t) 的因果性和稳定性也和这条链路强相关h(t) 0 对 t 0 成立系统是因果的∫|h(t)|dt 收敛系统是稳定有界输入有界输出的。这两条判断在选择题里出现的频率极高而且是纯送分题前提是你 h(t) 求对了。7. 常见问题排查与踩坑实录前面六节讲的是应该怎么做这一节讲的是实际做的时候会在哪儿翻车。这一节的内容大部分不是教材上写的是我做题和带人做题时反复看到的实录。7.1 高频错误速查表现象常见原因正确做法零输入响应系数算出来是复数特征根共轭硬用实系数方程解直接设 e^{αt}(C_1 cos βt C_2 sin βt)保持实数形式全响应比答案多一项零状态响应里忘了带 ε(t)检查所有指数项后面是否有 ε(t)卷积在 t 大于某个值后还在涨积分上限用了 t忽略了 f 的有限支撑取两个非零区间的交集谁窄用谁h(t) 里出现 δ(t) 但答案没有没判断 m 与 n 的大小关系m n 时 h(t) 不含 δm ≥ n 时含特解代入后同类项消不掉特解形式和激励撞车没乘 t检查激励指数是否为特征根是则升阶结果和拉氏变换法对不上卷积尺度变换漏了 1/|a|做尺度变换时主动补系数起始条件用了 0- 导致跳变丢失混淆了 0- 状态与 0 条件经典法待定系数一律用 0最后一行我单独强调0- 是物理事实0 是数学入口。题目给的永远是 0-你解的永远是 0中间那个跳变量需要你自己通过匹配法算出来。这个区分看着啰嗦但它是一大半错误的根源。7.2 计算细节上的几个死磕点第一个死磕点是积分限的分段点。分段点的位置来自两个约束的临界值比如 t 0 和 t 2这两个点必须单独验算确认左右两段的结果接得上。接不上就是分段错了别硬着头皮往下写。第二个死磕点是冲激函数的运算。三条常用规则抽样性质 ∫f(t)δ(t - t_0)dt f(t_0)乘积性质 f(t)δ(t) f(0)δ(t)尺度性质 δ(at) (1/|a|)δ(t)。这三条在简化卷积时经常用到尤其是当输入里含 δ 时直接用抽样性质比老老实实积分快得多。第三个死磕点是 ε(t) 与常数的关系。有界信号乘不乘 ε(t)在 t 0 时看不出差别但卷积结果里差一个整体平移考试会扣分。我的检查习惯是拿到任何时域结果先看 t 0 的部分是不是零因果系统不是零就说明漏了 ε(t)。第四个死磕点是初值定理式的自检。算完 h(t) 以后可以用 h(0) 反推右端最高阶项的系数两者应该对得上。这一步花二十秒能拦住大部分匹配法算错的情况。7.3 大题模板与时间分配最后说一下考场节奏。这一章的大题通常 15 到 20 分我建议按下面的模板来写抄下微分方程标注阶数 n 和右端最高阶 m30 秒写特征方程求特征根标注类型1 分钟求 h(t)算子法或匹配法写出并验算4 分钟求零输入响应套齐次解形式用 0- 条件定系数3 分钟求零状态响应写卷积式分段讨论积分限5 到 6 分钟相加得全响应必要时拆自由与强迫分量2 分钟检查ε(t) 是否齐全、t 0 是否为零、分段点是否连续1 分钟整个过程 15 分钟上下留 5 分钟复查。如果某一步卡住了先把后面能写的写上分步给分的情况下把 h(t) 求对往往就能拿到一半的分。还有一个小经验草稿纸上一定要画波形图。哪怕题目不要求作图把 f(τ) 和 h(t - τ) 的示意图画出来积分限一眼就能看出来比纯靠代数推快得多也不容易错。我见过太多人嫌画图麻烦结果在积分限上纠结了十分钟得不偿失。8. 一些不太写在教材上的体会零输入响应和零状态响应的分解本质上是在回答锅里的味道到底来自食材还是来自火候。这个类比可能不太严谨但它帮我记住了一件事一个系统的输出从来不只由输入决定初始储能也是一份实实在在的贡献只是它不声不响地待在那儿。我在带学生的过程中发现这一章真正的门槛不在计算量而在肯不肯把每一个符号的含义想清楚。0- 和 0 差的那个跳变、卷积里翻转平移的物理意义、h(t) 为什么是系统的身份证这些东西想通了题目只是它的练习题想不通做一百道也只是在练手速。所以如果你现在正卡在这一章别急着刷题先把 3.4 节的匹配法和 6.2 节的积分限确定各自独立地推三遍推到闭着眼睛都不会错为止。最后一个可扩展的方向这一章所有解方程的活儿到第 4 章之后基本都可以被拉普拉斯变换替代。零输入响应变成初始条件对应的自由项零状态响应变成 H(s)·F(s) 的反变换卷积变成乘法。那时候你回头看这一章会发现它不是被淘汰了而是被抽象了——卷积的那套分解叠加思想才是整门课真正的内核。所以这一章多花点时间是值得的后面几章会反复收它的利息。
返回列表