Trajectory Tracking between Josephson Junction and Classical Chaotic System via Iterative Learning Control

Abstract: This article addresses trajectory tracking between two non-identical systems with chaotic properties. To study trajectory tracking, we used the Rossler chaotic and resistive-capacitive-inductance shunted Josephson junction (RCLs-JJ) model in a similar phase space. In order to achieve goal tracking, two stages were required to approximate target tracking. The first stage utilizes the active control technique to transfer the output signal from the RCLs-JJ system into a quasi-Rossler system. Next, the RCLs-JJ system employs the proposed iterative learning control scheme in which the control signals are from the drive system to trace the trajectory of the Rossler system. The numerical results demonstrate the validity of the proposed method and the tracking system is asymptotically stable.


Introduction
Chaotic phenomena in semi-conductor devices found in the radio-frequency-base (rf-base) resistive-capacitive shunted Josephson Junction (RCs-JJ) were described by Kautz and Monaco [1] in the numerical study of three system parameters.Many studies have exhibited the chaotic behavior in superconducting resistive-capacitive-inductance Josephson Junction (RCLs-JJ) systems [2][3][4].The properties of RCLs-JJ including the homo-clinic, hetero-clinic, and super-harmonic bifurcations have been investigated to be excited by these varying parameters [5].The damped pendulum equation describes the junction behavior and demonstrates the chaotic strange attractor in the phase space [6].Synchronization is a significant topic, which was based on studying the tracking trajectory in chaotic systems.A non-linear controller by utilizing the back-stepping technique has been investigated to control bifurcation in the RCLs-JJ system [7].The chaos synchronization between two identical RCLs-JJ systems has been examined in which a number of different techniques to design the controller such as using active control [8], utilizing a common chaos to drive RCLs-JJ system approaching synchronization [9], applying the back-stepping [10][11][12][13], and using the time-delay feedback control [14], respectively.In other studies, the controller design or controller rule is directly determined by the Lyapunov function [11,12] and the RCLs-J junctions array synchronization [12].
Recently, the reconstruction of chaotic systems is similar to chaotic synchronization where the chaotic system was reconstructed from the unknown input [15][16][17].The chaotic systems synchronization utilized unknown inputs observed in the Takagi-Sugeno (T-S) fuzzy model [15].The Takagi-Sugeno (T-S) fuzzy model with unknown inputs applied the Lyapunov function and the Linear Matrix Inequalities (LMI) constraint approaching a zero synchronization error in secure communication [16].The evolutionary algorithms based on the measured data without any internal partial knowledge to reconstruct chaotic systems and demonstrates the numerical results in Reference [17].
In many studies, synchronization systems in the state space model for superconductor devices were described in identical RCLs-Josephson junction systems.Classical non-identical system synchronization is not meant to handle a superconducting system.Accordingly, many trajectory tracking studies have rarely been directly based on the RCLs-JJ system and classical chaotic systems together.This article focuses on the trajectory tracking between two non-identical systems in which one of the systems is the classical chaos system of the Rossler and the other is the RCLs-JJ model in state space representing a mesoscopic system.In this research, two stages tracking the trajectory method uses the RCLs-JJ system in the state space model to trace the trajectory of the Rossler system.The first tracking stage utilizes the active control technique [18] to transfer the RCLs-JJ system into the quasi-Rossler system.In the next stage, the proposed iterative learning control law approaches the trajectory of the Rossler system by correcting the repetition tolerance from the preceded output information between the quasi-Rossler system and the Rossler system [12].The RCLs-JJ system becomes the quasi-Rossler system and employs the iterative learning control (ILC) procedure to control signals from the Rossler system to trace the trajectory of the drive system.However, most ILC research update rules are designed based on the linear ILC law.The varying reference input [19] is very different from the chaotic signals, which are unexpected.Nevertheless, chaotic systems are unpredictable and can be synchronized [20,21] by designing a controller.The challenges involve extracting the message embedded within the chaotic signal from the normal signal of the regular system such as a digital encoding system.A chaotic signal cannot implement the code by the obvious method.The chaotic synchronization should lead to a secure communication system by using a synchronized chaotic electronic circuit [22].
In this study, the chaotic tracking trajectory that utilized the ILC method based on the Josephson junction chaos can also employ other methods such as a backstepping controller [7], an active control [8], a common chaos driving by Rossler [9], and LMI [16] and we examine the results.This article provides two main contributions.The first is a real-time feedforward procedure that uses iterative learning to modify the trajectory error between systems for tracking two non-identical systems.A few studies have been published on two different chaotic systems directly tracking or synchronizing.However, reports are especially rare for Josephson Junction systems and classical chaotic systems.Another contribution is the utility of the synchronization or minimum tracking error leads to the secure communication system by applying the Josephson Junction systems into quantum chaotic systems to deliver the message in the communication system.
The organization of this paper is explained.The next section describes the Rossler Chaotic and RCL-Shunted Josephson Junction system.Section 3 outlines the simulation results and provides an example to demonstrate the proposed learning control law in the RCLs-JJ system to trace the path of the Rossler system.Lastly, this paper highlights the future applications and conclusions in Section 4.

