6·02

信号分解 (signal decomposition):用阶跃、斜坡、impulse 写波形

About 33 minLecture 12

上一节 6·01 有了积木,这一节学两件事:一是看着一个分段的波形(题目给的图),把它写成 uu、rr(必要时还有 δ\delta)的加减;二是反过来,看到这样一个式子,画出波形、写出分段表达式。这样做的原因在最后一部分:电路是 linear、time-invariant 的,只要知道它对 u(t)u(t) 的反应,对 u(t−a)u(t - a) 的反应就是把它平移 aa,对整个波形的反应就是逐项相加。

① 方法:从左往右扫 课上第 12 讲

关键事实:u(t−a)u(t - a) 和 r(t−a)r(t - a) 在 t≤at \le a 都是 0。所以在一个式子里加上 Δ⋅u(t−a)\Delta \cdot u(t - a),aa 左边的波形一点不变,aa 右边整体抬高 Δ\Delta:这是一次跳变,斜率不变。加上 Δk⋅r(t−a)\Delta k \cdot r(t - a),左边也不变,右边的斜率多了 Δk\Delta k;它从 0 开始长,所以不产生跳变。于是可以从最左边出发,每走到一个"特殊时刻"就补一块积木,把波形修正成正确的样子,再继续往右走。

把波形写成 u 和 r 的四步
  1. 读图。写出分段定义,标出所有特殊时刻 aa:值跳变的点,和斜率改变的点。
  2. 在每个 aa 量两样东西。跳了多少:Δ\Delta = (aa 右边的值)−(aa 左边的值)。斜率变了多少:Δk\Delta k = (aa 右边的斜率)−(aa 左边的斜率)。
  3. 补积木。每个跳变加 Δ⋅u(t−a)\Delta \cdot u(t - a),每个斜率变化加 Δk⋅r(t−a)\Delta k \cdot r(t - a)。同一个 aa 两样都有,就两块都加。往下跳 Δ\Delta 为负;斜率变小 Δk\Delta k 为负。
  4. 相加并检查。每一段取一个 tt 代进去,和分段定义比较。

第 2 步量的是变化量,不是新的斜率:斜率从 +1 变成 −1,要加的是 (−1)−(+1)=−2(-1) - (+1) = -2 倍的 ramp,不是 −1 倍。

② 例题:读图写式子

例 1、2、3、6、8、10 的解答里都有一排小图:灰色虚线是要拼出的目标波形,青色是"到这一步为止已经写下的各项之和"。每补一块积木,青色线就在那个时刻之后被修正一次;最后一张图里青色和灰色完全重合,式子就写对了。

例 1

电源电压跌落 1 s

-1 1 2 3 4 1 3 0 t v(t)
0 到 1 为 3,1 到 2 跌到 1,之后恢复 3。

一个 3 V 的电源在 t=0t = 0 接通,在 1<t≤21 \lt t \le 2 s 期间跌落到 1 V,之后恢复 3 V 并一直保持。把 v(t)v(t) 写成 u 和 r 的组合。

读图。v=0v = 0(t≤0t \le 0),33(0<t≤10 \lt t \le 1),11(1<t≤21 \lt t \le 2),33(t>2t \gt 2)。各段都是水平的,斜率处处是 0,所以只有跳变:特殊时刻 0、1、2,跳变量依次是 3−0=+33 - 0 = +3、1−3=−21 - 3 = -2、3−1=+23 - 1 = +2。每个跳变补一个 Δ⋅u(t−a)\Delta \cdot u(t - a):

-1 1 2 3 4 1 3 0 t
① t=0t = 0 跳 +3:补 3u(t)3u(t)。
-1 1 2 3 4 1 3 0 t
② t=1t = 1 跳 −2:补 −2u(t−1)-2u(t - 1),右边整体降 2。
-1 1 2 3 4 1 3 0 t
③ t=2t = 2 跳 +2:补 2u(t−2)2u(t - 2),和目标重合 ✓
v(t)=3u(t)−2u(t−1)+2u(t−2)v(t) = 3u(t) - 2u(t - 1) + 2u(t - 2)

检查。t=−1t = -1:0 ✓。t=0.5t = 0.5:3 ✓。t=1.5t = 1.5:3−2=13 - 2 = 1 ✓。t=3t = 3:3−2+2=33 - 2 + 2 = 3 ✓。这个波形最后不回到 0,所以三个跳变加起来不是 0,而是终值 3−2+2=33 - 2 + 2 = 3。

例 2

斜升到 2,然后保持

