6·05

6.2 operator p 和 1/p:把求导、积分写成"乘法"

About 16 minBook · Davis 6.2

capacitor 的式子里一会儿求导、一会儿从 −∞-\infty 积分,写起来很长。Davis 用 Heaviside 的记号 pp 和 1/p1/p 把这两个运算写得像"乘一个因子",后面所有储能电路都用这套记号。这一节讲它们的定义,怎样作用在 f(t)u(t)f(t)u(t) 这种单边波形上,以及两个最容易丢东西的地方:积分下限带出的常数,和求导时冒出来的 impulse。

① 定义

differential operators / Heaviside operators p f(t)=dfdt,1p f(t)=∫−∞tf(α) dα\begin{gathered} p\,f(t) = \frac{df}{dt}, \\ \frac{1}{p}\,f(t) = \int_{-\infty}^{t} f(\alpha)\,d\alpha \end{gathered}

pp 表示"对时间求导",1/p1/p 表示"从 −∞-\infty 积到 tt"。它们单独写没有意义,必须从左边作用在一个波形上;方括号 p[ ⋅ ]p[\ \cdot\ ] 表示作用在括号里的整个波形上,不是乘法。1/p1/p 是定积分,所以结果里没有 +C。

capacitor 的两个式子于是变成 i=Cp vi = Cp\,v 和 v=1Cp iv = \frac{1}{Cp}\,i。1Cp\frac{1}{Cp} 读作"先用 1/p1/p 积分,再乘 1/C1/C",也就是 1C⋅p\frac{1}{C\cdot p},不是 1C⋅p\frac1C \cdot p。

它们是 linear operators(线性算子):p[af+bg]=a pf+b pgp[af + bg] = a\,pf + b\,pg,1p[af+bg]=a 1pf+b 1pg\frac1p[af + bg] = a\,\frac1p f + b\,\frac1p g,aa、bb 是常数。所以常数系数可以先提出来,和的每一项可以分开算。

② 1/p1/p 作用在 f(t)u(t)f(t)u(t) 上

单边波形 f(t)u(t)f(t)u(t) 在 t≤0t \le 0 时是 0。从 −∞-\infty 积到 tt:若 t≤0t \le 0,积出来是 0;若 t>0t \gt 0,从 −∞-\infty 到 0 那一段贡献 0,只剩从 0 到 tt 那一段。两种情况合起来,就是"从 0 积到 tt,再乘 u(t)u(t) 把负时间关掉":

技巧 B 1p[f(t)u(t)]=[∫0tf(α) dα] u(t)\frac{1}{p}\big[f(t)u(t)\big] = \Big[\int_{0}^{t} f(\alpha)\,d\alpha\Big]\,u(t)

关键是下限是 0。代入下限时原函数在 0 处的值一般不是 0,会留下一个常数。这个常数保证结果在 t=0t = 0 从 0 开始(面积从 0 开始累积)。

例题

例 1 常数

求 1p[6u(t)]\frac1p[6u(t)]。

f(t)=6f(t) = 6,∫0t6 dα=6t\int_0^t 6\,d\alpha = 6t,所以 1p[6u(t)]=6t u(t)\frac1p[6u(t)] = 6t\,u(t):阶跃积分一次是斜坡 (ramp)。从图上看,宽 tt、高 6 的矩形面积就是 6t6t。

例题

例 2 正弦(课本例 6.6 是同一类题)

求 1p[sin⁡(3t)u(t)]\frac1p[\sin(3t)u(t)]。

第 1 步,写成从 0 开始的定积分:1p[sin⁡(3t)u(t)]=[∫0tsin⁡3α dα]u(t)\frac1p[\sin(3t)u(t)] = \Big[\int_0^t \sin 3\alpha\,d\alpha\Big]u(t)。

第 2 步,找原函数并代入上下限。sin⁡3α\sin 3\alpha 的原函数是 −13cos⁡3α-\tfrac13\cos 3\alpha:

∫0tsin⁡3α dα=[−13cos⁡3α]0t=−13cos⁡3t+13\begin{aligned} \int_0^t \sin 3\alpha\,d\alpha &= \Big[-\tfrac13\cos 3\alpha\Big]_0^t \\ &= -\tfrac13\cos 3t + \tfrac13 \end{aligned}

第二项 +13+\tfrac13 来自下限:−(−13cos⁡0)=13-(-\tfrac13\cos 0) = \tfrac13。

第 3 步,乘 u(t)u(t):1p[sin⁡(3t)u(t)]=13(1−cos⁡3t) u(t)\frac1p[\sin(3t)u(t)] = \tfrac13(1 - \cos 3t)\,u(t)。

