C·08

和指数 convolution:例 15.14、例 13.55、闭式

About 8 minBook · A 15.14 · D 13.55

一阶电路(RC、RL)的 impulse response 是 e−atu(t)e^{-at}u(t),a=1/τa = 1/\tau,τ\tau 是 time constant。和它做 convolution 有一个固定的套路:把 e−ate^{-at} 提到积分外面。这一节先推出这个"指数闭式",再用它算指数输入、脉冲输入和一段正弦输入。

指数闭式 (closed form)

xx 是 one-sided 的时候,对 t>0t \gt 0:

x(t)∗e−atu(t)=e−at∫0tx(λ) eaλ dλ\begin{aligned} &x(t) * e^{-at}u(t) \\ &\quad= e^{-at}\int_0^{t} x(\lambda)\,e^{a\lambda}\,d\lambda \end{aligned}

t<0t \lt 0 时结果为 0。如果 xx 在某个时刻 TT 之后为 0,上限取 min⁡(t,T)\min(t, T)。

推导:

  1. 写出 h(t−λ)=e−a(t−λ) u(t−λ)h(t-\lambda) = e^{-a(t-\lambda)}\,u(t-\lambda)。
  2. xx one-sided、hh causal,所以积分区间是 0<λ<t0 \lt \lambda \lt t(C·03),在这个区间上 u(t−λ)=1u(t-\lambda) = 1。
  3. 拆指数:e−a(t−λ)=e−at eaλe^{-a(t-\lambda)} = e^{-at}\,e^{a\lambda}。
  4. e−ate^{-at} 不含积分变量 λ\lambda,是常数,提到积分外面。积分里只剩 x(λ) eaλx(\lambda)\,e^{a\lambda}。
例题

例 1 Davis 例 13.55:两个指数

求 x(t)=e−tu(t)x(t) = e^{-t}u(t) 与 y(t)=2e−2tu(t)y(t) = 2e^{-2t}u(t) 的 convolution。

-1 t 3 1 2 0 α y(t−α) x(α) = e-α
t=1.5t=1.5:y(t−α)y(t-\alpha) 是翻转后的指数,在 α=t\alpha = t 处最高(2),往左衰减。紫色是乘积。
1 2 3 4 5 .25 .5 0 t 2(e-t − e-2t) 峰值 t = ln 2
结果从 0 升起再衰减。

上下限。两个都是 one-sided 的,t≤0t \le 0 时不重叠,y=0y = 0;t>0t\gt 0 时重叠 0<α<t0\lt \alpha\lt t。乘积不是矩形或三角形,要真的积分。

用闭式(把 2e−2tu2e^{-2t}u 当作指数那一个,a=2a = 2,前面的 2 提出来):

∫0te−α 2e−2(t−α) dα=2e−2t∫0te−αe2α dα=2e−2t∫0teα dα=2e−2t (et−1)\begin{aligned} &\int_0^t e^{-\alpha}\,2e^{-2(t-\alpha)}\,d\alpha \\ &\quad= 2e^{-2t}\int_0^t e^{-\alpha}e^{2\alpha}\,d\alpha \\ &\quad= 2e^{-2t}\int_0^t e^{\alpha}\,d\alpha \\ &\quad= 2e^{-2t}\,(e^t - 1) \end{aligned}

化简并关掉负时间。2e−2t(et−1)=2(e−t−e−2t)2e^{-2t}(e^t - 1) = 2(e^{-t} - e^{-2t}):

e−tu(t)⊗2e−2tu(t)=2(e−t−e−2t)u(t)e^{-t}u(t)\otimes 2e^{-2t}u(t) = 2\big(e^{-t}-e^{-2t}\big)u(t)
检查