-1 1 2 3 4 1 2 0 t f(t)
0 到 2 以斜率 1 上升到 2,之后保持 2。

写成 u 和 r 的组合。

读图。f(t)=0f(t) = 0(t≤0t \le 0),tt(0<t≤20 \lt t \le 2),22(t>2t \gt 2)。特殊时刻 0 和 2,两处的值左右相等,都没有跳变,只有斜率变化。

-1 1 2 3 4 1 2 0 t
① t=0t = 0 斜率 0 → 1:补 r(t)r(t)。它会一直往上走(出了画面)。
-1 1 2 3 4 1 2 0 t
② t=2t = 2 斜率 1 → 0:补 −r(t−2)-r(t - 2),把继续的增长抵消掉 ✓

结果 f(t)=r(t)−r(t−2)f(t) = r(t) - r(t - 2)。检查:t=1t = 1:1−0=11 - 0 = 1 ✓。t=3t = 3:3−1=23 - 1 = 2 ✓。展开 t>2t \gt 2:t−(t−2)=2t - (t - 2) = 2,正是第 2 张图里"抵消"的意思。

例 3

冲过头再回落:0 → 2 → 1,然后保持

-1 1 2 3 4 1 2 0 t f(t)
0 到 1 以斜率 2 升到 2,1 到 2 以斜率 −1 落到 1,之后保持 1。

写成 u 和 r 的组合。(形状像一个会过冲的 step response 的折线近似。)

读图。f(t)=0f(t) = 0(t≤0t \le 0),2t2t(0<t≤10 \lt t \le 1),3−t3 - t(1<t≤21 \lt t \le 2),11(t>2t \gt 2)。在 t=1t = 1:2⋅1=22 \cdot 1 = 2,3−1=23 - 1 = 2,相等;在 t=2t = 2:3−2=13 - 2 = 1,相等。所以没有跳变,只有斜率变化。各段斜率依次是 0、+2、−1、0。

-1 1 2 3 4 1 2 0 t
① t=0t = 0:Δk=2−0=+2\Delta k = 2 - 0 = +2,补 2r(t)2r(t)。
-1 1 2 3 4 1 2 0 t
② t=1t = 1:Δk=−1−2=−3\Delta k = -1 - 2 = -3,补 −3r(t−1)-3r(t - 1),之后斜率是 −1,一直往下。
-1 1 2 3 4 1 2 0 t
③ t=2t = 2:Δk=0−(−1)=+1\Delta k = 0 - (-1) = +1,补 r(t−2)r(t - 2),停在 1 ✓
f(t)=2r(t)−3r(t−1)+r(t−2)f(t) = 2r(t) - 3r(t - 1) + r(t - 2)

第 2 张图是关键:斜率从 +2 变成 −1,要补的是变化量 −3-3 倍的 ramp,不是 −1-1 倍。检查:t=0.5t = 0.5:1 ✓。t=1.5t = 1.5:3−1.5=1.53 - 1.5 = 1.5 ✓。t=3t = 3:6−6+1=16 - 6 + 1 = 1 ✓。三个 Δk\Delta k 加起来 2−3+1=02 - 3 + 1 = 0,所以最后斜率是 0;但终值是 1,不是 0:斜率回到 0 不等于值回到 0。

例 4

阶梯:只有跳变,有一段在 0 以下

1.5 2.5 4 −1 1 2 0 t f(t)
0 到 1.5 为 2,1.5 到 2.5 为 −1,2.5 到 4 为 1,之后为 0。

写成阶跃的组合。

看解答
  1. 读图。各段都是水平的,所以只有跳变:t=0t = 0 从 0 到 2(Δ=+2\Delta = +2),t=1.5t = 1.5 从 2 到 −1(−3-3),t=2.5t = 2.5 从 −1 到 1(+2+2),t=4t = 4 从 1 到 0(−1-1)。
  2. 相加。f(t)=2u(t)−3u(t−1.5)+2u(t−2.5)−u(t−4)f(t) = 2u(t) - 3u(t - 1.5) + 2u(t - 2.5) - u(t - 4)。
  3. 检查。t=1t = 1:2 ✓。t=2t = 2:2−3=−12 - 3 = -1 ✓。t=3t = 3:2−3+2=12 - 3 + 2 = 1 ✓。t=5t = 5:2−3+2−1=02 - 3 + 2 - 1 = 0 ✓。

另一种写法:每一层乘一个"窗口"。每一段的值乘上该段的脉冲:

f=2 [u(t)−u(t−1.5)]−1 [u(t−1.5)−u(t−2.5)]+1 [u(t−2.5)−u(t−4)]\begin{aligned} f &= 2\,[u(t) - u(t-1.5)] \\ &\quad - 1\,[u(t-1.5) - u(t-2.5)] \\ &\quad + 1\,[u(t-2.5) - u(t-4)] \end{aligned}

合并同类项:u(t−1.5)u(t - 1.5) 的系数 −2−1=−3-2 - 1 = -3,u(t−2.5)u(t - 2.5) 的系数 1+1=21 + 1 = 2,u(t−4)u(t - 4) 的系数 −1-1。和上面一样。窗口写法不用量跳变,对曲线段也适用(见例 9)。

例 5

梯形:只有斜率变化

1 2 3 4 5 6 3 0 t f(t)
1 到 2 升到 3,2 到 4 保持 3,4 到 5 降回 0。

写成 ramp 的组合。

看解答
  1. 读图。没有跳变。各段斜率:t<1t \lt 1 为 0;11–22 为 3/1=33/1 = 3;22–44 为 0;44–55 为 −3-3;之后为 0。
  2. 斜率变化。t=1t = 1:3−0=+33 - 0 = +3;t=2t = 2:0−3=−30 - 3 = -3;t=4t = 4:−3−0=−3-3 - 0 = -3;t=5t = 5:0−(−3)=+30 - (-3) = +3。
  3. 相加。f(t)=3r(t−1)−3r(t−2)−3r(t−4)+3r(t−5)f(t) = 3r(t - 1) - 3r(t - 2) - 3r(t - 4) + 3r(t - 5)。
  4. 检查。t=1.5t = 1.5:3(0.5)=1.53(0.5) = 1.5 ✓。t=3t = 3:3(2)−3(1)=33(2) - 3(1) = 3 ✓。t=4.5t = 4.5:3(3.5)−3(2.5)−3(0.5)=10.5−7.5−1.5=1.53(3.5) - 3(2.5) - 3(0.5) = 10.5 - 7.5 - 1.5 = 1.5 ✓(从 3 降到 0 的一半)。t=6t = 6:15−12−6+3=015 - 12 - 6 + 3 = 0 ✓。
例 6

课上例题(第 12 讲):跳变和斜率变化同时出现

1 2 3 4 5 1 2 0 t f(t)
0 到 1 斜率 1,1 到 2 平在 1,t = 2 跳到 2,然后斜率 −1,t = 4 回到 0。

写成 u 和 r 的组合。

读图。f(t)=0f(t) = 0(t≤0t \le 0),tt(0<t≤10 \lt t \le 1),11(1<t≤21 \lt t \le 2),4−t4 - t(2<t≤42 \lt t \le 4),00(t>4t \gt 4)。特殊时刻 0、1、2、4;只有 t=2t = 2 有跳变(1 → 2),其余三处只是拐角(4−t4 - t 在 t=4t = 4 正好是 0)。从左往右,在每个时刻补积木:

1 2 3 4 5 1 2 0 t
① t=0t = 0:不跳,斜率 0 → 1,补 r(t)r(t)。
1 2 3 4 5 1 2 0 t
② t=1t = 1:不跳,斜率 1 → 0,补 −r(t−1)-r(t - 1)。
1 2 3 4 5 1 2 0 t
③ t=2t = 2:跳 +1 并且斜率 0 → −1,两块一起补:u(t−2)−r(t−2)u(t - 2) - r(t - 2)。
1 2 3 4 5 1 2 0 t
④ t=4t = 4:不跳,斜率 −1 → 0,补 r(t−4)r(t - 4),和目标重合 ✓

相加:

f(t)=r(t)−r(t−1)+u(t−2)−r(t−2)+r(t−4)\begin{aligned} f(t) &= r(t) - r(t - 1) + u(t - 2) \\ &\quad - r(t - 2) + r(t - 4) \end{aligned}

逐段检查(每段只留已经开始的项):0<t≤10 \lt t \le 1:tt ✓。1<t≤21 \lt t \le 2:t−(t−1)=1t - (t - 1) = 1 ✓。2<t≤42 \lt t \le 4:t−(t−1)+1−(t−2)=4−tt - (t - 1) + 1 - (t - 2) = 4 - t ✓,t=2+t = 2^+ 时是 2。t>4t \gt 4:(4−t)+(t−4)=0(4 - t) + (t - 4) = 0 ✓。在 t=2t = 2 这一点,u(0)=0u(0) = 0、r(0)=0r(0) = 0,所以 f(2)=1f(2) = 1,取的是左边的值,图上实心点在 1。

③ 同一个信号用 δ 和 p 再写一遍

6·01 说过,积木之间只差一次求导或积分:

p u(t−a)=δ(t−a)p r(t−a)=u(t−a)p2r(t−a)=δ(t−a)\begin{gathered} p\,u(t - a) = \delta(t - a) \\ p\,r(t - a) = u(t - a) \\ p^2 r(t - a) = \delta(t - a) \end{gathered}

反过来 u(t−a)=1pδ(t−a)u(t - a) = \frac1p\delta(t - a),r(t−a)=1p2δ(t−a)r(t - a) = \frac{1}{p^2}\delta(t - a)。pp 是 linear 的:对和逐项求导,常数提出来。所以一个写成 uu、rr 之和的信号,求导就是把每块积木"降一级":uu 变成 δ\delta,rr 变成 uu。这和技巧 J 是一回事:跳变 Δ\Delta 在导数里给 Δ δ\Delta\,\delta;斜率变化 Δk\Delta k 在导数里给一个高 Δk\Delta k 的台阶。

例 7

把例 1、2、3、6 的导数写成 δ 和 u 之和

对例 1、2、3、6 的结果求导(例 2、3 求到二阶),并反过来用 1/p 写出原信号。

(a) 例 1 的电压跌落 v=3u(t)−2u(t−1)+2u(t−2)v = 3u(t) - 2u(t - 1) + 2u(t - 2)。每个 uu 换成 δ\delta:

p v=3 δ(t)−2 δ(t−1)+2 δ(t−2)p\,v = 3\,\delta(t) - 2\,\delta(t - 1) + 2\,\delta(t - 2)
-1 1 2 3 4 −2 2 3 0 t 3δ(t) −2δ(t−1) 2δ(t−2)
三个 impulse,面积就是三次跳变:t=0t = 0 跳 +3,t=1t = 1 跳 −2,t=2t = 2 跳 +2。

反过来 v=1p[ ⋯ ]v = \frac1p[\,\cdots]:从左往右累加面积,t<0t \lt 0 是 0;0<t<10 \lt t \lt 1 已经加了 3;1<t<21 \lt t \lt 2 是 3−2=13 - 2 = 1;t>2t \gt 2 是 1+2=31 + 2 = 3。正好回到原波形 ✓。

(b) 例 2 f=r(t)−r(t−2)f = r(t) - r(t - 2)。每个 rr 换成 uu:p f=u(t)−u(t−2)p\,f = u(t) - u(t - 2),一个高 1 的脉冲,就是 ff 的斜率(00–22 为 1,其余为 0)。再求一次:p2f=δ(t)−δ(t−2)p^2 f = \delta(t) - \delta(t - 2)。所以 f=1p2[δ(t)−δ(t−2)]f = \frac{1}{p^2}[\delta(t) - \delta(t - 2)]。

(c) 例 3 的过冲波形 f=2r(t)−3r(t−1)+r(t−2)f = 2r(t) - 3r(t - 1) + r(t - 2)。每个 rr 直接换成 δ\delta,得 p2f=2 δ(t)−3 δ(t−1)+δ(t−2)p^2 f = 2\,\delta(t) - 3\,\delta(t - 1) + \delta(t - 2),即

f=1p2[2 δ(t)−3 δ(t−1)+δ(t−2)]f = \frac{1}{p^2}\big[2\,\delta(t) - 3\,\delta(t - 1) + \delta(t - 2)\big]

impulse 的面积 2、−3、1 正是三个拐角的 Δk\Delta k。中间那一级 p f=2u(t)−3u(t−1)+u(t−2)p\,f = 2u(t) - 3u(t - 1) + u(t - 2) 就是 ff 的斜率:00–11 为 2,11–22 为 2−3=−12 - 3 = -1,之后为 0。

(d) 例 6 的课上波形:r→ur \to u,u→δu \to \delta,逐项换:

p f=u(t)−u(t−1)+δ(t−2)−u(t−2)+u(t−4)\begin{aligned} p\,f &= u(t) - u(t - 1) + \delta(t - 2) \\ &\quad - u(t - 2) + u(t - 4) \end{aligned}

读这个结果:uu 那几项拼成"斜率函数"(00–11 为 1,11–22 为 0,22–44 为 1−1−1=−11 - 1 - 1 = -1,之后为 0),再加上 t=2t = 2 处那一跳给出的 δ(t−2)\delta(t - 2)。

1 2 3 4 5 −1 1 0 t 斜率部分 δ(t−2)
例 6 的 p fp\,f:青色是斜率部分,橙色箭头是 t=2t = 2 处面积 1 的 impulse。

④ 反过来:从式子画波形

从式子到分段表达式
  1. 列出每一项开始的时刻(u(t−a)u(t - a)、r(t−a)r(t - a) 里的 aa),排好序,把 tt 轴分成几段。
  2. 在每一段里,已经开始的项都在起作用:u(t−a)u(t - a) 贡献常数 1,r(t−a)r(t - a) 贡献 t−at - a;还没开始的项是 0。把起作用的项加起来、化简,就是这一段的表达式。
  3. 在分界点比较左右两段的值,判断有没有跳变,然后画图。
例 8

画出 2u(t) − u(2 − t) 和 r(t) − 2u(t − 1) − r(t − 3)

两个式子分别画出来,写成分段形式。(a) 里有一个 6·01 ③ 讲的"关掉"型阶跃 u(2−t)u(2 - t),它不是在 2 开始,而是在 2 结束,所以方框里第 1 步的"开始时刻"要改成"开关时刻"。

看解答

(a) f=2u(t)−u(2−t)f = 2u(t) - u(2 - t)。开关时刻:2u(t)2u(t) 在 0 打开;u(2−t)u(2 - t) 在 t<2t \lt 2 时为 1、在 2 关掉。所以 tt 轴分成三段:

  • t≤0t \le 0:2u(t)=02u(t) = 0,u(2−t)=1u(2 - t) = 1,f=−1f = -1。波形在负时间不是 0:它不是 one-sided waveform。
  • 0<t<20 \lt t \lt 2:2−1=12 - 1 = 1。
  • t≥2t \ge 2:u(2−t)=0u(2 - t) = 0(t=2t = 2 时是 u(0)=0u(0) = 0),f=2f = 2。
-2 -1 1 2 3 4 −1 1 2 0 t
① 2u(t)2u(t):在 0 打开,之后是 2。
-2 -1 1 2 3 4 −1 1 2 0 t
② −u(2−t)-u(2 - t):t<2t \lt 2 时是 −1,在 2 关掉,把左边整体降 1 ✓

端点:t=0t = 0 时 f(0)=2u(0)−u(2)=−1f(0) = 2u(0) - u(2) = -1,取左边的值;t=2t = 2 时 f(2)=2−u(0)=2f(2) = 2 - u(0) = 2,取右边的值。"关掉"型阶跃在跳变点取的是关掉以后的值,和 6·01 ③ 例 1(d) 一样。

-2 -1 1 2 3 4 −1 1 2 0 t f(t)
例 8(a):t = 0 处实心点在 −1,t = 2 处实心点在 2。

(b) g=r(t)−2u(t−1)−r(t−3)g = r(t) - 2u(t - 1) - r(t - 3)。开始时刻 0、1、3。

  • t≤0t \le 0:g=0g = 0。
  • 0<t≤10 \lt t \le 1:只有 r(t)=tr(t) = t,从 0 升到 1。
  • 1<t≤31 \lt t \le 3:t−2t - 2,从 −1 升到 1。
  • t>3t \gt 3:(t−2)−(t−3)=1(t - 2) - (t - 3) = 1。
1 2 3 4 −1 1 2 0 t
① r(t)r(t):从 0 开始斜率 1。
1 2 3 4 −1 1 2 0 t
② −2u(t−1)-2u(t - 1):1 以后整条线往下挪 2,斜率不变。
1 2 3 4 −1 1 2 0 t
③ −r(t−3)-r(t - 3):3 以后斜率归零,停在 1 ✓

在 t=1t = 1:左边的值 1,右边 1−2=−11 - 2 = -1,往下跳 2,但斜率没变(跳变前后都是 1):uu 项只把线整体平移,不改斜率;改斜率要靠 rr 项,这就是 −r(t−3)-r(t - 3) 的作用。少了它,gg 会一直以斜率 1 上升。

1 2 3 4 −1 1 2 0 t g(t)
例 8(b):t = 1 处实心点在 1(取左边的值),空心圈在 −1。

⑤ 函数乘阶跃:截断和平移

不是所有波形都由直线段组成。曲线段用"公式 × 窗口"来写:g(t) [u(t−a)−u(t−b)]g(t)\,[u(t - a) - u(t - b)] 在 a<t≤ba \lt t \le b 上等于 g(t)g(t)、其余为 0(窗口就是高 1 的脉冲);g(t) u(t−a)g(t)\,u(t - a) 是"从 aa 开始一直等于 g(t)g(t)"。这里有一个必须分清的区别:

截断和平移
  • 截断:g(t) u(t−a)g(t)\,u(t - a)。曲线本身不动,只是把 aa 以前的部分抹成 0。它在 t=a+t = a^+ 的值是 g(a)g(a)。
  • 平移:g(t−a) u(t−a)g(t - a)\,u(t - a)。把整条 g(t)u(t)g(t)u(t) 原样搬到右边 aa 处。它在 t=a+t = a^+ 的值是 g(0)g(0)。

r(t−a)=(t−a) u(t−a)r(t - a) = (t - a)\,u(t - a) 是平移(g(t)=tg(t) = t)。判断方法:看 uu 里面的式子和前面函数里的式子是不是同一个。

例 9

四个"函数乘阶跃"

画出并说明:(a) e−tu(t)e^{-t}u(t);(b) (t−1) u(t−1)(t - 1)\,u(t - 1) 和 t u(t−1)t\,u(t - 1);(c) e−(t−2)u(t−2)e^{-(t-2)}u(t - 2) 和 e−tu(t−2)e^{-t}u(t - 2);(d) sin⁡t [u(t)−u(t−π)]\sin t\,[u(t) - u(t - \pi)]。

(a) t≤0t \le 0 为 0;t>0t \gt 0 为 e−te^{-t}。在 0 处从 0 跳到 e0=1e^0 = 1,然后衰减。这是 capacitor 放电的形状(6·16)。

(b) (t−1)u(t−1)=r(t−1)(t - 1)u(t - 1) = r(t - 1):t=1t = 1 以前为 0,之后从 0 开始以斜率 1 上升,没有跳变。t u(t−1)t\,u(t - 1) 是截断:t≤1t \le 1 为 0,t>1t \gt 1 为 tt,在 t=1+t = 1^+ 已经是 1,所以多了一个高 1 的跳变。代数上 t u(t−1)=(t−1+1) u(t−1)=r(t−1)+u(t−1)t\,u(t - 1) = (t - 1 + 1)\,u(t - 1) = r(t - 1) + u(t - 1)。取 t=2t = 2:r(1)=1r(1) = 1,而 2 u(1)=22\,u(1) = 2,确实差 1。

(c) e−(t−2)u(t−2)e^{-(t-2)}u(t - 2) 是 (a) 平移 2:在 t=2t = 2 从 0 跳到 1 再衰减。e−tu(t−2)e^{-t}u(t - 2) 是截断:在 t=2t = 2 跳到 e−2≈0.135e^{-2} \approx 0.135。两者只差一个常数倍:e−tu(t−2)=e−2⋅e−(t−2)u(t−2)e^{-t}u(t - 2) = e^{-2} \cdot e^{-(t-2)}u(t - 2)。

(d) 窗口 u(t)−u(t−π)u(t) - u(t - \pi) 只留 0<t≤π0 \lt t \le \pi,在这一段等于 sin⁡t\sin t:一个半拱。两端 sin⁡0=sin⁡π=0\sin 0 = \sin\pi = 0,没有跳变。

-1 1 2 3 4 1 0 t e−tu(t)
(a) e−tu(t)e^{-t}u(t):t=0t = 0 处实心点在 0,空心圈在 1。
π/2 π 1 0 t f(t)
(d) f(t)=sin⁡t [u(t)−u(t−π)]f(t) = \sin t\,[u(t) - u(t - \pi)]。
1 2 3 1 2 3 0 t t·u(t−1) r(t−1)
(b) 青:r(t−1)r(t - 1)。橙:t u(t−1)t\,u(t - 1),在 t=1t = 1 跳了 1。
1 2 3 4 5 e⁻² 1 0 t e−(t−2)u(t−2) e−tu(t−2)
(c) 青:平移。橙:截断,起点只有 e−2e^{-2}。
最常见的错
  • t u(t−1)≠r(t−1)t\,u(t - 1) \ne r(t - 1):前者在 t=1t = 1 跳了 1。平移的 ramp 两个 tt 都要换。
  • e−tu(t−2)≠e−(t−2)u(t−2)e^{-t}u(t - 2) \ne e^{-(t-2)}u(t - 2):前者小了 e−2e^{-2} 倍。
  • 写完以后在每个开始时刻 aa 代 t=a+t = a^+ 检查一下:平移的 ramp 在那里应该是 0。

⑥ 为什么要拆:linear + time-invariant 课上第 11、12 讲

设一个系统(比如一个 RC 电路)满足两条:linear,即满足 superposition:输入乘常数、相加,输出也跟着乘、跟着加(定义和检验方法见 6·03);time-invariant,同样的输入晚 aa 秒进来,输出形状不变、也只是晚 aa 秒。再假设开始前电路里没有储能(zero state)。把它对 u(t)u(t) 的输出叫 step response,记作 s(t)s(t)。那么

