Next Article in Journal
Sound Wave Absorption Coefficient and Sound Velocity in Thermally Modified Wood
Previous Article in Journal
Compression Behavior of Rubber–Bentonite Mixture under Different Salinities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Function Projective Synchronization of Two Different Chaotic Systems with Uncertain Parameters

School of Environment and Architecture, University of Shanghai for Science and Technology, Shanghai 200093, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(14), 8135; https://doi.org/10.3390/app13148135
Submission received: 11 June 2023 / Revised: 7 July 2023 / Accepted: 10 July 2023 / Published: 12 July 2023

Abstract

:
This study proposes a new approach to realize generalized function projective synchronization (GFPS) between two different chaotic systems with uncertain parameters. The GFPS condition is derived by converting the differential equations describing the synchronization error systems into a series of Volterra integral equations on the basis of the Laplace transform method and convolution theorem, which are solved by applying the successive approximation method in the theory of integral equations. Compared with the results obtained by constructing Lyapunov functions or calculating the conditional Lyapunov exponents, the uncertain parameters and the scaling function factors considered in this paper have fewer restrictions, and the parameter update laws designed for the estimation of the uncertain parameters are simpler and easier to realize physically.

1. Introduction

Since Pecora and Carroll realized synchronization between two chaotic systems in 1990 [1,2], chaos synchronization has been widely explored and studied due to its potential applications in vast areas of physics and engineering science [3,4,5,6]. Synchronization is a process wherein two or more systems adjust a given property of their motion. Different types of synchronization phenomena have been numerically observed and experimentally verified in a variety of chaotic systems, such as complete synchronization [1,7], phase synchronization [8,9], anti-phase synchronization [10,11], lag synchronization [12,13], generalized synchronization [14,15] and projective synchronization [16,17].
Complete synchronization (CS) is characterized by a coincidence of states of two chaotic systems while evolving in time. CS only appears when two interacting systems are identical. However, it is difficult to construct two absolutely identical systems in the real world. Generalized synchronization (GS) refers to a functional relation of states of two different chaotic systems, which is an extension of CS. Projective synchronization (PS) is more complex than CS but simpler than GS. Additionally, PS is useful in extending binary digital to M-nary digital communication to achieve faster communication [18,19,20].
PS was first reported by Mainieri et al. [16] in a class of systems with partial linearity in which the drive and response vectors evolve in a proportional scale. Xu [21] showed that the scaling factor of PS in coupled partially linear systems is unpredictable and can be arbitrarily maneuvered by introducing a feedback control to master systems. The PS between two chaotic discrete dynamical systems was achieved by Xin and Wu [22] via linear state error feedback control. Wen and Xu [23] and Yan and Li [24] extended the projective synchronization feature to general nonlinear systems, including non-partially linear chaotic systems, by applying controllers to response systems, which is called generalized projective synchronization (GPS) and was studied later [25,26,27]. CS and anti-synchronization are particular cases of GPS. As an extension of GPS, modified projective synchronization (MPS) was introduced by Li et al. [28] and He et al. [29], where the responses of synchronized dynamical states can be synchronized up to a constant matrix. The combined projective synchronization (CPS) among different nonlinear systems has been reported by Feng et al. [30] by means of the active backstepping design method. Chen et al. [31,32] considered two chaotic systems synchronized up to a scaling function factor, which is called function projective synchronization (FPS). The FPS between two novel chaotic systems has been investigated by El Dessoky et al. [33] and Bekiros et al. [34] using the Lyapunov method of stability. Du et al. [35] and Srisuntorn et al. [36] discussed modified function projective synchronization (MFPS), in which two chaotic systems can be synchronized up to a desired scaling function matrix. As another extension of MPS and FPS, in generalized function projective synchronization (GFPS), the scaling function matrix is more generalized than that of MFPS, as has been discussed by Yu and Li [37], Li and Zhao [38], Wu et al. [39] and Al-Azzawi [40].
Most of the previous works have focused on the constant scaling factor and mainly on chaotic systems with exactly known parameters. However, in many practical situations, GFPS needs to be investigated while some parameters are uncertain. It is easy to understand that GFPS may be greatly affected by these uncertainties. Yu and Li [37] designed an adaptive controller to achieve GFPS between two different chaotic systems and gain some parameter update laws to estimate the unknown parameters of chaotic systems. Li and Zhao [38] firstly discussed the GFPS conditions of two different uncertain hyperchaotic systems when the scaling function factors were periodic functions and polynomial functions. However, only a few theoretical results about GFPS between two different uncertain chaotic systems have been obtained [41,42,43,44,45]. Furthermore, the aforementioned controllers all are designed through the construction of proper Lyapunov functions. The added controllers therefore are often very specific and may be too difficult to realize physically.
In this paper, an effective approach is proposed to design an adaptive controller to realize GFPS between two different uncertain chaotic systems. Moreover, the unknown parameters in the coupled system can be estimated by given parameter update laws. The key to realizing GFPS between two chaotic systems is the determination of the stability of the synchronization error system. In previous works, this problem usually was solved by using the Lyapunov method. However, there are two main difficulties in this method: the first one is that it is challenging to construct a proper Lyapunov function; the second is that the GFPS conditions derived using the Lyapunov method are usually sufficient and highly conserved. The controllers and parameter update laws designed based on such conditions are often very specific. In fact, GFPS can be viewed as a unique case of GS, which can be investigated using the auxiliary system approach [46]. To avoid the construction of Lyapunov functions, the approach presented in this paper converts the differential equations describing the synchronization error system between the auxiliary system and the response system into a series of Volterra integral equations by using the Laplace transform method. Then, the dynamical behavior of the error system near its origin can be analyzed by employing the successive approximation method [47]. In the theory of integral equations, the significance of the Laplace transform method is that the GFPS condition can be obtained without the construction of Lyapunov functions, even if the scaling function factors are periodic functions, polynomial functions or any other complex functions. Compared to previous works, our approach allows the easier design of controllers and simpler parameter update laws and removes some restrictions on the uncertain parameters. The feasibility and effectiveness of the approach are illustrated by several numerical simulations.
The rest of the paper is organized as follows. In Section 2, the GFPS scheme for two different uncertain chaotic systems is introduced. In Section 3, the GFPS condition is investigated using the Laplace transform method combined with the successive approximation method in the theory of Volterra integral equations. A numerical example is provided in Section 4 to demonstrate the effectiveness of our GFPS scheme. In Section 5, simple parameter update laws are designed according to the GFPS condition derived in Section 3 to estimate the uncertain parameters in chaotic systems, the validity of which are illustrated by numerical simulations. Finally, conclusions are drawn in Section 6.

2. GFPS Scheme between Two Different Uncertain Chaotic Systems