t=0t=0:0,连续地从 0 开始 ✓。t=1t=1:2(0.3679−0.1353)=0.4652(0.3679-0.1353) = 0.465。峰值:令导数 2(−e−t+2e−2t)=02(-e^{-t} + 2e^{-2t}) = 0,得 e−t=1/2e^{-t} = 1/2,t=ln⁡2t = \ln 2,峰值 2(12−14)=122(\tfrac12 - \tfrac14) = \tfrac12。面积:1×1=11 \times 1 = 1,∫0∞2(e−t−e−2t)dt=2(1−12)=1\int_0^\infty 2(e^{-t}-e^{-2t})dt = 2(1 - \tfrac12) = 1 ✓。

把例 1 的做法对一般的 aa、bb 做一遍,得到一个可以直接套用的通式:

两个指数的通式

a≠ba \ne b:

e−atu∗e−btu=e−at−e−btb−a u(t)e^{-at}u * e^{-bt}u = \frac{e^{-at}-e^{-bt}}{b-a}\,u(t)

a=ba = b:

e−atu∗e−atu=t e−atu(t)e^{-at}u * e^{-at}u = t\,e^{-at}u(t)

推导:用闭式,把 e−btue^{-bt}u 当指数那一个:e−bt∫0te−aλebλ dλ=e−bt⋅e(b−a)t−1b−a=e−at−e−btb−ae^{-bt}\int_0^t e^{-a\lambda}e^{b\lambda}\,d\lambda = e^{-bt}\cdot\dfrac{e^{(b-a)t}-1}{b-a} = \dfrac{e^{-at}-e^{-bt}}{b-a}。a=ba = b 时被积函数是 1,积分是 tt,结果 t e−att\,e^{-at}(C·03例 1 是 a=1a = 1 的情形)。用例 1 检查:a=1, b=2a=1,\ b=2,再乘前面的 2,得 2(e−t−e−2t)2(e^{-t}-e^{-2t}) ✓。

例题

例 2 Alexander 例 15.14:RL 电路的脉冲响应

is(t) 1 Ω 1 H io
电流源 is(t)i_s(t) 并联 1 Ω 和 1 H,求 inductor 电流 io(t)i_o(t)。输入 is=u(t)−u(t−2)i_s = u(t) - u(t-2)(高 1、0 到 2 s 的脉冲)。

书上用 Laplace 求出 h(t)=e−tu(t)h(t) = e^{-t}u(t)(还没学 Laplace 的话,直接把它当已知)。

解法 A:拆成两个阶跃(代数)。输入 =u(t)−u(t−2)= u(t) - u(t-2)。分配律:

io=u∗h⏟io′  −  u(t−2)∗h⏟io′′i_o = \underbrace{u * h}_{i_o'} \;-\; \underbrace{u(t-2) * h}_{i_o''}

io′i_o' 是 hh 的 running integral(C·07):io′=∫0te−λdλ=(1−e−t)u(t)i_o' = \int_0^t e^{-\lambda}d\lambda = (1 - e^{-t})u(t)。io′′i_o'' 的输入是晚 2 秒的阶跃,time-invariant 说输出也晚 2 秒:io′′=(1−e−(t−2))u(t−2)i_o'' = \big(1 - e^{-(t-2)}\big)u(t-2)。

0<t<20\lt t\lt 2:u(t−2)=0u(t-2) = 0,只有第一项,io=1−e−ti_o = 1 - e^{-t}。t>2t \gt 2:两项都在:

io=(1−e−t)−(1−e−(t−2))=e−(t−2)−e−t=(e2−1) e−t A\begin{aligned} i_o &= (1-e^{-t}) - \big(1-e^{-(t-2)}\big) \\ &= e^{-(t-2)} - e^{-t} = (e^2-1)\,e^{-t}\ \text{A} \end{aligned}

解法 B:翻转脉冲(图解)。固定 h(λ)=e−λh(\lambda) = e^{-\lambda},翻转平的 isi_s:is(t−λ)=1i_s(t-\lambda) = 1 在 [t−2, t][t-2,\ t] 上。临界时刻 0 和 2,每一段一张图:

-1 1 2 3 4 5 1 0 λ is(t−λ) h(λ)
① 0<t<20\lt t\lt 2:重叠 [0, t][0,\ t] → ∫0te−λdλ=1−e−t\int_0^t e^{-\lambda}d\lambda = 1-e^{-t}
-1 1 2 3 4 5 1 0 λ is(t−λ) h(λ)
② t>2t\gt 2:重叠 [t−2, t][t-2,\ t] → ∫t−2te−λdλ=(e2−1) e−t\int_{t-2}^t e^{-\lambda}d\lambda = (e^2-1)\,e^{-t}
1 2 3 4 5 6 .5 .865 1 0 t (s) 输入 is 输出 io io(2) = 1 − e-2
输入脉冲(灰虚线)和响应(青色):0 到 2 s 充上去,之后指数衰减。
查连续和面积

t=2t=2:1−e−2=0.86471 - e^{-2} = 0.8647;(e2−1)e−2=1−e−2=0.8647(e^2-1)e^{-2} = 1 - e^{-2} = 0.8647 ✓。t=3t=3:6.389×e−3=0.3186.389 \times e^{-3} = 0.318。面积:2×1=22 \times 1 = 2;∫02(1−e−t)dt+∫2∞(e2−1)e−tdt=1.1353+0.8647=2\int_0^2(1-e^{-t})dt + \int_2^\infty (e^2-1)e^{-t}dt = 1.1353 + 0.8647 = 2 ✓。

例题

例 3 半个正弦 * 指数(闭式 + 两次分部积分)

x(t)=sin⁡tx(t) = \sin t(0<t<π0 \lt t \lt \pi,半个正弦),其他为 0;h(t)=e−tu(t)h(t) = e^{-t}u(t)。求 y=x∗hy = x * h。

代入闭式。a=1a = 1。xx 在 T=πT = \pi 之后为 0,所以上限是 min⁡(t,π)\min(t, \pi):

y(t)=e−t∫0min⁡(t,π)eλsin⁡λ dλy(t) = e^{-t}\int_0^{\min(t,\pi)} e^{\lambda}\sin\lambda\,d\lambda
-1 π/2 π 4 5 1 0 λ x(λ) h(t−λ)
① 0<t<π0\lt t\lt \pi:重叠 [0, t][0,\ t],上限是 tt
-1 π/2 π 4 5 1 0 λ x(λ) h(t−λ)
② t>πt\gt \pi:半个正弦整个被盖住,重叠固定在 [0, π][0,\ \pi]

原函数(分部积分两次)。记 I=∫eλsin⁡λ dλI = \int e^{\lambda}\sin\lambda\,d\lambda。第一次让 sin⁡λ\sin\lambda 求导:I=eλsin⁡λ−∫eλcos⁡λ dλI = e^{\lambda}\sin\lambda - \int e^{\lambda}\cos\lambda\,d\lambda。第二次让 cos⁡λ\cos\lambda 求导:∫eλcos⁡λ dλ=eλcos⁡λ+I\int e^{\lambda}\cos\lambda\,d\lambda = e^{\lambda}\cos\lambda + I。代回去,II 出现在两边,移项:

I=eλsin⁡λ−eλcos⁡λ−II=12 eλ(sin⁡λ−cos⁡λ)\begin{aligned} I &= e^{\lambda}\sin\lambda - e^{\lambda}\cos\lambda - I \\ I &= \tfrac12\,e^{\lambda}(\sin\lambda - \cos\lambda) \end{aligned}

验证:对右边求导,12eλ(sin⁡λ−cos⁡λ)+12eλ(cos⁡λ+sin⁡λ)=eλsin⁡λ\tfrac12 e^{\lambda}(\sin\lambda - \cos\lambda) + \tfrac12 e^{\lambda}(\cos\lambda + \sin\lambda) = e^{\lambda}\sin\lambda ✓。它在 0 处的值是 −12-\tfrac12。

① 0<t<π0 \lt t \lt \pi。[I]0t=12et(sin⁡t−cos⁡t)+12[I]_0^t = \tfrac12 e^{t}(\sin t - \cos t) + \tfrac12,乘 e−te^{-t}:

y=12(sin⁡t−cos⁡t)+12e−ty = \tfrac12(\sin t - \cos t) + \tfrac12 e^{-t}

② t>πt \gt \pi。sin⁡π=0\sin\pi = 0、cos⁡π=−1\cos\pi = -1,[I]0π=12eπ+12[I]_0^\pi = \tfrac12 e^{\pi} + \tfrac12,是常数,乘 e−te^{-t}:

y=eπ+12 e−t≈12.07 e−ty = \frac{e^{\pi}+1}{2}\,e^{-t} \approx 12.07\,e^{-t}
π/2 π 4 5 6 .5 1 0 t 输入 sin t ½(sin t − cos t) y(t)
灰虚线:输入。青色:yy。紫色:12(sin⁡t−cos⁡t)\tfrac12(\sin t - \cos t),e−t/2e^{-t}/2 消失以后 yy 贴着它走;t>πt \gt \pi 输入关掉,yy 指数衰减。两个黑点:峰值 0.756(t≈2.28t \approx 2.28)和 y(π)=0.522y(\pi) = 0.522。
检查

起点 t=0t = 0:−12+12=0-\tfrac12 + \tfrac12 = 0 ✓。连续 t=πt = \pi:左段 12(0+1)+12e−π=0.5216\tfrac12(0 + 1) + \tfrac12 e^{-\pi} = 0.5216;右段 eπ+12e−π=12+12e−π=0.5216\tfrac{e^\pi + 1}{2}e^{-\pi} = \tfrac12 + \tfrac12 e^{-\pi} = 0.5216 ✓。面积:xx 面积 ∫0πsin⁡t dt=2\int_0^\pi \sin t\,dt = 2,hh 面积 1,yy 面积应为 2;∫0πy dt=1+12(1−e−π)\int_0^\pi y\,dt = 1 + \tfrac12(1 - e^{-\pi}),∫π∞y dt=12(1+e−π)\int_\pi^\infty y\,dt = \tfrac12(1 + e^{-\pi}),合计 2 ✓。

看结果的形状。12(sin⁡t−cos⁡t)=12sin⁡(t−π4)\tfrac12(\sin t - \cos t) = \tfrac{1}{\sqrt2}\sin(t - \tfrac{\pi}{4}):输出还是同频率的正弦,幅度缩小到 1/21/\sqrt2,晚了 π/4\pi/4。所以 yy 的峰值(0.756,在 t≈2.28t \approx 2.28)比输入的峰值(t=π/2t = \pi/2)晚。以后学 phasor 时会看到,这正是一阶低通在 ω=1\omega = 1 处的 1/21/\sqrt2 和 −45∘-45^\circ。

练习 Alexander 练习 15.13(书上的答案)

g(t)=1g(t) = 1(0<t<10\lt t\lt 1),f(t)=3e−tu(t)f(t) = 3e^{-t}u(t)。求 g∗fg * f。

提示

和例 2 一样,只是脉冲宽 1、指数前面有个 3。翻转矩形 gg,或者用闭式。

0<t<10\lt t\lt 1:∫0t3e−λdλ=3(1−e−t)\int_0^t 3e^{-\lambda}d\lambda = 3(1-e^{-t})。t>1t>1:∫t−1t3e−λdλ=3(e−(t−1)−e−t)=3(e−1)e−t\int_{t-1}^t 3e^{-\lambda}d\lambda = 3(e^{-(t-1)} - e^{-t}) = 3(e-1)e^{-t}。查连续 t=1t=1:3(1−e−1)=1.8963(1-e^{-1}) = 1.896,3(e−1)e−1=1.8963(e-1)e^{-1} = 1.896 ✓。面积 1×3=31 \times 3 = 3。

容易错
  • 提到外面的是 e−ate^{-at}(含 tt),留在积分里的是 e+aλe^{+a\lambda}(正号)。
  • 输入有限长时上限是 min⁡(t,T)\min(t, T):t>Tt \gt T 以后积分是常数,结果是"常数 × e−ate^{-at}",纯衰减。
  • 最后别忘了乘 u(t)u(t),或者写明 t<0t \lt 0 时为 0。