核对:t=0t = 0 时结果是 13(1−1)=0\tfrac13(1 - 1) = 0 ✓,面积从 0 开始。如果只写不定积分 −13cos⁡3t-\tfrac13\cos 3t,t=0t = 0 时是 −13-\tfrac13:面积还没开始累积就是负数,不可能。结果在 0 和 23\tfrac23 之间摆动、从不为负,因为 sin⁡3t\sin 3t 每个周期都是正的半周先来。

π/3 2π/3 π −1 2/3 1 0 t sin3t·u (1−cos3t)/3
例题

例 3 指数

求 1p[10e−5tu(t)]\frac1p[10e^{-5t}u(t)]。

∫0t10e−5α dα=[−2e−5α]0t=−2e−5t+2\int_0^t 10e^{-5\alpha}\,d\alpha = \big[-2e^{-5\alpha}\big]_0^t = -2e^{-5t} + 2,乘 u(t)u(t):2(1−e−5t) u(t)2(1 - e^{-5t})\,u(t)。从 0 开始,慢慢涨到 2。2 正好是整条指数曲线下的总面积 ∫0∞10e−5α dα=10/5\int_0^\infty 10e^{-5\alpha}\,d\alpha = 10/5。

③ pp 作用在 f(t)u(t)f(t)u(t) 上:乘积法则多出一个 impulse

把 f(t)u(t)f(t)u(t) 当成两个函数的乘积求导:p[f u]=(pf) u+f (pu)p[f\,u] = (pf)\,u + f\,(pu)。第一项是平常的导数。第二项里 pupu 是阶跃的导数:uu 在 t=0t = 0 跳了 1,所以它的导数是一根面积为 1 的 impulse,p u=δ(t)p\,u = \delta(t);再由 sampling property f(t)δ(t)=f(0)δ(t)f(t)\delta(t) = f(0)\delta(t)(这两条的来历见 6·01)。合起来:

技巧 C p[f(t)u(t)]=f′(t) u(t)+f(0) δ(t)p\big[f(t)u(t)\big] = f'(t)\,u(t) + f(0)\,\delta(t)

先看 f(0)f(0):f(0)=0f(0) = 0 就没有 impulse;f(0)≠0f(0) \ne 0 就多一根面积(strength)为 f(0)f(0) 的 impulse。

从图上看:f(0)≠0f(0) \ne 0 时,波形 f(t)u(t)f(t)u(t) 在 t=0t = 0 从 0 一下跳到 f(0)f(0),跳了多少,导数里那根 impulse 的面积就是多少。f(0)=0f(0) = 0 时波形在 0 处连续(最多是个拐角),没有 impulse。这正是 6·01 里"每个跳变给一根 impulse,拐角不给"的特例。

例题

例 4 f(0) = 0:没有 impulse

求 p[sin⁡(3t)u(t)]p[\sin(3t)u(t)]。

f(t)=sin⁡3tf(t) = \sin 3t,f′(t)=3cos⁡3tf'(t) = 3\cos 3t,f(0)=sin⁡0=0f(0) = \sin 0 = 0。所以 p[sin⁡(3t)u(t)]=3cos⁡(3t) u(t)p[\sin(3t)u(t)] = 3\cos(3t)\,u(t),没有 impulse。波形 sin⁡(3t)u(t)\sin(3t)u(t) 在 t=0t = 0 只是从平的变成斜的(拐角),值不跳。

例题

例 5 f(0) ≠ 0:有 impulse(课本例 6.17)

求 p[cos⁡(t)u(t)]p[\cos(t)u(t)]。

f(t)=cos⁡tf(t) = \cos t,f′(t)=−sin⁡tf'(t) = -\sin t,f(0)=cos⁡0=1f(0) = \cos 0 = 1。所以 p[cos⁡(t)u(t)]=−sin⁡(t) u(t)+δ(t)p[\cos(t)u(t)] = -\sin(t)\,u(t) + \delta(t)。波形 cos⁡(t)u(t)\cos(t)u(t) 在 t=0t = 0 从 0 跳到 1,导数里就有一根面积为 1 的 impulse。

例题

例 6 为什么 impulse 不能丢:检查一个漏掉 δ 的答案

有人算 p[2e−3tu(t)]p[2e^{-3t}u(t)],只写了 −6e−3tu(t)-6e^{-3t}u(t)。用 1/p1/p 把这个答案积回去,和 2e−3tu(t)2e^{-3t}u(t) 比一比,找出他漏掉的那一项。

积回去。技巧 B:∫0t−6e−3α dα=[2e−3α]0t=2e−3t−2\int_0^t -6e^{-3\alpha}\,d\alpha = \big[2e^{-3\alpha}\big]_0^t = 2e^{-3t} - 2,得 (2e−3t−2) u(t)(2e^{-3t} - 2)\,u(t)。

和原波形比。2e−3tu(t)−(2e−3t−2)u(t)=2u(t)2e^{-3t}u(t) - (2e^{-3t} - 2)u(t) = 2u(t):积回去少了一个 2u(t)2u(t)。tt 很大时他的结果趋于 −2,原波形趋于 0,明显不对。

找出漏掉的项。漏掉的项 m(t)m(t) 积分后必须正好补上 2u(t)2u(t),所以 m=p[2u(t)]=2δ(t)m = p[2u(t)] = 2\delta(t)。正确答案是 p[2e−3tu(t)]=−6e−3tu(t)+2δ(t)p[2e^{-3t}u(t)] = -6e^{-3t}u(t) + 2\delta(t),和技巧 C 一致:f(0)=2f(0) = 2。

这根 impulse 的面积 2 正好是波形在 t=0t = 0 从 0 跳到 2 的那一跳。1/p1/p 要把这一跳"积"回来,靠的就是它。

p 和 1/p 互为逆运算(Davis 定义 6.1、6.2)

对 one-sided(单边)波形,并且把跳变产生的 impulse 都留着,p 1pf=1p pf=fp\,\frac1p f = \frac1p\,p f = f。为什么要单边:1p pf=∫−∞tf′(α) dα=f(t)−f(−∞)\frac1p\,p f = \int_{-\infty}^t f'(\alpha)\,d\alpha = f(t) - f(-\infty),只有 f(−∞)=0f(-\infty) = 0 才等于 f(t)f(t)。例如常数 12:p[12]=0p[12] = 0,积回去还是 0,回不到 12;单边的 12u(t)12u(t):p[12u]=12δp[12u] = 12\delta,积回去正好是 12u12u。Davis 还规定这些运算是 causal(因果)的:输入是单边的,输出也是单边的,没有刺激就没有反应。

④ tneatt^n e^{at} 的积分:表格法

∫α e−α dα\int \alpha\,e^{-\alpha}\,d\alpha 这类积分要用分部积分 ∫u dv=uv−∫v du\int u\,dv = uv - \int v\,du。每用一次,多项式降一次幂,直到求导成 0。课本附录 C 的 tabular integration(表格法)把这个过程排成一张表,做几次都不容易乱:

  1. 左列写多项式部分,一直求导到 0:α→1→0\alpha \to 1 \to 0。
  2. 右列写指数部分,一直积分(每次都不加常数):e−α→−e−α→e−αe^{-\alpha} \to -e^{-\alpha} \to e^{-\alpha}。
  3. 左列第 kk 行乘右列第 k+1k+1 行(斜着往右下乘),符号依次取 +,−,+,…+, -, +, \dots,加起来就是原函数。
符号求导列积分列
+α\alphae−αe^{-\alpha}
−11−e−α-e^{-\alpha}
00e−αe^{-\alpha}
∫αe−α dα=(+) α⋅(−e−α)+(−) 1⋅e−α=−αe−α−e−α\begin{aligned} \int \alpha e^{-\alpha}\,d\alpha &= (+)\,\alpha\cdot(-e^{-\alpha}) + (-)\,1\cdot e^{-\alpha} \\ &= -\alpha e^{-\alpha} - e^{-\alpha} \end{aligned}

代入上下限 0 和 tt:[−αe−α−e−α]0t=(−te−t−e−t)−(0−1)=1−e−t−te−t\big[-\alpha e^{-\alpha} - e^{-\alpha}\big]_0^t = (-te^{-t} - e^{-t}) - (0 - 1) = 1 - e^{-t} - te^{-t}。下限又贡献了常数 1。求导核对:ddt(1−e−t−te−t)=e−t−e−t+te−t=te−t\frac{d}{dt}(1 - e^{-t} - te^{-t}) = e^{-t} - e^{-t} + te^{-t} = te^{-t} ✓。

例题

例 7 表格法做三行:1p[t2e−tu(t)]\frac1p\big[t^2e^{-t}u(t)\big]

求 h(t)=1p[t2e−tu(t)]h(t) = \frac1p\big[t^2e^{-t}u(t)\big],再用 pp 核对。

第 1 步,技巧 B 写成从 0 开始的定积分:h(t)=[∫0tα2e−α dα]u(t)h(t) = \Big[\int_0^t \alpha^2e^{-\alpha}\,d\alpha\Big]u(t)。

第 2 步,表格法。多项式 α2\alpha^2 要求导两次才到 0,所以比上面 αe−α\alpha e^{-\alpha} 那张表多一行:

符号求导列积分列
+α2\alpha^2e−αe^{-\alpha}
−2α2\alpha−e−α-e^{-\alpha}
+22e−αe^{-\alpha}
00−e−α-e^{-\alpha}
∫α2e−α dα=α2(−e−α)−2α e−α+2(−e−α)=−(α2+2α+2) e−α\begin{aligned} \int \alpha^2e^{-\alpha}\,d\alpha &= \alpha^2(-e^{-\alpha}) - 2\alpha\,e^{-\alpha} \\ &\quad + 2(-e^{-\alpha}) \\ &= -(\alpha^2 + 2\alpha + 2)\,e^{-\alpha} \end{aligned}

第 3 步,代上下限。上限给 −(t2+2t+2)e−t-(t^2 + 2t + 2)e^{-t}。下限 α=0\alpha = 0 处原函数等于 −2-2,减去它就是 +2+2:

h(t)=[ 2−(t2+2t+2) e−t ] u(t)h(t) = \big[\,2 - (t^2 + 2t + 2)\,e^{-t}\,\big]\,u(t)

第 4 步,用 pp 核对(技巧 C)。括号在 t=0t = 0 等于 2−2=02 - 2 = 0,没有 impulse。求导:−(2t+2)e−t+(t2+2t+2)e−t=t2e−t-(2t + 2)e^{-t} + (t^2 + 2t + 2)e^{-t} = t^2e^{-t} ✓,回到了被积函数。

看极限:t→∞t \to \infty 时 h→2h \to 2,正是整条曲线 t2e−tt^2e^{-t} 下的总面积。下限带出的常数 2,就是最后攒下的总量;丢了它,hh 会从 −2 开始、最后趋于 0,面积成了负的。

⑤ 1/p21/p^2:积两次

1p2=1p⋅1p\frac{1}{p^2} = \frac1p\cdot\frac1p:先积一次,对结果再积一次。每一次都用技巧 B(从 0 积到 tt,乘 u(t)u(t)),每一次的下限都可能带出常数,都不能丢。

例题

例 8 1/p² 作用在指数上:两次都有下限常数(先自己做)

求 1p2[4e−2tu(t)]\frac{1}{p^2}\big[4e^{-2t}u(t)\big]。

看解答

第一次积分:

∫0t4e−2α dα=[−2e−2α]0t=−2e−2t+2\begin{aligned} \int_0^t 4e^{-2\alpha}\,d\alpha &= \big[-2e^{-2\alpha}\big]_0^t \\ &= -2e^{-2t} + 2 \end{aligned}

得 2(1−e−2t) u(t)2(1 - e^{-2t})\,u(t)。下限贡献了 +2+2。

第二次积分(对第一次的结果再用一次技巧 B):

∫0t2(1−e−2α) dα=[ 2α+e−2α ]0t=2t+e−2t−1\begin{aligned} \int_0^t 2(1 - e^{-2\alpha})\,d\alpha &= \big[\,2\alpha + e^{-2\alpha}\,\big]_0^t \\ &= 2t + e^{-2t} - 1 \end{aligned}

下限这次贡献 −1-1(e0=1e^{0} = 1 被减掉)。结果是 (2t−1+e−2t) u(t)(2t - 1 + e^{-2t})\,u(t)。

核对:t=0t = 0 时 0−1+1=00 - 1 + 1 = 0 ✓。求导一次:(2−2e−2t) u(2 - 2e^{-2t})\,u(括号在 0 处为 0,无 impulse);再求导:4e−2tu4e^{-2t}u(2−2=02 - 2 = 0,无 impulse)✓。

形状:第一次积分趋于常数 2,所以第二次积分趋于斜率 2 的直线 2t−12t - 1。那个 −1-1 是因为第一次积分要花时间才涨到 2,比"一开始就是 2"少积了 1。如果漏掉第一次的 +2+2,第二次会得到 (e−2t−1) u(t)(e^{-2t} - 1)\,u(t):不涨反降,一看就错。

这一节最常见的三个错

(1) 积分只写原函数,忘了代下限 0 带出的常数(例 2、例 8)。
(2) 求导忘了看 f(0)f(0),丢掉 f(0)δ(t)f(0)\delta(t)(例 5、例 6)。
(3) 结果忘了乘 u(t)u(t):13(1−cos⁡3t)\tfrac13(1 - \cos 3t) 在 t<0t \lt 0 时并不是 0,乘上 u(t)u(t) 才是正确的单边波形。