不确定分数阶非线性系统自适应滑模同步控制

高 瑜，李 雄

（1.陕西铁路工程职业技术学院基础课部，陕西渭南714000；2.西安欧亚学院数理与信息技术应用中心，陕西西安710065）

近年来，分数阶混沌系统成为了最热门的非线性系统研究的领域，推动了分数阶非线性系统稳定性分析及同步控制方法的发展[1-3]。混沌系统同步控制由于在通信领域的广泛应用得到了研究者的重视，人们相继提出了很多种分数阶混沌系统的同步控制方法[4]，如滑模变结构控制法、自适应控制法、模糊控制法、脉冲控制法和Backstepping控制法等[5-7]。对于不确定分数阶非线性系统同步控制也有一些结果[8-9]，如文献[10]在系统不确定项满足有界的情况下利用滑模控制实现了不确定分数阶Duffing-Holmes系统的同步问题，文献[11]研究了不确定分数阶混沌系统的自适应模糊同步控制问题等。

自适应滑模控制方法常用来研究带不确定项的分数阶非线性系统，并且在稳定性分析中通常构造平方Lyapunov函数。随着文献[2]提出了分数阶系统的Lyapunov第二方法，对于分数阶非线性系统的控制及稳定性分析逐渐成为研究热点。但平方函数具有非常复杂的分数阶导数形式，这也使得分数阶非线性系统稳定性分析中无法应用平方Lyapunov函数。所以到目前为止几乎没有文献成功实现了分数阶混沌系统自适应滑模控制或同步。随着研究的深入，许多分数阶模型不仅需要满足渐近稳定，更需要在有限时间内稳定，这也给非线性系统稳定性分析带来了难度。文献[12]研究了分数阶非线性系统在有限时间内不存在稳定点的问题，推动了非线性系统有限时间稳定性理论的进一步发展。在非线性系统有限时间稳定理论中，文献[13]通过变量替换和函数构造提出了一个新的非线性系统有限时间稳定的充分条件，具有很强的推广性，但是只是针对一类整数阶非线性系统，对于分数阶非线性系统有限时间稳定性理论的研究还尚未深入。

文中主要研究了基于滑模控制的不确定分数阶非线性系统同步，首先针对二维分数阶混沌系统，通过构造分数阶滑模面及分数阶微分方程形式的自适应规则，设计了同步控制器，并利用分数阶Lyapunov第二方法证明了构造方法的合理性（需要指出的是本文系统中的不确定项可以是完全未知的）。以分数阶Duffing-Holmes系统和分数阶Van der Pol系统为实例，实现了驱动系统和响应系统的异结构有限时间同步控制（即在有限时间内误差系统趋于滑模面），验证了该方法和控制器的有效性。

1　预备知识

分数阶微积分概念提出了很多种定义，文中采用Caputo定义作为研究工具。

分数阶积分定义为

Caputo分数阶微分定义为

其中，Γ（∙）为Gamma函数。

当0<α<1时，Caputo分数阶微分的解等价于

定义1双参数Mittag-Leffler函数定义为

其中α,β>0,z为复数，Γ（∙）为 Gamma函数，其Laplace变化定义为

其中R（s）为s的实部，λ∈R，L（∙）为Laplace变换。

引理1如果那么x（t）在[0,+∞）上单调减少；同理则x（t）在 [0,+∞）上单调增加。

证明只需要证明引理1的前半部分，后半部分证明思路相同。因为,所以一定存在非负可积函数y（t）使得

对上式两边同时取α阶积分可得

若有t2>t1≥0时，有

这也就证明了x（t）在[0,+∞）上单调减少。

引理2设且具有连续的一阶导数，则

其中P为任意的n阶正定矩阵。

引理3（分数阶Lyapunov第二方法）设原点是如下分数阶非线性系统的平衡点：

其中x（t）∈Rn为系统变量，f（t,x（t）为满足局部Lipschitz条件的非线性函数。若存在Lyapunov函数V（t,x（t）和Κ类函数αi（i=1,2,3）使得

则系统渐近稳定。

引理4若满足以下等式：

其中x（t）和y（t）∈Rn具有连续的一阶导数 ，P,Q∈Rn×n为两个正定矩阵。若存在正定的矩阵M和正常数h使得

2　问题描述

考虑如下的二维不确定分数阶混沌系统

其中α∈（0,1），X（t）=[x1,x2]∈R2为系统输入变量，f（X,t）∈R为非线性函数，Δf（X）∈R为系统的不确定项，dx（t）∈R为随机扰动，u（t）∈R为控制变量。

考虑如下的响应系统

其中Y（t）=[y1,y2]∈R2为系统响应变量，g（Y,t）∈R为非线性函数，Δg（Y）∈R为系统的不确定项 ，dy（t）∈R为随机扰动。

定义如下的同步误差系统

定义 2.1如果存在正常数T=T（E（0）>0,使得当成立，则误差系统（3）关于T有限时间稳定。

假设 2.1系统不确定项 Δg（Y）,Δf（X）为有界变量，即存在正常数γ1，使得

假设 2.2系统随机扰动dx（t）,dy（t）为有界变量，即存在正常数γ2，使得

3　主要结果

3.1　分数阶滑模面设计

设计如下的分数阶滑模面

当系统发生滑模运动时，需满足如下条件

通过简单的证明推导，可以得出上式是渐近稳定的，即误差系统变量趋于零。

3.2　控制器设计

本文所要讨论的问题是如何设计同步控制器，使得误差系统能在有限时间内达到或趋近于滑模面

由式（4）可以得到，

为了实现同步误差系统能够在有限时间内稳定，本文设计如下的自适应规则：

因此，可以设计如下的控制器：

定理1给定初始条件下，设计如上的自适应滑模控制器和自适应规则的作用下，同步误差系统能够在有限时间内趋近滑模面，即实现了驱动系统和响应系统同步控制。

证明 构造 Lyapunov函数为V（t）=| |s（t），并对其两边同时取α阶导数得

将滑模面方程带入上式中得

由假设2.1和假设2.2可得

其中ξ=min{ξ1,ξ2}。

由上式可以进一步得

对上式两边同时取（0,t）上的积分得

4　数值仿真

仿真中驱动系统选取为分数阶Duffing-Holmes系统，响应系统选取为分数阶Van der Pol系统。系统不确定项与随机扰动分别选取如下：

选取系统初值

x1（0）=0.2,x2（0）=-0.2,y1（0）=-0.1,y2（0）=0.3.给定参数

由定理1设计如下滑模面与同步控制器：

仿真结果如图1，图2，图3所示。

图1　同步误差曲线图

5　结论

文中研究了不确定分数阶非线性系统自适应滑模同步控制，通过构造分数阶滑模面以及分数阶自适应规则，在满足系统所有变量有界的情况下，利用Lyapunov函数证明了定理的有效性和鲁棒性。基于该理论设计了自适应控制器，实现了二维分数阶Duffing-Holmes系统与分数阶Van der Pol系统异结构有限时间同步，通过合理选取初值与参数值进行数值仿真，可以得到误差系统能够在有限时间内趋于滑模面。该理论的研究有助于掌握分数阶非线性系统的相关性质，同步控制方法也具有良好的鲁棒性。本文所研究的方法仍需进一步改进，针对不同的阶次控制效果可能出现差异性，更严格的控制输入条件下（如带死区）实现自适应同步控制需要进一步的研究。

图2　驱动系统和响应系统同步曲线（x1-y1）

图3　驱动系统和响应系统同步曲线（x2-y2）