跳转到内容
新建笔记

信号与系统 01:数学、信号、系统与卷积

完整笔记入口 · 下一章:拉普拉斯与 RC

1. 函数不是“一个固定数”

跳转到“1. 函数不是“一个固定数””

x(t)x(t) 表示“给定时刻 tt,得到那个时刻的数值”。若电压满足 x(t)=2tx(t)=2t V,公式里的系数实际量纲为 V/s:

x(0)=0 V,x(0.5)=1 V,x(1)=2 V.x(0)=0\ \mathrm V,\qquad x(0.5)=1\ \mathrm V,\qquad x(1)=2\ \mathrm V.

x(t)x(t) 是整个变化规则,x(1)x(1) 是这个规则在 1 s 处的取值。离散信号 x[n]x[n] 也是规则,但 nn 只能取整数。均匀采样时:

tn=nTs,Ts=1fs,x[n]=xc(nTs).t_n=nT_s,\qquad T_s=\frac1{f_s},\qquad x[n]=x_{\mathrm c}(nT_s).

例如 fs=1000f_s=1000 Hz,Ts=0.001T_s=0.001 s,第 3 个编号的采样时刻是 t3=0.003t_3=0.003 s。这里从编号 0 开始,因此“编号 3”是第 4 次记录。

连续时间不等于幅值必须连续。阶跃在连续时间轴上定义,但在启动时刻跳变。离散时间不等于已经数字量化:理论上的 x[n]x[n] 仍可取任意实数,ADC 才会进一步将幅值编码为有限等级。

2. 导数是变化速度,积分是累积

跳转到“2. 导数是变化速度,积分是累积”

某段时间的平均速度:

ΔyΔt=y(t+Δt)−y(t)Δt.\frac{\Delta y}{\Delta t} =\frac{y(t+\Delta t)-y(t)}{\Delta t}.

时间间隔越来越小时,若极限存在:

y′(t)=dydt=lim⁡Δt→0y(t+Δt)−y(t)Δt.y'(t)=\frac{dy}{dt} =\lim_{\Delta t\to0}\frac{y(t+\Delta t)-y(t)}{\Delta t}.

若 yy 是电压,导数单位是 V/s;若 yy 是位移,导数是速度,单位 m/s。dy/dtdy/dt 在这里是一个整体记号。

电流定义为电荷的变化速度:

i(t)=dq(t)dt.i(t)=\frac{dq(t)}{dt}.

把时间分成很多小段,第 ℓ\ell 段电荷增量约为 i(tℓ)Δti(t_\ell)\Delta t。将所有小段相加,再取极限:

q(t)−q(0)=lim⁡Δt→0∑ℓi(tℓ)Δt=∫0ti(ξ) dξ.q(t)-q(0)=\lim_{\Delta t\to0}\sum_\ell i(t_\ell)\Delta t =\int_0^t i(\xi)\,d\xi.

ξ\xi 是积分内的临时变量,不是新物理量。积分的单位为 A·s=C。若 i=0.001i=0.001 A 持续 0.1 s,增加电荷 0.00010.0001 C。

两条后面会用到的公式

跳转到“两条后面会用到的公式”

乘积求导:

(uv)′=u′v+uv′.(uv)'=u'v+uv'.

对两边积分,把最后一项移过去:

∫uv′ dt=uv−∫u′v dt.\int u v'\,dt=uv-\int u'v\,dt.

这就是分部积分,它解释拉普拉斯中的导数规则。对 eate^{at}:

ddteat=aeat,∫eat dt=eata+常数,a≠0.\frac{d}{dt}e^{at}=ae^{at},\qquad \int e^{at}\,dt=\frac{e^{at}}a+\text{常数},\quad a\ne0.

3. 正弦中的幅值、频率与相位

跳转到“3. 正弦中的幅值、频率与相位”

统一采用余弦表示:

x(t)=Acos⁡(ωt+ϕ),ω=2πf,T=1f.x(t)=A\cos(\omega t+\phi),\qquad \omega=2\pi f,\qquad T=\frac1f.

AA 是峰值,ff 是每秒周期数,ϕ\phi 是在 t=0t=0 的相位角。弧度与角度换算为 2π2\pi rad=360∘360^\circ。

