无平衡点分数阶混沌系统全状态自适应控制
邵克勇, 陈丰, 王婷婷, 王季驰, 周立朋
东北石油大学 电气信息工程学院,黑龙江 大庆 163318

作者简介:邵克勇(1970-),男,教授,博士.研究方向:鲁棒控制,智能控制.E-mail:1127073951@qq.com

摘要

针对一类无平衡点的分数阶混沌系统,首先通过分数阶微分变换方法(FDTM)得到它的解序列。然后,研究了系统的Kaplan-Yorke维数和耗散性,基于系统的离散映射通过QR分解得到最大Lyapunov特征指数,通过该特征指数可以判断系统是否保持混沌。最后,给出一种全状态自适应控制方法,使系统的状态变量追踪期望轨迹,并通过数值模拟验证了本文算法的可行性。

关键词: 自动控制技术; 人工智能; 无平衡点的混沌系统; 分数阶微分变换方法; 全状态自适应控制; Lyapunov特征指数
中图分类号:TP273 文献标志码:A 文章编号:1671-5497(2017)04-1225-06
Full state based adaptive control of fractional order chaotic system without equilibrium point
SHAO Ke-yong, CHEN Feng, WANG Ting-ting, WANG Ji-chi, ZHOU Li-peng
College of Electrical Engineering and Information, Northeast Petroleum University, Daqing 163318,China
Abstract

For a class of chaotic system without equilibrium point, its solution series is obtained by the Fractional Differential Transformation Method (FDTM). Then, the Kaplan-Yorke dimension and the dissipativity of the system are investigated. The Lyapunov characteristic exponents are calculated based on the discrete map of the system through QR factorization algorithm, and the largest Lyapunov characteristic exponent is applied to judge whether the system keeps chaos or not. Finally, a full state based adaptive controller for the fractional order chaotic system is designed, which allows the states of the system to track the desired constant. The feasibility of the proposed algorithm is verified by numerical simulation.

Keyword: automatic control technology; artificial intelligence; chaotic system without equilibrium point; fractional differential transformation method(FDTM); full state based adaptive control; Lyapunov characteristic exponent
0 引 言

分数阶系统是包含分数阶积分与微分的动态系统, 它可以有效地应用于物理系统建模[1, 2, 3]。在整数阶微积分的意义下, Poincare-Bendixon定理指出:连续自治系统出现混沌现象的微分阶次必须大于等于3[4, 5]。但是, 人们在研究中发现, 在分数阶微积分的意义下, 在许多阶次低于3的分数阶非线性系统中, 仍然存在混沌现象[6]。由于分数阶混沌系统相比于整数阶具有更好的发展前景, 它的理论研究及实际应用受到国内外学者的普遍关注。然而该领域内的研究成果仍然偏少, 其中一个主要原因是分析非线性分数阶系统稳定性的理论工具具有局限性。文献[7]提出了基于线性分数阶系统稳定性理论的混沌控制方法, 然而这种方法只验证了局部的稳定性。文献[8]提出了主动控制的方法, 这种方法利用非线性控制律来消除控制系统中的非线性项, 再根据线性分数阶系统稳定性理论, 设计一个线性控制律来稳定此系统, 但当系统中的参数未知时, 不能用主动控制方法来设计控制器。文献[9]研究了一类整数阶新型金融系统的自适应控制。文献[10, 11]研究了分数阶混沌系统的自适应控制, 其中控制变量存在于局部系统方程中。

本文针对一类无平衡点的分数阶混沌系统, 设计了一种全状态自适应控制器, 在不消除非线性项的情况下来追踪期望轨迹。最后通过数值模拟, 验证了该算法的可行性。

1 预备知识

本节给出了分数阶微积分和微分变换方法的一些定义和性质[12]

1.1 分数阶微积分

t0Dtq表示分数阶微分积分算子。当 q为正值时, 它是一个微分算子; 当 q为负值时, 它是一个积分算子。分数阶微分的定义通常用Grunwald-Letnikov, Riemann-Liouville和Caputo来表示。

