Next Article in Journal
A New Scheme Using Cubic B-Spline to Solve Non-Linear Differential Equations Arising in Visco-Elastic Flows and Hydrodynamic Stability Problems
Previous Article in Journal
A New Optimal Family of Schröder’s Method for Multiple Zeros
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis, Circuit Simulation, and Control Synchronization of Fractional-Order Unified Chaotic System

School of Electronic Engineering, Xi’an University of Posts and Telecommunications, Xi’an, Shaanxi 710121, China
*
Authors to whom correspondence should be addressed.
Mathematics 2019, 7(11), 1077; https://doi.org/10.3390/math7111077
Submission received: 20 October 2019 / Revised: 3 November 2019 / Accepted: 5 November 2019 / Published: 8 November 2019
(This article belongs to the Section Engineering Mathematics)

Abstract

:
The traditional method of solving fractional chaotic system has the problem of low precision and is computationally cumbersome. In this paper, different fractional-order calculus solutions, the Adams prediction–correction method, the Adomian decomposition method and the improved Adomian decomposition method, are applied to the numerical analysis of the fractional-order unified chaotic system. The result shows that different methods have higher precision, smaller computational complexity, and shorter running time, in which the improved Adomian decomposition method works best. Then, based on the fractional-order chaotic circuit design theory, the circuit diagram of fractional-order unified chaotic system is designed. The result shows that the circuit simulation diagram of fractional-order unified chaotic system is basically consistent with the phase space diagram obtained from the numerical solution of the system, which verifies the existence of the fractional-order unified chaotic system of 0.9-order. Finally, the active control method is used to control and synchronize in the fractional-order unified chaotic system, and the experiment result shows that the method can achieve synchronization in a shorter time and has a better control performance.

1. Introduction

Since the definition of the Lorenz system by Lorenz in 1963, the study of modern chaos theory has made rapid progress. The typical ones are Chen system [1], Lü system [2], Rössler system [3], generalized Lü system [4], etc. In 2003, Liu et al. [5] discovered a chaos system that can produce four helical chaotic attractors in a three-dimensional chaos system. In 2004, Liu et al. [6] discovered a class of three-dimensional chaotic system with square terms. In 2005, Qi et al. [7] added a nonlinear term on the basis of Lorenz system and obtained a new chaotic system. In 2009, Bao et al. [8] proposed a three-dimensional continuous autonomous chaotic system with simple equations and complex attractor topologies.
The dynamic characteristics of fractional-order chaotic system is related to the corresponding order. In the same chaotic system, when the selection of order is different, the corresponding chaotic characteristics, such as bifurcation graph and Lyapunov exponent graph, are also different, so they have more abundant chaotic characteristics. Scholars at home and abroad have carried out this research. Chen et al. [9], Wang et al. [10] studied the chaotic circuit of fractional-order Liu system and other fractional-order chaotic system derived from fractional-order Liu system. Mao et al. [11], He et al. [12] studied fractional-order Lorenz system and circuit simulation. At the same time, the control and synchronization of fractional-order chaotic system also has made progress. Wu [13] proposed a new control method based on linear matrix inequality to solve the synchronization of two fractional-order chaotic system. Yan et al. [14] proposed a new fractional-order sliding mode controller to solve the synchronization problem of uncertain fractional-order chaotic system with modeling errors and external disturbances. Rabah et al. [15] designed a new robust fractional-order proportional integral-differential (FOPID) controller. Rajagopal et al. [16] derived and obtained a four-dimensional fractional-order discrete power model. Although some achievements have been made in these studies, the traditional fractional-order chaotic system solution methods for example calling function method have some problems, such as low accuracy and computationally cumbersome.
Based on these problems, in this paper, different methods are used to solve the numerical solution of the fractional-order unified chaotic system, and the phase space diagram of the corresponding order is obtained. Then, the circuit simulation diagram of the chaotic system is designed. Finally, three control laws are obtained by the active control method, and the control and synchronization of the fractional-order unified chaotic system are completed.

2. Intege-Order Unified Chaotic System

