【信号与系统】01 - 线性时不变系统

1. 信号与系统

1.1 序言

 　　实话说，我之前对信号系统一窍不通，但因为工作需要不得不学习它的基础知识。快餐式的学习又不是我的性格，于是就上手了奥本海姆的两本大部头经典，经过一番死磕，总算是把基本内容梳理完了。由于完全没有实战经验，我对信号系统的理解还很浅薄，加上所需数学知识还没有准备完善，学习过程中总是心里没底。但整理过程中，我仍然试图去了解问题的本质，即便有很多猜测和自以为的成分，却已是尽量做到逻辑完整、自洽。

　　信号学的内容比较多，以它为基础的延伸学科已然在各领域开枝散叶，比如语音信号处理和图形信号处理，都是如今AI某些分支的理论基础。万丈高楼平地起，本系列博客只打算学习经典信号学的基本概念和结论，更高阶的信号处理方法将放到《现代信号处理》专题里讨论。鉴于本人理论先行的习惯，打算还是先对连续信号建立较完整的理论，后面再集中讨论离散信号的不同、以及具体的处理策略，与奥本海姆的并行策略不同。

　　信号系统虽是一门应用学科，但就其理论体系而言，却是独立于具体的技术场景之外的。所以这里我也将不拿具体场景为例，而是把它当成一门纯粹的应用数学和算法设计课程来看待，注意力仍然集中在对讨论对象的定义、性质方法的推导、直至计算机的算法实现。教材虽厚，但大部分都是例题和习题，其实真正的理论部分内容有限。建议先粗略了解梗概，再细致推演细节，有数学基础的朋友甚至可以尝试自己“顺水推舟”，可能比直接看书更快、更深刻。

1.2 信号

　　课程名称就已经把它研究的对象概况清楚了：“信号”和“系统”。但我们肯定没法讨论任意的信号和系统，而只能针对某些广泛存在的、具有明确特征的信号与系统。具体而言，这里讨论的信号都可以抽象为一个实数域上的函数\(x(t)\)。由于大部分信号都体现出时序性，故这里的自变量也就用代表时间的\(t\)表示，但在讨论时还是仅作为一般的实（复）函数看待。

　　既然是一般实数域函数，那些我们已经熟知的函数性质和变换，这里就不赘述了。只是在信号系统中，它们有着更应景的名称，比如式（1）中的信号分别称为对信号\(x(t)\)的时移、时间反转、时间尺度变换。而更一般地，信号\(x(\alpha t+\beta)\)是对原信号时间轴的线性变换，想必你已经清楚它跟原信号之间的关系。

\[x(t-t_0);\;\;x(-t);\;\;x(kt)\tag{1}\]

　　还有类似偶信号、奇信号、周期信号的定义也都是简单直观的，其中周期信号满足\(x(t)=x(t+T)\)，它时移\(T\)后能与原信号重合。不难看出\(kT\)（\(k\)为整数）都是它的周期，而函数最小的周期被称为基波周期。但要注意，并不是所有的函数都有最小周期，比如常数函数、比如狄利克雷函数（有理数值为1、无理数处值为0）。教材里说除了常数函数都有基波周期，个人认为有误。

1.3 系统

　　“系统”被视为一个过程，它以一些信号作为输入，并以一定特征输出另一些信号。这样的描述和函数的定义并无本质差别，所以一个信号系统其实是一个广义的函数，只不过它的自变量和因变量都是实变量函数。相当一部分系统其实只有一个输入输出，故今后提到的系统都是指单输入输出系统\(x(t)\to y(t)\)。另外，简单的小系统可以通过综合运算，组合成更大的系统，这样的分拆思想是研究复杂系统的关键方法。

　　方框图可以为复杂系统提供直观的表示，图中的子系统一般用方框表示，像常用的积分、微分一般直接用\(\int,D\)表示。更简单的乘常数系统，可以直接在箭头上标出系数\(a\)。加法器\(\bigoplus\)是方框图中的常用部件，它一般用作多个系统的并联。相应地，复合系统\(S2(S1(x(t)))\)称为系统的串联或级联。右图中由两个子系统和一个加法器组成的系统也很常见（D系统的输出一般为负因子），它被称为反馈互联。

　　确定系统特性，是分析或设计系统之前必须做的事情。而且针对普遍特性的研究，才更容易有理论结果和实用价值，这里先介绍几个常见的系统特性。首先对一般的有界函数\(x(t)\)，总希望系统的输出\(y(t)\)也是有界的，这样的系统称为稳定的。以时间为自变量的系统中，往往\(y(t_0)\)的值与\(x(t)\)的\(t_0\)时刻之后的值无关，这样的系统叫因果系统。如果\(y(t_0)\)仅仅与\(x(t_0)\)有关，它也叫无记忆系统。以上几个系统性质比较符合经验特性，但一定程度上限制了系统的处理范围，对更广义的系统不一定成立，故今后不作为默认性质。真正对系统分析有用的，其实是下面要介绍的两个性质。

　　对输入信号最简单的变换，无非就是时移和叠加，很多重要系统对这些变换仍然有“不变”的输出。首先我们希望系统的输出不受输入信号时移的影响，即如果\(x(t)\to y(t)\)，那么\(x(t-t_0)\to y(t-t_0)\)。输出函数随输入函数一起时移，这样的系统称为时不变的。线性组合是信号的最简单复合，如果系统还能满足式（2），它就被称为线性系统。线性时不变系统普遍存在，而且它还有着非常完善的理论成果，因此本系列把线性时不变系统（LIT）作为默认系统，并透彻讨论它的特性。但要注意：在对陌生系统做分析时，一定要先确认它是LIT系统。

\[x_1(t)\to y_1(t);\;x_2(t)\to y_2(t)\;\;\Rightarrow\;\;ax_1(t)+bx_2(t)\to ay_1(t)+by_2(t)\tag{2}\]

2. LIT系统

2.1 单位冲激函数

　　为了得到LIT系统的一般表征，还是要结合“时不变”和“线性”两大特性，对任意输入\(x(t)\)作一个线性拆解。设\(\Delta(t-t_0)\)是仅在\(t_0\)处有非零值\(1\)的函数，则\(x(t)\)就是所有函数\(x(t_0)\Delta(t-t_0)\)的叠加（\(t_0\)取遍所有实数）。进而如果\(\Delta(t)\)的响应函数为\(Y(t)\)，那么\(x(t)\)的响应函数\(y(t)\)便是所有函数\(x(t_0)Y(t-t_0)\)的叠加。

　　由于\(t_0\)取遍所有实数，“累加”的确切意义便是积分，以下先把\(t_0\)视为连续变量并重新记为\(\tau\)。现在需要构造一个等价于\(\Delta(t-\tau)\)的微分\(\delta(t-\tau)\,\text{d}\tau\)，也就是说微分仅在\(\tau=t\)处有非零值\(1\)，那就要求\(\delta(t-\tau)\)在\(t\)处有无穷大值且积分（面积）为\(1\)。为此需要定义一个满足式（3）的特殊的函数\(\delta(t)\)，它被称为单位冲激函数，坐标系中往往以一个单位长度的向上向量表示（左图）。在此定义下，\(x(t)\)便有了如式（4）积分表达式。