t0DtqDtq=dqdtq, q> 01, q=0t0t()qI-q, q< 01

Riemann-Liouville定义[13], 表示如下:

t0Dtqft=1Γ1-qddt0tfτt-τq-n+12

式中: n-1< q< n; Γ·是Gamma函数。

根据方程(1)(2), 分数阶微分和积分算子的阶次为 0< q< 1, t0=0(本文始终考虑 t0=0的情况), 并且分别定义为:

0Dtqft=1Γ1-qddt0t(t-τ)-qfτ0Itqft=1Γqddt0t(t-τ)q-1fτ

由此, 可以得到 0Dtq=ddt(0It1-q)

选择Caputo求解分数阶微分方程的重要原因是它对常数求导结果为0, 然而Riemann-Liouville分数阶微分对常数求导结果不恒为0。

1.2 微分变换方法

函数 f(t)依据分数阶幂级数展开为:

ft=k=0Fktkq3

式中: 0< q< 1为分数阶微分的阶次; F(k)f(t)的分数阶微分变换, 表示如下:

Fk=1Γqk+10Dtqkf04

式中: k次可微Caputo分数阶微分。

近似函数 f(t)为:

ft=k=0NFktkq5

以下是Caputo分数阶微分变换的基本性质[14]:①如果f(t)=g(t)± h(t), 则F(k)=g(k)± H(k); ②如果f(t)=g(t)h(t), 则F(k)= m=0kG(m)H k-m; ③如果 f(t)=0Dtqg(t), F(k)=Γq(k+1)+1Γ(qk+1)G(k+1)

对于分数阶系统来说, 最常用的方法是用时间序列重构相空间, 计算最大Lyapunov指数[15], 这种做法仍有很多缺点, 如难以选择嵌入维数以及相空间重构的固有延迟参数[16]。本文利用FDTM得到一组分数阶微分方程的解, 用于计算分数阶系统中的Lyapunov特征指数。

引理1[17]x(t)Rn是一个可微函数, 对于任意 t0时, 存在以下关系:

0DtqxTtPxt/2xTtP0Dtqxt6q(0, 1)

式中: PRn×n为常数正定对称矩阵。

定理1[18, 19]x(t)=0为系统的一个平衡点, 且:

0Dtqxt=ft, xt7

假设存在一个Lyapunov函数 V(t, x)和一个k级函数γi(i=1, 2, 3)满足:

γ1xtVt, xγ2xt0DtqVt, x-γ3xt8

式中: q(0, 1), 则系统(7)是渐近稳定的。

2 主要结论

提出一个分数阶系统方程:

0Dtqx1=x20Dtqx2=x30Dtqx3=c1x2+c2x1x3+c3x2x3+c49

式中: 0Dtq是初始状态 t0=0系统的参数为 ci0(i=1, 2, 3, 4)阶次为 q(0, 1)的Caputo分数阶算子。

2.1 系统特性

2.1.1 平衡点

由于系统方程(9)在 c40时无平衡点, 因此不能通过Jacobi矩阵的特征值确定系统混沌阶数 q的范围。

2.1.2 分数阶Lyapunov维数和耗散性

选择参数:c1=1; c2=1; c3=1; c4=0.75; q=0.95。这个系统的差分变换为:

Γ(q(k+1)+1)Γ(qk+1)X1(k+1)=X2(k)Γ(q(k+1)+1)Γ(qk+1)X2(k+1)=X3(k)Γ(q(k+1)+1)Γ(qk+1)X3(k+1)=-X2(k)-  m=0kX1(m)X3(k-m)-  m=0kX2(m)X3(k-m)-0.7510

n阶逆变换的解为:

xt=m=0NXmt-t0yt=m=0NYmt-t0zt=m=0NZmt-t011

选取时间间隔Δ t=0.01、N=4, 混沌吸引子的相图如图1所示。

图1 混沌吸引子的相图Fig.1 Phase portraits of chaotic attractor

利用文献[20]中提到的方法, 令L1=0.012, L2=-0.0124, L3=-3.6386, L1+L2+L3=-3.6390< 0。根据Kaplan Yorke猜想, 定义系统(9)的分数阶Lyapunov维数[21]如下:

DKY=S+i=1sLiLs+1

式中: i=1sLi> 0, i=1s+1Li< 0; Li(i=1, 2, …, n)是系统(9)的Lyapunov指数; S为最大整数。显然, 如果其中一个Lyapunov指数是正定的, 并且它们的总和是负的, 则系统(9)是混沌和耗散的。另外, 该系统的Kaplan-Yorke维度是DKY=1.9677。如图2所示, 当q∈ [0.5, 1]时, 得出最大的Lyapunov指数, 且大于零, 则表示该系统是混沌的。

图2q变化的最大Lyapunov指数曲线Fig.2 Largest Lyapunov exponent curve versus q

2.2 分数阶系统的自适应控制

受控的分数阶混沌系统为:

0Dtqx1=x2+u1(x)0Dtqx2=x3+u2(x)0Dtqx3=c1x2+c2x1x3+c3x2x3+c4+u3(x)12

式中: xΤ=(x1, x2, x3)

假设:

u1(x)=K1xu2(x)=K2xu3(x)=K3x+U^3(x)13

式中: Ki=(ki, 1, ki, 2, ki, 3), i=1, 2, 3可以得出

0Dtqx1=x2+k1, 1x1+k1, 2x2+k1, 3x30Dtqx2=x3+k2, 1x1+k2, 2x2+k2, 3x30Dtqx3=c1x2+c2x1x3+c3x2x3+c4+ k3, 1x1+k3, 2x2+k3, 3x3+U^3(x)14

研究了使得系统的状态x1 、x2 、x3 达到常数r1 、r2 、r3 的控制器u1 (x)、u2 (x)、u3 (x)设计问题。

设控制误差为:

e1t=x1-r1e2t=x2-r2e3t=x3-r315

则动态误差为:

0Dtqe1t=e2t+r2+k1, 1e1t+r1+  k1, 2e2t+r2+k1, 3e3t+r30Dtqe2t=e3t+r3+k2, 1e1t+r1+  k2, 2e2t+r2+k2, 3e3t+r30Dtqe3t=c1e2t+r2+c2e1t+r1  e3t+r3+c3e3t+r3e3t+r3+  k3, 1e1t+r1+k3, 2e2t+r2+  k3, 3e3t+r3+c4+U^3x16

式中: ri(i=1, 2, 3)为任意常数。

为了简便, 假设 ri=0(i=1, 2, 3), 由式(16), 得:

0Dtqe1(t)=e2(t)+k1, 1e1(t)+k1, 2e2(t)+  k1, 3e3(t)0Dtqe2(t)=e3(t)+k2, 1e1(t)+k2, 2e2(t)+  k2, 3e3(t)0Dtqe3(t)=c1e2(t)+c2e1(t)e3(t)+  c3e2(t)e3(t)+k3, 1e1(t)+k3, 2e2(t)+  k3, 3e3(t)+c4+U^3(x)17

定理2 设计一个自适应控制器为:

u1(x)=K1xu2(x)=K2xu3(x)=K3x+U^3(x)18

式中: U^3(x)=-c^1(t)x2(t)-c^2(t)x1(t)x3(t)-c^3x(t)2(t)x3(t)-c^4(t), 选取参数自适应律为:

0Dtqc^1(t)=x2(t)x3(t)0Dtqc^2(t)=x2(t)x32(t)0Dtqc^3(t)=x2(t)x32t)0Dtqc^4(t)=x3(t)19

如果满足以下条件:

k1, 1k2, 1k2, 1k1, 2+1k2, 2k3, 2k1, 3k2, 3+1k3, 3< 0

则系统(12)是渐近稳定的。

证明 误差动态方程为:

0Dtqe1(t)=e2(t)+k1, 1e1(t)+k1, 2e2(t)+ k1, 3e3(t)0Dtqe2(t)=e3(t)+k2, 1e1(t)+k2, 2e2(t)+ k2, 3e3(t)0Dtqe3(t)=(c1-c^1(t))e2(t)+(c2-c^2(t))× e1(t)e3(t)+(c3-c^3(t))e2(t)e3(t)+ k3, 1e1(t)+k3, 2e2(t)+k3, 3e3(t)+(c4-c^4(t))20