延迟 td>0t_d>0:

y(t)=x(t−td)=Acos⁡(ω(t−td)+ϕ)=Acos⁡(ωt+ϕ−ωtd).\begin{aligned} y(t)&=x(t-t_d)\\ &=A\cos\bigl(\omega(t-t_d)+\phi\bigr)\\ &=A\cos(\omega t+\phi-\omega t_d). \end{aligned}

因此“输出相对输入”的相位差为 −ωtd-\omega t_d。延迟四分之一周期得到 −90∘-90^\circ;相位按 360∘360^\circ 周期等价,测量时还要考虑相位展开。

离散正弦:

x[n]=Acos⁡(Ωn+ϕ),Ω=ωTs=2πffs.x[n]=A\cos(\Omega n+\phi),\qquad \Omega=\omega T_s=\frac{2\pi f}{f_s}.

ω\omega 与 Ω\Omega 单位不同,不能直接把两者数字互换。

4. 复数为何能表示振荡

跳转到“4. 复数为何能表示振荡”

j2=−1j^2=-1,复数 a+jba+jb 可以看作平面上的横坐标 aa、纵坐标 bb。模与角度:

∣a+jb∣=a2+b2,arg⁡(a+jb)=atan2⁡(b,a).|a+jb|=\sqrt{a^2+b^2},\qquad \arg(a+jb)=\operatorname{atan2}(b,a).

用指数函数、正弦、余弦的幂级数,将奇偶幂分开:

ejθ=1+jθ+(jθ)22!+(jθ)33!+⋯=(1−θ22!+θ44!−⋯ )+j(θ−θ33!+θ55!−⋯ )=cos⁡θ+jsin⁡θ.\begin{aligned} e^{j\theta} &=1+j\theta+\frac{(j\theta)^2}{2!}+\frac{(j\theta)^3}{3!}+\cdots\\ &=\left(1-\frac{\theta^2}{2!}+\frac{\theta^4}{4!}-\cdots\right) +j\left(\theta-\frac{\theta^3}{3!}+\frac{\theta^5}{5!}-\cdots\right)\\ &=\cos\theta+j\sin\theta. \end{aligned}

这给出欧拉公式,以及:

cos⁡θ=ejθ+e−jθ2,sin⁡θ=ejθ−e−jθ2j.\cos\theta=\frac{e^{j\theta}+e^{-j\theta}}2,\qquad \sin\theta=\frac{e^{j\theta}-e^{-j\theta}}{2j}.

复指数同时组织两路正交振荡,真实电压可取实部。它是计算工具,真实测量不要求产生“虚数电压”。

∑n=02x[n]\sum_{n=0}^{2}x[n] 是 x[0]+x[1]+x[2]x[0]+x[1]+x[2] 的简写。令:

SN=1+q+q2+⋯+qN−1.S_N=1+q+q^2+\cdots+q^{N-1}.

乘以 qq 并相减:

qSN=q+q2+⋯+qN,qS_N=q+q^2+\cdots+q^N, SN−qSN=1−qN,SN=1−qN1−q,q≠1.S_N-qS_N=1-q^N,\qquad S_N=\frac{1-q^N}{1-q},\quad q\ne1.

若 ∣q∣<1|q|<1,N→∞N\to\infty 时 qN→0q^N\to0,于是:

∑n=0∞qn=11−q.\sum_{n=0}^{\infty}q^n=\frac1{1-q}.

q=1q=1 时有限和为 NN,不能代入含 1−q1-q 的分式。无限和的收敛条件不能删除;后面 Z 变换的收敛域就从这里来。

离散阶跃与离散冲激:

u[n]={1,n≥00,n<0,δ[n]={1,n=00,n≠0.u[n]= \begin{cases}1,&n\ge0\\0,&n<0,\end{cases} \qquad \delta[n]= \begin{cases}1,&n=0\\0,&n\ne0.\end{cases}

因此 δ[n]=u[n]−u[n−1]\delta[n]=u[n]-u[n-1]。

连续阶跃 u(t)u(t) 在 t<0t<0 为 0,在 t>0t>0 为 1;单点取值通常不改变普通积分。连续冲激 δ(t)\delta(t) 是广义函数,由单位面积和抽样性质定义:

∫−∞∞δ(t) dt=1,∫−∞∞f(t)δ(t−t0) dt=f(t0).\int_{-\infty}^{\infty}\delta(t)\,dt=1,\qquad \int_{-\infty}^{\infty}f(t)\delta(t-t_0)\,dt=f(t_0).

可以用宽 ε\varepsilon、高 1/ε1/\varepsilon 的矩形脉冲理解其极限,但不能把 δ(t)\delta(t) 当作“在 0 处取一个普通有限值的函数”。连续单位冲激量纲为每秒,离散 δ[n]\delta[n] 无量纲。

7. 信号能量、平均功率与 RMS

跳转到“7. 信号能量、平均功率与 RMS”

对电压信号,数学信号能量定义为:

E=∫−∞∞∣x(t)∣2 dt或Ed=∑n=−∞∞∣x[n]∣2.E=\int_{-\infty}^{\infty}|x(t)|^2\,dt \quad\text{或}\quad E_{\mathrm d}=\sum_{n=-\infty}^{\infty}|x[n]|^2.

前者单位 V²·s,后者按上述定义是 V²;离散和若近似连续积分,还须乘 TsT_s。这里的“信号能量”不自动等于焦耳,电阻消耗能量还需除以电阻。

平均功率型量:

P=lim⁡T→∞12T∫−TT∣x(t)∣2 dt.P=\lim_{T\to\infty}\frac1{2T}\int_{-T}^{T}|x(t)|^2\,dt.

有限记录的均方根:

xRMS=1N∑n=0N−1x[n]2.x_{\mathrm{RMS}}=\sqrt{\frac1N\sum_{n=0}^{N-1}x[n]^2}.

对零均值正弦,利用 cos⁡2θ=(1+cos⁡2θ)/2\cos^2\theta=(1+\cos2\theta)/2,在完整周期上第二项平均为零:

xRMS=A22=∣A∣2.x_{\mathrm{RMS}}=\sqrt{\frac{A^2}{2}}=\frac{|A|}{\sqrt2}.

若 x=d+Acos⁡θx=d+A\cos\theta,平方并取周期平均,交叉项为零:

xRMS2=d2+A22.x_{\mathrm{RMS}}^2=d^2+\frac{A^2}{2}.

所以含 DC 的总 RMS 与去均值后的交流 RMS 不同。

系统 T\mathcal T 是输入到输出的规则,y=T{x}y=\mathcal T\{x\}。

线性:同时满足叠加与缩放,即:

T{ax1+bx2}=aT{x1}+bT{x2}.\mathcal T\{a x_1+b x_2\}=a\mathcal T\{x_1\}+b\mathcal T\{x_2\}.

两点平均器满足此式;y=x2y=x^2 不满足,因为 (x1+x2)2(x_1+x_2)^2 多出 2x1x22x_1x_2;y=x+1y=x+1 也不线性,因为零输入输出非零。

时不变:若 x[n]→y[n]x[n]\to y[n],把输入改成 x[n−n0]x[n-n_0] 后,输出必须是 y[n−n0]y[n-n_0]。例如 y[n]=nx[n]y[n]=n x[n]:

T{x[n−n0]}=nx[n−n0],\mathcal T\{x[n-n_0]\}=n x[n-n_0],

而 y[n−n0]=(n−n0)x[n−n0]y[n-n_0]=(n-n_0)x[n-n_0],二者不同,所以时变。

因果:当前输出只依赖当前及过去输入。x[n−1]x[n-1] 因果,x[n+1]x[n+1] 需要未来输入;离线处理整段记录可以使用非因果方法。

无记忆:当前输出只依赖当前输入。2x[n]2x[n] 无记忆,两点平均器有记忆。

BIBO 稳定:任意有界输入,都产生有界输出。积分器输入常数 1 会得到不断增长的 tt,所以不满足此定义。

线性微分方程在固定非零初始状态下,完整输入输出映射通常是仿射映射。用于卷积和传递函数的线性映射是零状态响应;初始状态引起的响应单独相加。

9. 从冲激分解推导离散卷积

跳转到“9. 从冲激分解推导离散卷积”

任意序列都可写成:

x[n]=∑k=−∞∞x[k]δ[n−k].x[n]=\sum_{k=-\infty}^{\infty}x[k]\delta[n-k].

固定 nn 时,只有 k=nk=n 那项冲激为 1,其余为 0,因此右侧确实等于 x[n]x[n]。

定义零状态系统对 δ[n]\delta[n] 的响应为 h[n]h[n]。时不变性说明 δ[n−k]\delta[n-k] 的响应是 h[n−k]h[n-k];线性说明乘 x[k]x[k] 后输出乘同一系数,再把所有输出相加:

y[n]=T{∑kx[k]δ[n−k]}=∑kx[k]T{δ[n−k]}=∑kx[k]h[n−k].\begin{aligned} y[n] &=\mathcal T\left\{\sum_k x[k]\delta[n-k]\right\}\\ &=\sum_k x[k]\mathcal T\{\delta[n-k]\}\\ &=\sum_k x[k]h[n-k]. \end{aligned}

这就是卷积:

y[n]=(x∗h)[n]=∑kx[k]h[n−k].\boxed{y[n]=(x*h)[n]=\sum_k x[k]h[n-k]}.

无限展开需要相应的收敛与连续性条件;本章有限序列例子可以直接交换求和。这个分解逻辑与 MIT 卷积讲次 一致。

规则 y[n]=(x[n]+x[n−1])/2y[n]=(x[n]+x[n-1])/2。输入 δ[n]\delta[n] 得:

h[n]=12δ[n]+12δ[n−1].h[n]=\frac12\delta[n]+\frac12\delta[n-1].

输入 [1,2,3][1,2,3],区间外为零:

nn当前输入上一点当前输出
0100.5
1211.5
2322.5
3031.5

冲激响应只有两项非零,卷积因此也化成这条递推规则。

10. 连续卷积与零状态/零输入

跳转到“10. 连续卷积与零状态/零输入”

连续信号的冲激分解对应积分:

x(t)=∫−∞∞x(ξ)δ(t−ξ) dξ.x(t)=\int_{-\infty}^{\infty}x(\xi)\delta(t-\xi)\,d\xi.

同样利用线性与时不变性:

yzs(t)=∫−∞∞x(ξ)h(t−ξ) dξ.\boxed{y_{\mathrm{zs}}(t)=\int_{-\infty}^{\infty}x(\xi)h(t-\xi)\,d\xi}.

因果系统、输入从 0 开始时,积分缩为 0≤ξ≤t0\le\xi\le t。完整响应:

y(t)=yzs(t)+yzi(t).y(t)=y_{\mathrm{zs}}(t)+y_{\mathrm{zi}}(t).

“zs”表示零状态,由外部输入驱动;“zi”表示零输入,由原有电容电压等初始状态驱动。下一章会直接计算两部分。

对于普通函数型冲激响应:

∣y(t)∣≤sup⁡ξ∣x(ξ)∣∫−∞∞∣h(v)∣ dv.|y(t)|\le \sup_\xi|x(\xi)|\int_{-\infty}^{\infty}|h(v)|\,dv.

若积分有限,就能用一个固定常数界定所有有界输入的输出。离散情形把积分换成 ∑∣h[n]∣\sum|h[n]|。直接通路中的有限冲激项也可稳定,需要用广义函数对应的有限测度来理解,不能要求 δ(t)\delta(t) 是普通可积函数。

  1. 2 V 峰值、无 DC 的正弦 RMS?答案 2/2≈1.4142/\sqrt2\approx1.414 V。
  2. [1,2,3][1,2,3] 先放大 2 再延迟 1 点?答案 [0,2,4,6][0,2,4,6]。
  3. y[n]=x[n]+x[n−2]y[n]=x[n]+x[n-2] 是否线性、时不变、因果、有记忆?四项依次为是、是、是、是。
  4. 输入 [1,1][1,1]、冲激响应 [1,2][1,2] 的卷积?逐点为 [1,3,2][1,3,2]。
  5. 100 Hz 正弦延迟 1 ms 的相位?−2π100×0.001=−0.2π-2\pi100\times0.001=-0.2\pi rad,即 −36∘-36^\circ。

复习时先说出物理含义,再写公式;对所有涉及时间、频率和电压的结果检查单位。