\[\delta(t)=\left\{\begin{matrix}\infty,&(t=0)\\0,&(t\ne 0)\end{matrix}\right.,\;\;\int_{-\infty}^{\infty}\delta(t)\,\text{d}t=1\tag{3}\]

\[x(t)=\int_{-\infty}^{\infty}x(\tau)\delta(t-\tau)\,\text{d}\tau\tag{4}\]

　　教材上一般会用函数列的极限来解释\(\delta\)函数，但我们心里明白，式（3）的定义有严重的逻辑缺陷，一个值为无穷大的函数根本不能算作函数。这样的问题在工科教材上是无法解释清楚的，因为透彻地理解\(\delta(t)\)还需要《实变函数》和《泛函分析》的知识。在那里它叫狄拉克\(\delta\)函数，它是一个线性泛函，具体说是满足\(\delta(\text{d}t)=\text{d}H(t)\)的一个测度，其中\(H(t)\)是右图的Heviside阶梯函数。

　　\(\delta(t)\)的严格定义是满足式（4）的泛函，注意这里的函数已经不是一般意义的函数，甚至积分（勒贝格积分）也不是一般意义的积分（黎曼积分）。另外，还可以把\(\delta(t)\)看成一个广义的分布函数，这其实是对本课程非常适合的角度。相比一般函数，它们被称为奇异函数，数学上已经严格证明，奇异函数和常规函数之间可以正常进行四则运算以及微积分，而且计算的结果与直观的理解是一致的，以后可以放心使用。

2.2 LIT的分解

　　有了式（4）对输入信号\(x(t)\)分解，就可以继续讨论系统的输出。在微分意义下，把响应\(Y(t-\tau)\)写作\(h(t-\tau)\,\text{d}\tau\)，那么响应函数\(y(t)\)就是式（5），它就是LIT系统的表征函数。虽然一般把\(h(t)\)称为单位冲激响应，但形式\(h(t-\tau)\,\text{d}\tau\)才是真正意义上的响应函数。就好比把\(\delta(t-\tau)\,\text{d}\tau\)视为真正的输入信号，而\(\delta(t-\tau),\,h(t-\tau)\)只是响应系数而已。理清这个关系，在今后的推导中才不会困惑。

\[x(t)\,\to\, y(t)=\int_{-\infty}^{\infty}x(\tau)h(t-\tau)\text{d}\tau\tag{5}\]

　　刚才一直是把式（5）看成函数的叠加，如果在某个固定的时间点\(t=t_0\)观察这个式子，\(h(t_0-\tau)\)显然表示时间点\(\tau\)的单位冲激在\(t_0\)时刻的响应系数，而\(y(t_0)\)就是所有时间的响应叠加。在这样的理解下，我们就比较容易讨论因果性系统的输出函数。因果系统每个时刻的信号只能对之后的输出造成影响，在函数上的表现就是：当\(\tau>t_0\)时\(h(t_0-\tau)=0\)（也可以说，当\(t_0<0\)时\(h(t_0)=0\)）。使用换元法整理式（5），便得到因果LIT系统的表征函数如式（6）所示。特别地，无记忆系统的单位冲激响应可以计作\(h(t)=k\delta(t)\)，系统表征函数简化为\(y(t)=kx(t)\)，这就是我们直觉上认为“应该是”的那个线性系统。

\[x(t)\,\to\, y(t)=\int_0^{\infty}x(t-\tau)h(\tau)\text{d}\tau\tag{6}\]

　　对于一个稳定系统，要求\(x(t)\)有限时\(y(t)\)也有限，用反证法可知\(h(t)\)是绝对可积的（式（7））。反之如果式（7）成立，也不难证明\(y(t)\)有界（取\(|x(t)|\)的上限），从而我们知道：单位冲激响应绝对可积是系统稳定的充要条件。稳定系统的\(h(t)\)一般表现为渐进趋于\(0\)，如果它在有限时间后只有\(0\)值，这样的响应也叫有限冲激响应（FIR），这种系统也比较常见。另外，\(h(t)\)收敛于0的速度，一定程度上体现了系统响应的快慢。

\[\int_{-\infty}^{+\infty}|h(\tau)|\,\text{d}\tau<\infty\tag{7}\]

2.3 卷积

　　LIT的表征函数（5）由两个函数生成另一个函数，这样函数间的“运算”在实变函数里叫算子，而式（5）这样的算子被称为卷积。卷积就好像是把其中一个函数反转，然后移动着和另一个函数积分，明白式（5）的物理意义后，也能类推这种算子的作用了。一般把函数的卷积计作\(x(t)*y(t)\)，由式（5）中\(x(t),h(t)\)角色的对称性，不难得知卷积满足交换律（式（8））。

\[x(t)*y(t)=y(t)*x(t)\tag{8}\]

　　另外也不难证明（实变函数课程中有详细证明，也并不难），卷积还满足式（9）（10）的分配率和结合律（结合律本质上跟矩阵乘法的结合律相同）。既然LIT的表征函数是一个卷积\(x(t)*h(t)\)，再根据卷积的运算律，就能得到LIT系统的一些主要性质。比如系统\(h_1(t),h_2(t)\)并联后的系统满足：单位冲激响应为\(h_1(t)+h_2(t)\)；而串联（级联）后的系统满足：单位冲激响应为\(h_1(t)*h_2(t)\)，且串联的子系统可以任意调换顺序。反之如果系统的单位冲激响应可拆解为多个简单的响应，则可以方便地将系统拆解为子系统。

\[x(t)*[y(t)+z(t)]=x(t)*y(t)+x(t)*z(t)\tag{9}\]

\[x(t)*[y(t)*z(t)]=[x(t)*y(t)]*z(t)\tag{10}\]

　　不难验证，微分、积分本身就是LIT系统，根据卷积的结合律便有式（11）（12）成立，它们是卷积和微积分运算的互换法则。为了显式表示微积分的单位冲激响应，先考察恒定变换\(x(t)\to x(t)\)，根据式（4）可知它的单位冲激响应就是\(\delta(t)\)。利用结合律可知，微、积分系统的单位冲激响应就是\(\delta(t)\)的微、积分。把\(\delta(t)\)计作\(u_0(t)\)，它的（高阶）微积分计作\(u_n(t)\)（式（13），\(n\leqslant 0\)时都是奇异函数）。其中\(u_{-1}(t)\)就是上面提到的Heviside阶梯函数，它也被称为阶跃响应函数，简计作\(u(t)\)。再一次使用结合律，容易有式（14）的性质。

\[\dfrac{\text{d}}{\text{d}t}[x(t)*y(t)]=x(t)*\dfrac{\text{d}y(t)}{\text{d}t}\tag{11}\]

\[\int_{-\infty}^tx(\tau)*y(\tau)\,\text{d}\tau=x(t)*\int_{-\infty}^ty(\tau)\,\text{d}\tau\tag{12}\]

\[u_{n}(t)=\dfrac{\text{d}^n\delta(t)}{\text{d}t^n};\;\;u_{-n}(t)=\int_{-\infty}^t\dots\int_{-\infty}^\tau\delta(\tau)\,\text{d}\tau\tag{13}\]

\[u_m(t)*u_n(t)=u_{m+n}(t),\;\;(m,n\in\Bbb{Z})\tag{14}\]

本系列目录

　　　　01 - 线性时不变系统

　　　　02 - 傅里叶变换

　　　　03 - 系统函数的性质

　　　　04 - 离散时间信号

　　　　05 - 滤波、采样和通信

　　　　06 - 有理系统

 博客总目录在这里

---

【前序学科】微积分、傅里叶分析、常微分方程

【参考资料】

[1]  《信号与系统（2nd）》，奥本海姆，2013

　　本课程的主要参考书，内容详实，有丰富的例题、习题、以及现实中具体的问题。数学部分适合起点低的同学，但个人感觉缺少对本质的分析。过于周全的缺点便是，难以突显核心的结论，以及不能将关键信息立体地串联起来。