1. 引言
对于异结构系统的混沌同步或超混沌同步研究,如文献 [1] [2] 中所示,结构补偿还是相对常用的设计思路。结构补偿的作用是补偿响应系统与驱动系统在结构上的差异,以使两系统在同步控制器的作用下达到同步。如何设计出合适的结构补偿器值得研究。本文作者于文献 [3] 中提出中心平移同步法,该法通过将响应系统状态变量的中心平移至驱动系统需同步的状态变量,将结构补偿器的设计过程隐含于误差系统的建模中,将异结构系统的状态同步问题转化为非线性系统的稳定性控制问题。该方法为结构补偿器的设计提供了一种新思路。
对于具有不确定性系统的混沌同步或超混沌同步研究,如文献 [4] - [7] 中所示,自适应同步法仍是主流方法。文献 [8] 将扩张状态观测器和自抗扰控制器(ADRC)应用于混沌系统的观测和控制。考虑到扩张状态观测器不依赖于对象模型,因此自抗扰控制也可应用于具有不确定性系统的混沌同步研究中,如文献 [9] 中所示。在具有可变系数的三维混沌系统与五维超混沌系统的同步中,本文将三维混沌系统中的可变系数作为具有不确定性的参数,基于中心平移同步法设计同步控制器,使三维混沌系统的状态变量在可变参数的变化范围内均可与五维超混沌系统中的指定状态变量达到同步。该方法可应用于驱动系统模型具有不确定性但需同步状态及其导数可测和(或)响应系统模型确定但其包含变化范围已知的不确定参数的同步问题研究中。
2. 误差系统建模与同步控制器设计
2.1. 五维超混沌系统
五维超混沌系统的模型如下所示 [10]
(1)
其中,a = 23,b = 3,c = 18,m = 12,h = 4。
令五维超混沌系统(1)的初值为
,此时系统的Lyapunov指数分别为
,
,
,
,
,有两个大于0的Lyapunov指数,说明系统是超混沌的。
2.2. 三维混沌系统
三维混沌系统的模型如下所示 [11]
(2)
其中,
为可变系数。系统(2)随H变化的Lyapunov指数谱如图1中所示。当
,
时,系统(2)处于周期运动状态;当
,
,
时,系统(2)处于混沌运动状态。
2.3. 误差系统建模
令同步误差为
,根据文献 [3] ,误差系统的模型可表示为
(3)
未加入同步控制器
时,系统(3)等价于将三维混沌系统(2)中的状态变量
的中心平移至驱动系统(1)的指定状态变量
。现只需设计同步控制器
,使同步误差
全局渐近收敛于原点,即可使三维响应混沌系统(2)的状态变量在
内无论处于周期运动状态或是混沌运动状态,均可与五维驱动超混沌系统(1)中的指定状态变量达到同步。

Figure 1. Lyapunov exponent spectrum of the response system (2) versus increasing H
图1. 响应系统(2)随H变化的Lyapunov指数谱
2.4. 同步控制器设计
文献 [12] 中针对形如误差系统(3)的具有可变系数的三维混沌系统的稳定性控制问题,通过先对被控系统进行线性化,消去可变系数,求出使被控系统渐近稳定的控制参数范围,而后构造Lyapunov函数,基于线性化所得控制参数范围,采用试探法,求解可使被控系统在可变系数范围内全局渐近稳定的控制参数,有效的简化了控制参数的求解过程,且所得控制器可以是线性的。由于误差系统(3)与文献 [12] 中被控系统的模型完全相同,故可将同步控制器
设计为
。 (4)
将式(4)代入式(3),取径向无界的正定函数
,
作为误差系统(3)的Lyapunov函数,其导数为
, (5)
其中,矩阵Q的各阶顺序主子式分别为

