前言
物理系统天生都具有非线性,所以在某种程度上来说所有的系统都是非线性的。非线性系统可以用非线性差分方程表示。如果一些非线性系统在小范围内顺滑,通过一些假设,并且限定工作范围,它是可以被等效为线性系统并且用线性差分方程来描述的。在讲解非线性之前,我们先来回顾一下线性系统
线性系统
线性控制理论主要就是研究线性时不变(LTI)控制系统,一个具有如下形式的简单的LTI系统(x是状态向量,A是系统矩阵)
有以下几个性质:
- 如果A是非奇异矩阵的话那么该线性系统只有一个平衡点
- 如果无论任何初始条件下A的所有特征值都有负实部,那么系统的平衡点是稳定的
- 线性系统的瞬时响应是由系统的natural modes(不知道咋翻译)组成的,且它的一般解具有解析解
- 如果系统存在外部输入,公式变成
这样的系统有几个有趣的性质。1.满足叠加原理 2.具有有限输入有限输出稳定性 3.正弦输入会得到频率一样的正弦输出
非线性
非线性可以分成自然非线性和人为的非线性,自然非线性来自于系统硬件本身以及运动。例子有旋转运动的离心力以及接触平面的库仑摩擦力。通常这样的非线性特征是我们不希望遇到的,需要用控制系统补偿它们。而人为的非线性是由系统设计着引入的非线性。例如自适应控制以及启停式控制,都是传统的人为引入的非线性。
非线性也可以根据数学性质分为连续的和离散的。离散的非线性不能被局部等效为线性方程,因此离散非线性也被称为强非线性(滞后,回转间隙,精摩擦力等)。对于这样的强非线性系统,是否要在小范围内等效成线性取决于非线性的程度,这个我们后面5.2节会讲。
一个非线性系统特征举例
由于不满足线性和叠加原理,非线性系统的行为要比线性系统更复杂。举个例子
例1.1 一个简化的水下探测器的运动模型可以表示为:
v探测器的速度,u是控制输入。其中的非线性部分|v|v对应了一个经典的平方律拖拽现象,假设我们给一个正向单位阶跃输入,5秒后再给一个反向单位阶跃输入。系统的响应如下图(这里系统五秒后的输入图像貌似不对,不过不影响理解):
我们稍微解释一下这个输出,探测器在恒定推力的情况下速度上升,但受到的阻力也在上升。最终在4秒左右的时候达到平衡。我们发现系统对于正向输入的响应速度要比逆向的快,而且减速的时候高速区的减速要比低速区的快来不少,这是因为方程里明显的的阻尼项|v|在高速时的值要比低速时大
如果我们把输入变成10,系统响应变成下图。
可以发现系统是变快了,但是稳定时间不像线性系统并没有快十倍,最高速度也没有大十倍。我们算出两次不同输入得到的稳定速度
这是一个很传统的水下探测器的非线性特征,因此如果要实现对探测器良好的控制的话,需要仔细考虑其中的非线性部分
一些常见的非线性系统特征
现在我们总结一下非线性系统具有哪些值得注意的特征
1.多平衡点
不像线性系统只有一个平衡状态,非线性系统往往有多个平衡状态,即可以一直停在那里的状态,举个栗子
例1.2:一个一阶系统
如果假定初始状态为x(0)=xo,那么它的线性化为
可以求得它的解为(惭愧,不会算)
下图画出了不同初值时x(t)的值
可以发现无论初值怎么变,它只有一个平衡状态就是x=0
如果不采用线性化直接求导,得到解析解为
对于不同的初值得到的图像如下
可以发现有两个与初值相关的平衡点x=0与x=1。有了上面的例子我们现在可以简单讨论运动稳定性的问题。对于一个稳定的线性系统,它的稳定性不随初值变化,始终保持稳定,且有一个平衡点。但对非线性系统来说,由于初值的不同会导致系统去到不同的平衡点,甚至变得不稳定(通常在有限时间内变得不稳定,被称之为有限逃逸时间)。这告诉了我们非线性系统的稳定性与初值有关。
另外对于一个有界的外部输入,输入值也会影响非线性系统的稳定性产生影响,比如下面这个系统
如果输入u=-1,x收敛于0;如果u=1,|x|趋向于无穷大
2.极限环
非线性系统可能会呈现固定振幅固定周期无外部激励的自震荡现象,这些自震荡被称之为极限环或者自激励震荡。这个重要的现象由荷兰电子工程师Van der Pol在1920年代第一次研究发现
例1.3 范德玻尔方程
考虑如下二阶非线性差分方程:
m,c,k都是正的常量。该方程可以看成是一个带有与位置相关的阻尼项2c(x2-1)的质量-弹簧-阻尼 系统,如果x值大一些,那么阻尼项是正的且在消耗系统里的能量,这表示该系统有一个收敛的趋势。而当x小一些时,阻尼项为负并开始往系统里面增加能量(这里能量的减少与增加是怎么看出来的?),表示系统有一个发散的趋势。因此,由于非线性项随着x一直改变,系统的运动呈现一个持续不断,且与初值无关的震荡,既不会无限发散也不会收敛到0。如下图所示
这个“极限环”通过阻尼项一直从环境中吸收和释放能量来维持周期性的震荡,这与传统的质量弹簧系统是不一样的,传统的质量弹簧系统不会在震荡过程中与环境有能量交互。
线性系统里面也有周期震荡,但线性系统的周期震荡和这里的极限环是有区别的。首先,非线性系统的自激励震荡幅度与初值无关,而线性系统的周期震荡幅度由初值决定;另外,线性系统的对于系统参数的变化十分敏感,而非线性系统极限环受参数变化的影响很小。
极限环代表了一类非线性系统里面很重要的一个现象,在自然与工程中很常见。比如机翼的风抖,一种由结构振动和空气动力引发的极限环现象,这很常见也对飞行器很不安全。另外腿结构机器人的单足跳跃运动也是极限环的一个例子。可以发现,极限环有的时候对系统是有害的,有的时候是有用的。有害的时候作为一个控制工程师就要知道怎么消除它,有用的时候就要知道怎么激励并利用它。要做到这个需要理解极限环的性质,以及熟悉一些操控极限环的方法工具
3.分岔
当非线性动态系统的一些参数变化时,平衡点的稳定性和平衡点的数量有可能会随之改变,这就叫分岔现象,这些使系统发生质变的参数数值被称之为关键数值或者分岔数值(critical value or bifurcation value),研究分析现象就是分岔理论所要讨论的问题。
拿烟囱里的烟雾上升现象举个栗子,烟一开始会加速上升,因为它比周围空气轻。上升到一定程度后变成涡流。我们可以用无阻尼达芬方程来描述:
我们画出平衡点随着α变化的曲线
当α从正慢慢变成负的时候,平衡点从一个平衡点分裂成三个平衡点。这代表了α=0就是一个分岔数值,这样的分岔被称之为干草叉分岔。
还有另一种分岔,这种分岔往往包含了极限环的出现。在这种情况下,参数的变化使得一对复共轭特征值从左平面跨越到右平面,不稳定的系统分岔成了一个极限环,如下图
4.混沌
对于稳定线性系统,初值的小变化只能引起输出的小变化。而非线性系统,可能会呈现出混沌现象,意思是系统响应对初值十分敏感。混沌的一个重要特征就是输出结果的不可预测性,即便我们有了准确的非线性系统模型以及很强的计算机,系统的长期响应依旧不能很好地预测。
我们必须要能区分混沌与随机运动的区别,在随机运动中,系统模型或者输入有很强的不确定性,因此造成了输出的不可预测。在混沌中,系统与输入有着很小的不确定性,但依旧不可预测输出。让我们考虑这样一个简单的非线性系统
该方程代表了一类轻阻尼,受正弦力驱动并承受较大弹性形变的机械结构。下图展现了该系统两种不同初值下的输出
其中一个初值为(粗线)
另外一个初值为(细线)
可以发现由于x的五次方项带来的强非线性,系统在一段时间后呈现很大的差异,且基本看不出规律。
混沌现象在很多物理系统中都存在,最常见的就是流体力学中的扰动,大气动力学也呈现出明显的混沌行为,因此长期的气候预测是不可能的。
混沌主要在强非线性系统中出现,因此如果一个给定的系统初值或者外部输入使系统呈现强非线性区域,那么发生混沌的可能性也会相应增加。线性系统中不会有混沌,任意幅度的正弦输入得到的总是同样频率的正弦输出。
在反馈控制的内容中,研究混沌啥时候出现并且如何避免它是一个热点研究课题。
5.其他行为
除了上述行为还有其他非线性会出现的行为,都充分说明了非线性系统要比线性系统更复杂且难以预测