参数估计误差定义为

ec1(t)=c1-c^1(t)ec2(t)=c2-c^2(t)ec3(t)=c3-c^3(t)ec4(t)=c4-c^4(t)21

根据分数阶的性质, 可以得到:

0Dtqec1(t)=-0Dtqc^1(t)0Dtqec2(t)=-0Dtqc^2(t)0Dtqec3(t)=-0Dtqc^3(t)0Dtqec4(t)=-0Dtqc^4(t)22

因此, 构造Lyapunov函数为:

V(e1(t), e2(t), e3(t), ec1(t), ec2(t), ec3(t), ec4(t))=12(e1(t)2+e2(t)2+e3(t)2+ec1(t)2+ec2(t)2+ec3(t)2+ec4(t)2)(23)

根据引理1, 得:

0DtqV(e1(t), e2(t), e3(t), ec1(t), ec2(t), ec3(t), ec4(t))=120Dtq(e1(t)2+e2(t)2+e3(t)2+ec1(t)2+ec2(t)2+ec3(t)2+ec4(t)2)e1(t)0Dtqe1(t)+e2(t)0Dtqe2(t)+e3(t)0Dtqe3(t)+ec1(t)0Dtqec1(t)+ec2(t)0Dtqec2(t)+ec3(t)0Dtqec3(t)+ec4(t)0Dtqec4(t)=e1(t)e2(t)e3(t)Τk1, 1k2, 1k2, 1k1, 2+1k2, 2k3, 2k1, 3k2, 3+1k3, 3e1(t)e2(t)e3(t)+ec1(t)[e2(t)e3(t)-0Dtqc^1(t)]+ec2(t)[e1(t)e32(t)-0Dtqc^2(t)]+ec3(t)[e2(t)e32(t)-0Dtqc^3(t)]+ec4(t)[e3(t)-0Dtqc^4(t)](24)

由式(19)和(20)可以得到:

0DtqV(e1(t), e2(t), e3(t), ec1(t), ec2(t), ec3(t), ec4(t))< 0

根据定理1, 可得系统(12)渐近稳定。

需要指出的是, 在本文方法中, 控制器包含所有状态信息, 它使控制器更加灵活。

3 数值模拟

对于含有不确定参数ci 的系统(12), (i=1, 2, 3, 4), q=0.95, r1=r2=r3=0, 并且:

k1, 1k2, 1k3, 1k1, 2k2, 2k3, 2k1, 3k2, 3k3, 3=-2-1-1-3-3-1-1-2-4为负定。

其中, 初始状态 (x1(0), x2(0), x3(0))=(0.1, 2, -0.1), 其自适应控制的系统误差曲线如图3所示。

图3 自适应控制器的系统的误差曲线Fig.3 Error curves of system with adaptive controller

r1=2r2=1r3=0, 其他参数不变时, 可得到图4所示的自适应控制的状态响应曲线, 图4表明状态变量完全追踪期望轨迹。

图4 自适应控制的状态响应Fig.4 States response with adaptive controller

本文设计了一种自适应控制器, 使系统状态变量追踪期望轨迹, 但追踪变量轨迹仍有待研究。

4 结束语

研究了一类无平衡点的混沌系统, 通过FDTM得到它的解集, 并分析了该系统的特点。通过判断最大Lyapunov特征指数, 给出了保持系统混沌的阶次范围。针对分数阶混沌系统设计出一种全状态自适应控制器, 可以用于追踪期望轨迹。最后, 通过数值模拟验证了本文方法的可行性。

The authors have declared that no competing interests exist.

参考文献
[1] Meral F, Royston T, Magin R. Fractional calculus in viscoelasticity: an experimental study[J]. Communications in Nonlinear Science and Numerical Simulation, 2010, 15(4): 939-945. [本文引用:1]
[2] Magin R L. Fractional calculus models of complex dynamics in biological tissues[J]. Computers & Mathematics with Applications, 2010, 59(5): 1586-1593. [本文引用:1]
[3] Drapaca C, Sivaloganathan S. A fractional model of continuum mechanics[J]. Journal of Elasticity, 2012, 107(2) : 105-123. [本文引用:1]
[4] Wiggins S. Introduction to Applied Nonlinear Dynamical Systems and Chaos[M]//Marsden J E, Sirovich L, Antman S S. Texts in Applied Mathematics. Berlin: Springer, 2003. [本文引用:1]
[5] Hoppensteadt F C. Analysis and simulation of chaotic systems[J]. Journal of Applied Mathematics and Mechanics, 2002, 82(7): 472. [本文引用:1]
[6] Li C, Peng G. Chaos in Chen's system with a fractional order[J]. Chaos, Solitons & Fractals, 2004, 22(2): 443-450. [本文引用:1]
[7] Xu Y, Gu R, Zhang H, et al. Chaos in diffusionless lorenz system with a fractional order and its control[J]. International Journal of Bifurcation and Chaos, 2012, 22(4): 1250088. [本文引用:1]
[8] Danca M F, Garrappa R. Suppressing chaos in discontinuous systems of fractional order by active control[J]. Applied Mathematics and Computation, 2015, 257: 89-102. [本文引用:1]
[9] Tacha O I, Volos C K, Stouboulos I N, et al. Analysis, adaptive control and circuit simulation of a novel finance system with dissaving[J]. Archives of Control Sciences, 2016, 26(1): 95-115. [本文引用:1]
[10] Kuntanapreeda S. Adaptive control of fractional-order unified chaotic systems using a passivity-based control approach[J]. Nonlinear Dynamics, 2016: 84(4): 2505-2515. [本文引用:1]
[11] Wang Q, Zhang J, Ding D, et al. Adaptive mittag-leffler stabilization of a class of fractional order uncertain nonlinear systems[J]. Asian Journal of Control, 2016, 18(6): 2343-2351. [本文引用:1]
[12] Arikoglu A, Ozkol I. Solution of fractional differential equations by using differential transform method[J]. Chaos, Solitons & Fractals, 2007, 34(5): 1473-1481. [本文引用:1]
[13] Podlubny I. Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications[J]. Academic Press, 1999, 91(3): 427-436. [本文引用:1]
[14] Odibat Z, Momani S. A generalized differential transform method for linear partial differential equations of fractional order[J]. Applied Mathematics Letters, 2008, 21(2): 194-199. [本文引用:1]
[15] Eckmann J P, Kamphorst S O, Ruelle D, et al. Liapunov exponents from time series[J]. Physical Review A, 1986, 34(6): 4971-4979. [本文引用:1]
[16] Takens F. Detecting Strange Attractors in Turbulence[M]//Smale S. Dynamical Systems and Turbulence. Berlin: Springer, 1981: 366-381. [本文引用:1]
[17] Duarte-Mermoud M A, Aguila-Camacho N, Gallegos J A, et al. Using general quadratic Lyapunov functions to prove lyapunov uniform stability for fractional order systems[J]. Communications in Nonlinear Science and Numerical Simulation, 2015, 22(1): 650-659. [本文引用:1]
[18] Li Y, Chen Y, Podlubny I. Mittag-leffler stability of fractional order nonlinear dynamic systems[J]. Automatica, 2009, 45(8): 1965-1969. [本文引用:1]
[19] Li Y, Chen Y, Podlubny I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized mittag-leffler stability[J]. Computers & Mathematics with Applications, 2010, 59(5): 1810-1821. [本文引用:1]
[20] Caponetto R, Fazzino S. A semi-analytical method for the computation of the lyapunov exponents of fractional-order systems[J]. Communications in Nonlinear Science and Numerical Simulation, 2013, 18(1): 22-27. [本文引用:1]
[21] Frederickson P, Kaplan J L, Yorke E D, et al. The Liapunov dimension of strange attractors[J]. Journal of Differential Equations, 1983, 49(2): 185-207. [本文引用:1]