System Description and Transformation
The Rossler chaotic system has an initial condition X 0 , which is the drive system in a general form. .
is the state vector.The RCLs-JJ can be expressed by Equation ( 2) with initial conditions Y 0 = [0 0 0] T . .
where g(y 2 ) is given by g(y 2 ) = 0.366 as |y 2 | > 2.9 0.061 as |y 2 | ≤ 2.9, k in Equation ( 2) is the iterative number that uses the iterative learning control law (ILC) as shown in the equation below.
The ILC rule is a sequence of the control input signal for the response system as Equations ( 1) and ( 2) are nearly non-identical nonlinear systems, which is seen by their trajectory in Figure 1.
where ( ) is given by ( ) = 0.366 as | | > 2.9 0.061 as | | ≤ 2.9, k in Equation ( 2) is the iterative number that uses the iterative learning control law (ILC) as shown in the equation below.
The ILC rule is a sequence of the control input signal for the response system as ( ) , ,⋯ .Equations ( 1) and ( 2) are nearly non-identical nonlinear systems, which is seen by their trajectory in Figure 1. Figure 1 displays the time response of the state of two distinct systems known as the Rossler and the RCLs-JJ systems with different initial states.The trajectory error between the systems in each state should be enormous in Figure 1.The nonlinear system in Equation (2) should be transferred to the quasi-Rossler system to track the trajectory of the system in Equation (1).Therefore, the active control technique [23,24] is utilized in Equation (2).
According to the active control technique, the description of the system in Equation ( 2) and the dynamic transformation between the drive in Equation ( 1) and the response system in Equation ( 2) is shown in the equation below.2) should be transferred to the quasi-Rossler system to track the trajectory of the system in Equation (1).Therefore, the active control technique [23,24] is utilized in Equation ( 2).
According to the active control technique, the description of the system in Equation ( 2) and the dynamic transformation between the drive in Equation ( 1) and the response system in Equation ( 2) is shown in the equation below. .
where v a v b v c T is the active control function that eliminates the terms that have no z i for i = a, b, c.As a result, the active control function can be determined by the equation below.
where w a w b w c T is the error term in the active control procedure.Substituting Equation ( 8) into (7), Equation ( 7) was transformed into the formula below. .
The matrix form of Equation ( 8) can be rewritten as the formula below. .
The solution of w a w b w c T is shown below. .By substituting Equations ( 12) and ( 8) into the RCLs-Josephson Junction system in Equation ( 2) with an iterative learning control rule and changing the variable x to variable y, the system became the formula below. .
After the active control procedure, the RCLs-JJ system becomes a quasi-Rossler chaotic system so that the trajectory tracing between different systems traces the trajectory under identical systems.