The mathematical model of the three-dimensional continuous chaotic system [17] is as follows:
{ d x d t = ( 25 α + 10 ) ( y x ) d y d t = ( 28 35 α ) x x z + ( 29 α 1 ) y d z d t   = x y ( 8 + α ) z 3
Since the chaotic system can be transformed among Lorenz system, Chen system, and Lü system when α takes different value, the system is called unified chaotic system. When α ∈ [0,1], the system has a global chaotic characteristic. The attractor diagram for α with different values is shown in Figure 1. Time domain waveform is shown in Figure 2. Since the parameter α takes different value and the system corresponds to a different three-dimensional continuous chaotic system. When α ∈ [0.8,1], the Lyapunov exponent diagram is listed, as shown in Figure 3. The maximum Lyapunov exponent is greater than zero. It indicates that the system is in a chaotic state.

3. Fractional-Order Unified Chaotic System

3.1. Common Definitions of Fractional Calculus

(1) The Grünwald–Letnikov fractional differential is defined as [18]:
D t q a f ( t ) = d q f ( t ) d ( t a ) q = l i m N [ t a N ] q j = 0 N 1 ( 1 ) j ( q j ) f ( t j [ t a N ] )
(2) The Riemann–Liouville fractional differential is defined as [18]:
D t q a f ( t ) = { 1 Γ ( q ) a t ( t τ ) q 1 f ( τ ) d ( τ ) , q < 0 ( t ) , q = 0 D n [ D q n a f ( t ) ] , q > 0
where Γ(.) is a Gamma function, which is defined as:
Γ ( x ) = 0 t x 1 e t d t
(3) The Caputo fractional differential is defined as [18]:
D t q a f ( t ) = { 1 Γ ( m q ) 0 t f ( m ) ( τ ) d τ ( t τ ) q + 1 m , m 1 < q < m d m d t m f ( t ) , q = m
where the q-order differential and the power function respectively are:
D t 0 q t r = Γ ( r + 1 ) Γ ( r + 1 q ) ( t t 0 ) r q
D t 0 q C = 0

3.2. The Calling Function Method

Redefine the mathematical model of the fractional-order unified chaotic system:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) D q 2 x 2 = ( 28 35 α ) x 1 x 1 x 3 + ( 29 α 1 ) x 2 D q 3 x 3 = x 1 x 2 ( 8 + α ) x 3 / 3
When α = 0, construct the MATLAB function for the fractional-order unified chaotic system: Function [T Y] = FOLorenz (parameters, orders, TSim, Y0), where parameters are system parameters, orders are derivative order, TSim is simulation time, and Y0 is initial condition.
When q = 0.9, TSim = 100, Y0 = [0.1 0.1 0.1], the phase space diagram is shown in Figure 4.

3.3. The Adams–Bashforth–Moulton Prediction–Correction Method [19]

For the following differential equations:
{ d q x d t q = f ( t , x ) , 0 t T x k ( 0 ) = x 0 ( k ) , k = 0 , 1 , 2 , 3 , , q 1
Equivalent to the Volterra integral equation:
x ( t ) = k = 0 n 1 x 0 ( k ) t k k ! + 1 Γ ( q ) 0 t ( t τ ) q 1 f ( τ , x ( t ) ) d τ
The correction formula is:
x h ( t n + 1 ) = k = 0 q 1 x 0 ( k ) t n + 1 k k ! + h q Γ ( q + 2 ) f ( t n + 1 , x h p t ( n + 1 ) ) +   h q Γ ( q + 2 ) j = 0 n a j , n + 1 f ( t j , x h ( t , j ) )
where,
a j , n + 1 = { n q + 1 ( n q ) ( n + 1 ) q   , j = 0 ( n j 2 ) q + 1 + ( n j ) q + 1 2 ( n j + 1 ) q + 1 , 1 j n
The prediction formula is:
x h p ( t n + 1 ) = k = 0 n 1 x 0 ( k ) t n + 1 k k ! + 1 Γ ( q ) j = 0 n b j , n + 1 f ( t j , x h ( t j ) )
where,
b j , n + 1 = h q q ( ( n j + 1 ) q ( n j ) q ) ,   0 j n
The basic algorithm of Adams–Bashforth–Moulton prediction–correction method is shown in Equations (11) and (13), where the coefficients are shown in Equations (12) and (14). Although the calculation of the fractional time domain analysis method is more complicated than the frequency domain method, the order of the calculus operator can calculate the step change smaller, which is more conducive to analyze the dynamics of the fractional system.
Now redefine the mathematical model of the fractional-order unified chaotic system:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) D q 2 x 2 = ( 28 35 α ) x 1 x 1 x 3 + ( 29 α 1 ) x 2 D q 3 x 3 = x 1 x 2 ( 8 + α ) x 3 / 3  
The discretization of fractional-order unified chaotic system can be obtained:
{ x n + 1 = x 0 + ( 25 α + 10 ) h q Γ ( q + 2 ) ( y n + 1 x n + 1 + j = 0 n a j , n + 1 ( y j x j ) ) y n + 1 = y 0 + h q Γ ( q + 2 ) ( ( ( 28 35 α ) x n + 1 x n + 1 z n + 1 + ( 29 α 1 ) y n + 1 ) + j = 0 n a j , n + 1 ( ( 28 35 α ) x j x j z j + ( 29 α 1 ) y j ) ) z n + 1 = z 0 + h q Γ ( q + 2 ) ( x n + 1 y n + 1 ( 8 + α ) z n + 1 / 3 + j = 0 n a j , n + 1 ( x j y j ( 8 + α ) z j / 3 ) )
The steps to solve fractional differential equations by the Adams prediction–correction method are as follows:
Step 1:
Determine the step size and the interval of solution. Taking the fractional-order unified chaotic system as an example, the step size h = 0.01 is selected. The number of points is selected as N = 10,000, and the upper limit of integration is T = 100 (h = T/N).
Step 2:
Calculate the predictive factor.
Step 3:
Calculate the predictive value.
Step 4:
Calculate the correction factor.
Step 5:
Calculate the correction value.
Step 6:
Determine if the end condition is meet.
Step 7:
End the calculation if the condition is meet, or return to step 3.
The flow of numerical calculation of the Adams prediction–correction method is shown in Figure 5.
When q1 = q2 = q3 = 0.9, step length h = 0.001, and α = 0.8, the phase space diagram of each plane is shown in Figure 6. When q1 = q2 = q3 = 0.9, step length h = 0.001, and α = 1, the phase space diagram of each plane is shown in Figure 7.
In addition, according to the Adams prediction–correction method, the phase space of the fractional-order Duffing oscillator when q1 = q2 = q3 = 0.98 is solved, as shown in Figure 8. As a result, the Adams prediction–correction method is applicable to not only fractional-order unified chaotic system, but also other fractional-order chaotic system, such as fractional-order Duffing system.

3.4. The Adomian Decomposition Method

The Adomian decomposition method [20] is a time-domain approximation algorithm that is suitable for the solution of fractional-order system. It is a time-domain solution algorithm that can deal with linear and non-linear problems. It is not necessary to carry out the discretization and a large amount of computer memory, and provide a high-precision approximate numerical solution.
The mathematical model of the fractional-order unified chaotic system is:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) D q 2 x 2 = ( 28 35 α ) x 1 x 1 x 3 + ( 29 α 1 ) x 2 D q 3 x 3   = x 1 x 2 ( 8 + α ) x 3 / 3
By decomposing the linear and nonlinear terms in the system, we can obtain:
[ L x 1 L x 2 L x 3 ] = [ ( 25 α + 10 ) ( x 2 x 1 ) ( 28 35 α ) x 1 + ( 29 α 1 ) x 2 ( 8 + α ) 3 x 3 ]
[ N x 1 N x 2 N x 3 ] = [ 0 x 1 x 3 x 1 x 2 ]
The nonlinear terms are decomposed and the previous seven items are selected. The following equation can be obtained:
{ A 2 0 = x 1 0 x 3 0 A 2 1 = x 1 1 x 3 0 x 1 0 x 3 1 A 2 2 = x 1 2 x 3 0 x 1 1 x 3 1 x 1 0 x 3 2 A 2 3 = x 1 3 x 3 0 x 1 2 x 3 1 x 1 1 x 3 2 x 1 0 x 3 3 A 2 4 = x 1 4 x 3 0 x 1 3 x 3 1 x 1 2 x 3 2 x 1 1 x 3 3 x 1 0 x 3 4 A 2 5 = x 1 5 x 3 0 x 1 4 x 3 1 x 1 3 x 3 2 x 1 2 x 3 3 x 1 1 x 3 4 x 1 0 x 3   5
{ A 3 0 = x 1 0 x 2 0 A 3 1 = x 1 1 x 2 0 + x 1 0 x 2 1 A 3 2 = x 1 2 x 2 0 + x 1 1 x 2 1 + x 1 0 x 2 2 A 3 3 = x 1 3 x 2 0 + x 1 2 x 2 1 + x 1 1 x 2 2 + x 1 0 x 2 3 A 3 4 = x 1 4 x 2 0 + x 1 3 x 2 1 + x 1 2 x 2 2 + x 1 1 x 2 3 + x 1 0 x 2 4 A 3 5 = x 1 5 x 2 0 + x 1 4 x 2 1 + x 1 3 x 2 2 + x 1 2 x 2 3 + x 1 1 x 2 4 + x 1 0 x 2 5
According to the given conditions, the following equation can be obtained:
{ x 1 0 = x 1 ( t 0 ) x 2 0 = x 2 ( t 0 ) x 3 0 = x 3 ( t 0 )
Make c 1 0 = x 1 0 , c 2 0 = x 2 0 , c 3 0 = x 3 0 , the following equation can be obtained:
{ x 1 1 = [ ( 25 α + 10 ) ( c 2 0 c 1 0 ) ] ( t t 0 ) q Γ ( q + 1 ) x 2 1 = [ ( 28 35 α ) c 1 0 c 1 0 c 3 0 + ( 29 α 1 ) c 2 0 ] ( t t 0 ) q Γ ( q + 1 ) x 3 1 = [ c 1 0 c 2 0 ( 8 + α ) c 3 0 / 3 ] ( t t 0 ) q Γ ( q + 1 )
From the iteration relation, the other five coefficients of x can be deduced as follows:
{ c 1 2 = ( 25 α + 10 ) ( c 2 1 c 1 1 ) c 2 2 = ( 28 35 α ) c 1 1 + ( 29 α 1 ) c 2 1 c 1 1 c 3 0 c 1 0 c 3 1 c 3 2 = [ 8 + α 3 ] c 3 1 + c 1 1 c 2 0 + c 1 0 c 2 1
{ c 1 3 = ( 25 α + 10 ) ( c 2 2 c 1 2 ) c 2 3 = ( 28 35 α ) c 1 2 + ( 29 α 1 ) c 2 2 c 1 2 c 3 0 c 1 1 c 3 1 Γ ( 2 q + 1 ) Γ 2 ( q + 1 ) c 1 0 c 3 2 c 3 3 = [ 8 + α 3 ] c 3 2 + c 1 2 c 2 0 + c 1 1 c 2 1 Γ ( 2 q + 1 ) Γ 2 ( q + 1 ) + c 1 0 c 2 2
{ c 1 4 = ( 25 α + 10 ) ( c 2 3 c 1 3 ) c 2 4 = ( 28 35 α ) c 1 3 + ( 29 α 1 ) c 2 3 c 1 3 c 3 0 ( c 1 1 c 3 1 + c 1 1 c 3 2 ) Γ ( 3 q + 1 ) Γ ( q + 1 ) Γ ( 2 q + 1 ) c 1 0 c 3 3 c 3 4 = [ 8 + α 3 ] c 3 3 + c 1 3 c 2 0 + ( c 1 2 c 2 1 + c 1 1 c 2 2 ) Γ ( 3 q + 1 ) Γ ( q + 1 ) Γ ( 2 q + 1 ) + c 1 0 c 2 3
{ c 1 5 = ( 25 α + 10 ) ( c 2 4 c 1 4 )   c 2 5 = ( 28 35 α ) c 1 4 + ( 29 α 1 ) c 2 4 c 1 4 c 3 0 ( c 1 3 c 3 1 + c 1 1 c 3 3 ) Γ ( 4 q + 1 ) Γ ( q + 1 ) Γ ( 3 q + 1 ) c 1 2 c 3 2 Γ ( 4 q + 1 ) Γ 2 ( 2 q + 1 ) c 1 0 c 3 4 c 3 5 = [ 8 + α 3 ] c 3 4 + c 1 4 c 2 0 + ( c 1 3 c 2 1 + c 1 1 c 2 3 ) Γ ( 4 q + 1 ) Γ ( q + 1 ) Γ ( 3 q + 1 ) + c 1 2 c 2 2 Γ ( 4 q + 1 ) Γ 2 ( 2 q + 1 ) + c 1 0 c 2 4
{ c 1 6 = ( 25 α + 10 ) ( c 2 5 c 1 5 ) c 2 6 = ( 28 35 α ) c 1 5 + ( 29 α 1 ) c 2 5 c 1 5 c 3 0 ( c 1 4 c 3 1 + c 1 1 c 3 4 ) Γ ( 5 q + 1 ) Γ ( q + 1 ) Γ ( 4 q + 1 )   ( c 1 3 c 3 2 + c 1 2 c 3 3 )   Γ ( 5 q + 1 ) Γ ( 2 q + 1 ) Γ ( 3 q + 1 ) c 1 0 c 3 5   c 3 6 = [ 8 + α 3 ] c 3 5 + c 1 5 c 2 0 + ( c 1 4 c 2 1 + c 1 1 c 2 4 ) Γ ( 5 q + 1 ) Γ ( q + 1 ) Γ ( 4 q + 1 ) + ( c 1 3 c 2 2 + c 1 2 c 2 3 ) Γ ( 5 q + 1 ) Γ ( 2 q + 1 ) Γ ( 3 q + 1 ) + c 1 0 c 2 5
The solution to the system equation can be expressed as:
x j ( t ) = c j 0 + c j 1 ( t t 0 ) q Γ ( q + 1 ) + c j 2 ( t t 0 ) 2 q Γ ( 2 q + 1 ) + c j 3 ( t t 0 ) 3 q Γ ( 3 q + 1 ) + c j 4 ( t t 0 ) 4 q Γ ( 4 q + 1 ) + c j 5 ( t t 0 ) 5 q Γ ( 5 q + 1 ) + c j 6 ( t t 0 ) 6 q Γ ( 6 q + 1 )
where j = 1, 2, 3. Thus, when q1 = q2 = q3 = 0.96, α = 0, h = 0.001, the phase space diagram of each plane is shown in Figure 9.

3.5. The Improved Adomian Decomposition Method

The algorithm is further optimized according to the Adomian decomposition method [21].
The fractional-order unified chaotic system is taken as an example to decompose the nonlinear part. According to the given conditions, the following equation can be obtained:
{ x 1 0 = x 1 ( t 0 ) x 2 0 = x 2 ( t 0 ) x 3 0 = x 3 ( t 0 )
Make c 1 0 = x 1 0 , c 2 0 = x 2 0 , c 3 0 = x 3 0 , the following equation can be obtained:
{ x 1 1 = [ ( 25 α + 10 ) ( c 2 0 c 1 0 ) ] ( t t 0 ) q q x 2 1 = [ ( 28 35 α ) c 1 0 c 1 0 c 3 0 + ( 29 α 1 ) c 2 0 ] ( t t 0 ) q q x 3 1 = [ c 1 0 c 2 0 ( 8 + α ) c 3 0 / 3 ] ( t t 0 ) q q
From the iteration relation, the other five coefficients of x can be deduced as follows:
{ c 1 2 = ( 25 α + 10 ) ( c 2 1 c 1 1 ) c 2 2 = ( 28 35 α ) c 1 1 + ( 29 α 1 ) c 2 1 c 1 1 c 3 0 c 1 0 c 3 1 c 3 2 = [ 8 + α 3 ] c 3 1 + c 1 1 c 2 0 + c 1 0 c 2 1
{ c 1 3 = ( 25 α + 10 ) ( c 2 2 c 1 2 ) c 2 3 = ( 28 35 α ) c 1 2 + ( 29 α 1 ) c 2 2 c 1 2 c 3 0 2 c 1 1 c 3 1 c 1 0 c 3 2 c 3 3 = [ 8 + α 3 ] c 3 2 + c 1 2 c 2 0 + 2 c 1 1 c 2 1 + c 1 0 c 2 2
{ c 1 4 = ( 25 α + 10 ) ( c 2 3 c 1 3 ) c 2 4 = ( 28 35 α ) c 1 3 + ( 29 α 1 ) c 2 3 c 1 3 c 3 0 3 ( c 1 2 c 3 1 + c 1 1 c 3 2 ) c 1 0 c 3 3 c 3 4 = [ 8 + α 3 ] c 3 3 + c 1 3 c 2 0 + 3 ( c 1 2 c 2 1 + c 1 1 c 2 2 ) + c 1 0 c 2 3
{ c 1 5 = ( 25 α + 10 ) ( c 2 4 c 1 4 ) c 2 5 = ( 28 35 α ) c 1 4 + ( 29 α 1 ) c 2 4 c 1 4 c 3 0 4 ( c 1 3 c 3 1 + c 1 1 c 3 3 ) 6 c 1 2 c 3 2 c 1 0 c 3 4 c 3 5 = [ 8 + α 3 ] c 3 4 + c 1 4 c 2 0 + 4 ( c 1 3 c 2 1 + c 1 1 c 2 3 ) + 6 c 1 2 c 2 2 + c 1 0 c 2 4
{ c 1 6 = ( 25 α + 10 ) ( c 2 5 c 1 5 ) c 2 6 = ( 28 35 α ) c 1 5 + ( 29 α 1 ) c 2 5 c 1 5 c 3 0 5 ( c 1 4 c 3 1 + c 1 1 c 3 4 ) 10   ( c 1 3 c 3 2 + c 1 2 c 3 3 ) c 1 0 c 3 5 c 3 6 = [ 8 + α 3 ] c 3 5 + c 1 5 c 2 0 + 5 ( c 1 4 c 2 1 + c 1 1 c 2 4 ) + 10 ( c 1 3 c 2 2 + c 1 2 c 2 3 ) + c 1 0 c 2 5
The solution to the system equation can be expressed as:
x j ( t ) = c j 0 + c j 1 ( t t 0 ) q 1 ! q + c j 2 ( t t 0 ) 2 q 2 ! q 2 + c j 3 ( t t 0 ) 3 q 3 ! q 3 + c j 4 ( t t 0 ) 4 q 4 ! q 4 + c j 5 ( t t 0 ) 5 q 5 ! q 5 + c j 6 ( t t 0 ) 6 q 6 ! q 6
It can be seen from the above calculation process that the improved Adomian decomposition method has no Gamma function in the process of decomposition, so the decomposition process is simpler than the Adomian decomposition method, which it is more convenient in the programming of the program and has little calculation amount.
When q1 = q2 = q3 = 0.96, α = 0, h = 0.001, the phase space diagram is shown in Figure 10.
For the convenience of comparison, Table 1 can obtain the calling function method, the Adams prediction–correction method, the Adomian decomposition method, and the improved Adomian decomposition method performance comparison results by running the program and timing. From the results, the improved Adomian decomposition method works best.
The above three methods are directly calculated in the time domain, avoiding the approximation error caused by transforming time domain into frequency domain and improving the accuracy compared with the calling function method.

4. Circuit Simulation of Fractional-Order Unified Chaotic System

In practical engineering application, there are many algorithms for solving fractional differential equations. The most commonly used is time domain and complex frequency domain transformation method. By solving transfer function 1 s q , its expanded form is obtained and then transformed to time domain for solution. In reference [22], the expansion of 1 s q is determined by Bode graphic approximation method, and the 1 s q expansion of q from 0.1 to 0.9 is derived in reference [23]. Generally, when designing a fractional-order chaotic circuit, the corresponding integer order chaotic circuit should be designed firstly, and the capacitance in the integer circuit is replaced by a fractional-order module to form a fractional-order chaotic circuit [24].
The key point of designing fractional-order chaotic circuit is the circuit design of the integral operator in the circuit diagram. Because the integral operator of q can be described in the frequency domain by the transfer function F(s)= 1 s q , the following unit circuit can be used to achieve the fractional order equivalent circuit [25], as shown in Figure 11.
The system function [26] of the equivalent circuit between A and B can be expressed as:
H ( s ) = R 1 R 1 C 1 + R 2 R 2 C 2 + R 3 R 3 C 3 + + R n R n C n
Use Table 2 and Table 3 to calculate the cell resistance and capacitance values for the fractional 1 s q property.
The circuit diagram of the fractional-order unified chaotic system is shown in Figure 12.
For fractional-order unified chaotic system, the linear part is composed of basic operational amplifier circuit, such as reverse proportional amplifier, reverse adder, subtracter, in-phase proportional amplifier, reverse integrator, etc., and the nonlinear term is represented by multiplier. The parameters of each component are as follows:
{ R 1 = R 8 = R 20 = R 21 = 10 K Ω R 2 = R 9 = 35 K Ω R 4 = R 15 = R 19 = R 22 = R 23 = R 27 = 100 K Ω R 10 = 10.15 K Ω R 11 = 29 K Ω R 14 = 10.35 K Ω R 29 = 30 K Ω
And when α = 0.5,
{ R 3 = 22.5 K Ω R 12 = 20 K Ω R 13 = 58 K Ω R 28 = 8.5 K Ω
When q1 = q2 = q3 = 0.9, use Table 2 and Table 3 to calculate:
{ R 5 = R 16 = R 24 = 1.55 M Ω R 6 = R 17 = R 25 = 6154 M Ω R 7 = R 18 = R 26 = 2.5 M Ω C 1 = C 2 = C 7 = 0.73 μ F C 3 = C 5 = C 8 = 0.52 μ F C 4 = C 6 = C 9 = 1.1 μ F
The simulation circuit diagram of the fractional-order unified chaotic system is built in the Multisim 14.1 version. The simulation result is shown in Figure 13.
It can be seen from Figure 13 that the circuit simulation diagram of the fractional-order unified chaotic system is basically consistent with the phase space diagram obtained by the numerical solution of the system, and it is verified that the 0.9-order fractional-order unified chaotic system is exists.

5. Control and Synchronization of Fractional-Order Unified Chaotic System

5.1. Control of Fractional-Order Unified Chaotic System

The controlled model of the fractional-order unified chaotic system is:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) + u 1 D q 2 x 2 = ( 28 35 α ) x 1 x 1 x 3 + ( 29 α 1 ) x 2 + u 2 D q 3 x 3 = x 1 x 2 ( 8 + α ) x 3 3 + u 3
Select:
{ u 1 = v 1 u 2 = x 1 x 3 + v 2 u 3 = x 1 x 2 + v 3
where, u 1 , u 2 , u 3 are the control input, v 1 , v 2 , v 3 are the linear functions of x 1 , x 2 , x 3 .
Get the following formula:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) + v 1 D q 2 x 2 = ( 28 35 α ) x 1 + ( 29 α 1 ) x 2 + v 2 D q 3 x 3 = ( 8 + α ) x 3 3 + v 3
To ensure stability, select v 1 , v 2 , v 3 , respectively:
{ v 1 = k 1 x 1 v 2 = k 2 x 2 v 3 = 0
where k 1 , k 2 is the gain, and the system becomes:
{ D q 1 x 1 = ( 25 α + 10 ) ( x 2 x 1 ) k 1 x 1 D q 2 x 2 = ( 28 35 α ) x 1 + ( 29 α 1 ) x 2   k 2 x 2 D q 3 x 3 = ( 8 + α ) x 3 3
The system is written in matrix form as:
D q x = A x
A = [ 25 α 10 k 1 25 α + 10 0 28 35 α 29 α 1 k 2 0 0 0 ( α + 8 ) / 3 ]
The characteristic polynomial is:
Det ( λ E A )   =   [ ( λ   +   25 α   +   10   +   k 1 ) ( λ     29 α   +   1   +   k 2 )     ( 25 α   +   10 ) ( 28     35 α ) ] = λ 2     29 λ α   +   λ   +   λ k 2   +   25 λ     725 α 2 +   25 α   +   25 α k 2   +   10 λ     290 α   +   10   +   10 k 2   +   λ k 1   29 α k 1   +   k 1   +   k 1 k 2     700 α   +   875 α 2   280   +   350 α
Collation can be obtained as follows:
Det ( λ E     A )   =   ( λ   +   ( α   +   8 ) / 3 )   ×   [ λ 2   +   ( k 1   +   k 2     4 α   +   11 ) λ   +   150 α 2   +   ( 29 k 1   +   25 k 2     615 ) α +   ( k 1   +   10 ) (   k 2   +   1 )     280 ]
According to the lemma [27] d q x d t q = Ax, 0 < q < 1, if any eigenvalues of the matrix A meet | a r g ( e i g ( A ) ) | > q π 2 , and the system is asymptotically stable. When the eigenvalues of the matrix are negative, the above research system is stable. So, the solution to k 1 , k 2 is worth noting that the eigenvalue is negative. The correlation properties of the matrix eigenvalues can be used to the following inequalities:
{ k 1 + k 2 4 α 11 29 k 1 25 k 2 615 ( k 1 + 10 ) ( k 2 + 1 ) 280
{ k 1 + k 2 4 α 11 29 k 1 25 k 2 315 k 1 k 2 28 k 1 + 35 k 2 735
To solve the inequality system, we can get:
{ k 1 = 0 k 2 = 27
and
{ k 1 = 10 k 2 = 23
When q1 = q2 = q3 = 0.9 and h = 0.001, the control effect is shown in Figure 14. Figure 14 shows that the three control laws can make the state of fractional-order unified chaotic system asymptotically approach zero.

5.2. Synchronization of Fractional-Order Unified Chaotic System

The response system is:
{ D q 1 y 1 = ( 25 β + 10 ) ( y 2 y 1 ) + u 1 D q 2 y 2 = ( 28 35 β ) y 1 y 1 y 3 + ( 29 β 1 ) y 2 + u 2 D q 3 y 3 = y 1 y 2 ( 8 + β ) y 3 3 + u 3
The error between the response system and the drive system is:
{ e 1 = y 1 x 1 e 2 = y 2 x 2 e 3 = y 3 x 3
The dynamic equation of the error system is obtained as follows:
{ D q 1 e 1 = ( 25 β + 10 ) ( e 2 e 1 ) + 25 ( β α ) ( x 2 x 1 ) + u 1 D q 2 e 2 = ( 28 35 β ) e 1 + ( 29 β 1 ) e 2 35 ( β α ) x 1 + e 1 e 3 e 1 y 3 e 3 y 1 + 29 ( β α ) x 2 + u 2 D q 3 e 3 = ( 8 + β ) 3 e 3 ( β α ) 3 x 3 e 1 e 2 + e 1 y 2 + e 2 y 1 + u 3
Since it is the same fractional chaotic system, the selected gain k 1 and k 2 can be used. When α = 1, β = 0.5, k 1 = 10 ,   k 2 = 23 , the solid line represents the variation curve of the drive system state variable x(t), and the broken line represents the variation curve of the response system y(t). The synchronization effect diagram is shown in Figure 15. The control input diagram and the synchronization error diagram are shown in Figure 16.
Figure 15 and Figure 16 show that under the same control law, the response system of the fractional-order unified chaotic system is synchronized with the drive system, and has good control performance.

6. Conclusions

In this paper, different numerical methods of fractional calculus are used to analyze the fractional-order unified chaotic system, that is the Adams prediction–correction method, the Adomian decomposition method and the improved Adomian decomposition method. Secondly, 0.9-order fractional-order unified chaotic circuit simulation diagram is designed. Finally, the control and synchronization of the fractional-order unified chaotic system is completed by the active control method.

Author Contributions

G.L., X.Z., and H.Y. conceived the idea and research theme. X.Z. designed and performed the experiments. G.L., X.Z., and H.Y. analyzed the experimental results. G.L., X.Z., and H.Y. wrote and revised the paper.

Funding

This work was supported by the National Natural Science Foundation of China (no. 51709228).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, G.R.; Ueta, T. Yet another chaotic attractor. Int. J. Bifurc. Chaos 1999, 9, 1465–1466. [Google Scholar] [CrossRef]
  2. Lü, J.H.; Chen, G.R.; Zhang, S.C. Dynamical analysis of a new chaotic attractor. Int. J. Bifurc. Chaos 2002, 12, 1001–1015. [Google Scholar] [CrossRef]
  3. Rössler, O.E. An equation for continuous chaos. Phys. Lett. A 1976, 57, 397–398. [Google Scholar] [CrossRef]
  4. Lü, J.H.; Chen, G.R.; Cheng, D.Z. A new chaotic system and beyond: the generalized Lorenz-like system. Int. J. Bifurc. Chaos 2004, 14, 1507–1537. [Google Scholar] [CrossRef]
  5. Liu, W.B.; Chen, G.R. A new chaotic system and its generation. Int. J. Bifurc. Chaos 2003, 13, 261–267. [Google Scholar] [CrossRef]
  6. Liu, C.X.; Liu, T.; Liu, L.; Liu, K. A new chaotic attractor. Chaos Solitons Fractals 2004, 22, 1031–1038. [Google Scholar] [CrossRef]
  7. Qi, G.Y.; Chen, G.R.; Du, S.Z.; Chen, Z.Q.; Yuan, Z.Z. Analysis of a new chaotic system. Phys. A 2005, 352, 295–308. [Google Scholar] [CrossRef]
  8. Bao, B.C.; Liu, Z.; Xu, J.P. New chaotic system and its hyperchaos generation. J. Syst. Eng. Electron. 2009, 20, 1179–1187. [Google Scholar]
  9. Chen, X.R.; Liu, C.X.; Wang, F.Q.; Li, Y.X. Study on the fractional-order Liu chaotic system with circuit experiment and its control. Acta Phys. Sin. 2008, 57, 1416–1422. [Google Scholar]
  10. Wang, F.Q.; Liu, C.X. Study on the critical chaotic system with fractional order and circuit experiment. Acta Phys. Sin. 2006, 55, 3922–3927. [Google Scholar]
  11. Mao, X.X.; Sun, K.H.; Liu, W.H. Image encryption algorithm based on fractional order unified chaotic system. Transducer Microsyst. Technol. 2017, 36, 138–141. [Google Scholar]
  12. He, S.B.; Sun, K.H.; Wang, H.H. Solution of the fractional-order chaotic system based on Adomian decomposition algorithm and its complexity analysis. Acta Phys. Sin. 2014, 63, 030502. [Google Scholar]
  13. Wu, C.M. Linear feedback synchronization control method for fractional order chaotic systems. Control Eng. China 2019, 26, 898–902. [Google Scholar]
  14. Yan, X.M.; Shang, T.; Zhao, X.G. Synchronization of uncertain fractional-order chaotic systems based on the fractional-order sliding mode controller. Acta Math. Appl. Sin. 2018, 41, 765–776. [Google Scholar]
  15. Rabah, K.; Ladaci, S.; Lashab, M. Bifurcation-based fractional-order PIλDμ controller design approach for nonlinear chaotic systems. Front. Inf. Technol. Electron. Eng. 2018, 19, 180–191. [Google Scholar] [CrossRef]
  16. Rajagopal, K.; Karthikeyan, A.; Srinivasan, A. Bifurcation and chaos in time delayed fractional order chaotic memfractor oscillator and its sliding mode synchronization with uncertainties. Chaos Solitons Fractals 2017, 103, 347–356. [Google Scholar] [CrossRef]
  17. Tao, C.H.; Lu, J.A.; Lü, J.H. The feedback synchronization of a unified chaotic system. Acta Phys. Sin. 2002, 51, 1497–1501. [Google Scholar]
  18. Wu, Q.; Huang, J.H. Fractional Calculus; Tsinghua University Press: Beijing, China, 2016. [Google Scholar]
  19. Psihoyios, G.; Simos, T.E. Trigonometrically fitted Adams-Bashforth-Moulton methods for periodic initial value problems. Proceedings Second MIT Conference on Compurational Fluid and Solid Mechanics, Cambridge, MA, USA, 17–20 June 2003. [Google Scholar]
  20. Zhu, Y.G.; Chang, Q.S.; Wu, S.C. A new algorithm for calculating Adomian polynomials. Appl. Math. Comput. 2005, 169, 402–416. [Google Scholar] [CrossRef]
  21. Lei, T.F.; Zhang, X.; Zang, H.Y.; Xia, X.X. Experimental simulation of fractional Rabinovich hyperchaotic system based on improved Adomian decomposition method. Exp. Technol. Manag. 2018, 35, 137–140. [Google Scholar]
  22. Charef, A.; Sun, H.H.; Tsao, Y.Y.; Onaral, B. Fractal system as represented by singularity function. IEEE Trans. Autom. Control 1992, 37, 1465–1470. [Google Scholar] [CrossRef]
  23. Ahmad, W.M.; Sprott, J.C. Chaos in fractional-order autonomous nonlinear systems. Chaos Solitons Fractals 2003, 16, 339–351. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, X.H.; Sun, Q. Construction of heterogeneous fractional-order chaotic system and its multiple circuits simulation. J. Syst. Simul. 2014, 26, 1460–1466. [Google Scholar]
  25. Lei, T.F.; Fu, H.Y.; Zhang, X.; Wang, Z.L. Dynamics analysis and circuit implementation of a fractional-order chaotic systems based on adomian decomposition. Chin. J. Electron Devices 2019, 42, 193–200. [Google Scholar]
  26. Sun, K.H.; Yang, J.L.; Qiu, S.S. Circuit simulation and implement for fractional-order chaotic systems. Comput. Simul. 2011, 28, 117–119. [Google Scholar]
  27. Geng, Y.F.; Wang, L.Z. Function projective synchronization of fractional-order united chaotic system based on sliding mode control. J. Tianjin Norm. Univ. Nat. Sci. Ed. 2019, 39, 23–26. [Google Scholar]
Figure 1. Phase space diagram of x-y plane. (a) α = 0; (b) α = 0.5; (c) α = 0.8; (d) α = 1.
Figure 1. Phase space diagram of x-y plane. (a) α = 0; (b) α = 0.5; (c) α = 0.8; (d) α = 1.
Mathematics 07 01077 g001aMathematics 07 01077 g001b
Figure 2. Time domain waveform of integer-order unified chaotic system with α = 1. (a) x(t); (b) y(t).
Figure 2. Time domain waveform of integer-order unified chaotic system with α = 1. (a) x(t); (b) y(t).
Mathematics 07 01077 g002
Figure 3. Lyapunov exponent of integral-order unified chaotic system with α ∈ [0.8,1].
Figure 3. Lyapunov exponent of integral-order unified chaotic system with α ∈ [0.8,1].
Mathematics 07 01077 g003
Figure 4. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Figure 4. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Mathematics 07 01077 g004aMathematics 07 01077 g004b
Figure 5. The flow of numerical calculation of the Adams prediction–correction method.
Figure 5. The flow of numerical calculation of the Adams prediction–correction method.
Mathematics 07 01077 g005
Figure 6. Phase space diagram of each plane when α = 0.8. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Figure 6. Phase space diagram of each plane when α = 0.8. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Mathematics 07 01077 g006
Figure 7. Phase space diagram of each plane when α = 1. (a) x-y-z plane; (b) x-y plane; (c) x-z plane, (d) y-z plane.
Figure 7. Phase space diagram of each plane when α = 1. (a) x-y-z plane; (b) x-y plane; (c) x-z plane, (d) y-z plane.
Mathematics 07 01077 g007
Figure 8. Duffing oscillator phase space diagram when q1 = q2 = q3 = 0.98.
Figure 8. Duffing oscillator phase space diagram when q1 = q2 = q3 = 0.98.
Mathematics 07 01077 g008
Figure 9. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Figure 9. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Mathematics 07 01077 g009
Figure 10. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Figure 10. Phase space diagram of each plane. (a) x-y-z plane; (b) x-y plane; (c) x-z plane; (d) y-z plane.
Mathematics 07 01077 g010
Figure 11. Equivalent circuit of fractional integral operator 1 s q .
Figure 11. Equivalent circuit of fractional integral operator 1 s q .
Mathematics 07 01077 g011
Figure 12. Fractional-order unified chaotic system circuit diagram.
Figure 12. Fractional-order unified chaotic system circuit diagram.
Mathematics 07 01077 g012
Figure 13. Circuit simulation diagram. (a) x-y plane; (b) x-z plane; (c) y-z plane.
Figure 13. Circuit simulation diagram. (a) x-y plane; (b) x-z plane; (c) y-z plane.
Mathematics 07 01077 g013
Figure 14. Control effect diagram. (a) α = 1, k 1 = 10 , k 2 = 23 ; (b) α = 0, k 1 = 10 , k 2 = 23 .
Figure 14. Control effect diagram. (a) α = 1, k 1 = 10 , k 2 = 23 ; (b) α = 0, k 1 = 10 , k 2 = 23 .
Mathematics 07 01077 g014
Figure 15. Synchronization effect diagram.
Figure 15. Synchronization effect diagram.
Mathematics 07 01077 g015
Figure 16. Control input and error diagram. (a) Control input diagram; (b) Error diagram.
Figure 16. Control input and error diagram. (a) Control input diagram; (b) Error diagram.
Mathematics 07 01077 g016
Table 1. Performance comparison of different algorithms.
Table 1. Performance comparison of different algorithms.
Calling Function MethodAdams Prediction–Correction MethodAdomian Decomposition MethodImproved Adomian Decomposition Method
Time ComplexityO( n 2 ), SlowlyO( n 2 ), SlowlyO(n), QuicklyO(n), Quickly
Space ComplexityO(n), ComplexO(n), ComplexO(1), SimpleO(1), Simple
N = 10001.823s1.298s0.487s0.306s
N = 20005.032s4.680s0.510s0.403s
N = 500017.897s15.362s0.858s0.623s
N = 10,00031.259s28.191s1.458s0.919s
Table 2. Resistance value in q-order fractional unit circuit.
Table 2. Resistance value in q-order fractional unit circuit.
1 s q R 1 / M Ω R 2 / M Ω R 3 / M Ω R 4 / M Ω R 5 / M Ω R 6 / M Ω
1 s 0.1 0.6360.38150.5672
1 s 0.2 1.1300.60700.35000.4250.2498
1 s 0.3 2.0500.93800.48300.2620.14300.106
1 s 0.4 3.7441.39200.63100.2940.123.0.068
1 s 0.5 6.8241.94400.74400.2960.07540.030
1 s 0.6 12.3302.444800.73800.2330.07540.006
1 s 0.7 21.9002.60000.52600.1130.0246
1 s 0.8 37.8501.75400.17000.0170.0018
1 s 0.9 62.8400.25000.0025
Table 3. Capacitance value in q-order fractional unit circuit.
Table 3. Capacitance value in q-order fractional unit circuit.
1 s q C 1 / μ F C 2 / μ F C 3 / μ F C 4 / μ F C 5 / μ F C 6 / μ F
1 s 0.1 15.7200.15720.0006335
1 s 0.2 27.9902.93000.2850.0132
1 s 0.3 22.6405.52001.2000.24600.029
1 s 0.4 15.0205.92601.9200.60500.1830.036
1 s 0.5 9.2465.14502.1290.84800.3240.925
1 s 0.6 5.5274.08501.9900.92600.4200.156
1 s 0.7 3.2843.13901.7000.88600.4540.207
1 s 0.8 1.9802.40001.3900.78000.420
1 s 0.9 1.2321.84001.100

Share and Cite

MDPI and ACS Style

Li, G.; Zhang, X.; Yang, H. Numerical Analysis, Circuit Simulation, and Control Synchronization of Fractional-Order Unified Chaotic System. Mathematics 2019, 7, 1077. https://doi.org/10.3390/math7111077

AMA Style

Li G, Zhang X, Yang H. Numerical Analysis, Circuit Simulation, and Control Synchronization of Fractional-Order Unified Chaotic System. Mathematics. 2019; 7(11):1077. https://doi.org/10.3390/math7111077

Chicago/Turabian Style

Li, Guohui, Xiangyu Zhang, and Hong Yang. 2019. "Numerical Analysis, Circuit Simulation, and Control Synchronization of Fractional-Order Unified Chaotic System" Mathematics 7, no. 11: 1077. https://doi.org/10.3390/math7111077

APA Style

Li, G., Zhang, X., & Yang, H. (2019). Numerical Analysis, Circuit Simulation, and Control Synchronization of Fractional-Order Unified Chaotic System. Mathematics, 7(11), 1077. https://doi.org/10.3390/math7111077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop