ISSN 
下载:
-
由于实际工程应用的需要以及研究的深入，实际问题常常面对复杂的工作环境和不确定因素，对系统精确的建模就变得尤为重要，分数阶微分模型为解决这一问题提供了新的思路，得到了越来越多的关注。相比于整数阶系统的控制发展，分数阶系统的控制和同步控制理论还很不完善，在某些方面甚至是空白，如对于线性定常的分数阶系统，目前还没有关于其稳定性的劳斯型判据，也很难找到合适的李雅普诺夫函数来进行系统分析和控制器设计。
目前，发展最为成熟的可用于指导控制器设计并进行稳定性分析的方法是针对线性定常同元次分数阶系统基于线性矩阵不等式 (LMI) 的分析与综合方法[1-2]，学者们也尝试了其他的控制方法，如根轨迹法[3]、H∞控制[4]、自抗扰控制[5]及最优控制[6]等。滑模控制瞬态反应快、对不确定性具有鲁棒性等特点，在分数阶混沌系统的同步控制设计中也得到了一些应用[7-14]。文献[15]设计了一种分数阶滑模面，开拓了分数阶混沌系统的滑模控制的新道路。文献[16]设计了一种新型四维自治超混沌系统，并找到其分数阶状态，设计了整数阶滑模控制器实现了异结构同步。文献[17]将自适应控制和模糊控制相结合，解决了带有不确定、参数未知和外部扰动的分数阶混沌系统的异同步问题。文献[18]针对分数阶混沌系统存在未知参数和未知外界扰动的情况，提出了一种分数阶自适应更新律，通过Mittag-Leffler稳定性定理证明了同步误差系统收敛。
上述及现有文献大多一是直接利用系统变量并基于Lyapunov稳定性理论进行分析，但是由分数阶定义可知，分数阶系统具有对历史特性的记忆性和全局特性，状态变量实际是对之前状态积累的一个反映，且积累的影响随时间增长而减弱，直接对状态变量进行分析显然是不准确的；二是所设计的自适应律随时间增长而增大，可能引起控制量无界的问题。
稳定性作为系统的基本性能要求，是控制器设计中最重要的环节。关于分数阶系统的稳定性理论仍处于探索阶段，目前比较常用的有分数阶稳定性理论、Mittag-Leffler定理[19]、分数阶拓展稳定性理论[20]、基于K类函数的分数阶Lyapunov直接法[21]。文献[22]介绍了阶次小于1的分数阶系统的稳定性理论；文献[23]从另一角度证明了分数阶稳定性理论不仅适用于分数阶自治系统，也同样适用于分数阶非自治系统。文献[24-26]对几种常用的分数阶系统稳定性理论方法进行了总结分析，这些方法虽然比最原始的稳定性判据更方便，但在指导控制器设计方面还存在较大困难。
综上所述，分数阶系统与控制研究既具有重要的工程意义和较广的应用前景，又充满困难和挑战。因此本文选取了一类存在参数摄动、未知函数及外部扰动等不确定因素的分数阶混沌系统进行研究，在传统滑模控制中引入分数阶微积分算子，既保持了传统滑模控制的理论特性，又增加了新的应用特性；引入频率分布模型[27]对系统模型，直接分析系统真实状态变量，设计一种新型分数阶滑模同步控制器，旨在丰富分数阶混沌系统的同步控制研究。
-
分数阶微积分实际上是指任意阶次的微积分，在其发展的历史中，学者们多次从不同角度对微积分进行了定义，下面简单介绍几种常用的定义。
定义1 分数阶积分由Riemann-Liouville积分定义[28-29]。函数f(t) 的n阶积分 (n为正实数) 为：
$${I^n}(f(t)) = \frac{1}{{\Gamma (n)}}\int_{{\rm{ }}0}^{{\rm{ }}t} {(t - \tau )} f(\tau ){\rm{d}}\tau $$ (1) 式中，$\Gamma (n) = \int_{{\rm{ }}0}^{{\rm{ }}\infty } {{t^{n - 1}}{{\rm{e}}^t}{\rm{d}}t} $，为Euler's Gamma函数。
分数阶算子是非整数阶基本运算${}_aD_t^\alpha $的微分和积分运算，其中a和t是运算的起止边界。3种常用的分数阶定义分别为Grünwald-Letnikov型定义、Riemann-Liouvill型定义和Caputo型定义[29-30]。
定义2 Riemann-Liouvill型分数阶导数定义为：
$${}_aD_t^\alpha f(t) = \frac{1}{{\Gamma (n - \alpha )}}\frac{{{{\rm{d}}^n}}}{{{\rm{d}}{t^n}}}\int_{{\rm{ }}a}^{{\rm{ }}t} {\frac{{f(\tau ){\rm{d}}\tau }}{{{{(t - \tau )}^{\alpha - n + 1}}}}} $$ (2) 式中，$n - 1 < \alpha < n$。
Riemann-Liouvill型分数阶导数有超奇异性，20世纪60年代，意大利的Caputo提出了弱奇异的分数阶微分定义。
定义3 Caputo分数阶导数定义为：
$$_aD_t^\alpha f(t) = \frac{1}{{\Gamma (n - \alpha )}}\int_a^t {\frac{{f(\tau ){\rm{d}}\tau }}{{{{(t - \tau )}^{\alpha - n + 1}}}}} \;\;\;n - 1 < \alpha < n$$ (3) 由定义可以看出，输入函数的初值以衰减的形式进行输出，零初值下分数阶微分是卷积分的形式。因此，分数阶微分实际上是一个积分，且随时间衰减记忆，这就是分数阶微分算子最大的特点。
本文采用工程应用较多的Caputo型定义，在处理具体问题时，可以选取不同类型的定义。为方便表达，用${D^\alpha }$代替${}_aD_t^\alpha $。
-
将分数阶混沌系统写为如下的形式：
$${D^\alpha }\boldsymbol{x}(t) = \boldsymbol{Ax}(t) + \boldsymbol{f}(x, t)$$ (4) 式中，$\boldsymbol{x}(t) = {[{x_1}(t), {x_2}(t), \cdots, {x_n}(t)]^{\rm{T}}}$为混沌系统的状态向量；$\boldsymbol{f}(x, t) \in {R^n}$为已知的光滑非线性函数；$\boldsymbol{A} \in {R^{n \times n}}$为系统的已知参数矩阵，$\boldsymbol{Ax}(t)$和$\boldsymbol{f}(x, t)$分别表示混沌系统中的线性部分和非线性部分。
由于式 (4) 所表示的系统既包括线性部分又包括非线性部分，具有非常广泛的代表性，能够包含大多数混沌系统，将现有文献中的分数阶混沌系统总结如表 1所示。
表 1 4种不同类型的分数阶混沌系统
类型 模型 A $\boldsymbol{f}(x, t)$ Chen系统 Lorenz系统 Lü系统 Liu系统 将式 (4) 表示的系统作为驱动系统，则被控响应系统表达式为：
$$\begin{array}{*{20}{c}} {{D^\alpha }\boldsymbol{y}(t) = (\boldsymbol{A} + \Delta {\boldsymbol{A}_y})\boldsymbol{y}(t) + \boldsymbol{g}(y, t) + }\\ {\Delta \boldsymbol{g}(y, t) + \boldsymbol{d}(t) + \boldsymbol{u}(t)} \end{array}$$ (5) 式中，$\boldsymbol{y}(t) \in {R^n}$为响应系统的状态向量；$\Delta {\boldsymbol{A}_y} \in {R^{n \times n}}$未知参数矩阵；$\boldsymbol{g}(y, t) \in {R^n}$为响应系统的非线性部分；$\Delta \boldsymbol{g}(y, t) \in {R^n}$表示不确定非线性函数向量，$\Delta {\boldsymbol{A}_y}$和$\Delta \boldsymbol{g}(y, t)$分别表示系统中的线性和非线性的不确定性；${\boldsymbol{u}}(t) \in {R^n}$为需要设计的同步控制器。
定义同步误差$\boldsymbol{e} = \boldsymbol{y}-\boldsymbol{ x}$。控制目标是设计合适的控制器${\boldsymbol{u}}(t)$，使得$\mathop {\lim }\limits_{t \to \infty } ||\boldsymbol{e}(t)|| \to 0$，实现驱动系统与响应系统的同步。
由式 (4) 和式 (5) 可得同步误差系统方程为：
$$\begin{array}{*{20}{c}} {{D^\alpha }\boldsymbol{e}(t) = \boldsymbol{A}e + \Delta \boldsymbol{A}y(t) + \boldsymbol{g}(y, t)-\boldsymbol{f}(x, t) + }\\ {\Delta \boldsymbol{g}(y, t) + \boldsymbol{d}(t) + \boldsymbol{u}(t)} \end{array}$$ (6) 假设1 不确定项$\Delta {\boldsymbol{A}_y} \in {R^{n \times n}}$，$\Delta \boldsymbol{g}(y) \in {R^n}$和外部扰动$\boldsymbol{d}(t) \in {R^n}$均有界，且满足：
$$||\Delta \boldsymbol{A}|{|_\infty } \le {\varepsilon _1}, \;\;||\Delta \boldsymbol{g}(y) + \boldsymbol{d}|{|_\infty } \le {\varepsilon _2}$$ (7) 式中，${\varepsilon _1}$和${\varepsilon _2}$为正常数。
-
以简单的分数阶线性滑模面为例，简要说明分数阶滑模控制相比于传统整数阶滑模控制的优势。分数阶滑模面的表达式为：
$${D^\alpha }x = - kx$$ (8) 当$\alpha = 1$时，式 (8) 为整数阶滑模面，收敛速度为${e^{ - kt}}$；当$\alpha \ne 1$时，式 (8) 为分数阶滑模面，收敛速度为${t^{ - \alpha }}$。图 1为分数阶滑模的基本原理[31]，由图可以看出，分数阶滑模运动相对于整数阶滑模运动产生的抖震较小，且能取得较高的控制精度。
图 1 整数阶滑模面与分数阶滑模面收敛特性对比图
$${D^\alpha }\boldsymbol{x}(t) = {\boldsymbol{A}_D}\boldsymbol{x}(t), \;\;\;x(0) = {x_0}$$ (9) 式中，$0<\alpha<1$，$\boldsymbol{x}(t) \in {R^n}$为系统状态向量；${\boldsymbol{A}_D} \in {R^{n \times n}}$为系统参数矩阵。当$\left| {\arg (\lambda ({\boldsymbol{A}_D}))} \right| \ge \alpha \frac{{\rm{\pi }}}{2}$时分数阶系统稳定；当$\left| {\arg (\lambda ({\boldsymbol{A}_D}))} \right| \ge \alpha \frac{{\rm{\pi }}}{2}$时分数阶系统渐近稳定，此时，状态向量以${t^{ - \alpha }}$速度收敛至零。
图 2 分数阶线性系统稳定性理论原理图
引理2[30] 分数阶系统${D^\alpha }y(t) = \upsilon (t)$，$0 < \alpha < 1$，$y(t) \in R$，$\upsilon (t) \in R$可以被看作如下一个线性频率分布模型：
$$\left\{ \begin{array}{l} \frac{{\partial z(\omega, t)}}{{\partial t}} = - \omega z(\omega, t) + \upsilon (t)\\ y(t) = \int_{{\rm{ }}0}^{{\rm{ }}\infty } {{\mu _\alpha }(\omega )z(\omega, t){\rm{d}}\omega } \end{array} \right.$$ (10) 式中，权重函数${\mu _\alpha }(\omega ) = \sin (\alpha \pi ){\rm{/\pi }}{\omega ^\alpha }$；系统状态$z(\omega, t) \in R$。
定义滑模函数为：
$$\boldsymbol{s}(t) = \boldsymbol{e}(t) + \boldsymbol{C}{I^\alpha }\boldsymbol{e}(t)$$ (11) 式中，$\boldsymbol{C} = {\rm{diag[}}{c_1}, {c_2}, \cdots, {c_n}{\rm{]}} \in {R^{n \times n}}$；$\boldsymbol{s}(t) \in {R^{n \times n}}$；$\boldsymbol{e}(t) \in {R^{n \times n}}$；${I^\alpha }$为由Riemann-Liouville积分定义的分数阶积分。
对式两边关于时间求$\alpha $阶次导，有：
$${D^\alpha }s(t) = {D^\alpha }\boldsymbol{e}(t) + \boldsymbol{Ce}(t)$$ (12) 为方便表达，记$\boldsymbol{s}(t)$为s。当系统发生滑模运动时，须满足$\boldsymbol{s} = 0$和${D^\alpha }\boldsymbol{s} = 0$，那么有：
$${D^\alpha }\boldsymbol{e}(t) = - \boldsymbol{Ce}(t)$$ (13) 根据引理1，选择合适的矩阵C，保证分数阶同步误差系统稳定，且C将决定$||\boldsymbol{e}|{|_1} \to 0$的速度。
根据引理2，可以得到滑模面的等价频率分布模型为：
$$\left\{ {\begin{array}{*{20}{l}} {\frac{{\partial \boldsymbol{z}(\omega, t)}}{{\partial t}} =-\omega \boldsymbol{z}(\omega, t) + {D^\alpha }\boldsymbol{s}}\\ {\boldsymbol{s} = \int_0^\infty {\mu (\omega )\boldsymbol{z}(\omega ){\rm{d}}\omega } } \end{array}} \right.$$ (14) 式中，权重函数为$\mu (\omega ) = \sin (\alpha {\rm{\pi }}){\rm{/\pi }}{\omega ^\alpha }$；$\boldsymbol{z}(\omega, t) \in R$为实际误差变量。选择Lyapunov函数为：
$${V_s}(t) = \frac{1}{2}\int_{ 0}^{ \infty } {\mu (\omega ){\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega $$ 对其求导得：
$$\begin{array}{c}{\dot V_s}(t) = \int_{ 0}^{ \infty } {\mu (\omega ){\boldsymbol{z}^{\rm{T}}}(\omega, t)\frac{{\partial \boldsymbol{z}(\omega, t)}}{{\partial t}}{\rm{d}}\omega } = \\ \int_{ 0}^{ \infty } {\mu (\omega ){\boldsymbol{z}^{\rm{T}}}(\omega, t)(-\omega \boldsymbol{z}(\omega, t) + {D^\alpha }\boldsymbol{s}){\rm{d}}\omega } = \\-\int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega + {\boldsymbol{s}^{\rm{T}}}{D^\alpha }\boldsymbol{s} = \\-\int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega + {\boldsymbol{s}^{\rm{T}}}({D^\alpha }e + Ce) = \\ - \int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega + \\ {\boldsymbol{s}^{\rm{T}}}(\boldsymbol{Ae} + \Delta \boldsymbol{Ay}(t) + \boldsymbol{g}(y, t) - \boldsymbol{f}(x, t) + \Delta \boldsymbol{g}(y, t) + \\ \boldsymbol{d}(t) + \boldsymbol{u}(t) + \boldsymbol{Ce})\end{array}$$ (15) -
设计控制器如下：
$$\begin{array}{c} \boldsymbol{u} =-\boldsymbol{Ae}-{{\hat \varepsilon }_1}{\mathop{\rm sgn}} (\boldsymbol{s})||y(t)||-{{\hat \varepsilon }_2}{\mathop{\rm sgn}} (\boldsymbol{s}) - \\ \boldsymbol{g}(y, t) + \boldsymbol{f}(x, t) - \boldsymbol{Ce} - \boldsymbol{Ks} \end{array}$$ (16) 式中，${\hat \varepsilon _1}, {\hat \varepsilon _2}$是${\varepsilon _1}, {\varepsilon _2}$估计值。将式 (16) 代入式 (15) 可得：
$$\begin{array}{c} {{\dot V}_s}(t) \le-\int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega + \\ {\boldsymbol{s}^{\rm{T}}}(\Delta \boldsymbol{A}y(t) + \Delta \boldsymbol{g}(y, t)-{{\hat \varepsilon }_1}{\mathop{\rm sgn}} (\boldsymbol{s})||\boldsymbol{y}(t)||-\\ {{\hat \varepsilon }_2}{\mathop{\rm sgn}} (\boldsymbol{s})) - {\boldsymbol{Ks}}{^{\rm{T}}}\boldsymbol{s} \le \\ - \int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)\, } {\rm{d}}\omega + \\ ||\boldsymbol{s}||{\varepsilon _1}||\boldsymbol{y}(t)|| + ||\boldsymbol{s}||{\varepsilon _2} - {{\hat \varepsilon }_1}||\boldsymbol{s}||||\boldsymbol{y}(t)|| - \\ {{\hat \varepsilon }_2}||\boldsymbol{s}|| - {\boldsymbol{Ks}^{\rm{T}}}\boldsymbol{s} = \\ - \int_{ 0}^{ \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega - \\ {{\tilde \varepsilon }_1}||\boldsymbol{s}||||\boldsymbol{y}(t)|| - {{\tilde \varepsilon }_2}||\boldsymbol{s}|| - {\boldsymbol{Ks}^{\rm{T}}}\boldsymbol{s} \end{array}$$ (17) 式中，${\tilde \varepsilon _1} = {\hat \varepsilon _1}-{\varepsilon _1}$和${\tilde \varepsilon _2} = {\hat \varepsilon _2}-{\varepsilon _2}$为估计误差。设计分数阶参数自适应律为：
$${D^\alpha }{\hat \varepsilon _1} = {q_1}(||\boldsymbol{s}||||\boldsymbol{y}(t)|| - {\gamma _1}{\hat \varepsilon _1})$$ (18) $${D^\alpha }{\hat \varepsilon _2} = {q_2}(||\boldsymbol{s}|| - {\gamma _2}{\hat \varepsilon _2})$$ (19) 式中，${q_1}, {q_2}, {\gamma _1}, {\gamma _2} > 0$为调节参数。显然，
$$\left\{ \begin{array}{l} {D^\alpha }{{\tilde \varepsilon }_1} = {D^\alpha }{{\hat \varepsilon }_1}-{D^\alpha }{\varepsilon _1} = {D^\alpha }{{\hat \varepsilon }_1}\\ {D^\alpha }{{\tilde \varepsilon }_2} = {D^\alpha }{{\hat \varepsilon }_2}-{D^\alpha }{\varepsilon _2} = {D^\alpha }{{\hat \varepsilon }_2} \end{array} \right.$$ (20) 根据引理2，得到参数估计误差的频率分布模型为：
$$\left\{ \begin{array}{l} \frac{{\partial {z_{{\varepsilon _1}}}(\omega, t)}}{{\partial t}} = - \omega {z_{{\varepsilon _1}}}(\omega, t) + {D^\alpha }{{\tilde \varepsilon }_1}\\ {{\tilde \varepsilon }_1} = \int_{ 0}^{ \infty } {\mu (\omega ){z_{{\varepsilon _1}}}(\omega, t)\, } {\rm{d}}\omega \end{array} \right.$$ (21) $$\left\{ \begin{array}{l} \frac{{\partial {z_{{\varepsilon _2}}}(\omega, t)}}{{\partial t}} = - \omega {z_{{\varepsilon _2}}}(\omega, t) + {D^\alpha }{{\tilde \varepsilon }_2}\\ {{\tilde \varepsilon }_2} = \int_{ 0}^{ \infty } {\mu (\omega ){z_{{\varepsilon _2}}}(\omega, t)\, } {\rm{d}}\omega \end{array} \right.$$ (22) 选择Lyapunov函数为：
$$\begin{array}{c} {V_\varepsilon }(t) = \frac{1}{{2{q_1}}}\int_{ 0}^{ \infty } {\mu (\omega )z_{{\varepsilon _1}}^2(\omega, t)\, } {\rm{d}}\omega + \\ \frac{1}{{2{q_2}}}\int_{ 0}^{ \infty } {\mu (\omega )z_{{\varepsilon _2}}^2(\omega, t)\, } {\rm{d}}\omega \end{array}$$ (23) 对其求导得：
$$\begin{array}{l} {{\dot V}_\varepsilon }(t) =-\frac{1}{{{q_1}}}\int_{ 0}^{ \infty } {\mu (\omega )\omega z_{{\varepsilon _1}}^2(\omega, t)} \omega + \frac{1}{{{q_1}}}{{\tilde \varepsilon }_1}{D^\alpha }{{\tilde \varepsilon }_1}-\\ \;\;\;\;\;\;\;\;\frac{1}{{{q_2}}}\int_{ 0}^{ \infty } {\mu (\omega )\omega z_{{\varepsilon _2}}^2(\omega, t)} {\rm{d}}\omega + \frac{1}{{{q_2}}}{{\tilde \varepsilon }_2}{D^\alpha }{{\tilde \varepsilon }_2} \end{array}$$ (24) 定理1 考虑形如式 (4) 的分数阶混沌系统作为驱动系统，形如式 (5) 的分数阶混沌系统作为响应系统，选取分数阶滑模面如式 (11) 的形式，选择式 (16) 表示的控制律，选择式 (18) 和式 (19) 表示的分数阶自适应律，那么同步误差系统 (6) 渐近稳定，且收敛至零的邻域。
证明：选择全局Lyapunov函数为：
$$V(t){\rm{ = }}{V_s}(t) + {V_\varepsilon }(t)$$ (25) 对式 (25) 关于时间t求导得，并利用不等式$2{\boldsymbol{\tilde \theta }^{\rm{T}}}\boldsymbol{\hat \theta } \ge {\boldsymbol{\tilde \theta }^2}-{\boldsymbol{\theta }^2}$，可得：
$$\begin{array}{*{20}{c}} {\dot V(t){\rm{ = }}{{\dot V}_s}(t) + {{\dot V}_\varepsilon }(t) \le }\\ {-\int_0^\infty {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega-}\\ {\frac{1}{{{q_1}}}\int_0^\infty {\mu (\omega )\omega \boldsymbol{z}_{{\varepsilon _1}}^2(\omega, t)} {\rm{d}}\omega-}\\ {\frac{1}{{{q_2}}}\int_0^\infty {\mu (\omega )\omega \boldsymbol{z}_{{\varepsilon _2}}^2(\omega, t)} {\rm{d}}\omega - }\\ {\frac{{{\gamma _1}}}{2}\tilde \varepsilon _1^2 + \frac{{{\gamma _1}}}{2}\varepsilon _1^2 - \frac{{{\gamma _2}}}{2}\tilde \varepsilon _2^2 + \frac{{{\gamma _2}}}{2}\varepsilon _2^2 - \boldsymbol{K}{\boldsymbol{s}^{\rm{T}}}\boldsymbol{s} \le }\\ { - \int_0^\infty {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} {\rm{d}}\omega - }\\ {\frac{1}{{{q_1}}}\int_0^\infty {\mu (\omega )\omega \boldsymbol{z}_{{\varepsilon _1}}^2(\omega, t)} {\rm{d}}\omega - }\\ {\frac{1}{{{q_2}}}\int_0^\infty {\mu (\omega )\omega \boldsymbol{z}_{{\varepsilon _2}}^2(\omega, t)} {\rm{d}}\omega + \frac{{{\gamma _1}}}{2}\varepsilon _1^2 + \frac{{{\gamma _2}}}{2}\varepsilon _2^2} \end{array}$$ (26) 将$\mu (\omega )$代入，根据定积分中值定理可知存在$\xi \in [0, \infty )$，使得：
$$\begin{array}{c} {V_s}(t) = \\ \frac{1}{2}\frac{{\sin (\alpha {\rm{\pi }})}}{{\rm{\pi }}}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {{\omega ^{-\alpha }}{\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega = \\ \frac{1}{2}\frac{{\sin (\alpha {\rm{\pi }})}}{{\rm{\pi }}}{\xi ^{-\alpha }}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {{\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega \end{array}$$ (27) 类似地，存在$\zeta \in [0, \infty )$，使得：
$$\begin{array}{c} \int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {\mu (\omega )\omega {\boldsymbol{z}^{\rm{T}}}(\omega, t)\boldsymbol{z}(\omega, t)} \, {\rm{d}}\omega = \\ \frac{{\sin (\alpha {\rm{\pi }})}}{{\rm{\pi }}}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {{\omega ^{1- \alpha }}{\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega = \\ \frac{{\sin (\alpha {\rm{\pi }})}}{{\rm{\pi }}}{\zeta ^{1- \alpha }}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {{\boldsymbol{z}^{\rm{T}}}(\omega, t)} \boldsymbol{z}(\omega, t){\rm{d}}\omega = \\ 2\frac{{{\zeta ^{1- \alpha }}}}{{{\xi ^{ - \alpha }}}}{V_s} = {\mu _1}{V_s} \end{array}$$ (28) 式中，${\mu _1} = 2\frac{{{\zeta ^{1 - \alpha }}}}{{{\xi ^{ - \alpha }}}}$为正常数。同理可知，存在常数${\mu _2}, {\mu _3} > 0$，使得：
$$\begin{array}{l} \frac{1}{{{q_1}}}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {\mu (\omega )\omega z_{{\varepsilon _1}}^2(\omega, t)} \, {\rm{d}}\omega = \\ {\mu _2}\frac{1}{{2{q_1}}}\int_0^\infty {\mu (\omega )z_{{\varepsilon _1}}^2(\omega, t)} \, {\rm{d}}\omega \end{array}$$ (29) $$\begin{array}{l} \frac{1}{{{q_2}}}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {\mu (\omega )\omega z_{{\varepsilon _2}}^2(\omega, t)} \, {\rm{d}}\omega = \\ {\mu _3}\frac{1}{{2{q_2}}}\int_{{\kern 1pt} 0}^{{\kern 1pt} \infty } {\mu (\omega )z_{{\varepsilon _2}}^2(\omega, t)\, } {\rm{d}}\omega \end{array}$$ (30) 因此由不等式 (26) 可知：
$$\dot V(t) \le - \rho V + Q$$ (31) 式中，$\rho = \min \{ {\mu _1}, {\mu _2}, {\mu _3}\} $；$Q = \frac{{{\gamma _1}}}{2}\varepsilon _1^2 + \frac{{{\gamma _2}}}{2}\varepsilon _2^2$为一有界的常数。
$$V(t) \le V(0){{\rm{e}}^{ - \rho t}} + \frac{Q}{\rho }(1 - {{\rm{e}}^{ - \rho t}})$$ (32) 因此，当时间趋向于无穷时，$V(t) \le Q/\rho $，系统同步误差及参数估计误差收敛到原点的邻域$\Omega = \{ V \le Q/\rho \} $内，系统渐近收敛。
-
选取分数阶超Chen系统进行仿真验证，驱动系统为：
$${D^\alpha }\boldsymbol{x} = \left[{\begin{array}{*{20}{c}} {-35}&{35}&0&1\\ 7&{12}&0&0\\ 0&0&{-8}&0\\ 0&0&0&{0.3} \end{array}} \right]\boldsymbol{x} + \left[\begin{array}{c} 0\\ -{x_1}{x_3}\\ {x_1}{x_2}\\ {x_2}{x_3} \end{array} \right]$$ (33) 初始值选取为$\boldsymbol{x}(0) = {(3, 2, 1, - 1)^{\rm{T}}}$。驱动系统出现混沌吸引子，如图 3所示。
图 3 分数阶超Chen混沌吸引子 ($\alpha = 0.98$)
被控响应系统模型为：
$$\begin{array}{c} {D^\alpha }\boldsymbol{y} = \left[{\begin{array}{*{20}{c}} {-35}&{35}&0&1\\ 7&{12}&0&0\\ 0&0&{-8}&0\\ 0&0&0&{0.3} \end{array}} \right]\boldsymbol{y} + \left[\begin{array}{c} 0\\ -{y_1}{y_3}\\ {y_1}{y_2}\\ {y_2}{y_3} \end{array} \right] + \\ \Delta \boldsymbol{g}(y) + \boldsymbol{d} \end{array}$$ (34) 式中，函数不确定项和外部干扰项为：
$$\left\{ \begin{array}{l} \Delta {g_1}(y) = 0.2\sin (4t){y_1}\\ \Delta {g_2}(y) = - 0.25\cos (5t){y_2}\\ \Delta {g_3}(y) = 0.15\cos (2t){y_3}\\ \Delta {g_4}(y) = - 0.2\sin (3t){y_4} \end{array} \right., \;\;\;\left\{ \begin{array}{l} {d_1} = 0.12\cos (2t)\\ {d_2} = - 0.2\sin (3t)\\ {d_3} = 0.1\cos (5t)\\ {d_4} = - 0.15\sin (t) \end{array} \right.$$ (35) 自适应律增益选择为${q_1} = 0.01$，${q_2} = 0.01$，响应系统阶次选为$\alpha = 0.98$，$\boldsymbol{C} = {\rm{diag(}}2, 2.15, 3, 1.6{\rm{)}}$，$\boldsymbol{K} = {\rm{diag(}}2, 3, 3, 2{\rm{)}}$。仿真结果如图 4所示。
图 4 同步误差曲线
从图中可以看到，驱动系统与响应系统的同步误差渐近收敛到零，意味着两个系统实现了同步。参数曲线如图 5所示，从图中也可以看出，使用分数阶非单增型自适应律可以使估计误差渐近收敛，避免参数估计随时间无限增大的问题，及在时间足够长的情况下，可能引起的控制量无界问题。
图 5 参数估计曲线
-
本文对分数阶混沌系统滑模控制研究现状进行了简要介绍，分析了目前研究存在的不足和难点，研究了一类带有不确定性和外部扰动的分数阶系统的滑模同步问题，充分利用已知条件，将系统不确定性分为系统名义矩阵部分和不确定部分。设计了一类新型分数阶PI滑模面和非单增型自适应鲁棒同步控制律，与现有文献中的滑模控制设计相比，能够有效减小抖震并缩短收敛时间，避免了随时间增长可能引起的控制量无界的问题。最后，选择了频率分布表示的Lyapunov函数证明该控制律能够控制系统误差状态收敛到滑模面上，避免了大多文献中直接用伪状态变量对同步误差系统进行分析的错误。
A Novel Sliding Mode Synchronization Method of Uncertain Fractional-Order Chaotic Systems
-
摘要: 针对一类存在参数摄动、未知函数及外部扰动等不确定因素的分数阶混沌系统的同步控制问题，设计了一类具有新颖的分数型积分滑模面的同步控制器。所设计的新型分数阶滑模面抖震更小、收敛速率更快。提出了一种改进的分数阶非增长型自适应律，有效避免了随时间增长可能引起的控制量无界的问题。引入频率分布模型分析系统模型，并基于Lyapunov稳定性定理证明同步误差收敛，避免了直接用伪状态变量对同步误差系统进行分析的错误，形成了分数阶运算和整数阶同步控制方法有机结合的新方法。仿真结果证明了该方法的有效性。Abstract: A novel adaptive sliding mode synchronization controller for uncertain fractional order chaotic systems is proposed in this paper. The proposed controller can realize the synchronization of systems with parameters perturbation, function uncertainties, and exterior disturbance. The new sliding mode surface with fractional order has advantages of week chattering and high convergent rate. A novel fractional order adaptive updating law is proposed to prevent the estimations from increasing infinitely. The sliding mode surface and the parameter estimation error of the controller are modelled by utilizing frequency distribution model, and the convergence of the error system is verified by Lyapunov functions, avoiding the mistake caused by directly applying pseudo state variables for the analysis of system synchronization error. Simulation results demonstrate the effectiveness of the proposed scheme.
-
Key words:
- adaptive control /
- chaos /
- fractional-order /
- sliding mode control /
- synchronization
-
[1] SABATIER J, MOZE M, FARGE C. LMI stability conditions for fractional order systems[J]. Computers & Mathematics with Applications, 2010, 59(5): 1594-1609. https://www.researchgate.net/publication/223596550_LMI_Stability_conditions_for_fractional_order_systems [2] FARGES C, MOZE M, SABATIER J. Pseudo-state feedback stabilization of commensurate fractional order systems[J]. Automatica, 2010, 46(10): 1730-1734. doi: 10.1016/j.automatica.2010.06.038 [3] MERRIKH-BAYAT F, AFSHAR M, KARIMIGHARTEMANI M. Extension of the root-locus method to a certain class of fractional-order systems[J]. ISA Transations, 2009, 48(1): 48-53. doi: 10.1016/j.isatra.2008.08.001 [4] PADULA F, ALCANTARA S, VILANOVA R, et al. H∞ control of fractional linear systems[J]. Automatica, 2013, 49(7): 2276-2280. doi: 10.1016/j.automatica.2013.04.012 [5] LI M, LI D, WANG J, et al. Active disturbance rejection control for fractional-order system[J]. ISA Transactions, 2013, 52(3): 365-374. doi: 10.1016/j.isatra.2013.01.001 [6] DJENNOUNE S, BETTAYEB M. Optimal synergetic control for fractional-order systems[J]. Automatic, 2013, 49(7): 2243-2249. doi: 10.1016/j.automatica.2013.04.007 [7] 吴学礼, 刘杰, 张建华, 等.基于不确定性变时滞分数阶超混沌系统的滑模自适应鲁棒的同步控制[J].物理学报, 2014, 63(16): 160507. doi: 10.7498/aps.63.160507 WU Xue-li, LIU Jie, ZHANG Jian-hua, et al. Synchronizing a class of uncertain and variable time-delay fractional-order hyperchaotic systems by adaptive sliding robust mode control[J]. Acta Phys Sin, 2014, 63(16): 160507. doi: 10.7498/aps.63.160507 [8] GAO L. Modified sliding mode synchronization of typical three-dimensional fractional-order chaotic systems[J]. Neurocomputing, 2015, 166: 53-58. doi: 10.1016/j.neucom.2015.04.031 [9] YANG C C. Synchronization of second-order chaotic systems via adaptive terminal sliding mode control with input nonlinearity[J]. Journal of the Franklin Institute, 2012, 349: 2019-2032. doi: 10.1016/j.jfranklin.2012.03.012 [10] XIANG W, HUANGPU Y G. Second-order terminal sliding mode controller for a class of chaotic systems with unmatched uncertainties[J]. Communications in Nonlinear Science and Numerical Simulation, 2010, 15(11): 3241-3247. doi: 10.1016/j.cnsns.2009.12.012 [11] WANG C C, NS P, YAU H T. Chaos control in AFM system using sliding mode control by backstepping design[J]. Commun Nonlinear Sci Numer Simulat, 2010, 15(3): 741-751. doi: 10.1016/j.cnsns.2009.04.031 [12] HAERI M, EMADZADEH A A. Synchronization different chaotic systems using active sliding mode control[J]. Chaos Solitons & Fractal, 2007, 31(1): 119-129. https://www.researchgate.net/publication/254010020_Modified_Function_Projective_Synchronization_of_Two_Different_Chaotic_Systems_Using_Active_Sliding_Mode_Control [13] 潘光, 魏静.一种分数阶混沌系统同步的自适应滑模控制器设计[J].物理学报, 2015, 64(4): 040505. http://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201504007.htm PAN Guang, WEI Jing. Design of an adaptive sliding mode controller for synchronization of fractional order chaotic systems[J]. Acta Phys Sin, 2015, 64(4): 040505. http://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201504007.htm [14] 张洪, 陈天麒.参数未知的时延混沌系统滑模变结构同步控制[J].电子与信息学报, 2005, 27(12): 1937-1941. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX200512019.htm ZHANG Hong, CHEN Tian-qi. Synchronization of Time-delay chaotic systems with uncertainty based on sliding mode controller[J]. Journal of Electronics & Information Technology, 2005, 27(12): 1937-1941. http://www.cnki.com.cn/Article/CJFDTOTAL-DZYX200512019.htm [15] MOHAMMAD P A. Robust stabilization and synchronization of a class of fractional-order chatic systems via a novel fractional sliding mode controller[J]. Commun Nonlinear Sci Numer Simulat, 2012(17): 2670-2681. https://www.researchgate.net/publication/241105766_Robust_stabilization_and_synchronization_of_a_class_of_fractional-order_chaotic_systems_via_a_novel_fractional_sliding_mode_controller [16] YUAN G, LIANG C H, WU Q Q, et al. A new fractional-order hyperchaotic system and its modified projective synchronzation[J]. Chaos, Solitons & Fractals, 2015(76): 190-204. https://www.researchgate.net/publication/275673296_A_new_fractional-order_hyperchaotic_system_and_its_modified_projective_synchronization [17] BOULKROUNE A, BOUZERIBA A, BOUDEN T. Fuzzy generalized projective synchronization of incommensurate fractional-order chaotic systems[J]. Neurocomputing, 2015(44): 791-801. https://www.researchgate.net/publication/283894915_Fuzzy_generalized_projective_synchronization_of_incommensurate_fractional-order_chaotic_systems [18] WANG Q, DING D S, QI D L. Mittag-Leffler synchronization of fractional order uncertain chaotic systems[J]. Chinese Physics B, 2015, 24(6): 060508. doi: 10.1088/1674-1056/24/6/060508 [19] LI Y, CHEN Y, PODLUBNY I. Mittag-Leffler stability of fractional order nonlinear dynamic systems[J]. Automatica, 2009, 45(8): 1965-1969. doi: 10.1016/j.automatica.2009.04.003 [20] DELACARI H, BALEANU D, SADATI J. Stability analysis of Caputo fractional-order nonlinear systems revised[J]. Nonlinear Dynamics, 2012, 67(4): 2433-2439. doi: 10.1007/s11071-011-0157-5 [21] NUSSBAUM R D. Some remarks on a conjecture in parameter adaptive control[J]. Systems & Control Letters, 1983, 3(5): 243-246. [22] MATIGNON D. Stability results for fractional differential equations with applications to control processing[J]. Computational Engineering in Systems Applications, 1996(2): 963-968. https://www.researchgate.net/publication/2581881_Stability_Results_For_Fractional_Differential_Equations_With_Applications_To_Control_Processing [23] 胡建兵, 赵灵冬.分数阶系统稳定性理论与控制研究[J].物理学报, 2013, 62(24): 240504. doi: 10.7498/aps.62.240504 HU Jian-bing, ZHAO Ling-dong. Stability theorem and control of fractional systems[J]. Acta Phys Sin, 2013, 62(24): 240504. doi: 10.7498/aps.62.240504 [24] MACHADO J T, KIRYAKOVA V, MAINARDI F. Recent history of fractional calculus[J]. Communications in Nonlinear Science and Numerical Simulation, 2011, 16(3): 1140-1153. doi: 10.1016/j.cnsns.2010.05.027 [25] LI C P, ZHANG F R. A survey on the stability of fractional differential equations[J]. European Physical Journal, 2011, 193(1): 27-47. https://www.researchgate.net/publication/225144966_A_survey_on_the_stability_of_fractional_differential_equations [26] CHEN Y Q, PETRÁŠ I, XUE D. Fractional order control-tutorial[C]//Proceedings of the American Control. St, Louis, Mo, USA: [s.n.], 2009: 1397-1141. [27] TRIGEASSOU J C, MAAMRI N, SABATIER OUTSTALOUP J A. A Lyapunov approach to the stability of fractional differential equations[J]. Signal Processing, 2011, 91: 437-445. doi: 10.1016/j.sigpro.2010.04.024 [28] OUSTALOUP A. La dérivation non-entière: Théorie, synthèse et applications[M]. Paris: Hermès Editeur, 1995. [29] PODLUBNY I. Fractional differential equations: an introduction tractional derivatives, fractional differentials equations, to methods of their solution and some of their applications[M]. San Diego, Calif, USA: Academic Press, 1998. [30] MONJE C A, CHEN Y, VINAGRE B M, et al. Fractional-order systems and controls[M]//Fundamentals and Applications. London, UK: Springer, 2010. [31] ZHANG B T, PI Y G, LUO Y. Fractional order sliding-mode control based on parameters auto-tuning for velocity control of permanent magnet synchronous motor[J]. ISA Transactions, 2012, 51(5): 649-656. doi: 10.1016/j.isatra.2012.04.006 -
点击查看大图
图(5) / 表(1)
计量
- 文章访问数: 4561
- HTML全文浏览量: 1299
- PDF下载量: 469
- 被引次数: 0