x(t)=∑kΔk u(t−ak)⇓y(t)=∑kΔk s(t−ak)\begin{gathered} x(t) = \sum_k \Delta_k\,u(t - a_k) \\ \Downarrow \\ y(t) = \sum_k \Delta_k\,s(t - a_k) \end{gathered}

理由:time-invariant 给出"u(t−ak)u(t - a_k) 的输出是 s(t−ak)s(t - a_k)",linear 给出"乘系数、相加"。所以只要算出一个 s(t)s(t),阶梯、脉冲输入的输出全是平移和加法。输入是一串 impulse 时同理:δ(t−a)\delta(t - a) 的输出是 impulse response 平移 aa。这两条性质在 C·01 还会细讲,并由此推出 convolution。

例 10

反过来用:已知 step response 和输出,求输入

某个 zero state 的 linear、time-invariant 系统,step response 是 s(t)=r(t)−r(t−1)s(t) = r(t) - r(t - 1):输入一个单位阶跃,输出用 1 s 线性升到 1,然后保持。给它加一个只有跳变的输入(阶梯)x(t)x(t),测得的输出是下图的折线:

y(t)={0,t≤02t,0<t≤12,1<t≤28−3t,2<t≤3−1,3<t≤3.5t−4.5,3.5<t≤4.50,t>4.5y(t) = \begin{cases} 0, & t \le 0 \\ 2t, & 0 \lt t \le 1 \\ 2, & 1 \lt t \le 2 \\ 8 - 3t, & 2 \lt t \le 3 \\ -1, & 3 \lt t \le 3.5 \\ t - 4.5, & 3.5 \lt t \le 4.5 \\ 0, & t \gt 4.5 \end{cases}
1 2 3 3.5 4.5 −1 1 2 0 t y(t)
测得的输出:三段斜坡,每段 1 s,中间是平台。

求输入 x(t)x(t)。

  1. 想清楚要找什么。由上面的结论,输入 x=∑kΔku(t−ak)x = \sum_k \Delta_k u(t - a_k) 的输出是 y=∑kΔks(t−ak)y = \sum_k \Delta_k s(t - a_k)。所以只要把 yy 写成若干个平移、缩放的 ss 之和,平移量就是输入跳变的时刻 aka_k,系数就是跳变量 Δk\Delta_k。
  2. 先用 ① 的方法把 yy 写成 ramp。yy 处处连续,没有跳变;斜率依次是 0、2、0、−3、0、1、0,变化发生在 0、1、2、3、3.5、4.5,Δk\Delta k 依次是 +2,−2,−3,+3,+1,−1+2, -2, -3, +3, +1, -1: y=2r(t)−2r(t−1)−3r(t−2)+3r(t−3)+r(t−3.5)−r(t−4.5)\begin{aligned} y &= 2r(t) - 2r(t - 1) - 3r(t - 2) \\ &\quad + 3r(t - 3) + r(t - 3.5) - r(t - 4.5) \end{aligned}
  3. 两两配成 ss。s(t−a)=r(t−a)−r(t−a−1)s(t - a) = r(t - a) - r(t - a - 1),所以相差 1 s、系数相反的一对 ramp 就是一个平移的 ss: y=2s(t)−3s(t−2)+s(t−3.5)y = 2s(t) - 3s(t - 2) + s(t - 3.5) 下面三张小图把这三个 ss 逐个加上(灰色虚线是测得的 yy):
2 3.5 −1 2 0 t
① 2s(t)2s(t):0 到 1 s 的斜坡,升 2。
2 3.5 −1 2 0 t
② −3s(t−2)-3s(t - 2):2 到 3 s 的斜坡,降 3。
2 3.5 −1 2 0 t
③ +s(t−3.5)+s(t - 3.5):3.5 到 4.5 s 的斜坡,升 1,和 y 重合 ✓
  1. 读出输入。s(t−a)s(t - a) 来自 u(t−a)u(t - a),系数照搬: x(t)=2u(t)−3u(t−2)+u(t−3.5)x(t) = 2u(t) - 3u(t - 2) + u(t - 3.5) 即 0<t≤20 \lt t \le 2 时为 2,2<t≤3.52 \lt t \le 3.5 时为 −1,之后为 0。
  2. 正向核对。把 xx 代回去算两个点:y(2.5)=2s(2.5)−3s(0.5)+s(−1)=2−1.5+0=0.5y(2.5) = 2s(2.5) - 3s(0.5) + s(-1) = 2 - 1.5 + 0 = 0.5,题目给的 8−3×2.5=0.58 - 3 \times 2.5 = 0.5 ✓。y(4)=2−3+s(0.5)=−0.5y(4) = 2 - 3 + s(0.5) = -0.5,题目给的 4−4.5=−0.54 - 4.5 = -0.5 ✓。
1 2 3 3.5 4.5 −1 1 2 0 t x(t) y(t)
橙色:求出的输入。灰色虚线:输出。输入每跳一次 Δ,输出就用 1 s 走完同样的 Δ。

没有做任何积分:只用了"拆成积木、平移、相加",而且是倒着用的。这个系统的 ss 是折线,所以能直接读图;换成 RC 电路的 s(t)=(1−e−t/τ)u(t)s(t) = (1 - e^{-t/\tau})u(t),第 3 步的配对思路完全一样。

练习 1:读图写式子并求导

g(t)g(t):00 到 22 以斜率 2 上升(升到 4);t=2t = 2 跳到 1,22 到 33 平在 1;33 到 44 以斜率 −1 降回 0;之后为 0。(a) 写成 uu 和 rr 的组合。(b) 求 p gp\,g。

提示

t=2t = 2 同时有跳变(Δ=1−4\Delta = 1 - 4)和斜率变化(2→02 \to 0),两块积木都要补,和例 6 的 t=2t = 2 一样。

(a) t=0t = 0:Δk=+2\Delta k = +2,补 2r(t)2r(t)。t=2t = 2:Δ=1−4=−3\Delta = 1 - 4 = -3,补 −3u(t−2)-3u(t - 2);Δk=0−2=−2\Delta k = 0 - 2 = -2,补 −2r(t−2)-2r(t - 2)。t=3t = 3:Δk=−1\Delta k = -1,补 −r(t−3)-r(t - 3)。t=4t = 4:Δk=+1\Delta k = +1,补 r(t−4)r(t - 4)。

g=2r(t)−2r(t−2)−3u(t−2)−r(t−3)+r(t−4)\begin{aligned} g &= 2r(t) - 2r(t - 2) - 3u(t - 2) \\ &\quad - r(t - 3) + r(t - 4) \end{aligned}

检查 2<t≤32 \lt t \le 3:2t−2(t−2)−3=12t - 2(t - 2) - 3 = 1 ✓;3<t≤43 \lt t \le 4:1−(t−3)=4−t1 - (t - 3) = 4 - t ✓;t>4t \gt 4:4−t+(t−4)=04 - t + (t - 4) = 0 ✓。

(b) r→ur \to u,u→δu \to \delta:p g=2u(t)−2u(t−2)−3 δ(t−2)−u(t−3)+u(t−4)p\,g = 2u(t) - 2u(t - 2) - 3\,\delta(t - 2) - u(t - 3) + u(t - 4)。

练习 2:从式子画波形

画出 q(t)=3u(t)−3r(t−1)+3r(t−2)q(t) = 3u(t) - 3r(t - 1) + 3r(t - 2),写成分段形式。

提示

开始时刻 0、1、2,每一段只留已经开始的项。

t≤0t \le 0:0。0<t≤10 \lt t \le 1:3。1<t≤21 \lt t \le 2:3−3(t−1)=6−3t3 - 3(t - 1) = 6 - 3t,从 3 降到 0。t>2t \gt 2:6−3t+3(t−2)=06 - 3t + 3(t - 2) = 0。形状:0 处跳到 3,平一段,再用 1 秒斜降到 0。

练习 3:有曲线段的波形

w(t)w(t):0<t≤10 \lt t \le 1 时为 2;t>1t \gt 1 时为 2e−(t−1)2e^{-(t-1)}(从 2 开始衰减);t≤0t \le 0 时为 0。写成"公式 × 窗口"的形式,并说明 t=1t = 1 处有没有跳变。

提示

第一段用窗口 u(t)−u(t−1)u(t) - u(t - 1),第二段用 u(t−1)u(t - 1)。第二段的函数是不是"平移"的形式?

w=2[u(t)−u(t−1)]+2e−(t−1)u(t−1)w = 2[u(t) - u(t - 1)] + 2e^{-(t-1)}u(t - 1),也可以合并成 2u(t)+2(e−(t−1)−1)u(t−1)2u(t) + 2\big(e^{-(t-1)} - 1\big)u(t - 1)。t=1+t = 1^+ 时第二段是 2e0=22e^0 = 2,和左边的 2 相等,所以没有跳变(第二段正是平移形式,起点等于 2e02e^0)。