要使矩阵Q是负定的,即
,只需
。文献 [12] 中已证明,当
时,在
内,
,此时矩阵Q是负定的,则
。又因为
是径向无界的,则误差系统(3)在原点处是全局渐近稳定的。现取
,则同步控制器
可表示为
(6)
3. 数值仿真
中心平移同步法的实质是将不易直观设计出的结构补偿器
, (7)
及同步控制器
加入到三维混沌系统(2)中。令五维超混沌系统(1)作为驱动系统,可得响应系统的模型为
(8)
式(8)与式(3)等价。
定义1 加入同步控制器与结构补偿器后,响应系统(8)的Lyapunov指数称为条件Lyapunov指数 [13] [14] 。
定理1(同步定理) 只有当响应系统(8)的所有条件Lyapunov指数都是负值时,才能达到响应系统(8)与驱动系统(1)的同步 [13] [14] 。
令响应系统(8)的初值为
,将式(6)及式(7)代入式(8)中,可得响应系统(8)随H变化的条件Lyapunov指数谱如图2中所示。从图2中可以看出,不同H值下系统(8)的条件Lyapunov指数均为负值,说明在
内,响应系统(8)与驱动系统(1)中的对应状态变量可以达到同步。
4. 算例
令
。未加入同步控制器
和结构补偿器
时,响应系统(8)的Lyapunov指数分别为
,
,
,说明系统是混沌的。此时响应系统(8)与驱动系统(1)各对应状态变量的误差及同步曲线分别如图3(a)和图3(b)中所示。图中横轴t表示式(8)的求解区间,为一无量纲量,后文各图中横轴t的定义与此处相同。从图3中可以看出,状态变量
、
、
的差异很大。
现将式(6)及式(7)代入式(8)中,可得响应系统(8)的条件Lyapunov指数分别为
,
,
,均为负值,说明响应系统(8)与驱动系统(1)中的对应状态变量可以达到同步。
加入同步控制器
和结构补偿器
后,响应系统(8)与驱动系统(1)各对应状态变量的误差及同步曲线分别如图4(a)和图4(b)中所示。从图4中可以看出,误差
、
、
分别渐近收敛于0,且收敛速度很快,
、
、
均可以很好的达到同步。这说明将中心平移同步法应用于对异结构系统指定状态变量间的同步控制是可行的与有效的。

Figure 2. Sub-Lyapunov exponent spectrum of the response system (8) versus increasing H
图2. 响应系统(8)随H变化的条件Lyapunov指数谱

(a) (b)
Figure 3. Error and synchronization curves of the corresponding state variables of the response and the drive systems before adding uc and us to the response system: (a) Error curves; (b) Synchronization curves
图3. 未加入uc和us时响应系统与驱动系统对应状态变量的误差及同步曲线:(a) 误差曲线;(b) 同步曲线

(a) (b)
Figure 4. Error and synchronization curves of the corresponding state variables of the response and the drive systems after adding uc and us to the response system: (a) Error curves; (b) Synchronization curves
图4. 加入uc和us后响应系统与驱动系统对应状态变量的误差及同步曲线:(a) 误差曲线;(b) 同步曲线
5. 中心平移同步法的优点
中心平移同步法的实质是将结构补偿器的设计过程隐含在误差系统的建模过程中,将混沌同步问题转化为非线性系统的稳定性控制问题。因此无需确定驱动系统的模型,只需驱动系统中各需同步状态变量及其导数已知或可通过状态观测器估算,即可通过搭建结构已知的响应系统模型完成对同步控制器的设计以及对驱动系统中需同步状态变量的同步跟踪。这有效的消除了驱动系统中的不确定性对同步控制器设计的影响。
6. 结论
本文基于中心平移同步法,通过将响应系统中待同步状态变量的中心平移至驱动系统中的指定状态变量,构造误差系统,使误差系统与响应系统具有相同的模型结构,这样便可将异结构系统的状态同步问题转化为非线性系统的稳定性控制问题。中心平移同步法可有效简化同步控制器中控制参数的求解过程。本文通过数值仿真验证了该同步方法的可行性与有效性。该方法为具有不确定性的异结构混沌、超混沌系统的同步问题提供了一种新的解决思路。但该方法的电路可实现性需进一步研究,即如何在不构建结构补偿器的基础上将中心平移同步法直接应用于待同步系统的电路实现中有待解决。