Chaos control of a multi-dimensional chaotic mapping system by modified stability transformation method

Jilei Zhou1 , Canchang Liu2 , Chuanbo Ren3

1, 2, 3School of Transportation and Vehicle Engineering, Shandong University of Technology, Zibo, China

3Corresponding author

Journal of Vibroengineering, Vol. 19, Issue 2, 2017, p. 1103-1115. https://doi.org/10.21595/jve.2016.17317
Received 23 June 2016; received in revised form 21 October 2016; accepted 11 November 2016; published 31 March 2017

Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

For the multi-dimensional chaotic mapping system with various forms of bifurcations, it is found that the spectral radius around the static bifurcation point is approximately equal to one and the convergence speed for traditional stability transformation method (STM) which is used to stabilize the unstable fixed points is fairly slow. In this paper, a modified STM is proposed to efficiently stabilize a 3D chaotic mapping system to stable fixed points. Firstly, according to the information of the fixed point, the stability matrix C is derived, demonstrating that it is unnecessarily an involutory matrix for STM that can also stabilize chaotic system to the fixed point. Then, the critical parameter qcrit that satisfies convergence condition and the optimal parameter qopt corresponding to best performance for STM are determined respectively. Moreover, STM is combined with Newton method (NT) to overcome the disadvantage of slow convergence around the static bifurcation point without requiring a priori information of fixed point. It is indicated that the number of iterations, the absolute and relative errors between the convergent value and analytical fixed point for combined STM-NT decrease enormously comparing with that of traditional STM. Finally, numerical analysis verifies the high efficiency of modified STM proposed in this paper.

Keywords: chaotic mapping system, modified stability transformation method, relaxes parameter, chaos control.

1. Introduction

Chaos is an intrinsically complex phenomenon for deterministic nonlinear dynamical system. Generally, there are two complementary attributes to define chaos in both the temporal and spatial aspects: exponentially sensitive dependence on initial conditions, and pattern structure of strange attractors [1]. Due to its theoretical significance and engineering applicability, the chaos control has been paid extensive attention in diverse fields. Chaos control can catch the desired fixed points or periodic orbits embedded in the chaotic attractor of nonlinear dynamical system through implementing the target guidance and position. A key issue for chaos control is how to determine the control variable varying with time in terms of the past or present state, and incorporate it into a feedback term to achieve the purpose of controlling unstable periodic orbits (UPOs). Therefore, it is necessary to detect the UPOs of chaotic system in advance. However, detecting the UPOs of nonlinear systems with both numerical calculation and experimental implementation is one of the most challenging problems for nonlinear dynamics. Owing to the inherent unstable nature of UPOs, the analytic expressions for every periodic orbit cannot be obtained, only if the chaotic system is a simple nonlinear polynomial with low degree and the period of concern is low.

The most popular algorithm for detecting UPOs is Newton method (NT) [2], which has the second order convergence speed, but it is heavily dependent on the initial condition, as well as on the behavior of the partial derivatives. Further, the extended algorithms based on Newton method including relaxed Newton method [2, 3], parameter adjustment hybrid method [4] etc. have also local convergence. To overcome this limitation, the NT is combined with a global convergent method named stability transformation method (STM) to improve the convergent efficiency in this paper. From another perspective, Pyragas, et al. proposed and subsequently developed time-delayed feedback control (DFC) in a series of papers [5-7] to detect UPOs by stabilizing continuous chaotic system to period orbit. In addition, Parsopoulos et al. [8, 9] applied the stochastic optimization method (say, particles warm optimization) to locate chaotic system to high periodic orbits. Because this algorithm does not need the detailed information of initial condition and partial derivatives, it is an efficient alternative.

In order to locate chaotic system (including both discrete and continuous system) to an arbitrary UPO, Schmelcher and Diakonos [10-12] proposed the stability transformation method (STM), which does not require the specific information of target UPOs and has a global convergence. Meanwhile, it is worth noting that Kawai and Tsubone [18] developed a new control algorithm based on STM and DFC to control the UPOs by a simple circuit implementation. The author and partner [13] applied STM to stabilize a novel 3D autonomous chaotic system to multi-equilibrium. It is found that STM can be regarded as a special form of speed feedback control method, which facilitates the practical implementation of STM. In the subsequent paper [14], we further evaluated and compared the generality, efficiency and applicability of several popular chaos control approaches, and revealed the essential connections among these approaches. Recently, STM has been utilized to control the nonlinear iterative solution of First Order Reliability Method (FORM) [15-17]. Meng and Li [18] further expanded STM to the reliability-based design optimization which indicates that STM can greatly improve the optimization efficiency. Although STM has been widely used in chaos control, structure reliability and optimization design, the convergence speed is very slow when stabilized fixed points are around the static bifurcation points of chaotic system. In addition, the choices of stability matrix and relax parameter have direct influence on the performance of STM, and how to determine the optimal relax parameter and stability matrix that make STM be the fastest convergence speed remains to be further studied. Therefore, in this paper a novel modified STM is proposed to efficiently stabilize a 3D chaotic mapping system with various forms of bifurcations to stable fixed points.

Focusing on these problems aforementioned, this paper provides some new insights on stabilizing and controlling the unstable periodic fixed points of a 3D chaotic mapping system with various forms of bifurcations by modified STM. According to the information of unstable periodic fixed orbits, the stability matrix C, critical relax parameter qcrit and optimal relax parameter qopt are derived respectively when STM satisfies given convergence condition. In order to overcome the disadvantage of slow convergence around the static bifurcation point, this paper combines STM with NT to stabilize chaotic system to any unstable periodic fixed point. Finally, the efficiency of improved STM is verified by numerical simulation. The aim of this work is to deepen the understanding for STM and promote its practical application.

2. Review of stability transformation method

In order to detect any UPO of chaotic system, Schmelcher and Diakonos [12] suggested an appropriate linear transformation to modify the eigenvalues of Jacobian matrix of dynamical system, and stabilize the original system to the unstable periodic fixed points without changing their value and location. The detailed introduction about STM can be referred to the literature [12].

The expression of STM for controlled system is written as:

(1)
X k + 1 = X k + q C F X k - X k ,

in which, q is a relaxation parameter that has closely relationship with the eigenvalue of Jacobi matrix of original chaotic system and 0 <q< 1, and the stability matrix C here is taken as n×n dimensional involutory matrix (namely, only one element in each row and each column in this matrix is 1 or –1, and the others is 0). Nevertheless, the later analysis in Section 5 indicates that non-involutory matrix can also achieve the objective of detecting unstable periodic fixed point. According to feedback control strategy, the corresponding block diagram for STM is displayed in Fig. 1.

From Eq. (1), it is easy to know that the STM neither requires the detail information of the target periodic fixed points nor calculates the Jacobian matrix of the system, namely the local information on the derivative of periodic fixed points. Therefore, it is a global convergent method and has low requirement for initial condition.

Fig. 1. Control block diagram of STM

Control block diagram of STM

3. Stabilize unstable periodic fixed points of chaotic mapping system by STM

Consider the following 3D mapping system [19]:

(2)
x k + 1 = a 1 x k + a 2 y k + a 3 y k 2 , y k + 1 = b 1 - b 2 z k , z k + 1 = c x k .

where a1= 0.05, a2= 0.25, a3= 0.12, b1= 4.0, c= 2.15, 0.4 b2 1.2. Practically, the mapping system Eq. (2) can be viewed as a kind of nonlinear iterative procedure expressed as Xk+1=F(Xk), which is adopted in many iterative algorithms for structural reliability assessment [15-17] and optimal deign [18]. The mapping system Eq. (2) does not always converge to some fixed point, and yield the phenomena of periodic oscillation, bifurcation and chaos with the control parameter b2 varying in the interval [0.4, 1.2], which is indicated in Fig. 2(a). In order to clearly verify the type of the attractor, the three Lyapunov exponents (named L1, L2 and L3) of the mapping system are also displayed in Fig. 2(b). It is demonstrated from Fig. 2(b) that when the control parameter b2< 0.544, L1, L2 and L3 are less than zero, then the mapping system is stable and convergent to a fixed point. When the parameter b2 [0.544, 0.683), the Lyapunov exponents L1= 0, L2< 0 and L3< 0, the attractor of the mapping system will change from a fixed point to a series limit cycles and then emerge a two-dimensional torus. However, the attractor will become a chaos while L1> 0, L2< 0 and L3< 0, especially become a super chaos with more than two Lyapunov exponents greater than zero.

Fig. 2. a) Bifurcation plot of component y versus parameter b2; b) Lyapunov exponents versus b2 for three-dimensional mapping system

a) Bifurcation plot of component y versus parameter b2;  b) Lyapunov exponents versus b2 for three-dimensional mapping system

a)

a) Bifurcation plot of component y versus parameter b2;  b) Lyapunov exponents versus b2 for three-dimensional mapping system

b)

3.1. Periodic-p fixed points and bifurcation forms

Let:

(3)
X k + p = x k + p y k + p z k + p ,           F X k = a 1 x k + a 2 y k + a 3 y k 2 b 1 - b 2 z k c x k .

The periodic-p iterative mapping system can be written as:

(4)
X k + p = F p X k .

If Eq. (4) has fixed points, let X=Xk+p=Xk, and substitute it to Eq. (3) and Eq. (4) respectively, and then all the analytical solutions of periodic-p fixed points are obtained. Consequently, the explicit expression of the periodic-1fixed point with respect to parameter b2 can be derived:

y ( 1 ) = a 1 - 1 - a 2 b 2 c + b 2 c ( a 2 - a 1 - 1 b 2 c ) 2 - 4 a 3 b 1 ( a 1 - 1 ) b 2 c 2 a 3 b 2 c ,
x ( 1 ) = - y 1 - b 1 b 2 c ,           z 1 = - y 1 - b 1 b 2 .

Similarly, the implicit expression of the periodic-2 fixed points with respect to b2 is also achieved:

(5)
A b 2 ( y 2 ) 4 + B b 2 ( y 2 ) 3 + C b 2 ( y 2 ) 2 + D b 2 y 2 + E b 2 = 0 ,
x ( 2 ) = b 1 - y ( 2 ) b 2 c ,           z ( 2 ) = a 3 c ( y 2 ) 2 + a 2 c - a 1 b 2 y 2 + a 1 b 1 b 2 ,

in which:

A b 2 = c 3 a 3 3 b 2 3 , B b 2 = 2 c 3 a 3 2 a 1 b 2 2 ,
C ( b 2 ) = c 3 a 3 a 2 2 b 2 3 - 2 c 2 a 3 2 b 1 b 2 2 - 2 c 2 a 3 a 2 a 1 b 2 2 + 2 c 2 a 3 2 a 1 b 1 b 2 2
            - c 2 a 3 a 2 b 2 2 + c a 3 a 1 2 + c a 3 a 1 b 2 ,
D b 2 = 2 c 2 a 3 a 2 b 1 a 1 - 1 b 2 2 - c 2 a 2 2 b 2 2 + 2 c a 1 a 2 b 2 - 2 c a 3 a 1 2 b 1 b 2 + 2 c a 3 a 1 b 1 b 2 + 1
            - a 1 2 E b 2 = c a 3 b 1 2 b 2 - 2 c a 3 a 1 b 1 2 b 2 + c a 3 a 1 2 b 1 2 b 2 - c a 2 a 1 b 1 b 2 + c a 2 b 1 b 2 + a 1 2 b 1 - b 1 .

The Jacobi matrix of mapping system at the fixed point Xk is expressed as:

(6)
J = F ( X k ) X k = x k + 1 x k x k + 1 y k x k + 1 z k y k + 1 x k y k + 1 y k y k + 1 z k z k + 1 x k z k + 1 y k z k + 1 z k = a 1 a 2 + 2 a 3 y k 0 0 0 - b 2 c 0 0 .

The eigenvalues λ1, λ2, λ3 of Jacobi matrix for mapping system with parameter b2 can be obtained according to Eq. (6), which is displayed in Fig. 3. It is indicated from Fig. 3 that the module of eigenvalues are respectively equal or greater than one at the parameter b2= 0.544 and 0.605, and then the mapping system will emerge bifurcations. When b2= 0.544 (see Fig. 2(a)), the eigenvalues of Jacobi matrix at periodic-1 fixed point are gained:

λ 1 = 0.5113 + 0.8563 i ,           λ 2 = 0.5113 - 0.8563 i ,           λ 3 = - 0.9727 .

Then it is easy to get λ1λ2=0.99781, and the two conjugate complex roots of eigenvalues will pass through unit circle as shown in Fig. 4, which implies that the mapping system produces Naimark-Sacker quasi-period bifurcation (see Fig. 2). Similarly, the eigenvalues of Jacobi matrix at periodic-1fixed point are obtained when b2= 0.605:

λ 1 = 0.5249 + 0.8799 i ,           λ 2 = 0.5249 - 0.8799 i ,           λ 3 = - 0.9999 .

Hence, the negative real eigenvalue λ3 will pass through unit circle, demonstrating that the mapping system presents period-doubling bifurcation. Meanwhile, the solid line in Fig. 5 displays analytical solutions of periodic-2 fixed points for mapping system, which indicates that b2= 0.605 is the static bifurcation point of the mapping system.

Fig. 3. Eigenvalue curves of mapping system

 Eigenvalue curves  of mapping system

Fig. 4. The forms of eigenvalues passing through unit circle

 The forms of eigenvalues passing  through unit circle

3.2. Stabilize unstable periodic fixed points by STM

With regard to periodic-p fixed points of original mapping system, the control form of STM can be formulated as:

(7)
X k + p = X k + q C F p X k - X k .

Taking periodic-2 fixed points as example, all the analytical solutions of periodic-2 fixed points for mapping system are obtained from Section 3.1, as solid lines shown in Fig. 5. It is indicated that the number of analytical solutions for periodic-2 fixed point change from one to three with the parameter b2, i.e. static bifurcation emerges. In order to detect all the periodic-2 points with STM, it is required to choose different involutory matrix C (i.e. the different stability transformation systems) and initial points in the basin of chaotic attractor. The system parameter b2 changes in the interval [0.4, 1.2]. As for every stability transformation system with appropriate relaxation parameter q, an initial point is suitably selected in the chaotic attractor and the STM is performed. Then the Nth iterative solution is compared with the exact periodic-2 point, as dashed lines displayed in Fig. 5. From Fig. 5, some interesting conclusions are obtained as follows:

1) The required iterative number N will decrease with relax parameter q, but less unstable fixed points can be detectable by STM. Taking the stability matrix C1=I and initial condition X02(–3,1,0), when q= 0.1, the N1= 2000th iterative value of y component always converge to the upper branch of analytical fixed points (see red dashed line in Fig. 5(a)). However, while q= 0.6, the N2= 500th iterative value of y component jump from the upper branch to the lower one with system parameter b2 increasing in the interval [0.4, 1.2]. Especially, when b2> 1.144, the iterative value of y component will generate oscillation and could not converge to the corresponding analytical fixed points (see red dashed line in Fig. 5(b)).

2) Keeping the initial points unchanged, STM can stabilize the mapping system to the same or distinct analytical fixed points by choosing different stability matrix. Taking the initial point X02(–3, 1, 0) and q= 0.6, when the stability matrix C1=I, the iterative value of y component converges to the upper branch of analytical fixed points with b2 varying in the interval [0.4, 0.98], but to the lower branch with b2 varying in the interval [0.98, 1.144] (see the red dashed line in Fig. 5(b)). However, when the stability matrix C3= [0 1 0; –1 0 0; 0 0 1], the iterative value of y component converges to the lower branch of analytical fixed points with b2 varying in the whole interval [0.4, 1.2] (see the blue dot dashed line in Fig. 5(b)).

3) Keeping stability matrix C unchanged, STM can stabilize the mapping system to distinct analytical fixed points by choosing different initial points. Adopting stability matrix C1=I and relax parameter q= 0.1, when the initial points are taken as X01(0,0,0) and X02(–3,1,0), the iterative value of y component respectively converges to the lower and upper branch of analytical fixed points with b2 varying in the interval [0.4, 1.2] (see blue dotted and red dashed lines in Fig. 5(a)).

Fig. 5. STM control of period-2 fixed points for mapping system when b2 varying in the interval [0.4, 1.2] a) C1=I, C2= [0 0 0;0 1 0;1 0 0], X01= (0,0,0); b) C3= [0 1 0; –1 0 0;0 0 1], X02= (–3,1,0)

 STM control of period-2 fixed points for mapping system when b2 varying in the interval [0.4, 1.2] a) C1=I, C2= [0 0 0;0 1 0;1 0 0], X01= (0,0,0); b) C3= [0 1 0; –1 0 0;0 0 1], X02= (–3,1,0)
 STM control of period-2 fixed points for mapping system when b2 varying in the interval [0.4, 1.2] a) C1=I, C2= [0 0 0;0 1 0;1 0 0], X01= (0,0,0); b) C3= [0 1 0; –1 0 0;0 0 1], X02= (–3,1,0)

4. Convergence analysis for STM

The essence of STM is to introduce an appropriate linear transformation to modify the eigenvalues of Jacobian matrix for original chaotic system, and stabilize to the unstable periodic fixed points without changing their value and location. The Jacobi matrix J* of STM for transformed system is expressed as:

(8)
J * = I + q C J I ,

in which J is the Jacobi matrix of the original mapping system. As for periodic-2 mapping system, the Jacobi matrix J can be attained:

(9)
J = F ( 2 ) ( X k ) X k = a 1 2 a 1 a 2 + 2 a 1 a 3 y k - a 2 b 2 - 2 a 3 b 1 b 2 + 2 a 3 b 2 2 z k - b 2 c 0 0 c a 1 c a 2 + 2 c a 3 y k 0 .

It is known from Eq. (8) that, if the spectral radius of matrix J* is less than 1 (i.e., ρJ*< 1), STM will converge to periodic fixed point of original system. However, in contrast if the spectral radius of matrix J* is more than 1 (i.e., ρJ*> 1), STM does not converge and will generate undesired oscillation. Fig. 6 demonstrates three eigenvalue curves of Jacobian matrix for different stability transformation systems which coincides well with the convergence property of the curves in Fig. 5. In addition, it is found that when the parameter b2= 0.605, the spectral radius of the Jacobi matrix for STM is always close to 1 (i.e., ρJ* 1) no matter how the parameter q and matrix C are selected. This implies that b2= 0.605 is a critical parameter that corresponds to the static bifurcation point of original mapping system.

Fig. 6. Eigenvalue curves of stability transformation system when a) C=C1, q=q3; b) C=C2, q=q2; c) C=C3, q=q3

 Eigenvalue curves of stability transformation system  when a) C=C1, q=q3; b) C=C2, q=q2; c) C=C3, q=q3

a)

 Eigenvalue curves of stability transformation system  when a) C=C1, q=q3; b) C=C2, q=q2; c) C=C3, q=q3

b)

 Eigenvalue curves of stability transformation system  when a) C=C1, q=q3; b) C=C2, q=q2; c) C=C3, q=q3

c)

When the system parameter b2= 0.605 (i.e. the static bifurcation point of original system), the spectral radius of stability transformation system ρ(J(X*))ρ(J*(X*))1 and therefore the convergence speed of STM is very slow around the static bifurcation point. It is supposed that the iteration of STM stops when the absolute error of sequential iterative values satisfies Xk+1-Xk 10-5. Fig. 7(a) displays the required iterative number when STM satisfies convergence condition with different relax parameter q. It is indicated from Fig. 7(a) that the maximum iterative number by STM occurs around the static bifurcation point b2= 0.605, which is consistent with the theoretical analysis. Furthermore, the smaller the parameter q is, the more iterative number required by STM is. When q 0, the required iterative number N. In addition, Fig. 7(b) and Fig. 7(c) respectively demonstrate the absolute errors and relative errors between iterative solution and periodic fixed point when STM (initial point X02 and stability matrix C3) satisfies convergence condition (≤ 10-5) with relax parameter q= 0.01. It is observed that the maximum absolute and relative errors happen at the bifurcation point b2= 0.605. Specifically, when the relax parameter q= 0.01, the maximum absolute error Xk+1-X*= 0.2337×10-5 and the maximum relative errorXk+1-X*X*×100 %= 10.17 %. Therefore, it is worth noting that detecting unstable periodic fixed points of chaotic system requires careful concern on the bifurcation point owing to the slow convergence of STM.

Fig. 7. a) Iterative numbers; b) absolute error between iterative solution and fixed point; c) relative error between iterative solution and fixed point for STM with varying q

 a) Iterative numbers; b) absolute error between iterative solution and fixed point;  c) relative error between iterative solution and fixed point for STM with varying q

a)

 a) Iterative numbers; b) absolute error between iterative solution and fixed point;  c) relative error between iterative solution and fixed point for STM with varying q

b)

 a) Iterative numbers; b) absolute error between iterative solution and fixed point;  c) relative error between iterative solution and fixed point for STM with varying q

c)

5. Determination of stability matrix C and relaxation parameter q

5.1. Choice of stability matrix C

With regarding to the chaotic mapping system, the spectral radius ρJ of the eigenvalues of the Jacobi matrix J is larger than 1. The purpose of introducing the stability matrix for STM is to make the spectral radius ρJ* of the transformation system less than 1 under the condition of small enough parameter q without changing the location of unstable periodic fixed point. Pingel et al. [12] have determined the minimal set of stability matrix for a two- and three-dimensional system separately by using geometrical and topological arguments. Herein, we try to derive the stability matrix in light of the information of the eigenvalues at periodic fixed point.

Let J be an arbitrary n×n invertible matrix which has n distinct eigenvalues λ1, λ2,…, λn and the corresponding eigenvectors are p1, p2, …, pn. Then, there exists an invertible matrix P=p1,p2,,pn satisfying P-1JP=Λ:=diagλ1,λ2,,λn.

Thus, perform the multiplicative operation of Jocobian matrix J* Eq. (8) of stability transformation system:

(10)
P - 1 J * P = P - 1 [ I + q C ( J - I ) ] P = P - 1 I P + q P - 1 C ( J - I ) P
            = P - 1 I P + q P - 1 C P P - 1 J - I P = I + q P - 1 C P Λ - I = I + q C D Λ - I ,

where CD=P-1CP, and the stability matrix C can be obtained from the following three circumstances:

1) If all the eigenvalues of (JI) have negative real parts, i.e. Reλi1< 0 (i= 1, 2,…, n), let CD=I, and then C=PCDP-1=PIP-1=I.

2) If all the eigenvalues of (JI) have positive real part, i.e. Reλi1> 0 (i= 1, 2,…, n), let CD=I, then C=PCDP-1=PIP-1=I.

3) If (JI) has the eigenvalues with both negative and positive real part, let CD=diagδ1,δ2,,δn be a diagonal matrix, in which:

(11)
δ i = 1 ,                 R e ( λ i - 1 ) < 0 , - 1 ,           R e ( λ i - 1 ) > 0 ,

then the stability matrix C=PCDP˗1. However, in some situation the stability matrix C is unnecessarily taken as an involutory matrix depending on the Jacobian matrix of original dynamical system, which sheds a new light on the development of STM.

From the above processing, the matrix J* has become a diagonal matrix, and the next step is to choose a suitable parameter q, such that all the eigenvalues of the matrix J* have absolute value less than one i.e. λi*=1+qδi(λi-1)<1 (i= 1, 2,…, n). According to the stability theory, Eq. (7) will converge to some periodic fixed point.

5.2. Determination of critical parameter qcrit and optimal parameter qopt

The prominent advantages of STM as compared with other chaos control algorithms are the global convergence property, the low requirement for initial point and fast speed convergence to linear domain of periodic fixed point. However, the rate of convergence will be slow when iterative values are located in linear domain of periodic fixed point. To overcome this limitation, in terms of whether the information of periodic fixed points are known or not in advance, this paper will improve the performance of STM from two aspects as followings

For the information of periodic fixed points to be stabilized is known in advance, increase relax parameter q and the convergence speed will be faster. However, STM will not converge sometimes if q is too large. Hence, we should determine an optimal qoptthat makes STM converge fastest to periodic fixed point.

For the unknown information of periodic fixed points, we can combine STM with NT to improve the computational efficiency of STM.

As aforementioned, the relax parameter q of STM is closely related to the stability property of the periodic fixed points. As for a periodic fixed point, there exists a critical parameter qcrit which make the controlled system be at the boundary of stability and instability, namely the module of the maximum eigenvalue of the Jacobi matrix J* be equal to 1. It is obtained from Eq. (9):

(12)
λ i * = 1 + q δ i ( λ i - 1 ) 1 ,           i = 1,2 , . , n .

Let λi=ai+ibi and substitute it into Eq. (12) and the following module is achieved:

(13)
λ i * = 1 + q δ i ( a i - 1 ) + i δ i q b i 1 [ 1 + q δ i ( a i - 1 ) ] 2 + q 2 b i 2 1 ,           i = 1,2 , , n ,

in which 0 <q< 1, δi is presented in Eq. (11). From Eq. (13), the formula of parameter q is further attained:

q i n f i = 1,2 , n - 2 δ i a i - 1 ( a i - 1 ) 2 + b i 2 = i n f i = 1,2 , n 2 a i - 1 ( a i - 1 ) 2 + b i 2 .

Accordingly, the critical value qcrit is written as:

(14)
q c r i t = i n f i = 1,2 , n 2 a i - 1 ( a i - 1 ) 2 + b i 2 .

Again, take periodic-2 fixed points of controlled mapping system as an example. It is seen from Fig. 5(b) and Fig. 6(a) that when C1=I, q= 0.6, the module of the maximum eigenvalue of Jacobi matrix J* of transformed system Eq. (7) with system parameter b2= 1.144 is approximate equal to 1. That is to say, q= 0.6 is the critical value of STM at the system parameter b2= 1.144. On the other hand, when the mapping system parameter b2= 1.144, the analytical periodic-2 fixed point for mapping system can be acquired form implicit expression Eq. (5) and then substituted into Eq. (9), so the corresponding eigenvalues of Jacobi matrix J* are obtained:

λ 1 * = 0.5984 + 1.0846 i ,           λ 2 * = 0.5984 - 1.0846 i ,           λ 3 * = - 1.1944 .

In terms of Eq. (14) the critical value qcrit is calculated:

q c r i t = i n f i = 1,2 , 3 2 a i - 1 ( a i - 1 ) 2 + b i 2 = 0.6005 .

Which is in agreement with predefined value qcrit= 0.6.

Clearly, the convergence speed of STM becomes very slow when q is close to the critical parameter qcrit, and the transformed system is at the boundary of stability and instability. Fig. 8 illustrates the required iterative number when STM (in which b2= 1.144, C=I and X0(–3,1,0)) satisfies convergence condition with relax parameter q varying in the interval (0, qcrit). It is demonstrated from Fig. 8 that the iterative number versus q presents approximately an U-type curve, and when the relaxation parameter q is close to zero or critical value qcrit respectively, the required iterative number tends to infinity, i.e. theoretically ρJ*= 1. Therefore, there exists an optimal value qopt for the parameter q in the interval 0 <q<qcrit which makes the required iterations be the least. It is seen that the U-type curve changes slightly nearby q=qcrit/2 = 0.3, thus here the optimal parameter qopt is set as:

(15)
q o p t q c r i t 2 = i n f i = 1,2 , n a i - 1 ( a i - 1 ) 2 + b i 2 .

Fig. 8. Iterative numbers for STM when q varying in the interval (0, qcrit)

 Iterative numbers for STM when q varying in the interval (0, qcrit)

Nevertheless, in most instances the information of periodic fixed points are unknown in advance, and the exact optimal parameter qopt could not be determined. For this situation, this paper combines STM with NT to stabilize periodic fixed points. Firstly, a small enough parameter q is taken and a set of initial points and involutory matrix Ci are chosen. Next, STM is conducted until the absolute error of sequential iterative points satisfies a given value Xk+1-Xkε (ε 0.01~0.1). In this way, the iterative value has been in linear domain of periodic fixed point. Given the fact that Newton method has the second order convergence speed at the linear domain of initial condition, the iterative value Xk of STM can be taken as the initial point of Newton method and continue to conduct the iteration until satisfying the convergent condition. It is pointed out that in order to enlarge the linear domain, a relaxation factor can be introduced for NT, i.e. the relaxed Newton method can be taken as an alternative [14]. In practice, the NT is written as:

(16)
X k + p = X k - D G - 1 ( X k ) G ( X k ) = X k - ( J - - I ) - 1 F p X k - X k ,

where, Xk is the kth iterative point of STM which satisfies the given convergent condition, and GXk=FpXk˗Xk, J-=Fp(Xk)/Xk, DG(Xk)=J--I. The iterative value will converge to the desired periodic fixed point in a high precision with a small number of iterations. Fig. 9(a) shows the required iterative number versus system parameter b2 by the combination of STM and NT. Compared with Fig. 7(a), the efficiency is improved considerably. Meanwhile, the corresponding absolute and relative errors between iterative solution and fixed point are displayed in Fig. 9(b) and Fig. 9(c). When the relax parameter q= 0.01, the maximum absolute error Xk+1-X*= 0.1298×10-5 and the maximum relative error Xk+1-X*X*×100 %= 0.56 % which are remarkably decreased compared with those in Fig. 7(b) and Fig. 7(c). However, the most iterative number, the maximum absolute and relative errors still happen at the static bifurcation point of original mapping system.

Fig. 9. a) Iterative numbers; b) absolute error between iterative solution and fixed point; c) relative error between iterative solution and fixed point when converging to fixed point using STM-NT

 a) Iterative numbers; b) absolute error between iterative solution and fixed point; c) relative error between iterative solution and fixed point when converging to fixed point using STM-NT

a)

 a) Iterative numbers; b) absolute error between iterative solution and fixed point; c) relative error between iterative solution and fixed point when converging to fixed point using STM-NT

b)

 a) Iterative numbers; b) absolute error between iterative solution and fixed point; c) relative error between iterative solution and fixed point when converging to fixed point using STM-NT

c)

6. Conclusions

For deeply understanding the numerical performance of STM for chaos control of nonlinear dynamical system and convergence control of engineering iterative systems, several associated issues of STM is thoroughly studied.

Firstly, a 3D chaotic mapping system with various forms of bifurcations is taken as an example in this paper and the effects of three factors of STM such as initial condition, stability matrix C and relaxation parameter q on the capability of control algorithm are investigated in detail. It is found that the required iterative number N will decrease with relax parameter q, but less unstable fixed points can be detectable by STM. Especially, the stabilized periodic fixed point by STM will jump from one branch to another of analytical fixed point even if selecting the same relax parameter, initial condition and stability matrix (as red dashed line shown in Fig. 5(b)). The other detailed conclusions can be seen in Section 3.2.

Next, the eigenvalue curves of Jacobian matrix for different stability transformation systems are investigated. It is revealed that the maximum iterative number for STM occurs around the static bifurcation point, and thus the convergence speed of STM is very slow and hard to converge. Meanwhile, the absolute and relative errors between iterative solution and periodic fixed point with different relax parameter q are demonstrated respectively when STM satisfies convergence condition (≤ 10-5). It is observed that the maximum absolute and relative errors happen at the bifurcation point.

At last, according to the information of unstable periodic fixed point, the stability matrix of STM is derived for three different cases. However, in some situation the matrix C is unnecessarily taken as involutory matrix depending on the Jacobian matrix of original dynamical system, which sheds a new light on the development of STM. Meanwhile, the critical parameter qcrit on the basis of stability matrix C in STM are determined. It is found that the iterative number versus q presents approximately a U-type curve, and the smaller or the closer to the critical parameter qcrit the relax parameter q is, the slower the rate of convergence of STM is. Moreover, the optimal parameter qopt in the interval (0, qcrit) that makes the required iterations by STM be the least is obtained. When a priori information of the periodic fixed point is unknown in advance, STM is combined with NT to overcome the disadvantage of slow convergence around the static bifurcation point.

Acknowledgements

The supports of the National Natural Science Foundation of China (Grant No. 51275280) and Shandong Province Natural Science Foundation of China (Grant No. ZR2015EL048) are much appreciated.

References

  1. Ott E. Chaos in Dynamical Systems. 2nd Edition, Cambridge University Press, Cambridge, 2002. [Publisher]
  2. Miller J. R., Yorke J. A. Finding all periodic orbits of maps using Newton methods: sizes of basins. Physica D: Nonlinear Phenomena, Vol. 135, Issues 3-4, 2000, p. 195-211. [Publisher]
  3. Pérez R., Lopes V. L. R. Recent applications and numerical implementation of Quasi- Newton methods for solving nonlinear systems of equations. Numerical Algorithms, Vol. 35, Issue 2, 2004, p. 261-285. [Publisher]
  4. Bu S. L., Wang S. Q., Ye H. Q. Stabilizing unstable discrete systems. Physical Review E, Vol. 64, Issue 4, 2001, p. 046209. [Publisher]
  5. Pyragas K. Continuous control of chaos by self-controlling feedback. Physics Letters A, Vol. 170, Issue 6, 1992, p. 421-428. [Publisher]
  6. Pyragas K. Control of chaos via extended delay feedback. Physics Letters A, Vol. 206, Issues 5-6, 1995, p. 323-330. [Publisher]
  7. Pyragas V., Pyragas K. Adaptive modification of the delayed feedback control algorithm with a continuously varying time delay. Physics Letters A, Vol. 375, Issue 44, 2011, p. 3866-3871. [Publisher]
  8. Petalas Y. G., Parsopoulos K. E., Vrahatis M. N. Stochastic optimization for detecting periodic orbits of nonlinear mappings nonlinear phenomena in complex systems, Vol. 11, Issue 2, 2008, p. 285-291. [Search CrossRef]
  9. Parsopoulos K. E., Vrahatis M. N. On the computation of all global minimizers through particle swarm optimization. IEEE Transactions on Evolutionary Computation, Vol. 8, Issue 3, 2004, p. 211-224. [Publisher]
  10. Schmelcher P., Diakonos F. K. Detecting unstable periodic orbits of chaotic dynamical systems. Physical Review Letters, Vol. 78, 1997, p. 4733-4736. [Publisher]
  11. Pingel D., Schmelcher P., Diakonos F. K. Detecting unstable periodic orbits in chaotic continuous-time dynamical systems. Physical Review E, Vol. 64, Issue 2, 2001, p. 1-10. [Publisher]
  12. Pingel D., Schmelcher P., Diakonos F. K. Stability transformation: a tool to solve nonlinear problems. Physics Reports, Vol. 400, Issue 2, 2004, p. 67-148. [Publisher]
  13. Zhou J. L., Yang D. X. Chaos control of a new 3D autonomous system by stability transformation method. Nonlinear Dynamics, Vol. 73, Issues 1-2, 2013, p. 1-13. [Publisher]
  14. Yang D. X., Zhou J. L. Connections among several chaos feedback control approaches and chaotic vibration control of mechanical systems. Communications in Nonlinear Science and Numerical Simulation, Vol. 19, Issue 11, 2014, p. 3954-3968. [Publisher]
  15. Yang D. X., Li G., Cheng G. D. Convergence analysis of first order reliability method using chaos theory. Computers and Structures, Vol. 84, Issues 8-9, 2006, p. 563-571. [Publisher]
  16. Yang D. X. Chaos control for numerical instability of first order reliability method. Communications in Nonlinear Science and Numerical Simulation, Vol. 15, Issue 10, 2010, p. 3131-3141. [Publisher]
  17. Keshtegar B. Stability iterative method for structural reliability analysis using a chaotic conjugate map. Nonlinear Dynamics, Vol. 84, Issue 4, 2016, p. 2161-2174. [Publisher]
  18. Meng Z., Li G., Wang B. P., Hao P. A hybrid chaos control approach of the performance measure functions for reliability-based design optimization. Computers and Structures, Vol. 146, 2015, p. 32-43. [Publisher]
  19. Wang X. M., Li Y. F. Study of controlling discrete high-dimensional chaotic systems by prediction feedback. Journal of Jilin University (Information Science Edition), Vol. 29, Issue 1, 2011, p. 36-40. [Search CrossRef]

Cited By

Shock and Vibration
Fa Che, Chao Yin, Jilei Zhou, Zhinan Hu, Xingkui Zhao, Lu Sheng, Qian Zhang, Xiaoyu Cheng, Gangting Tang, Han Zhang, Dong Liu, Said Elias
2021