Robust sliding mode control for fractional-order chaotic economical system with parameter uncertainty and external disturbance*
Zhou Kea), Wang Zhi-Huib)†, Gao Li-Kea), Sun Yueb), Ma Tie-Dongb)
Electric Power Research Institute, Guangxi Power Grid Corporation, Nanning 530023, China
College of Automation, Chongqing University, Chongqing 400044, China

Corresponding author. E-mail: cquwzh@163.com

*Project supported by the National Natural Science Foundation of China (Grant Nos. 51207173 and 51277192).

Abstract

This paper presents a modified sliding mode control for fractional-order chaotic economical systems with parameter uncertainty and external disturbance. By constructing the suitable sliding mode surface with fractional-order integral, the effective sliding mode controller is designed to realize the asymptotical stability of fractional-order chaotic economical systems. Comparing with the existing results, the main results in this paper are more practical and rigorous. Simulation results show the effectiveness and feasibility of the proposed sliding mode control method.

Keyword: 05.45.Xt; 05.45.Jn; 05.45.Pq; fractional-order chaotic systems; economical system; sliding mode control; parameter uncertainty
1. Introduction

Fractional calculus has a history of more than 300 years. Although it has a long history, the applications of fractional calculus to physics and engineering are only a recent focus of interest.[1, 2] Nowadays, it has been found that many systems in interdisciplinary fields, such as viscoelastic materials, [3] electrode-electrolyte polarization, [4] thermoelectric systems, [5] dielectric polarization, [6] etc., can be elegantly described by fractional differential derivatives.

In the past decade, econophysics have become an alternative scientific methodology to understand the highly complex dynamics of real financial or economical systems. These economists present a model to explain the central features of economic data (such as irregular growth, structural changes, erratic macroeconomic fluctuations, and so on). The model is often affected by exogenous shock variables and displays irregular and random behavior. Typical random shocks coming from external influence include weather variables, political events, and other human factors. It is worth mentioning that chaos can give a reasonable explanation of the above mentioned complexity.[7] Nowadays, the importance of chaos in economics has increased tremendously and many chaotic models have been used to study the complex economical systems.[811]

On the other hand, research interest has been growing in fractional calculus not only from physicists and engineers but also from researchers in life science and economics. In fact, financial variables possessing long memories make fractional models more appropriate for dynamic behaviors in finance. The magnitude of the financial variables (such as foreign exchange rates, gross domestic product, interest rates, production and stock market prices) can have a long memory. Research on fractional-order economical system has a wider range of applications. Chen[7] studied nonlinear dynamics and chaos in a fractional-order financial system. Wang et al.[12] proposed a delayed fractional-order financial system and discussed its complex dynamical behaviors by numerical simulations. Abd-Elouahab et al.[13] extended the nonlinear feedback control in ODE systems to fractional-order financial systems in order to eliminate the chaotic behavior. Wang et al.[14] addressed the synchronization problem of chaotic fractional-order economical system via active control. Chen et al.[15] studied the control and synchronization of fractional-order financial systems by the Routh– Hurwitz stability condition and linear control method.

Moreover, the control and synchronization of fractional-order dynamical systems has attracted a great deal of interest.[1622] The most common control methods include the active control method, [2325] sliding mode control method, [2629] adaptive control method, [3032] observer-based control method, [3335] impulsive control method, [3638] etc. Among these methods, sliding mode control (SMC) is an efficient method to deal with the robust control scheme. The main feature of SMC is to switch the control law to force the states of the system from the initial states onto some predefined sliding mode surface. The system on the surface has desired performance, such as stability, disturbance rejection capability, and tracking ability.

Recently, sliding mode control schemes for the fractional-order economical system are becoming a matter of concern for more and more researchers.[39, 40] In Ref.  [39], the sliding mode control just includes one single feedback controller, which is simple and easy to implement. However, the dynamics analysis of the designed sliding surface is not given, which cannot ensure desired performance. In Ref.  [40], the designed sliding surface satisfies exponential reaching law, but the full state feedback controller is necessary. Moreover, the controlled system is a fractional-order linear system with bound input, and its stability does not conform to the stability theorem for fractional-order linear system. Motivated by the aforementioned discussion, in this paper we address the modified sliding mode control scheme for a class of fractional-order chaotic economical systems. Compared with Refs.  [39] and [40], the proposed sliding mode surface meets the existence condition with desired dynamics, and the designed single controller is easy to implement to realize the robust control of the fractional-order economical system.

The rest of this paper is organized as follows. In Section 2, the problem formulation for the sliding mode control is provided. Main results for the sliding mode control of fractional-order chaotic economical system with or without parameter uncertainty and external disturbance are derived in Section 3. In Section 4, two illustrative examples for fractional-order economical system are presented to demonstrate the effectiveness and feasibility of the proposed method. Finally, concluding remarks are given in Section 5.

2. Problem formulation

In Refs.  [10] and [11], a dynamic model of finance, composed of three first-order differential equations, is presented. The model describes the time-variation of three state variables, i.e., the interest x, the investment demand y, and the price index z, which can be expressed by

where a> 0 is the saving amount, b> 0 is the cost per investment, and c> 0 is the elasticity of demand of the commercial markets. Recently, Chen generalized system  (1) to a fractional-order version[7]

where Dα (· ) denotes the Caputo fractional derivative of α ∈ (0, 1), defined as follows:[1]

where m is the first integer larger than α , and Γ (α ) is the Euler gamma function given by

By using the linear fractional-order stability theorem (system Dα X= AX is asymptotical stable if |arg(eig(A))|> α π /2 holds), it is easy to see that if y= 0, the following two-dimensional subsystem of system  (2):

is asymptotically stable about the origin x= 0 and z= 0 for all x and z (where has two eigenvalues with negative real parts such that |arg(eig(A))|> α π /2 holds, and assume that α 1= α 3= α in the following). Thus, if the proposed controller can make the state y converge to zero in finite time T, then it can further yield the asymptotical stability of subsystem  (3).

In the following section, we will propose the novel single state sliding mode controller to realize the finite time stability of state y, which is added to the second state equation,

Remark 1 The controlled system  (4) just contains one single controller u(t), which is more energy-efficient than full state feedback controller. On the other hand, the controller u(t) is added into the second equation (to realize the state y(t)= 0 in finite time) rather than the first and third ones in system   (4). This is due to the fact that the linear subsystem  (3) is asymptotically stable according to the linear fractional-order stability theorem. If the controller is added into the first or third one to make x(t)= 0 or z(t)= 0, then the subsystems

and

are not linear, and cannot satisfy the stability condition for linear fractional-order system.

3. Main results

The main idea of sliding mode control is to design a discontinuous control to force the system state trajectories to some predefined sliding mode surface, which is composed of two steps: 1) constructing a sliding mode surface that has desired properties; 2) developing a switch control law to make the trajectories reach the surface and stay on it evermore.

We choose the following nonsingular fractional-order sliding mode surface as:

where sign(· ) is the sign function, k> 0, and γ ∈ (0, 1) are constants.

Once the system trajectory reaches the sliding surface  (7), it yields the following equations:

From Eqs.    (7) and  (8), the dynamics of sliding mode can be achieved as follows:

it further yields

Before we give the main results, the fractional-order Lyapunov stability theorem can be stated as follows.

Lemma 1[41] Let x= 0 be an equilibrium point for the nonautonomous fractional-order system,

where f(t, x) satisfies the Lipschitz condition with Lipschitz constant l> 0 and α ∈ (0, 1). Assume that there exists a Lyapunov function V(t, x) satisfying

where β 1, β 2, β 3, and β are positive constants. Then the equilibrium point of the system  (11) is Miattag– Leffler (asymptotic) stable.

In the following theorem, the finite-time stability of the sliding surface  (7) is proved.

Theorem 1 The sliding mode dynamic  (10) is asymptotic stable and its trajectory converges to the equilibrium y= 0 in a finite time.

Proof Considering the following Lyapunov function:

and its derivative along the trajectory of system  (10) is

Since |y| and |y|γ are positive, it yields that

From Lemma  1, the Lyapunov function V(t) satisfies condition  (12) (if let β 1= 0.5, β 2= 2, and β = 1) and condition  (13) (with β 3= k> 0), thus the system state y(t) will converge to zero asymptotically. In the following, the finite-time convergence of y(t) will be further investigated.

From Eq.  (15), we have

It is obvious that

By taking the integral of both sides of Eq.  (18) from t1 to t2 (t1 is the reaching time for state y(t) to the sliding surface s(t) and t2 is the convergence time for y(t) to zero), and considering y(t2)= 0, thus we have

Therefore, the state trajectory y(t) can converge to zero in finite time

We complete the proof.

After constructing a fractional-order sliding mode surface, the second step is to design a nonlinear feedback controller u(t) to force the error state y(t) to reach the sliding surface and stay on it evermore. In what follows, we propose the suitable controller to ensure the existence of the sliding motion.

Theorem 2 If system  (2) is controlled with the following controller:

where k> 0, ξ > 0, and γ ∈ (0, 1) are constants, then the system state y(t) will converge to the sliding surface s(t)= 0 in finite time.

Proof Selecting a Lyapunov function in the form of V(t)= |s(t)| and the time derivative of V(t) is obtained as follows:

From systems   (4),   (21), and (22), we have

Obviously, the system state trajectories will converge to the sliding surface s(t)= 0 asymptotically. Similar to the proof of Theorem  1, we can also obtain the results that the sliding motion occurs in finite time. Referring to Eqs.  (17)– (19), the reaching time t1 can be derived as

where t0 is the initial time. The detailed proof is similar to that of Theorem 1, which is omitted here. We complete the proof.

Remark 2 From Theorems 1 and 2, we can obtain that the state trajectories of system  (2) will first converge to s(t)= 0 in finite time t1, and then the state y(t) will converge to zero in finite time t2. Based on the above discussion, once y(t) converges to zero in finite time, the subsystem  (3) will be asymptotically stable about the origin x= 0 and z= 0.

Remark 3 In Ref.  [39], the same control scheme for fractional-order economical system  (2) is presented. However, when the states are driven to the designed sliding surface

it satisfies s= 0 and ṡ = 0, i.e.,

Obviously, typical fractional-order stability theorem does not apply to nonlinear differential equation  (26), and the given sliding surface  (25) does not have a desired system dynamics. By comparison, in this paper, the proposed sliding surface  (7) can satisfy the finite time stability of system state y(t), which can further obtain the asymptotical stability of system states x(t) and z(t). Therefore, the main results of this paper are more reasonable and rigorous than the existing results.[39]

Considering that system  (4) is perturbed by parameter uncertainty and external disturbance, which comes from the introduction of the single controller u(t), the systems can be expressed as follows:

where Δ f(x, y, z) and d(t) are parameter uncertainty and external disturbance, respectively, which are assumed to be bounded as |Δ f(x, y, z)|< η 1 and |d(t)|< η 21,   η 2> 0). The following theorem will give sufficient conditions of finite time convergence to sliding surface s(t) for system  (27) with parameter uncertainty and external disturbance.

Theorem 3 If system  (27) is controlled with the following controller:

where k> 0, ξ > 0, and γ ∈ (0, 1) are constants, s(t) has the same form as (7), then the system state y(t) will converge to the sliding surface s(t)= 0 in finite time.

Proof Selecting a Lyapunov function in the form of V(t)= |s(t)| and the time derivative of V(t) is obtained as follows:

From Eqs.  (27)– (29), we obtain

where the inequalities sign(s)Δ f(x, y, z)-η 1 ≤ |Δ f(x, y, z)|-η 1≤ 0 and sign(s)d(t)-η 2≤ |d(t)|-η 2≤ 0 are used. Obviously, the system state trajectories will converge to the sliding surface s(t)= 0 asymptotically. Similar to the proof of Theorem  1, we can also obtain the results that the sliding motion occurs in finite time. Referring to Eqs.  (17)– (19), the reaching time t1 can be derived as

where t0 is the initial time. The detailed proof is similar to that of Theorem  1, which is omitted here. We complete the proof.

Once the states of system  (27) reach the sliding surface  (7), the desired convergence characteristic (state y(t) converges to zero in finite time) will be subsequently achieved, and then the asymptotical stability of states x(t) and z(t) can be obtained by the designed sliding surface  (7).

4. Simulation results

In the simulation, two illustrative examples are given to verify and demonstrate the effectiveness of the proposed method, and the time-domain numerical algorithm based on predictor– corrector scheme[39] is used for simulations.

4.1. Control of nominal fractional-order economical system

According to Ref.  [39], the system is chaotic when α 1= α 2= α 3= 0.9 and (a, b, c)= (1, 0.1, 1). The chaotic trajectories of the fractional-order economical system without the controller u(t) are shown in Fig.  1. By Theorems  1 and 2, the time response of system   (4) and the sliding surface are shown in Figs.  2 and 3, respectively. Throughout the simulation, the initial condition is taken as [1, 1, − 1]T, and the control parameters in Eqs.  (7) and (21) are chosen as k= 0.1, γ = 0.95, and ξ = 0.5.

Fig.  1. The chaotic trajectories of the fractional-order economical system  (2).

4.2. Control of fractional-order economical system with parameter uncertainty and external disturbance

Let the parameter uncertainty be Δ f(x, y, z)= 0.5sin(π x)cos(π y)sin(2π z) and the external disturbance be d(t)= 0.1cos(2t), where |Δ f(x, y, z)|≤ η 1= 0.5 and |d(t)|≤ η 2= 0.1. Based on Theorems  1 and 3, the time response of system  (27) and the sliding surface are shown in Figs.  4 and 5, respectively. Throughout the simulation, the initial condition is taken as [1, 1, − 1]T, and the control parameters in Eqs.  (7) and (28) are chosen as k= 0.1, γ = 0.95, and ξ = 0.5. Note that the introduction of parameter uncertainty and external disturbance may change the chaotic behavior of the original fractional-order chaotic economical system, and this change cannot affect the control performance. From the simulation results (see Fig.  4), the proposed controller  (28) can still realize the robust control of the original fractional-order economical system.

Fig.  2. Time response of system  (4). (a) x(t), (b) y(t), (c) z(t).

Fig.  3. Time response of sliding mode surface  (7).

Fig.  4. Time response of system  (27). (a) x(t), (b) y(t), (c) z(t).

Fig.  5. Time response of sliding mode surface  (7).

5. Conclusions

In this paper, we present a modified sliding mode control method to stabilize typical fractional-order chaotic economical systems. By constructing an efficient sliding mode surface and control law, some effective and practical sufficient conditions are presented to guarantee the system dynamics with or without uncertainty and external disturbance can converge to the origin asymptotically. Finally, two numerical examples for nominal and uncertain fractional-order economical systems are given to demonstrate the effectiveness of the method.

