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

Posted 2021-03-11 edward-bian

1. 信号与系统

1.1 序言

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

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

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

1.2 信号

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

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

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

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

1.3 系统

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

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

   

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

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

[x_1(t) o y_1(t);;x_2(t) o y_2(t);;Rightarrow;;ax_1(t)+bx_2(t) o ay_1(t)+by_2(t) ag{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)视为连续变量并重新记为( au)。现在需要构造一个等价于(Delta(t- au))的微分(delta(t- au), ext{d} au)，也就是说微分仅在( au=t)处有非零值(1)，那就要求(delta(t- au))在(t)处有无穷大值且积分（面积）为(1)。为此需要定义一个满足式（3）的特殊的函数(delta(t))，它被称为单位冲激函数，坐标系中往往以一个单位长度的向上向量表示（左图）。在此定义下，(x(t))便有了如式（4）积分表达式。

  

[delta(t)=left{egin{matrix}infty,&(t=0),&(t e 0)end{matrix} ight.,;;int_{-infty}^{infty}delta(t), ext{d}t=1 ag{3}]

[x(t)=int_{-infty}^{infty}x( au)delta(t- au), ext{d} au ag{4}]

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

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

2.2 LIT的分解

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

[x(t), o, y(t)=int_{-infty}^{infty}x( au)h(t- au) ext{d} au ag{5}]

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

[x(t), o, y(t)=int_0^{infty}x(t- au)h( au) ext{d} au ag{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( au)|, ext{d} au<infty ag{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) ag{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) ag{9}]

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

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

[dfrac{ ext{d}}{ ext{d}t}[x(t)*y(t)]=x(t)*dfrac{ ext{d}y(t)}{ ext{d}t} ag{11}]

[int_{-infty}^tx( au)*y( au), ext{d} au=x(t)*int_{-infty}^ty( au), ext{d} au ag{12}]

[u_{n}(t)=dfrac{ ext{d}^ndelta(t)}{ ext{d}t^n};;;u_{-n}(t)=int_{-infty}^tdotsint_{-infty}^ audelta( au), ext{d} au ag{13}]

[u_m(t)*u_n(t)=u_{m+n}(t),;;(m,ninBbb{Z}) ag{14}]

 

本系列目录

　　　　01 - 线性时不变系统

　　　　02 - 傅里叶变换

　　　　03 - 系统函数的性质

　　　　04 - 离散时间信号

　　　　05 - 滤波、采样和通信

　　　　06 - 有理系统

 博客总目录在这里

 

---

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

【参考资料】

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

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