capacitor 的式子里一会儿求导、一会儿从 −∞ 积分,写起来很长。Davis 用 Heaviside 的记号 p 和 1/p 把这两个运算写得像"乘一个因子",后面所有储能电路都用这套记号。这一节讲它们的定义,怎样作用在 f(t)u(t) 这种单边波形上,以及两个最容易丢东西的地方:积分下限带出的常数,和求导时冒出来的 impulse。
① 定义
differential operators / Heaviside operators
pf(t)=dtdf,p1f(t)=∫−∞tf(α)dα
p 表示"对时间求导",1/p 表示"从 −∞ 积到 t"。它们单独写没有意义,必须从左边作用在一个波形上;方括号 p[ ⋅ ] 表示作用在括号里的整个波形上,不是乘法。1/p 是定积分,所以结果里没有 +C。
capacitor 的两个式子于是变成 i=Cpv 和 v=Cp1i。Cp1 读作"先用 1/p 积分,再乘 1/C",也就是 C⋅p1,不是 C1⋅p。
它们是 linear operators(线性算子):p[af+bg]=apf+bpg,p1[af+bg]=ap1f+bp1g,a、b 是常数。所以常数系数可以先提出来,和的每一项可以分开算。
② 1/p 作用在 f(t)u(t) 上
单边波形 f(t)u(t) 在 t≤0 时是 0。从 −∞ 积到 t:若 t≤0,积出来是 0;若 t>0,从 −∞ 到 0 那一段贡献 0,只剩从 0 到 t 那一段。两种情况合起来,就是"从 0 积到 t,再乘 u(t) 把负时间关掉":
技巧 B
p1[f(t)u(t)]=[∫0tf(α)dα]u(t)关键是下限是 0。代入下限时原函数在 0 处的值一般不是 0,会留下一个常数。这个常数保证结果在 t=0 从 0 开始(面积从 0 开始累积)。
例题例 1 常数
求 p1[6u(t)]。
f(t)=6,∫0t6dα=6t,所以 p1[6u(t)]=6tu(t):阶跃积分一次是斜坡 (ramp)。从图上看,宽 t、高 6 的矩形面积就是 6t。
例题例 2 正弦(课本例 6.6 是同一类题)
求 p1[sin(3t)u(t)]。
第 1 步,写成从 0 开始的定积分:p1[sin(3t)u(t)]=[∫0tsin3αdα]u(t)。
第 2 步,找原函数并代入上下限。sin3α 的原函数是 −31cos3α:
∫0tsin3αdα=[−31cos3α]0t=−31cos3t+31第二项 +31 来自下限:−(−31cos0)=31。
第 3 步,乘 u(t):p1[sin(3t)u(t)]=31(1−cos3t)u(t)。
核对:t=0 时结果是 31(1−1)=0 ✓,面积从 0 开始。如果只写不定积分 −31cos3t,t=0 时是 −31:面积还没开始累积就是负数,不可能。结果在 0 和 32 之间摆动、从不为负,因为 sin3t 每个周期都是正的半周先来。
例题例 3 指数
求 p1[10e−5tu(t)]。
∫0t10e−5αdα=[−2e−5α]0t=−2e−5t+2,乘 u(t):2(1−e−5t)u(t)。从 0 开始,慢慢涨到 2。2 正好是整条指数曲线下的总面积 ∫0∞10e−5αdα=10/5。
③ p 作用在 f(t)u(t) 上:乘积法则多出一个 impulse
把 f(t)u(t) 当成两个函数的乘积求导:p[fu]=(pf)u+f(pu)。第一项是平常的导数。第二项里 pu 是阶跃的导数:u 在 t=0 跳了 1,所以它的导数是一根面积为 1 的 impulse,pu=δ(t);再由 sampling property f(t)δ(t)=f(0)δ(t)(这两条的来历见 6·01)。合起来:
技巧 C
p[f(t)u(t)]=f′(t)u(t)+f(0)δ(t)先看 f(0):f(0)=0 就没有 impulse;f(0)=0 就多一根面积(strength)为 f(0) 的 impulse。
从图上看:f(0)=0 时,波形 f(t)u(t) 在 t=0 从 0 一下跳到 f(0),跳了多少,导数里那根 impulse 的面积就是多少。f(0)=0 时波形在 0 处连续(最多是个拐角),没有 impulse。这正是 6·01 里"每个跳变给一根 impulse,拐角不给"的特例。
例题例 4 f(0) = 0:没有 impulse
求 p[sin(3t)u(t)]。
f(t)=sin3t,f′(t)=3cos3t,f(0)=sin0=0。所以 p[sin(3t)u(t)]=3cos(3t)u(t),没有 impulse。波形 sin(3t)u(t) 在 t=0 只是从平的变成斜的(拐角),值不跳。
例题例 5 f(0) ≠ 0:有 impulse(课本例 6.17)
求 p[cos(t)u(t)]。
f(t)=cost,f′(t)=−sint,f(0)=cos0=1。所以 p[cos(t)u(t)]=−sin(t)u(t)+δ(t)。波形 cos(t)u(t) 在 t=0 从 0 跳到 1,导数里就有一根面积为 1 的 impulse。
例题例 6 为什么 impulse 不能丢:检查一个漏掉 δ 的答案
有人算 p[2e−3tu(t)],只写了 −6e−3tu(t)。用 1/p 把这个答案积回去,和 2e−3tu(t) 比一比,找出他漏掉的那一项。
积回去。技巧 B:∫0t−6e−3αdα=[2e−3α]0t=2e−3t−2,得 (2e−3t−2)u(t)。
和原波形比。2e−3tu(t)−(2e−3t−2)u(t)=2u(t):积回去少了一个 2u(t)。t 很大时他的结果趋于 −2,原波形趋于 0,明显不对。
找出漏掉的项。漏掉的项 m(t) 积分后必须正好补上 2u(t),所以 m=p[2u(t)]=2δ(t)。正确答案是 p[2e−3tu(t)]=−6e−3tu(t)+2δ(t),和技巧 C 一致:f(0)=2。
这根 impulse 的面积 2 正好是波形在 t=0 从 0 跳到 2 的那一跳。1/p 要把这一跳"积"回来,靠的就是它。
p 和 1/p 互为逆运算(Davis 定义 6.1、6.2)
对 one-sided(单边)波形,并且把跳变产生的 impulse 都留着,pp1f=p1pf=f。为什么要单边:p1pf=∫−∞tf′(α)dα=f(t)−f(−∞),只有 f(−∞)=0 才等于 f(t)。例如常数 12:p[12]=0,积回去还是 0,回不到 12;单边的 12u(t):p[12u]=12δ,积回去正好是 12u。Davis 还规定这些运算是 causal(因果)的:输入是单边的,输出也是单边的,没有刺激就没有反应。
④ tneat 的积分:表格法
∫αe−αdα 这类积分要用分部积分 ∫udv=uv−∫vdu。每用一次,多项式降一次幂,直到求导成 0。课本附录 C 的 tabular integration(表格法)把这个过程排成一张表,做几次都不容易乱:
- 左列写多项式部分,一直求导到 0:α→1→0。
- 右列写指数部分,一直积分(每次都不加常数):e−α→−e−α→e−α。
- 左列第 k 行乘右列第 k+1 行(斜着往右下乘),符号依次取 +,−,+,…,加起来就是原函数。
| 符号 | 求导列 | 积分列 |
| + | α | e−α |
| − | 1 | −e−α |
| 0 | e−α |
∫αe−αdα=(+)α⋅(−e−α)+(−)1⋅e−α=−αe−α−e−α
代入上下限 0 和 t:[−αe−α−e−α]0t=(−te−t−e−t)−(0−1)=1−e−t−te−t。下限又贡献了常数 1。求导核对:dtd(1−e−t−te−t)=e−t−e−t+te−t=te−t ✓。
例题例 7 表格法做三行:p1[t2e−tu(t)]
求 h(t)=p1[t2e−tu(t)],再用 p 核对。
第 1 步,技巧 B 写成从 0 开始的定积分:h(t)=[∫0tα2e−αdα]u(t)。
第 2 步,表格法。多项式 α2 要求导两次才到 0,所以比上面 αe−α 那张表多一行:
| 符号 | 求导列 | 积分列 |
| + | α2 | e−α |
| − | 2α | −e−α |
| + | 2 | e−α |
| 0 | −e−α |
∫α2e−αdα=α2(−e−α)−2αe−α+2(−e−α)=−(α2+2α+2)e−α第 3 步,代上下限。上限给 −(t2+2t+2)e−t。下限 α=0 处原函数等于 −2,减去它就是 +2:
h(t)=[2−(t2+2t+2)e−t]u(t)第 4 步,用 p 核对(技巧 C)。括号在 t=0 等于 2−2=0,没有 impulse。求导:−(2t+2)e−t+(t2+2t+2)e−t=t2e−t ✓,回到了被积函数。
看极限:t→∞ 时 h→2,正是整条曲线 t2e−t 下的总面积。下限带出的常数 2,就是最后攒下的总量;丢了它,h 会从 −2 开始、最后趋于 0,面积成了负的。
⑤ 1/p2:积两次
p21=p1⋅p1:先积一次,对结果再积一次。每一次都用技巧 B(从 0 积到 t,乘 u(t)),每一次的下限都可能带出常数,都不能丢。
例题例 8 1/p² 作用在指数上:两次都有下限常数(先自己做)
求 p21[4e−2tu(t)]。
看解答
第一次积分:
∫0t4e−2αdα=[−2e−2α]0t=−2e−2t+2得 2(1−e−2t)u(t)。下限贡献了 +2。
第二次积分(对第一次的结果再用一次技巧 B):
∫0t2(1−e−2α)dα=[2α+e−2α]0t=2t+e−2t−1下限这次贡献 −1(e0=1 被减掉)。结果是 (2t−1+e−2t)u(t)。
核对:t=0 时 0−1+1=0 ✓。求导一次:(2−2e−2t)u(括号在 0 处为 0,无 impulse);再求导:4e−2tu(2−2=0,无 impulse)✓。
形状:第一次积分趋于常数 2,所以第二次积分趋于斜率 2 的直线 2t−1。那个 −1 是因为第一次积分要花时间才涨到 2,比"一开始就是 2"少积了 1。如果漏掉第一次的 +2,第二次会得到 (e−2t−1)u(t):不涨反降,一看就错。
这一节最常见的三个错
(1) 积分只写原函数,忘了代下限 0 带出的常数(例 2、例 8)。
(2) 求导忘了看 f(0),丢掉 f(0)δ(t)(例 5、例 6)。
(3) 结果忘了乘 u(t):31(1−cos3t) 在 t<0 时并不是 0,乘上 u(t) 才是正确的单边波形。