Trajectory Tracking between Systems via the Iterative Learning Control
The RCLs-JJ system has now became the quasi-Rossler chaotical system.The ILC procedure and control signal from the drive system utilizes the response system to track the drive system.When an is found and the iteration number k is sufficient, the tracked error dynamic system should be equal to zero, that is .
The tracking trajectory has changed to two similar systems.In many studies, the synchronization between identical systems employ the control signal from the drive system [25][26][27].Accordingly, The RCLs-JJ system in Equation ( 2) uses the control signals from the drive system, x 1 and x 3 , which can be rewritten as the equation below. .
The control signals from the Rossler system of Equation ( 14) are similar to the quasi-Ross system in Equation ( 13) and the iterative learning control law U (k) , which is defined by error dynamics.The dynamic error system between the Rossler system in Equation ( 1) and the quasi-Rossler system in Equation ( 13) can be exhibited by the formula below.
where the matrix T in Equation (15) and realeig(B 2 ) is the real part of eigenvalue of B 2 .
When the z i=a,b,c = e j=1,2,3 , the term M x i , y j G(e) in Equation ( 14) can absorb the z a z b z c T to choose the appropriate matrix M. By induction, the expansion of u (k) in Equation ( 16) can be written as the equation below.

Lyapunov Stability of Systems
Equation ( 7) may be the dynamical error system in the active control procedure.Therefore, the Lyapunov function is defined by the equation below.
Theorem 1.The Lyapunov function in the active control procedure, which is used to transfer the RCLs-JJ system in Equation ( 2) to the quasi-Rossler system in Equation ( 13), can be defined by Equation (18).
Proof.Equation ( 18) should prove the first derivative is negative and the dynamic system is stable at the equilibrium (0, 0, 0).The first derivative of the Lyapunov function is shown below.
Substituting Equation ( 12) into Equation ( 9 Theorem 2. Based on the u (k) in Equation ( 16), the Lyapunov function is defined in the iterative control stage to trace the trajectory of the Rossler system below.
Proof.Let u (k−1) be defined as the formula below.
Applying −u (k−1) to Equation (15), we obtain the error dynamics as This implies that Equation ( 15) employs the iterative learning control law, which is asymptotically stable at equilibrium.The Lyapunov function is generally a function of the trajectory error or the synchronization error.The error norm can be extended to be a cost function having the lowest cost when the trajectory between systems is close to each other [28][29][30].The cost function in this work is defined by the formula below.

Results and Discussion
To verify the proposed iterative learning control law, we utilize an example to demonstrate the tracing error and trajectory between the Rossler dynamical system in Equation ( 1) with an initial state (x 10 , x 20 , x 30 ) = (0.2, 0.4, 0.1) and the RCLs-JJ system in Equation ( 2) with an initial state (y 10 , y 20 , y 30 ) = (0, 0, 0), respectively.

Deciding Iterative Control Learning Law
The Rossler system in Equation ( 1) is given by the formula below. .
The RCLs-JJ model in Equation ( 2) is given by the formula below. .
where the parameters in Equation ( 25) have been defined in Equations ( 3) and ( 4) in which the values of entries in matrix B are β L = 2.6, β c = 0.707, i N = 1.132, and ψ(y 1 ) = 1.132 − sin(y 1 ).The U (k) of the system in Equation ( 25) is the ILC rule in Equation ( 5).The u (k) defined in Equation ( 16) minimizes the tracking error between the Rossler system and the RCLs-JJ system.The matrix M x i , y j in Equation ( 15) and B 2 of Equation ( 16) can alternate respectively as the equation below.
where the x 1 is from Rossler system and matrix B 2 is the decomposition of matrix A in Equation (1).The time interval of the simulation range is from 0 to 300 s and the minimum time step is 0.01 s.This article utilized the program of the Euler method and figures in MATLAB to investigate trajectory tracking by using the ILC law in Equation ( 17).In the active control procedure, we use the Simulink in MATLAB to transfer the RCLs-JJ system to the Rossler system.

Exhibiting Simulation Results and Discussion
Figure 2 displays two distinct phase portraits of the two systems known as the Rossler and RCLs-JJ systems with different initial states.The chaotic behavior of the Rossler system for each phase portrait is displayed in Figure 2a-c.The phase portrait behaviors of the RCLs-JJ system is shown in Figure 2d-f.The chaotic behavior of the RCLs-JJ system is clear only in Figure 2e and the others are not chaos.The challenge is how to find the controller approach in the minimum trajectory error.
The first stage overcomes the non-identical trajectory between the two systems and employs the active control to change the RCLs-JJ system into a quasi-Rossler system from Equation (6) to Equation (13).The phases of two systems after utilizing the active control technique are exhibited in Figure 3a-c.
The new phase portraits of the RCLs-JJ system are almost completely different from the original phase portraits and are now closer to the Rossler system.Therefore, the new system is called a quasi-Rossler system.The time response of each component of the two systems indicated in Figure 3d-f in the paths of Rossler and RCLs-JJ systems, respectively, are distant from each other.
The phases of two systems after utilizing the active control technique are exhibited in Figure 3a-c.
The new phase portraits of the RCLs-JJ system are almost completely different from the original phase portraits and are now closer to the Rossler system.Therefore, the new system is called a quasi-Rossler system.The time response of each component of the two systems indicated in Figure 3d-f  Original phase portrait of the Rossoler system with [ x 10 x 20 x 30 ] T = [ 0.2 0.4 0.1 ] T in the left column, (a) the Rossoler system x 1 , x 2 , (b) the Rossoler system x 2 , x 3 , and (c) the Rossoler system x 1 , x 3 , and the original phase portrait of the RCL-shunted Josephson Junction system with an initial condition at the original in the right column, (d) the RCLs-JJ system y 1 , y 2 , (e) the RCLs-JJ system y 2 , y 3 , and (f) the RCLs-JJ system y 1 , y 3 .( Figure 4 demonstrates the tracking error between the Rossler and the quasi-Rossler systems, which transfers from the RCLs-JJ system.The vibration of the tracking error has many large amplitudes in the second (y 2 -x 2 ) and third (y 3 -x 3 ) components at the specific moment.
Figure 5 demonstrates the phase portrait of the x 1 (y 1 ) and x 2 (y 2 ) by utilizing the ILC update law to track the trajectory.The two trajectories almost overlap, which validates the tracking error phenomena in Figure 6.
The tracking errors oscillation in Figure 6 are between 0.1408 and −0.1959 for the first tracking error, between 0.1434 and −0.4217 for the second, and between 0.4784 and −0.8344 for the third tracking error, respectively.
Figure 5 demonstrates the phase portrait of the x1 (y1) and x2 (y2) by utilizing the ILC update law to track the trajectory.The two trajectories almost overlap, which validates the tracking error phenomena in Figure 6.
The tracking errors oscillation in Figure 6 are between 0.1408 and −0.1959 for the first tracking error, between 0.1434 and −0.4217 for the second, and between 0.4784 and −0.8344 for the third tracking error, respectively.Figure 7 exhibits the tracking error after using the ILC law.The third ingredient varies the most when comparing to subfigures in Figure 1.The lager vibration will always occur at a particular Figure 7 exhibits the tracking error after using the ILC law.The third ingredient varies the most when comparing to subfigures in Figure 1.The lager vibration will always occur at a particular moment t i because the tracking error between two systems became larger at same moments.Comparing Figures 4 and 6, the tracking errors are successfully suppressed between 0.2 and −0.2 for the first two components and the fault of the third component between 0.5 and −0.8 by applying the proposed ILC update law.The error rate is asymptotically stable.
Figure 8 shows the cost function after utilizing the ILC procedure defined by the tracking error norm, which is bounded and less than one.The cost function is created by Figure 4 after applying the divergent active control procedure.The trajectories of the two systems is near each other.
Comparing Figures 4 and 6, the tracking errors are successfully suppressed between 0.2 and −0.2 for the first two components and the fault of the third component between 0.5 and −0.8 by applying the proposed ILC update law.The error rate is asymptotically stable.
Figure 8 shows the cost function after utilizing the ILC procedure defined by the tracking error norm, which is bounded and less than one.The cost function is created by Figure 4 after applying the divergent active control procedure.The trajectories of the two systems is near each other.moment ti because the tracking error between two systems became larger at same moments.Comparing Figures 4 and 6, the tracking errors are successfully suppressed between 0.2 and −0.2 for the first two components and the fault of the third component between 0.5 and −0.8 by applying the proposed ILC update law.The error rate is asymptotically stable.
Figure 8 shows the cost function after utilizing the ILC procedure defined by the tracking error norm, which is bounded and less than one.The cost function is created by Figure 4 after applying the divergent active control procedure.The trajectories of the two systems is near each other.

Conclusions
This article proposed a learning control law to trace the trajectory between two non-identical nonlinear systems and successfully used a two-stage approach of combining an active control technique and an iterative learning control update law to significantly inhibit and improve the tracking errors in the numerical results.The simulation example helped infer the trajectory tracking process and assert the proposed ILC rule.In this research, the ILC tracking method has been applied to the Josephson junction chaos.Many studies providing other methods such as the backstepping controller [7], the active control [8], common chaos driving by Rossler [9], and LMI [16] should show similar demonstrations.
Future studies should provide a feedforward procedure in a real-time iterative learning process to modify the trajectory error between non-identical systems and apply synchronization or the minimum tracking error procedure into a secure system in quantum communication.The real-time iterative learning procedure would use signal tracking in the deep learning procedure.

Figure 1 .
Figure 1.The time response of the Rossler system and the resistive-capacitive-inductance shunted the Josephson Junction (RCLs-JJ) system.

Figure 1 .
Figure 1.The time response of the Rossler system and the resistive-capacitive-inductance shunted the Josephson Junction (RCLs-JJ) system.

Figure 1
Figure1displays the time response of the state of two distinct systems known as the Rossler and the RCLs-JJ systems with different initial states.The trajectory error between the systems in each state should be enormous in Figure1.The nonlinear system in Equation (2) should be transferred to the quasi-Rossler system to track the trajectory of the system in Equation(1).Therefore, the active control technique[23,24] is utilized in Equation (2).

=
−z a −z b −z c T a −z b −z c T , and taking s 1 = s 2 = s 3 = 1, the following formula is found..

Figure 2 .
Figure 2. Original phase portrait of the Rossoler system with = 0.2 0.4 0.1 in the left column, (a) the Rossoler system x1, x2, (b) the Rossoler system x2, x3, and (c) the Rossoler system x1, x3, and the original phase portrait of the RCL-shunted Josephson Junction system with an initial condition at the original in the right column, (d) the RCLs-JJ system y1, y2, (e) the RCLs-JJ system y2, y3, and (f) the RCLs-JJ system y1, y3.

Figure 3 .
Figure 3.After applying the active control procedure to (a-c), the phase portrait of the Rossler (x i -blue) and RCLs-JJ system (y j -red) and (c-f) time response of the system is stated in the Rossler (x i -blue) system and the RCLs-JJ system (y j -black).(a) Phase portrait x 1 (y 1 ), x 2 (y 2 ); (b) Phase portrait x 2 (y 2 ), x 3 (y 3 ); (c) Phase portrait x 1 (y 1 ), x 3 (y 3 ); (d)Time response of state x 1 -blue, y 1 -black; (e) Time response of state x 2 -blue, y 2 -black; (f) Time response of state x 3 -blue, y 3 -black.

Figure 4 .Figure 4 .
Figure 4.The tracking error between the Rossler and RCLs-JJ systems through an active control procedure.

Figure 5 .
Figure5.The phase portrait of x 1 (y 1 ) and x 2 ( y 2 ) using the iterative learning control update law.

Figure 6 .
Figure 6.The tracking error between the Rossler and the RCLs-JJ systems utilizing the ILC update rule.

Figure 6 .
Figure 6.The tracking error between the Rossler and the RCLs-JJ systems utilizing the ILC update rule.

Figure 7 .
Figure 7.The time response of x3 for the Rossler system and y3 for RCLs-JJ using the ILC (iterative learning control) update rule.

Figure 8 .Figure 7 .
Figure 8.The cost function from the error norm after applying the ILC (iterative learning control) update law.

Figure 7 .
Figure 7.The time response of x3 for the Rossler system and y3 for RCLs-JJ using the ILC (iterative learning control) update rule.

Figure 8 .Figure 8 .
Figure 8.The cost function from the error norm after applying the ILC (iterative learning control) update law.