Consider two chaotic systems (drive–response systems) with uncertain parameters given in the following form:
x ˙ = f 1 ( x , ϕ ) , y ˙ = f 2 ( y , ψ ) + u ,
where x , y R n ; ϕ R m , ψ R l are uncertain parameter vectors, and u is the controller to be designed. For simplicity, assume that f 1 ( 0 , ϕ ) = f 2 ( 0 , ψ ) = 0 hold for any values of ϕ and ψ . If there exists a scaling function matrix h = d i a g { h 1 ( x ) , h 2 ( x ) , , h n ( x ) } , such that
lim t | | y h ( x ) x | | = 0 ,
then the two chaotic systems in system (1) are said to be GFPS with respect to the scaling function matrix h [37,40]. The two chaotic systems are said to be FPS if
h 1 ( x ) = h 2 ( x ) = = h n ( x )
in Equation (2). Moreover, we find MPS between the two chaotic systems if h = d i a g { h 1 , h 2 , , h n } ,   h i ( i = 1 , 2 , , n ) are real constants in Equation (2). Obviously, both FPS and MPS are specific cases of GFPS.
The controller u and the parameter update laws in system (1) are designed as below:
u = J ( x ) f 1 ( x , ϕ ^ ) f 2 ( y , ψ ^ ) + k δ 1 , ϕ ^ ˙ = g 1 ( ϕ ^ ϕ , δ 1 ) , ψ ^ ˙ = g 2 ( ψ ^ ψ , δ 1 ) ,
where J ( x ) = d ( h ( x ) x ) / d x , k is the control parameter, and δ 1 = y h ( x ) x is the synchronization error; ϕ ^ and ψ ^ are the estimated values of uncertain parameter vectors ϕ and ψ , respectively. In addition, g 1 ( ϕ ^ ϕ , 0 ) = g 2 ( ψ ^ ψ , 0 ) = 0 hold if
ϕ ^ = ϕ a n d ψ ^ = ψ .
The GFPS occurring in system (1) can be regarded as a specific case of GS between the two chaotic systems, which can be detected using the auxiliary system approach [46]. Consider the following identical copy of the response system driven by the same driving signal x,
z ˙ = f 2 ( z , ψ ) + u ¯ , u ¯ = J ( x ) f 1 ( x , ϕ ¯ ) f 2 ( z , ψ ¯ ) + k δ 2 , ϕ ¯ ˙ = g 1 ( ϕ ¯ ϕ , δ 2 ) , ψ ¯ ˙ = g 2 ( ψ ¯ ψ , δ 2 ) ,
in which z R n and δ 2 = z h ( x ) x . ϕ ¯ and ψ ¯ are the estimated values of uncertain parameter vectors ϕ and ψ , respectively. Additionally, g 1 ( ϕ ¯ ϕ , 0 ) = g 2 ( ψ ¯ ψ , 0 ) = 0 hold for
ϕ ¯ = ϕ a n d ψ ¯ = ψ .
Then, GFPS conditions (2) and (4) lead to
lim t | | y z | | = 0 , lim t | | ϕ ^ ϕ ¯ | | = 0 , lim t | | ψ ^ ψ ¯ | | = 0 .

3. Detection of GFPS Condition Based on the Laplace Transform

By letting
e 1 = y z 2 , e 2 = ϕ ^ ϕ ¯ 2 , e 3 = ψ ^ ψ ¯ 2 e 4 = y + z 2 , e 5 = ϕ ^ + ϕ ¯ 2 , e 6 = ψ ^ + ψ ¯ 2
systems (1), (3) and (5) become
e ˙ 1 = k e 1 + J ( x ) 2 [ f 1 ( x , e 2 + e 5 ) f 1 ( x , e 5 e 2 ) ] + 1 2 [ f 2 ( e 1 + e 4 , ψ ) f 2 ( e 4 e 1 , ψ ) ] + 1 2 [ f 2 ( e 4 e 1 , e 6 e 3 ) f 2 ( e 4 + e 1 , e 6 + e 3 ) ] , e ˙ 2 = 1 2 [ g 1 ( e 2 + e 5 ϕ , e 1 + e 4 h ( x ) x ) g 1 ( e 5 e 2 ϕ , e 4 e 1 h ( x ) x ) ] , e ˙ 3 = 1 2 [ g 2 ( e 3 + e 6 ψ , e 1 + e 4 h ( x ) x ) g 2 ( e 6 e 3 ψ , e 4 e 1 h ( x ) x ) ] , e ˙ 4 = k e 4 + J 2 [ f 1 ( x , e 2 + e 5 ) + f 1 ( x , e 5 e 2 ) ] + 1 2 [ f 2 ( e 1 + e 4 , ψ ) + f 2 ( e 4 e 1 , ψ ) ] 1 2 [ f 2 ( e 4 e 1 , e 6 e 3 ) + f 2 ( e 4 + e 1 , e 6 + e 3 ) ] , e ˙ 5 = 1 2 [ g 1 ( e 2 + e 5 ϕ , e 1 + e 4 h ( x ) x ) + g 1 ( e 5 e 2 ϕ , e 4 e 1 h ( x ) x ) ] , e ˙ 6 = 1 2 [ g 2 ( e 3 + e 6 ψ , e 1 + e 4 h ( x ) x ) + g 2 ( e 6 e 3 ψ , e 4 e 1 h ( x ) x ) ] , x ˙ = f 1 ( x , ϕ ) .
GFPS condition (7) can be written as
lim t | | e i | | = 0 , i = 1 , 2 , 3 .
For sufficiently small e 1 , 2 , 3 near ( 0 , 0 , 0 ) , the right-hand side of the first three equations in system (8) can be expanded as
e ˙ 1 = k e 1 + F 1 ( x , e 2 , e 3 ) , e ˙ 2 = α 1 e 1 + α 2 e 2 + F 2 ( e 1 , e 2 ) , e ˙ 3 = α 3 e 1 + α 4 e 3 + F 3 ( e 1 , e 3 ) .
where
F 1 ( x , e 2 , e 3 ) = J ( x ) n = 1 1 ( 2 n 1 ) ! ( 2 n 1 ) f 1 ( x , ϕ ) ϕ ( 2 n 1 ) e 2 ( 2 n 1 ) n = 1 1 ( 2 n 1 ) ! ( 2 n 1 ) f 2 ( h ( x ) x , ψ ) ψ ( 2 n 1 ) e 3 ( 2 n 1 ) , F 2 ( e 1 , e 2 ) = n = 1 1 ( 2 n + 1 ) ! [ ( 2 n + 1 ) g 1 ( 0 , 0 ) δ 1 ( 2 n + 1 ) e 1 ( 2 n + 1 ) + ( 2 n + 1 ) g 1 ( 0 , 0 ) ( ϕ ^ ϕ ) ( 2 n + 1 ) e 2 ( 2 n + 1 ) ] , F 3 ( e 1 , e 3 ) = n = 1 1 ( 2 n + 1 ) ! [ ( 2 n + 1 ) g 2 ( 0 , 0 ) δ 1 ( 2 n + 1 ) e 1 ( 2 n + 1 ) + ( 2 n + 1 ) g 2 ( 0 , 0 ) ( ψ ^ ψ ) ( 2 n + 1 ) e 3 ( 2 n + 1 ) ] , α 1 = g 1 ( 0 , 0 ) δ 1 , α 2 = g 1 ( 0 , 0 ) ( ϕ ^ ϕ ) , α 3 = g 2 ( 0 , 0 ) δ 1 , α 4 = g 2 ( 0 , 0 ) ( ψ ^ ψ ) .
Consider the Laplace transform, defined as follows:
E j ( s ) = L [ e j ] = 0 + e j ( t ) e s t d t , e j ( t ) = L 1 [ E j ] = 1 2 π i σ i σ + i E j ( s ) e s t d s , j = 1 , 2 , 3 .
Both e j and E j are vectors with the same dimensions for j = 1 , 2 , 3 . Taking the Laplace transform of e 1 , 2 , 3 on both sides of system (10) yields
( S M ) [ E 1 , E 2 , E 3 ] T = [ e 10 , e 20 , e 30 ] T + [ F ^ 1 , F ^ 2 , F ^ 3 ] T ,
where
S = s I n 0 0 0 s I m 0 0 0 s I l , M = k I n 0 0 α 1 α 2 0 α 3 0 α 4 ,
I n , I m , I l represent the n × n , m × m and l × l real identity matrices, respectively. T denotes the transpose of vectors. e j 0 , j = 1 , 2 , 3 , are given initial values of system (10) and F ^ 1 , 2 , 3 are the Laplace transforms of F 1 , 2 , 3 , respectively, i.e.,
F ^ j = L [ F j ] = 0 + F j ( t ) e s t d t , j = 1 , 2 , 3 .
The solution to Equation (12) can be derived by using Cramer’s rule [48]. Take, for example, the expression of E 1 , which is an n-dimensional vector and always can be given by the form below:
E 1 = E 11 E 12 E 1 n = 1 D ( s ) w 11 ( s ) w 12 ( s ) w 1 n ( s ) + 1 D ( s ) w 21 ( s ) F ^ 11 ( s ) w 22 ( s ) F ^ 12 ( s ) w 2 n ( s ) F ^ 1 n ( s ) ,
where ( E 11 , E 12 , , E 1 n ) and ( F ^ 11 , F ^ 12 , , F ^ 1 n ) are the n components of the vectors E 1 and F ^ 1 , respectively. D ( s ) is the determinant of matrix ( S M ) . w 1 j , w 2 j , j = 1 , , n , are polynomials of s, and their highest powers are less than that of D ( s ) . Taking the inverse Laplace transform on Equation (14) and applying the convolution theorem, one has
e 1 j = L 1 [ E 1 j ] = L 1 [ w 1 j ( s ) D ( s ) ] + L 1 [ w 2 j ( s ) D ( s ) ] F 1 j ( t ) , j = 1 , , n .
where ( e 11 , e 12 , , e 1 n ) and ( F 11 , F 12 , , F 1 n ) are the n components of the vectors e 1 and F 1 , respectively. ∗ denotes the convolution operation. First, we introduce the following theorem.
Theorem 1.
The necessary condition for e 1 j 0 with t , j = 1 , 2 , , n , in Equation (15) is that all eigenvalues of matrix M defined in Equation (13) have negative real parts.
Proof. 
Without loss of generality, for any i = 1 , 2 and j = 1 , 2 , , n , w i j ( s ) / D ( s ) have the following form
w i j ( s ) D ( s ) = p = 0 v a p s p p = 0 n + m + l b p s p ,
where v < n + m + l , a p and b p are constants. Consider the following three cases of the inverse Laplace transforms for Equation (16):
  • D ( s ) has n + m + l single roots s 1 , s 2 , , s ( n + m + l )
    w i j ( s ) D ( s ) = q 1 s s 1 + q 2 s s 2 + + q ( n + m + l ) s s ( n + m + l ) ,
    where q p = w i j ( s ) D ( s ) ( s s p ) | s = s p , p = 1 , 2 , , n + m + l .
    L 1 [ w i j ( s ) D ( s ) ] = q 1 e s 1 t + q 2 e s 2 t + + q ( n + m + l ) e s ( n + m + l ) t
  • D ( s ) has a pair of conjugate complex roots s 1 , 2 = β ± j ω
    w i j ( s ) D ( s ) = w i j ( s ) ( s β j ω ) ( s β + j ω ) Q ( s ) = q 1 s β j ω + q 2 s β + j ω + P ( s ) Q ( s ) ,
    where q 1 = w i j ( s ) D ( s ) | s = β + j ω , q 2 = w i j ( s ) D ( s ) | s = β j ω ,
    L 1 [ w i j ( s ) D ( s ) ] = q 1 e ( β + j ω ) t + q 2 e ( β j ω ) t + L 1 [ P ( s ) Q ( s ) ]
  • D ( s ) has ( n + m + l ) -repeated real roots s = s 0
    w i j ( s ) D ( s ) = q 1 s s 0 + q 2 ( s s 0 ) 2 + + q ( n + m + l ) ( s s 0 ) ( n + m + l ) ,
    where q ( n + m + l p ) = 1 p ! d p d s p [ ( s s 0 ) ( n + m + l ) w i j ( s ) D ( s ) ] | s = s 0 , p = 1 , 2 , , n + m + l 1 ,
    q ( n + m + l ) = [ ( s s 0 ) ( n + m + l ) w i j ( s ) D ( s ) ] | s = s 0 ,
    L 1 [ w i j ( s ) D ( s ) ] = ( q 1 + q 2 t + + q ( n + m + l ) t ( n + m + l 1 ) ) e s 0 t
From the analysis of the inverse Laplace transform of Equation (16), it is clear that the inverse Laplace transform of w i j ( s ) / D ( s ) , i = 1 , 2 , in Equation (15) must be a sum of exponential functions. When t , all roots of D ( s ) must have negative real parts if
L 1 [ w 1 j ( s ) D ( s ) ] , L 1 [ w 2 j ( s ) D ( s ) ] 0 , j = 1 , 2 , , n
Assume that matrix M only has eigenvalues with negative real parts; from Equation (15), when t + , e 1 j can be written in the following form:
e 1 j = 0 t R 1 j ( t τ ) F 1 j ( τ ) d τ , j = 1 , 2 , , n
where
R 1 j = L 1 [ w 2 j ( s ) D ( s ) ] , j = 1 , 2 , , n
One can easily confirm that the condition that matrix M only has eigenvalues with negative real parts also is necessary for e 2 j 0 , j = 1 , 2 , , m and e 3 j 0 , j = 1 , 2 , , l . Similarly, under such conditions, when t + , e 2 j and e 3 j have the following form:
e 2 j = 0 t R 2 j ( t τ ) F 1 j ( τ ) d τ + 0 t R 3 j ( t τ ) F 2 j ( τ ) d τ , j = 1 , 2 , , m e 3 j = 0 t R 4 j ( t τ ) F 1 j ( τ ) d τ + 0 t R 5 j ( t τ ) F 2 j ( τ ) d τ + 0 t R 6 j ( t τ ) F 3 j ( τ ) d τ , j = 1 , 2 , , l
in which R i j 0 for t + , i = 2 , 3 , , 6 . ( F 21 , F 22 , , F 2 m ) and ( F 31 , F 32 , , F 3 l ) are the m and l components of the vectors F 2 and F 3 , respectively, which are defined in Equation (10).
Theorem 2.
e 1 j = 0 , j = 1 , 2 , , n , e 2 j = 0 , j = 1 , 2 , , m , e 3 j = 0 , j = 1 , 2 , , l , are the unique continuous solutions of Equations (17) and (18).
Proof. 
In fact, Equations (17) and (18) are a series of Volterra integral equations, which can be solved by using the successive approximation method [47] in the theory of integral equations. Consider the Volterra integral equation as below:
A ( t ) = B ( t ) + 0 t T ( t τ ) H ( τ , A ( τ ) ) d τ ,
where A ( t ) R n and | A ( t ) | < , T ( t ) is a n × n matrix, and | T ( t ) | L [ 0 , t 0 ] for any 0 < t 0 < . Both B ( t ) and H ( t , A ( t ) ) are vectors with n continuous components. Furthermore, for any η > 0 , there must exist a constant ρ ( η ) > 0 such that
| H ( t , A 1 ) H ( t , A 2 ) | ρ ( η ) | A 1 A 2 | , ( | A 1 | , | A 2 | η )
From the results given by Nohel [47], if the above conditions are satisfied, the successive approximations
ω 0 ( t ) = 0 , ω n + 1 ( t ) = B ( t ) + 0 t T ( t τ ) H ( τ , ω n ( τ ) ) d τ , n = 0 , 1 , 2 ,
will uniformly converge to the unique continuous solution A ( t ) = ω ( t ) of Equation (19).
Comparing Equations (17)–(19), it is easy to check that F 1 , 2 , 3 , which are defined in Equation (10), are smooth enough to satisfy condition (20). From Equation (21), e 1 j = 0 , j = 1 , 2 , , n , e 2 j = 0 , j = 1 , 2 , , m , e 3 j = 0 , j = 1 , 2 , , l , are the unique continuous solutions of Equations (17) and (18). □
According to Theorems 1 and 2, we have the following result.
Theorem 3.
GFPS in systems (1) and (3) can be achieved if the matrix M defined in Equation (13) only has eigenvalues with negative real parts.

4. Discussion of GFPS Scheme

To illustrate the effectiveness of the GFPS scheme proposed in the previous sections on the basis of the Laplace transform method, the following chaotic systems with uncertain parameters, which have been considered in [37,38], are discussed:
x ˙ = f ( x ) + Φ ( x ) Θ , y ˙ = g ( y ) + Ψ ( y ) Ω + U ,
where x , y R n are the state vectors; Φ ( x ) : R n R n × p , Ψ ( y ) : R n R n × q , Θ R p and Ω R q represent the vectors of uncertain parameters; and U R n is a controller. For simplicity, f ( 0 ) + Φ ( 0 ) Θ = g ( 0 ) + Ψ ( 0 ) Ω = 0 hold for any values of Θ and Ω . Define the error vector of GFPS as
e = y h ( x ) x ,
where h ( x ) = d i a g { h 1 ( x ) , h 2 ( x ) , , h n ( x ) } , h i ( x ) : R n R ( i = 1 , 2 , , n ) are scaling function factors that compose the scaling function matrix h ( x ) . According to [38], GFPS in (22) can be achieved by the controller and the parameter update laws as follows:
U = J ( f ( x ) + Φ ( x ) Θ ^ ) g ( y ) Ψ ( y ) Ω ^ + K e , Θ ^ ˙ = ( J Φ ( x ) ) T e ( Θ ^ Θ ) , Ω ^ ˙ = Ψ T ( y ) e ( Ω ^ Ω ) ,
where J = d ( h ( x ) x ) / d x , Θ ^ and Ω ^ are estimated values of the uncertain parameters in Equation (22); K is a n × n matrix, which can be decomposed into K = K 1 + K 2 , in which K 1 and K 2 satisfy the following assumptions:
( i ) K 1 T = K 1 , ( i i ) K 2 = d i a g { b 1 , b 2 , , b n } , b i < 0 , ( i = 1 , 2 , , n ) .
From the first condition in condition (24), all real parts of the eigenvalues of K 1 are zero since it is a real antisymmetric matrix. Then, all eigenvalues of matrix K ( K = K 1 + K 2 ) must have negative real parts if condition (24) is satisfied.
From Theorem 3 in Section 3, GFPS in system (22) appears due to the controller and parameter update laws presented in Equation (23) if all eigenvalues of the matrix
K 0 0 0 I p 0 0 0 I q ,
have negative real parts, in which I p and I q are the p × p and q × q real identity matrices, respectively. Obviously, the above GFPS condition implies that all eigenvalues of matrix K should have negative real parts.
Remark 1.
The GFPS condition derived in [38] is a special condition of Theorem 3 in Section 3.
An example is given below to demonstrate that the condition (24) is not necessary for GFPS occurring in system (22) under the controller and parameter update laws in (23). Consider the hyperchaotic Chen system and hyperchaotic Lorenz system as the drive and response systems, respectively [38].
x ˙ 1 = a x 1 + a x 2 + x 4 , x ˙ 2 = d x 1 + c x 2 x 1 x 3 , x ˙ 3 = b x 3 + x 1 x 2 , x ˙ 4 = r x 4 + x 2 x 3 ,
and
y ˙ 1 = a 1 y 1 + a 1 y 2 + y 4 + u 1 , y ˙ 2 = b 1 y 1 y 2 y 1 y 3 + u 2 , y ˙ 3 = c 1 y 3 + y 1 y 2 + u 3 , y ˙ 4 = d 1 y 4 y 1 y 3 + u 4 ,
where a , b , c , d , r , a 1 , b 1 , c 1 and d 1 are uncertain parameters to be identified. The true values of the uncertain parameters of systems (26) and (27) are chosen as a = 35 , b = 3 , c = 12 , d = 7 , r = 0.08 , a 1 = 10 , b 1 = 28 , c 1 = 8 / 3 and d 1 = 1.3 . In addition, we choose the scaling function matrix as h ( x ) = d i a g { h 1 ( x 1 ) , h 2 ( x 2 ) , h 3 ( x 3 ) , h 4 ( x 4 ) } , where
h 1 ( x 1 ) = 2 sin ( x 1 + 1 ) + 1 , h 2 ( x 2 ) = 3 sin ( 2 x 2 + 3 ) + 1 , h 3 ( x 3 ) = sin ( 3 x 3 1 ) + 1 , h 4 ( x 4 ) = 2 sin ( x 4 + 4 ) + 1 ,
Then, one can find that J = d ( h ( x ) x ) / d x = d i a g { J 1 , J 2 , J 3 , J 4 } , where
J 1 = 2 cos ( x 1 + 1 ) x 1 + 2 sin ( x 1 + 1 ) + 1 , J 2 = 6 cos ( 2 x 2 3 ) x 2 3 sin ( 2 x 2 3 ) + 1 , J 3 = 3 cos ( 3 x 3 1 ) x 3 sin ( 3 x 3 1 ) + 1 , J 4 = 2 cos ( x 4 4 ) x 4 + 2 sin ( x 4 4 ) + 1 ,
According to the theoretical analysis, the controller and parameter update laws presented in (23) have the following form:
u 1 = J 1 x 4 y 4 + a ^ 1 ( y 1 y 2 ) a ^ J 1 ( x 1 x 2 ) + i = 1 4 k 1 i e i , u 2 = y 2 + y 1 y 3 J 2 x 1 x 3 b ^ 1 y 1 + d ^ J 2 x 1 + c ^ J 2 x 2 + i = 1 4 k 2 i e i , u 3 = J 3 x 1 x 2 y 1 y 2 + c ^ 1 y 3 b ^ J 3 x 3 + i = 1 4 k 3 i e i , u 4 = y 1 y 3 + J 4 x 2 x 3 d ^ 1 y 4 + r ^ J 4 x 4 + i = 1 4 k 4 i e i ,
a ^ ˙ = ( x 1 x 2 ) J 1 e 1 + ( a a ^ ) , b ^ ˙ = x 3 J 3 e 3 + ( b b ^ ) , c ^ ˙ = x 2 J 2 e 2 + ( c c ^ ) , d ^ ˙ = x 1 J 2 e 2 + ( d d ^ ) , r ^ ˙ = x 4 J 4 e 4 + ( r r ^ ) ,
a ^ ˙ 1 = ( y 2 y 1 ) e 1 + ( a 1 a ^ 1 ) , b ^ ˙ 1 = y 1 e 2 + ( b 1 b ^ 1 ) , c ^ ˙ 1 = y 3 e 3 + ( c 1 c ^ 1 ) , d ^ ˙ 1 = y 4 e 4 + ( d 1 d ^ 1 ) ,
where k i j ( i = 1 , 2 , 3 , 4 ; j = 1 , 2 , 3 , 4 ) are real constants that compose the matrix K; e i = ( y i h i x i ) , i = 1 , 2 , 3 , 4 ; a ^ , b ^ , c ^ , d ^ , r ^ , a ^ 1 , b ^ 1 , c ^ 1 and d ^ 1 are estimated values of the unknown parameters in systems (26) and (27). From the analytical result given in [38], GFPS between systems (26) and (27) will be achieved if the K matrix is chosen as
K = 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 ,
which satisfies condition (24). However, the GFPS conditions obtained in this paper for systems (26) and (27) only require that all eigenvalues of K matrix have negative real parts. Numerical simulations for systems (26), (27) are carried out below by selecting
K = 1 1 0 0 2 1 1 0 0 1 1 1 0 0 1 1 .
Remark 2.
The K matrix (34) is not antisymmetric, which does not satisfy the first condition in (24).
The four eigenvalues of K matrix (34) are 1 ± 1.85 i and 1 ± 0.77 i . The numerical results are shown in Figure 1 and Figure 2, in which the initial values of the systems (26) and (27) are taken as x 1 ( 0 ) = 3 , x 2 ( 0 ) = 5 , x 3 ( 0 ) = 7 , x 4 ( 0 ) = 8 , y 1 ( 0 ) = 12 , y 2 ( 0 ) = 2 , y 3 ( 0 ) = 3 , y 4 ( 0 ) = 5 and the estimated parameters have initial conditions a ^ ( 0 ) = 4 , b ^ ( 0 ) = 6 , c ^ ( 0 ) = 2 , d ^ ( 0 ) = 1 , r ^ ( 0 ) = 5 , a ^ 1 ( 0 ) = 1 , b ^ 1 ( 0 ) = 3 , c ^ 1 ( 0 ) = 0.3 and d ^ 1 ( 0 ) = 0.1 . Although K matrix (34) does not satisfy condition (24), Figure 1 shows that GFPS between systems (26) and (27) is achieved with the scaling function factors (28). Figure 2 illustrates that the estimated parameters a ^ , b ^ , c ^ , d ^ , r ^ , a ^ 1 , b ^ 1 , c ^ 1 , d ^ 1 in systems (26) and (27) adapt themselves to the true values a = 35 , b = 3 , c = 12 , d = 7 , r = 0.08 , a 1 = 10 , b 1 = 28 , c 1 = 8 / 3 and d 1 = 1.3 , respectively. The numerical results in Figure 1 and Figure 2 demonstrate that GFPS between systems (26) and (27) can be achieved with the scaling function factors (28) and matrix (34), and the parameter estimate errors converge to zero.
Remark 3.
The GFPS condition derived using the Laplace transform method is more generalized than that proposed in [38], which is obtained using the Lyapunov method.
To further verify the theory’s parts in more detail, another K matrix,
K = 1 2 3 4 2 1 2 3 2 3 4 4 2 3 4 4 ,
is chosen to perform the numerical simulations for systems (26) and (27) with scaling function factors (28), in which all the initial conditions remain the same.
Remark 4.
For matrix (35), neither of the two conditions in (24) is satisfied.
Matrix (35) has four eigenvalues with negative real parts, 3.4925 ± 0.6286 i , 0.5075 ± 0.2452 i , which satisfies the condition given in Theorem 1 in Section 3. The numerical results are presented in Figure 3 and Figure 4. Figure 3 illustrates that GFPS between systems (26) and (27) occurs with the scaling function factors (28). Figure 4 shows that the estimated parameters a ^ , b ^ , c ^ , d ^ , r ^ , a ^ 1 , b ^ 1 , c ^ 1 , d ^ 1 in systems (26) and (27) also can converge to their true values a = 35 , b = 3 , c = 12 , d = 7 , r = 0.08 , a 1 = 10 , b 1 = 28 , c 1 = 8 / 3 and d 1 = 1.3 , respectively. The numerical results in Figure 1, Figure 2, Figure 3 and Figure 4 demonstrate that GFPS between systems (26) and (27) can be achieved with the scaling function factors (28), and the parameter estimate errors converge to zero as long as all eigenvalues of K matrix have negative real parts, which verifies the correctness of our theoretical results in Section 3.

5. The Design of Simple Parameter Update Laws

According to the GFPS condition given in Theorem 1 in Section 3, simple parameter update laws can be designed for the estimation of the uncertain parameters in system (1). From M matrix defined in Equation (13), g 1 ( ϕ ^ ϕ , e 1 ) and g 2 ( ψ ^ ψ , e 1 ) in Equation (3) can be selected so that the eigenvalues of matrix M are independent of ϕ and ψ . Consider that g 1 ( ϕ ^ ϕ , e 1 ) and g 2 ( ψ ^ ψ , e 1 ) in Equation (3) have the following form:
g 1 ( ϕ ^ ϕ , δ 1 ) = γ 11 ( ϕ ^ ϕ ) + γ 12 Γ 1 ( δ 1 ) , g 2 ( ψ ^ ψ , e 1 ) = γ 21 ( ψ ^ ψ ) + γ 22 Γ 2 ( δ 1 ) ,
where γ 11 , γ 12 , γ 21 , γ 22 are control parameters; Γ 1 , 2 ( δ 1 ) are the functions of δ 1 , Γ 1 , 2 ( δ 1 ) δ 1 | δ 1 = 0 = 0 . With parameter update laws (36), matrix M in Equation (13) becomes
M = k I n 0 0 0 γ 11 I m 0 0 0 γ 21 I l ,
where I n , I m and I l represent n × n , m × m and l × l real identity matrices, respectively. All eigenvalues of matrix (37) have negative real parts as long as matrices k, γ 11 and γ 21 have no eigenvalues with non-negative real parts. For the sake of simplicity, one can choose all k, γ 11 and γ 21 to be scalars. Therefore, GFPS in system (1) can be achieved by the controller and parameter update laws presented in Equations (3) and (36) if k , γ 11 , γ 21 < 0 .
To illustrate the effectiveness of the parameter update laws (36), drive–response systems constructed with the Lorenz and Chen systems are selected to perform numerical simulations. The Lorenz chaotic system [49], as the drive system, can be described by the following equations:
x ˙ 1 = σ ( x 2 x 1 ) , x ˙ 2 = r x 1 x 2 x 1 x 3 , x ˙ 3 = x 1 x 2 β x 3 ,
where x 1 , x 2 , x 3 are state variables, and σ , r , β are uncertain parameters to be identified. In particular, when σ = 10 , r = 28 and β = 8 / 3 , system (38) displays a chaotic attractor.
The Chen chaotic system [50], as the response system, is given as below:
y ˙ 1 = a ( y 2 y 1 ) + u 1 , y ˙ 2 = ( c a ) y 1 y 1 y 3 + c y 2 + u 2 , y ˙ 3 = y 1 y 2 b y 3 + u 3 ,
where y 1 , y 2 , y 3 are state variables, a , b , c are uncertain parameters to be estimated and u 1 , u 2 , u 3 are the control laws to be designed. System (39) exhibits chaotic dynamics when a = 35 , b = 3 and c = 28 . The scaling function matrix between systems (38) and (39) is chosen as h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 ) } , and one has
J = d ( h ( x ) x ) d x = d i a g { 2 cos ( x 1 + 1 ) x 1 + 2 sin ( x 1 + 1 ) + 1 , 6 cos ( 2 x 2 3 ) x 2 3 sin ( 2 x 2 3 ) + 1 , 3 cos ( 3 x 3 1 ) x 3 sin ( 3 x 3 1 ) + 1 } .
According to Equations (3) and (36), the controller u 1 , u 2 , u 3 and the parameter update laws can be given by
u 1 = ( 2 cos ( x 1 + 1 ) x 1 + 2 sin ( x 1 + 1 ) + 1 ) σ ^ ( x 2 x 1 ) + k ( y 1 ( 2 sin ( x 1 + 1 ) + 1 ) x 1 ) + a ^ ( y 1 y 2 ) , u 2 = ( 6 cos ( 2 x 2 3 ) x 2 3 sin ( 2 x 2 3 ) + 1 ) ( r ^ x 1 x 2 x 1 x 3 ) + k ( y 2 ( 3 sin ( 2 x 2 + 3 ) + 1 ) x 2 ) ( c ^ a ^ ) y 1 + y 1 y 3 c ^ y 2 , u 3 = ( 3 cos ( 3 x 3 1 ) x 3 sin ( 3 x 3 1 ) + 1 ) ( x 1 x 2 β ^ x 3 ) + k ( y 3 ( sin ( 3 x 3 1 ) + 1 ) x 3 ) y 1 y 2 + b ^ y 3 ,
σ ^ ˙ = γ 11 ( σ ^ σ ) + γ 12 ( 1 c o s ( y 1 ( 2 s i n ( x 1 + 1 ) + 1 ) x 1 ) ) , r ^ ˙ = γ 11 ( r ^ r ) + γ 12 ( 1 c o s ( y 2 ( 3 s i n ( 2 x 2 + 3 ) + 1 ) x 2 ) ) , β ^ ˙ = γ 11 ( β ^ β ) + γ 12 ( 1 c o s ( y 3 ( s i n ( 3 x 3 1 ) + 1 ) x 3 ) ) ,
a ^ ˙ = γ 21 ( a ^ a ) + γ 22 ( 1 c o s ( y 1 ( 2 s i n ( x 1 + 1 ) + 1 ) x 1 ) ) , b ^ ˙ = γ 21 ( b ^ b ) + γ 22 ( 1 c o s ( y 2 ( 3 s i n ( 2 x 2 + 3 ) + 1 ) x 2 ) ) , c ^ ˙ = γ 21 ( c ^ c ) + γ 22 ( 1 c o s ( y 3 ( s i n ( 3 x 3 1 ) + 1 ) x 3 ) ) ,
where σ ^ , r ^ , β ^ , a ^ , b ^ and c ^ are estimated values of the uncertain parameters in systems (38) and (39). The true values of the uncertain parameters of systems (38) and (39) are chosen as σ = 10 , r = 28 , β = 8 / 3 and a = 35 , b = 3 , c = 28 . According to the previous analysis results in the section, if k , γ 11 , γ 21 , γ 12 , γ 22 are chosen as k = γ 11 = γ 21 = 1 and γ 12 = γ 22 = 0.01 , GFPS between uncertain Lorenz and Chen chaotic systems will be achieved under the control of the controller (40), and the uncertain parameters will be estimated using the parameter update laws (41) and (42). The initial conditions are given by x 1 ( 0 ) = 0.1 , x 2 ( 0 ) = 0.2 , x 3 ( 0 ) = 0.3 , y 1 ( 0 ) = 0.4 , y 2 ( 0 ) = 0.5 , y 3 ( 0 ) = 0.6 , σ ^ = 9 , r ^ = 25 , β ^ = 3 , a ^ = 33 , b ^ = 2 and c ^ = 29 to perform the numerical simulations. As shown in Figure 5 and Figure 6, GFPS between systems (38) and (39) is achieved, and the uncertain parameters are also identified.
If the scaling function matrix is taken as h ( x ) = d i a g { 0.02 x 1 2 + 0.2 x 1 + 1 , 0.01 x 2 2 + 0.1 x 2 1 , 0.03 x 3 2 0.1 x 3 + 0.5 } and all other parameter conditions in systems (38) and (39) remain the same, GFPS between systems (38) and (39) still can be achieved and the uncertain parameters are also identified under the controller and parameter update laws given by Equations (3) and (36), respectively. The simulation results are shown in Figure 7 and Figure 8, which illustrate the validity of the parameter update laws (36) proposed in the section. It is worth pointing out that the parameter update laws (36) are still valid even if the scaling function factors are other complex functions.

6. Conclusions

In this paper, an adaptive controller is proposed to realize GFPS between two different chaotic systems (drive–response systems) with uncertain parameters based on the Laplace transform method. Compared to previous research, first, the uncertain chaotic systems considered in this paper are more generalized and have fewer restrictions on uncertain parameters; second, the GFPS scheme proposed in this paper is more generalized, and it can be used even when the scaling function factors are periodic functions, polynomial functions or other complex functions. Third, the parameter update laws designed for the estimation of uncertain parameters based on the GFPS condition derived using the Laplace transform method are simpler.
The greatest advantage of the GFPS scheme based on the Laplace transform method is that it can be used to design simple controllers and parameter update laws without constructing Lyapunov functions. The added controllers designed by this approach are more generalized and easier to realize physically than those designed through the Lyapunov function method.
It should be noted that we only consider the local stability of the origin in the GFPS error system between two chaotic systems in this paper. Moreover, to obtain the GFPS conditions, we require the smoothness of the trajectories in two chaotic systems to guarantee the existence, uniqueness and continuousness of solutions in a class of Volterra integral equations. According to the works in the literature [51,52], the global stability of the GFPS error system between two non-smooth chaotic systems is worthy of discussion, and this will be one of our next research points.

Author Contributions

Conceptualization, B.Z.; methodology, B.Z. and Y.Z.; software, Y.Z.; validation, Y.Z.; writing—original draft preparation, B.Z. and Y.Z.; writing—review and editing, B.Z.; funding acquisition, B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 11672185.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Acknowledgments

The authors greatly appreciate the anonymous referees for their valuable suggestions and questions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pecora, L.M.; Carroll, T.L. Synchronization in chaotic systems. Phys. Rev. Lett. 1990, 64, 821–824. [Google Scholar] [CrossRef] [PubMed]
  2. Carroll, T.L.; Pecora, L.M. Synchronizing chaotic circuits. IEEE Trans. Circ. Syst. I 1991, 38, 453–456. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, G.; Dong, X. From Chaos to Order: Methodologies, Perspectives and Applications; World Scientific: Singapore, 1998. [Google Scholar]
  4. Klug, F. Synchronization transitions in supply chain networks. Int. J. Syst.-Sci.-Oper. Logist. 2023, 10, 2224104. [Google Scholar]
  5. Zhong, X.H.; Jia, Z.W.; Li, Q.T.; Wang, L.S.; Guo, Y.Y.; Sun, Y.H.; Wang, Y.C.; Wang, A.B. Multi-wavelength broadband chaos generation and synchronization using long-cavity FP lasers. IEEE J. Sel. Top. Quantum Electron. 2023, 29, 0600207. [Google Scholar] [CrossRef]
  6. Khezeli, M.; Najafi, E.; Molana, M.H.; Seidi, M. A sustainable and resilient supply chain (RS-SCM) by using synchronisation and load-sharing approach: Application in the oil and gas refinery. Int. J. Syst.-Sci.-Oper. Logist. 2023, 10, 2198055. [Google Scholar] [CrossRef]
  7. Li, J.; Li, H.L.; Yang, J.P.; Yang, J.K.; Zhang, L. Hybrid control-based synchronization of fractional-order delayed complex-valued fuzzy neural networks. Comput. Appl. Math. 2023, 42, 154. [Google Scholar] [CrossRef]
  8. Rosenblum, M.G.; Pikovsky, A.S.; Kurths, J. Phase synchronization of chaotic oscillators. Phys. Rev. Lett. 1996, 76, 1804–1807. [Google Scholar] [CrossRef] [Green Version]
  9. Haerter, P.; Viana, R.L. Synchronization of phase oscillators due to nonlocal coupling mediated by the slow diffusion of a substance. Braz. J. Phys. 2023, 53, 114. [Google Scholar] [CrossRef]
  10. Cao, L.; Lai, Y. Antiphase synchronism in chaotic systems. Phys. Rev. E 1998, 58, 382–386. [Google Scholar] [CrossRef] [Green Version]
  11. Thazhathethil, R.; Pallimanhiyil, A.S. In-phase and anti-phase bursting dynamics and synchronisation scenario in neural network by varying coupling phase. J. Biol. Phys. 2023. [Google Scholar] [CrossRef]
  12. Rosenblum, M.G.; Pikovsky, A.S.; Kurths, J. From phase to lag synchronization in coupled chaotic oscillators. Phys. Rev. Lett. 1997, 78, 4193–4296. [Google Scholar] [CrossRef]
  13. Kumar, V.; Heiland, J.; Benner, P. Exponential lag synchronization of Cohen-Grossberg neural networks with discrete and distributed delays on time scales. Neural Process. Lett. 2023. [Google Scholar] [CrossRef]
  14. Kocarev, L.; Parlitz, U. Generalized synchronization, predictability, and equivalence of unidirectionally coupled dynamical systems. Phys. Rev. Lett. 1996, 76, 181–1819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Xing, Y.A.; Dong, W.J.; Zeng, J.; Guo, P.T.; Zhang, J.; Ding, Q. Study of generalized chaotic synchronization method incorporating error-feedback coefficients. Entropy 2023, 25, 818. [Google Scholar] [CrossRef]
  16. Mainieri, R.; Rehacek, J. Projective synchronization in the threedimensional chaotic systems. Phys. Rev. Lett. 1999, 82, 3042–3045. [Google Scholar] [CrossRef]
  17. Han, J.; Chen, G.C.; Wang, L.M.; Zhang, G.D.; Hu, J.H. Direct approach on fixed-time stabilization and projective synchronization of inertial neural networks with mixed delays. Neurocomputing 2023, 535, 97–106. [Google Scholar] [CrossRef]
  18. Chee, C.Y.; Xu, D.L. Chaos-based M-nary digital communication technique using controlled projective synchronization. IEEE Proc. Circuits Devices Syst. 2006, 153, 357–360. [Google Scholar] [CrossRef]
  19. Zhang, F.F.; Leng, S.; Li, Z.F.; Jiang, C.M. Time delay complex Chen chaotic system and secure communication scheme for wireless body area networks. Entropy 2020, 22, 1420. [Google Scholar] [CrossRef]
  20. Baluni, S.; Yadav, V.K.; Das, S. Quasi projective synchronization of time varying delayed complex valued Cohen-Grossberg neural networks. Inf. Sci. 2022, 612, 231–240. [Google Scholar] [CrossRef]
  21. Xu, D.L. Control of projective synchronization in chaotic systems. Phys. Rev. E. 2001, 63, 027201. [Google Scholar] [CrossRef]
  22. Xin, B.G.; Wu, Z.H. Projective synchronization of chaotic discrete dynamical systems via linear state error feedback control. Entropy 2015, 17, 2677–2687. [Google Scholar] [CrossRef] [Green Version]
  23. Wen, G.; Xu, D. Nonlinear observer control for full-state projective synchronization in chaotic continuous-time system. Chaos Solitons Fractals 2005, 26, 71–77. [Google Scholar] [CrossRef]
  24. Yan, J.; Li, C. Generalized projective synchronization of a unified chaotic system. Chaos Solitons Fractals 2005, 26, 1119–1124. [Google Scholar] [CrossRef]
  25. Zhou, Y.; Wang, H.; Liu, H. Generalized function projective synchronization of incommensurate fractional-order chaotic systems with inputs saturation. Int. J. Fuzzy Syst. 2019, 21, 823–836. [Google Scholar] [CrossRef]
  26. Zhang, L.; Wang, Y.; Wang, Q.; Qiao, S.; Wang, F. Generalized projective synchronization for networks with one crucial node and different dimensional nodes via a single controller. Asian J. Control 2020, 24, 1471–1483. [Google Scholar] [CrossRef]
  27. Li, X.; Zhao, X.; Lu, Y.; Xu, T.; Hou, P.; Cui, X. Generalized projection synchronization and generalized dislocation synchronization of fractional-order financial chaotic systems with different orders. In Proceedings of the International Symposium on Electrical, Electronics and Information Engineering, Soul, Republic of Korea, 19–21 February 2021; pp. 635–638. [Google Scholar]
  28. Li, G.H. Modified projective synchronization of chaotic system. Chaos Solitons Fractals 2007, 32, 1786–1790. [Google Scholar] [CrossRef]
  29. He, J.M.; Lei, T.F.; Chen, F.Q. Global matrix projective synchronization of delayed fractional-order neural networks. Soft Comput. 2023, 27, 8991–9000. [Google Scholar] [CrossRef]
  30. Feng, C.F.; Tan, Y.R.; Wang, Y.H.; Yang, H.J. Active backstepping control of combined projective synchronization among different nonlinear systems. Automatika 2017, 58, 295–301. [Google Scholar] [CrossRef]
  31. Chen, Y.; Li, X. Function projective synchronization between two identical chaotic systems. Int. J. Mod. Phys. C 2007, 18, 883–888. [Google Scholar] [CrossRef]
  32. Chen, Y.; An, H.L. The function cascade synchronization approach with uncertain parameters or not for hyperchaotic systems. Appl. Math. Comput. 2008, 197, 96–110. [Google Scholar] [CrossRef] [Green Version]
  33. El-Dessoky, M.M.; Alzahrani, E.O.; Almohammadi, N.A. Chaos control and function projective synchronization of noval chaotic dynamical system. J. Comput. Anal. Appl. 2019, 26, 1. [Google Scholar]
  34. Bekiros, S.; Yao, Q.J.; Mou, J.; Alkhateeb, A.F.; Jahanshahi, H. Adaptive fixed-time robust control for function projective 332 synchronization of hyperchaotic economic systems with external perturbations. Chaos Solitons Fractals 2023, 172, 113609. [Google Scholar] [CrossRef]
  35. Du, H.Y.; Zeng, Q.S.; Wang, C.H. Modified function projective synchronization of chaotic system. Chaos Solitons Fractals 2009, 197, 2399–2404. [Google Scholar] [CrossRef]
  36. Srisuntorn, R.; Weera, W.; Botmart, T. Modified function projective synchronization of master-slave neural networks with mixed interval time-varying delays via intermittent feedback control. AIMS Math. 2022, 7, 18632–18661. [Google Scholar] [CrossRef]
  37. Yu, Y.; Li, H.X. Adaptive generalized function projective synchronization of uncertain chaotic systems. Nonlinear Anal. RWA 2010, 11, 2456–2464. [Google Scholar] [CrossRef]
  38. Li, Z.B.; Zhao, X.S. Generalized function projective synchronization of two different hyperchaotic systems with unknown parameters. Nonlinear Anal. RWA 2011, 12, 2607–2615. [Google Scholar] [CrossRef]
  39. Wu, X.J.; Fu, Z.Y.; Kurths, J. A secure communication scheme based generalized function projective synchronization of a new 5D hyperchaotic system. Phys. Scr. 2015, 90, 045210. [Google Scholar] [CrossRef] [Green Version]
  40. Al-Azzawi, S.F.; Al-Talib, Z.S. Generalized function projective synchronization via nonlinear controller strategy. J. Interdiscip. Math. 2022, 25, 1753–1765. [Google Scholar] [CrossRef]
  41. Zhou, X.B.; Xiong, L.L.; Cai, X.M. Adaptive Switched Generalized Function Projective Synchronization between Two Hyperchaotic Systems with Unknown Parameters. Entropy 2013, 16, 377–388. [Google Scholar] [CrossRef]
  42. Sun, J.W.; Guo, J.C.; Zhang, X.C. Adaptive generalized hybrid function projective dislocated synchronization of new four- dimensional uncertain chaotic systems. Appl. Math. Comput. 2015, 252, 304–314. [Google Scholar] [CrossRef]
  43. Hamel, S.; Boulkroune, A. A generalized function projective synchronization sheme for uncertain chaotic systems subject to input nonlinearities. Int. J. Gen. Syst. 2016, 45, 689–710. [Google Scholar] [CrossRef]
  44. Li, Z.B.; Tang, J.S. Generalized binary function projective synchronization of chaotic systems with unknown parameters. Optik 2017, 137, 101–107. [Google Scholar] [CrossRef]
  45. Li, Q.P.; Liu, S.Y.; Chen, Y.G. Finite-time adaptive modified function projective multi-lag generalized compound synchronization for multiple uncertain chaotic systems. Int. J. Appl. Math. Comput. Sci. 2018, 28, 613–624. [Google Scholar] [CrossRef] [Green Version]
  46. Abarbanel, H.D.I.; Rulkov, N.F.; Sushchik, M.M. Gneralized synchronization of chaos: The auxiliary system approach. Phys. Rev. E 1996, 53, 4528–4535. [Google Scholar] [CrossRef] [Green Version]
  47. Nohel, J.A. Some problems in nonlinear Volterra integral equations. Bull. Am. Math. Soc. 1962, 68, 323–330. [Google Scholar] [CrossRef] [Green Version]
  48. Kyrchei, I. Cramer’s rule for generalized inverse solutions. In Advances in Linear Algebra Research; Kyrchei, I., Ed.; Nova Science Publishers: Hauppauge, NY, USA, 2015; pp. 79–132. [Google Scholar]
  49. Lorenz, E.N. Deterministic nonperiodic flow. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  50. Chen, G.; Ueta, T. Yet another chaotic attrator. Int. J. Bifurc. Chaos 1999, 9, 1465–1466. [Google Scholar] [CrossRef]
  51. Zhang, T.W.; Xiong, L.L. Periodic motion for impulsive fractional functional differential equations with piecewise Caputo derivative. Appl. Math. Lett. 2020, 101, 106072. [Google Scholar] [CrossRef]
  52. Liu, Y.T. Global exponential stability and synchronization of discrete-time fuzzy bidirectional associative memory neural networks via Mittag-Leffler difference approach. Int. J. Fuzzy Syst. 2023, 25, 1922–1934. [Google Scholar] [CrossRef]
Figure 1. GFPS between systems (26) and (27) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (34).
Figure 1. GFPS between systems (26) and (27) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (34).
Applsci 13 08135 g001
Figure 2. The estimates of the uncertain parameters of systems (26) and (27) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (34). (a) uncertain parameters in system (26) (b) uncertain parameters in system (27).
Figure 2. The estimates of the uncertain parameters of systems (26) and (27) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (34). (a) uncertain parameters in system (26) (b) uncertain parameters in system (27).
Applsci 13 08135 g002
Figure 3. GFPS between systems (26) and (27) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (35).
Figure 3. GFPS between systems (26) and (27) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (35).
Applsci 13 08135 g003
Figure 4. The estimates of the uncertain parameters of systems (26) and (27) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (35). (a) uncertain parameters in system (26) (b) uncertain parameters in system (27).
Figure 4. The estimates of the uncertain parameters of systems (26) and (27) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 , 2 s i n ( x 4 + 4 ) + 1 } and K matrix (35). (a) uncertain parameters in system (26) (b) uncertain parameters in system (27).
Applsci 13 08135 g004
Figure 5. GFPS between systems (38) and (39) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 ) } . (a) ( 2 s i n ( x 1 + 1 ) + 1 ) x 1 vs. y 1 (b) ( 3 s i n ( 2 x 2 + 3 ) + 1 ) x 2 vs. y 2 (c) ( s i n ( 3 x 3 1 ) + 1 ) ) x 3 vs. y 3 .
Figure 5. GFPS between systems (38) and (39) is achieved with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 ) } . (a) ( 2 s i n ( x 1 + 1 ) + 1 ) x 1 vs. y 1 (b) ( 3 s i n ( 2 x 2 + 3 ) + 1 ) x 2 vs. y 2 (c) ( s i n ( 3 x 3 1 ) + 1 ) ) x 3 vs. y 3 .
Applsci 13 08135 g005
Figure 6. The estimates of the uncertain parameters of systems (38) and (39) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 ) } . (a) uncertain parameters in system (38) (b) uncertain parameters in system (39).
Figure 6. The estimates of the uncertain parameters of systems (38) and (39) for GFPS with the scaling function matrix h ( x ) = d i a g { 2 s i n ( x 1 + 1 ) + 1 , 3 s i n ( 2 x 2 + 3 ) + 1 , s i n ( 3 x 3 1 ) + 1 ) } . (a) uncertain parameters in system (38) (b) uncertain parameters in system (39).
Applsci 13 08135 g006
Figure 7. GFPS between systems (38) and (39) is achieved with the scaling function matrix h ( x ) = d i a g { 0.02 x 1 2 + 0.2 x 1 + 1 , 0.01 x 2 2 + 0.1 x 2 1 , 0.03 x 3 2 0.1 x 3 + 0.5 ) } . (a) ( 0.02 x 1 2 + 0.2 x 1 + 1 ) x 1 vs. y 1 (b) ( 0.01 x 2 2 + 0.1 x 2 1 ) x 2 vs. y 2 (c) ( 0.03 x 3 2 0.1 x 3 + 0.5 ) x 3 vs. y 3 .
Figure 7. GFPS between systems (38) and (39) is achieved with the scaling function matrix h ( x ) = d i a g { 0.02 x 1 2 + 0.2 x 1 + 1 , 0.01 x 2 2 + 0.1 x 2 1 , 0.03 x 3 2 0.1 x 3 + 0.5 ) } . (a) ( 0.02 x 1 2 + 0.2 x 1 + 1 ) x 1 vs. y 1 (b) ( 0.01 x 2 2 + 0.1 x 2 1 ) x 2 vs. y 2 (c) ( 0.03 x 3 2 0.1 x 3 + 0.5 ) x 3 vs. y 3 .
Applsci 13 08135 g007
Figure 8. The estimates of the uncertain parameters of systems (38) and (39) for GFPS with the scaling function matrix h ( x ) = d i a g { 0.02 x 1 2 + 0.2 x 1 + 1 , 0.01 x 2 2 + 0.1 x 2 1 , 0.03 x 3 2 0.1 x 3 + 0.5 } . (a) uncertain parameters in system (38) (b) uncertain parameters in system (39).
Figure 8. The estimates of the uncertain parameters of systems (38) and (39) for GFPS with the scaling function matrix h ( x ) = d i a g { 0.02 x 1 2 + 0.2 x 1 + 1 , 0.01 x 2 2 + 0.1 x 2 1 , 0.03 x 3 2 0.1 x 3 + 0.5 } . (a) uncertain parameters in system (38) (b) uncertain parameters in system (39).
Applsci 13 08135 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhen, B.; Zhang, Y. Generalized Function Projective Synchronization of Two Different Chaotic Systems with Uncertain Parameters. Appl. Sci. 2023, 13, 8135. https://doi.org/10.3390/app13148135

AMA Style

Zhen B, Zhang Y. Generalized Function Projective Synchronization of Two Different Chaotic Systems with Uncertain Parameters. Applied Sciences. 2023; 13(14):8135. https://doi.org/10.3390/app13148135

Chicago/Turabian Style

Zhen, Bin, and Yu Zhang. 2023. "Generalized Function Projective Synchronization of Two Different Chaotic Systems with Uncertain Parameters" Applied Sciences 13, no. 14: 8135. https://doi.org/10.3390/app13148135

APA Style

Zhen, B., & Zhang, Y. (2023). Generalized Function Projective Synchronization of Two Different Chaotic Systems with Uncertain Parameters. Applied Sciences, 13(14), 8135. https://doi.org/10.3390/app13148135

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