Reference
1 Podlubny I 1999 Fractional Differential Equations New York Academic [Cited within:2]
2 Hilfer R 2001 Applications of Fractional Calculus in Physics Singapore World Scientific [Cited within:1]
3 Hu K X and Zhu K Q 2009 Chin. Phys. Lett. 26 108301 DOI:10.1088/0256-307X/26/10/108301 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541]
4 Ichise M, Nagayanagi Y and Kojima T 1971 J. Electroanal. Chem. Interfacial Electrochem. 33 253 DOI:10.1016/S0022-0728(71)80115-8 [Cited within:1]
5 Ezzat M A 2011 Appl. Math. Modelling 35 4965 DOI:10.1016/j.apm.2011.04.004 [Cited within:1] [JCR: 1.706]
6 Sun H, Abdelwahab A and Onaral B 1984 IEEE Trans. Autom. Control 29 441 DOI:10.1109/TAC.1984.1103551 [Cited within:1]
7 Chen W C 2008 Chaos Solitons Fract. 36 1305 DOI:10.1016/j.chaos.2006.07.051 [Cited within:3] [JCR: 1.246]
8 Chian A L, Rempel E L and Rogers C 2006 Chaos Soliton. Fract. 29 1194 DOI:10.1016/j.chaos.2005.08.218 [Cited within:1] [JCR: 1.246]
9 Cesare L D and Sportelli M 2005 Chaos Soliton. Fract. 25 233 DOI:10.1016/j.chaos.2004.11.044 [Cited within:1] [JCR: 1.246]
10 Ma J H and Chen Y S 2001 Appl. Math. Mech. 22 1240 http://118.145.16.217/magsci/article/article?id=11472009 [Cited within:1] [JCR: 0.647]
11 Ma J H and Chen Y S 2001 Appl. Math. Mech. 22 1375 http://118.145.16.217/magsci/article/article?id=11472066 [Cited within:2] [JCR: 0.647]
12 Wang Z, Huang X and Shi G D 2011 Comput. Math. Appl. 62 1531 DOI:10.1016/j.camwa.2011.04.057 [Cited within:1] [JCR: 0.75]
13 Abd-Elouahab M S, Hamri N E and Wang J W 2010 Math. Probl. Eng. 2010 270646 DOI:10.1155/2010/270646 [Cited within:1] [JCR: 1.383]
14 Wang Z and Huang X 2011 Procedia Engineering 15 516 DOI:10.1016/j.proeng.2011.08.098 [Cited within:1]
15 Chen L P, Chai Y and Wu R C 2011 Discrete Dyn. Nat. Soc. 2011 958393 DOI:10.1155/2011/958393 [Cited within:1] [JCR: 0.82]
16 Wang X Y and Zhang Y L 2011 Chin. Phys. B 20 100506 DOI:10.1088/1674-1056/20/10/100506 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
17 Wang M J, Wang X Y and Niu Y J 2011 Chin. Phys. B 20 010508 DOI:10.1088/1674-1056/20/1/010508 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
18 Fu J, Yu M and Ma T D 2011 Chin. Phys. B 20 120508 DOI:10.1088/1674-1056/20/12/120508 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
19 Wang X Y, Zhao G B and Yang Y H 2013 Int. J. Mod. Phys. B 27 1350034 DOI:10.1142/S0217979213500343 [Cited within:1] [JCR: 0.358]
20 Luo C and Wang X Y 2013 Nonlinear Dyn. 71 241 DOI:10.1007/s11071-012-0656-z [Cited within:1]
21 Guan Z H, Huang F J and Guan W J 2005 Phys. Lett. A 346 153 DOI:10.1016/j.physleta.2005.08.006 [Cited within:1] [JCR: 1.11]
22 Li C G and Chen G R 2004 Chaos Soliton. Fract. 22 549 DOI:10.1016/j.chaos.2004.02.035 [Cited within:1] [JCR: 1.246]
23 Bhalekar S and Daftardar-Gejji V 2010 Commun. Nonlinear Sci. Numer. Simulat. 15 3536 DOI:10.1016/j.cnsns.2009.12.016 [Cited within:1]
24 Taghvafard H and Erjaee G H 2011 Commun. Nonlinear Sci. Numer. Simulat. 16 4079 DOI:10.1016/j.cnsns.2011.02.015 [Cited within:1]
25 Wang X Y, Zhang X P and Ma C 2012 Nonlinear Dyn. 69 511 DOI:10.1007/s11071-011-0282-1 [Cited within:1]
26 Yin C, Dadras S, Zhong S M and Chen Y Q 2013 Appl. Math. Model. 37 2469 DOI:10.1016/j.apm.2012.06.002 [Cited within:1] [JCR: 1.706]
27 Sun N, Zhang H G and Wang Z L 2011 Acta Phys. Sin. 60 050511(in Chinese) http://118.145.16.217/magsci/article/article?id=17393125 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
28 Liu D and Yan X M 2009 Acta Phys. Sin. 58 3747(in Chinese) http://118.145.16.217/magsci/article/article?id=17390389 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
29 Liu F C, Li J Y and Zang X F 2011 Acta Phys. Sin. 60 141012(in Chinese) http://118.145.16.217/magsci/article/article?id=17392858 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
30 Zhao L D, Hu J B and Liu X H 2010 Acta Phys. Sin. 59 2305(in Chinese) http://118.145.16.217/magsci/article/article?id=17391582 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
31 Zhang R X and Yang S P 2010 Chin. Phys. B 19 020510 DOI:10.1088/1674-1056/19/2/020510 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
32 Ma T D, Jiang W B, Fu J, Chai Y, Chen L P and Xue F Z 2012 Acta Phys. Sin. 61 160506(in Chinese) http://118.145.16.217/magsci/article/article?id=17394965 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
33 Senejohnny D M and Delavari H 2012 Commun. Nonlinear Sci. Numer. Simulat. 17 4373 DOI:10.1016/j.cnsns.2012.03.004 [Cited within:1]
34 Liu L, Liang D L, Liu C X and Zhang Q 2012 Int. J. Mod. Phys. B 26 1250166 DOI:10.1142/S0217979212501664 [Cited within:1] [JCR: 0.358]
35 Wang X Y and Hu Z W 2012 Int. J. Mod. Phys. B 26 1250176 DOI:10.1142/S0217979212501767 [Cited within:1] [JCR: 0.358]
36 Wang X Y, Zhang Y L, Li D and Zhang N 2011 Chin. Phys. B 20 030506 DOI:10.1088/1674-1056/20/3/030506 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
37 Ma T D, Jiang W B and Fu J 2012 Acta Phys. Sin. 61 090503(in Chinese) http://118.145.16.217/magsci/article/article?id=17394284 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
38 Ma T D, Jiang W B, Fu J and Xue F Z 2012 Acta Phys. Sin. 61 100507(in Chinese) http://118.145.16.217/magsci/article/article?id=17394445 [Cited within:1] [JCR: 1.016] [CJCR: 1.691]
39 Dadras S and Momeni H R 2010 Physica A 389 2434 DOI:10.1016/j.physa.2010.02.025 [Cited within:7] [JCR: 1.676]
40 Wang Z, Huang X and Shen H 2012 Neurocomputing 83 83 DOI:10.1016/j.neucom.2011.11.018 [Cited within:3] [JCR: 1.634]
41 Li Y, Chen Y Q and Podlubny I 2009 Automatica 45 1965 DOI:10.1016/j.automatica.2009.04.003 [Cited within:1] [JCR: 2.919]