下面这些性质让很多 convolution 不用真的积分就能写出来,也能用来检查答案。Alexander 列了七条(编号沿用他的);后面再补一条平移规则和两条检查规则(起止时刻、面积)。每一条都说明为什么成立,并给一个最小的例子。
① 代数性质:交换律、分配律、结合律
- 性质 1,交换律:x∗h=h∗x。证明见C·03。
- 性质 2,分配律 (distributive):f∗(x+y)=f∗x+f∗y。积分对被积函数是线性的,所以和可以拆开积;常数也可以提出来:(kx)∗h=k(x∗h)。用途:把分段信号拆成阶跃或矩形之和,例如脉冲 =u(t)−u(t−2)(C·08例 2 解法 A)。
- 性质 3,结合律 (associative):f∗(x∗y)=(f∗x)∗y。用途:两个系统串联,输入 x 先过 h1 得 w=x∗h1,再过 h2 得 y=w∗h2=x∗(h1∗h2):总的 impulse response 是 h1∗h2。
② 和 δ 卷积:原样复制、平移
- 性质 4:f∗δ=f。写出积分 ∫f(λ)δ(t−λ)dλ:作为 λ 的函数,δ(t−λ) 是位于 λ=t、面积 1 的 impulse;由 sampling(sifting)property(6·01),积分把 f 在 λ=t 的值挑出来,得 f(t)。
- 性质 5:f∗δ(t−t0)=f(t−t0)。同理,δ(t−t0−λ) 位于 λ=t−t0,挑出 f(t−t0)。带面积:f∗Aδ(t−t0)=Af(t−t0)。和 impulse 卷积 = 平移并缩放;这正是C·01例 1:3δ(t−2)∗te−tu(t)=3(t−2)e−(t−2)u(t−2)。
- 性质 6:f∗δ′=f′。
f 是 0 到 1 的三角形。f∗δ(t−2.5):形状不变,整个往右挪 2.5。
有 impulse 时结果可以不连续
C·06边界表的"结果连续"只在两个函数都不含 impulse 时成立。比如 δ∗u=u 在 0 处跳变,f∗δ(t−t0) 会把 f 的跳变原样搬过去。和 δ 卷积时直接用性质 4、5,不要画边界表。
③ 和 u 卷积:running integral
- 性质 7:f∗u=∫−∞tf(λ)dλ,即 C·07的 running integral(Davis:f⊗u=p1f)。小检查:δ∗u=∫−∞tδ=u;[u(t)−u(t−1)]∗u=r(t)−r(t−1)(脉冲从头积:先升 1 秒,然后平在 1;r=tu(t) 是斜坡)。
④ 平移规则
平移规则
如果 x∗h=y,那么
x(t−t1)∗h(t−t2)=y(t−t1−t2)两个都晚,结果晚"两个延迟之和"。
为什么:x(t−t1)=x∗δ(t−t1),h(t−t2)=h∗δ(t−t2)(性质 5)。用交换律、结合律重新组合:x∗δ(t−t1)∗h∗δ(t−t2)=(x∗h)∗[δ(t−t1)∗δ(t−t2)]。最后一项用性质 5(f=δ(t−t1))是 δ(t−t1−t2),再用一次性质 5 就得 y(t−t1−t2)。用途:先用从 0 开始的信号算,最后整体平移,比直接对平移过的信号积分省事。
⑤ 两条检查规则:起止时刻和面积
起止时刻和面积
(a) x 只在 (a,b) 上不为 0,h 只在 (c,d) 上不为 0,则 y=x∗h 只在 (a+c, b+d) 上不为 0:起点相加,终点相加,宽度相加。
(b) 两个面积都有限时:
∫−∞∞y(t)dt=(∫−∞∞xdt)(∫−∞∞hdt)
(a) 的理由在C·06:第一次重叠在 t−c=a,最后一次在 t−d=b。(b) 的理由:交换积分顺序,
∫ydt=∫∫x(λ)h(t−λ)dλdt=∫x(λ)[∫h(t−λ)dt]dλ=∫x(λ)Ahdλ=AxAh
中括号里是把 h 平移 λ 以后的面积,平移不改变面积,所以等于 h 的面积 Ah。例:第 04 节例 1,Ax=2、Ah=2,梯形面积 4 ✓。注意 u 的面积是无穷大,和 u 卷积时这条用不了。
⑥ Davis 的求导性质
Davis 用 p 表示 d/dt。(13.8-33) 说 convolution 的导数可以随意分给两个因子(j 取 0,1,…,k 中任意一个):
pk{x⊗y}=pjx⊗pk−jy
理由(13.8-34):两边的 Laplace 变换都是 skX(s)Y(s)(Davis 约定 L{px}=sX(s),没有初值项)。小检查:k=1, j=1,x=y=u。左边 p{tu(t)}:按 Davis 的规则 p[fu]=f′u+f(0)δ,f=t 时 f(0)=0,得 u(t)。右边 pu⊗u=δ⊗u=u。相等 ✓。Davis 用它把含 impulse 的 convolution 变成普通函数的 convolution 再求导;
⑦ 例题
例题例 1 只用性质,一行一题
r(t)=tu(t) 是斜坡。不积分,求:(a) [δ(t)−δ(t−2)]∗r(t);(b) u(t+1)∗u(t−3);(c) 反过来:找 x(t),使 x∗u(t)=u(t−1)−u(t−3);(d) [δ(t)+δ(t−1)]∗[δ(t)−δ(t−1)]∗u(t)。
看解答
(a) 分配律 + 性质 5:一项是 r 原样,一项是 r 晚 2:
[δ−δ(t−2)]∗r=r(t)−r(t−2)
分段:t<0 为 0;0<t<2 为 t;t>2 为 t−(t−2)=2。斜坡被"削平"在 2。
(b) 平移规则:u(t+1) 是 u 提前 1(t1=−1),u(t−3) 是 u 晚 3(t2=3)。u∗u=tu(t),整体平移 t1+t2=2:
u(t+1)∗u(t−3)=(t−2)u(t−2)
延迟可以是负的(提前),规则照样成立。
(c) 性质 7:x∗u 是 x 的 running integral,所以 x 是右边的导数。u(t−1)−u(t−3) 在 1 处跳上 1、在 3 处跳下 1,导数是两个 impulse:
x(t)=δ(t−1)−δ(t−3)
检查:从头积 δ(t−1)−δ(t−3),过 1 时加 1、过 3 时减 1,正好是 1 到 3 的矩形 ✓。
(d) 先把两组 impulse 卷开(分配律 + 性质 5:δ(t−a)∗δ(t−b)=δ(t−a−b))。和 (1+z)(1−z)=1−z2 一样,中间两项相消:
[δ+δ(t−1)]∗[δ−δ(t−1)]=δ−δ(t−1)+δ(t−1)−δ(t−2)=δ(t)−δ(t−2)
再和 u 卷(结合律,先算哪两个都行):u(t)−u(t−2),0 到 2 的矩形。
例题例 2 中间有缺口的 h:算一块,另一块是平移
x(t)=2(0<t<2);h(t)=1(0<t<1 和 2<t<3),中间 1 到 2 是缺口,其他为 0。求 y=x∗h。
x:高 2,0 到 2。h:0 到 1 和 2 到 3 高 1,中间是缺口。
拆 h。记 Rab 为在 (a,b) 上高 1 的矩形:h=R01+R23,而 R23 就是 R01 晚 2。
只算一块。x∗R01:宽 2 高 2 和宽 1 高 1 的矩形,用C·06的通式:0<t<1 为 2t,1<t<2 为 2,2<t<3 为 2(3−t)。
另一块 = 平移。平移规则:x∗R23 = x∗R01 晚 2:2<t<3 为 2(t−2),3<t<4 为 2,4<t<5 为 2(5−t)。
相加。只有 2<t<3 两块重叠:2(3−t)+2(t−2)=2。所以 y=2t(0–1),2(1–4),10−2t(4–5),其他为 0。
灰虚线:x∗R01 和它晚 2 的拷贝;青色:两者之和。中间平的一段正好填住了缺口。
为什么缺口"不见了":看翻转的 x。x(t−λ) 是宽 2 的窗口 [t−2, t],y = 2 × 窗口里 h 的面积:
① 0<t<1:窗口只盖住 R01 的 [0, t] → y=2t② 1<t<4:窗口宽 2,总盖住 1 个单位的 h(图里 0.5 + 0.5)→ y=2③ 4<t<5:只剩 R23 的 [t−2, 3] → y=2(5−t)
检查
起止:0+0=0 到 2+3=5 ✓。连续:t=1,4 两边都是 2,t=5 为 0 ✓。面积:Ax=4,Ah=2,y 面积 1+6+1=8 ✓。也可以拆成 h=R03−R12(大矩形减小矩形),结果相同。
如果输入只有 1 宽(高还是 2,比缺口窄),同样的做法给出 2t、2(2−t)、2(t−2)、2(4−t):在 t=2 掉到 0。窄的输入"看得见"缺口,宽的看不见。
| # | 性质 | 一句话 |
| 1 | x∗h=h∗x | 翻哪个都行 |
| 2 | f∗(x+y)=f∗x+f∗y | 拆开分别卷,常数提出来 |
| 3 | f∗(x∗y)=(f∗x)∗y | 串联系统:h1∗h2 |
| 4 | f∗δ=f | 原样复制 |
| 5 | f∗Aδ(t−t0)=Af(t−t0) | 平移并缩放 |
| 6 | f∗δ′=f′ | 求导 |
| 7 | f∗u=∫−∞tf | running integral |
| 平移 | x(t−t1)∗h(t−t2)=y(t−t1−t2) | 延迟相加 |
| 检查 | 起点、终点相加;Ay=AxAh | 无 impulse 时还要连续 |