Adaptive Generalized Synchronization between Circuit and Computer Implementations of the Rössler System

: The synchronization between chaotic systems implemented in similar ways—e.g., computer models or circuits—is a well-investigated topic. Nevertheless, in many practical applications, such as communication, identification, machine sensing, etc., synchronization between chaotic systems of different implementation types—e.g., between an analog circuit and computer model—might produce fruitful results. In this research, we study the synchronization between a circuit modeling the Rössler chaotic system and a computer model using the same system. The theoretical possibility of this kind of synchronization is proved, and experimental evidence of this phenomenon is given with special attention paid to the numerical methods for computer model simulation. We show that synchronization between a circuit with uncertain parameters and a computer model is possible, and the parameters obtained from the synchronized computer model are in high correspondence with the circuit element specification. The obtained results establish the possibility of using adaptive generalized synchronization for the parameter identification of real systems. It was also found that Heun’s method yielded the most accurate results in synchronization among second-order numerical integration methods. The best among the first-order methods appears to be the Euler–Cromer method, which can be of interest in embedded applications. method ( Heun ), implicit midpoint method ( IMP ), explicit midpoint method ( EMP ) and trapezoidal rule ( Trapezoidal ). The best results are achieved using second-order Heun and trapezoidal methods.


Introduction
The synchronization of dynamical systems is a widely used principle in communication, measurement and control. The possibility of synchronizing chaotic systems gained attention after work by Pecora and Carroll [1], while the phenomenon of chaotic synchronization was independently discovered earlier by several research groups [2,3]. That primarily discovered type of synchronization was later called complete synchronization (CS) because it implies the fully correlated behavior of coupled systems. A decade after the first works on CS were published, it was found that nonidentical chaotic systems can also synchronize [4]. This type of synchronization was called generalized synchronization (GS), and a simple approach for detecting generalized synchronization was proposed based on an auxiliary system that maps the case of GS to the case of CS [5]. GS may occur in both slightly dissimilar systems as well as in totally different systems [6,7]. The work presented in [8] introduces the concept of the consistency of a nonlinear system's output signal, which is closely related to synchronization and can be treated as the property of the signal produced by a system in GS with the input signal.
While CS can occur in identical systems, many applications imply that the master (drive) and slave (response) systems may have different parameter values. To synchronize these systems, it was proposed to equip a slave system with a special controller aimed at fitting the mismatched parameters. Typically, the parameter adaptation law is derived from the Lyapunov function derivative [9,10]. This approach can be naturally extended to GS [11,12].
Interest in adaptive synchronization, especially its generalized case, is due to the possibility of using this principle in multiple applications. The most prominent example is parameter identification, because adjusted parameters of the slave system can then often be directly measured and related to uncertain parameters of the master system [10,[13][14][15]. Particular examples in which this approach has been used are the design of a secure communication system [9], parameter identification of a gene regulatory network [16] and investigation of circadian rhythms in mammals [17]. It also may be used in chaotic sensors, typically consisting of a master system which serves as a measurement unit and a slave system which serves as a reference [18,19].
While researchers usually report synchronization between systems implemented in a similar way-e.g., between two computer models or between two analog circuits [20]synchronization between similar systems implemented in different ways-e.g., the synchronization between a circuit and a computer model is also an important but poorly elaborated problem. Meanwhile, experimental evidence of synchronization between a real system and a model has been reported in several studies; e.g., [21], which shows that a different physical nature is not an obstacle for synchronization. Computer implementation is often more reliable, stable and easy to investigate, while a circuit with uncertain parameters may play the role of a sensor, a plant in a control task or provide some real data for identification. In our previous works, we found that the CS between digital and analog implementations of the same chaotic system meets notable restrictions due to truncation and round-off errors in a computer model, noise and imperfect operation of the circuit and the existence of signal-corrupting barriers such as analog-to-digital converters (ADCs) and digital-to-analog converters (DAC) in a coupling channel [22,23]. These results forced us to conduct the current research, in which we focus on the adaptive generalized synchronization between the analog and numeric implementation of the same system.
In this study, we consider a variant of adaptive synchronization between Rössler oscillators. We experimentally confirm the following hypotheses: 1. Synchronization between analog and numerical models of the same system is possible despite some parameters of the analog system being uncertain. 2. It is possible to perform analog system parameter identification using the controller derived from the Lyapunov function. 3. Numerical methods used for slave system implementation substantially affect the synchronization accuracy.
This paper is organized as follows. Section 2 describes the main theoretical highlights. Section 3 describes our experiments and reports their results. Finally, we present conclusions.
Data obtained from a designed circuit and some program codes in MATLAB can be found in the authors' repository [24].

Adaptive and Generalized Synchronization
This section presents the theoretical background and shows how the different natures of synchronized systems affect their behavior after coupling.

Terms and Definitions
Let us consider a master systemẋ and a slave systemẏ (t) = G(x(t), y(t)), where x ∈ R n , y ∈ R m and t ∈ R 1 . For n = m the systems (1) and (2) are said to be completely synchronized if lim t→∞ x(t) − y(t) = 0.
Equation (3) shows that the dynamics of the master and slave systems become identical over time irrespective of the different initial conditions they had. One way to completely synchronize two similar systems is to introduce a control law C F and design the slave system as follows:ẏ (t) = F(y(t)) + C F (x(t), y(t)).
If there is a vector of uncertain parameters g in the master system, a case of adaptive synchronization is considered. A parameter update law for the slave system C g can be designed, leading to an overall master-slave configuration of The control law C g is designed to ensure that the vector of the slave parameters g s matches the vector of the master parameters g m over time: The technique to obtain C F and C g is as follows. Denote the error between system phase variables and parameters as e(t) = (y(t) − x(t); g s (t) − g m ) . Then, consider the Lyapunov candidate function as where P is a positive definite matrix; e.g., an identity matrix. The Lyapunov candidate function derivativeV(t) is then algebraically expanded, and C F and C g are found analytically from equations, providingV(t) < 0 for all t > 0. If it is impossible to find such functions, a weaker condition is considered,V(t) ≤ 0. In this case, it should be ensured that If this condition is satisfied, then, according to Barbalat's lemma, lim t→∞ e 2 = 0. Section 3 gives the particular examples of C F and C g for the Rössler system with uncertain parameters.
Two different systems (1) and (2) are said to be in generalized synchronization if a nonlinear function φ exists such that In the case of a generalized synchronization between systems with similar dimensions, the error in dynamics ∆(t) = x(t) − y(t) persists. The error trajectory ∆(t) lies within an attractor A e ⊂ R n with non-zero amplitude.
There is a simple approach to detect generalized synchronization between systems (1) and (2) using an auxiliary system:ż where G is similar in both slave and auxiliary systems. If GS occurs, then, due to the existence and uniqueness theorem, the error between y(t) and z(t) vanishes over time.
In fact, in this case, we deal with two identical slave systems driven by one master system, and we know nothing about function φ, which is convenient in many applications.

Synchronization between Systems Implemented in Different Ways
In any experimental study, we are restricted in terms of implementation accuracy since each method of implementation imposes several errors on the system. Therefore, any experimental study deals not with an ideal continuous system F but with its distorted counterpartF. Figure 1 shows the main error sources in analog-to-numeric synchronization. Each kind of error in chaotic system implementation leads to the non-predictable disruption of system behavior due to the sensitive dependence of system dynamics on parameter values and nonlinearity types.
Circuit implementation suffers from non-ideal components and analog noise. For example, typical values of resistor tight tolerance are within 0.5-5%, which substantially limits the accuracy of the parameter setting. The analog integrators usually implemented with operational amplifiers are not strictly linear. For instance, the high-quality operational amplifier OPA2134 has a total harmonic distortion at 1 kHz equal to 8 × 10 −5 % [25], which is relatively small but restricts the integration accuracy to values much lower than achievable by a digital computer. Temperature drift and noise of various natures also affect the circuit implementation accuracy.
A computer model can be implemented much more accurately but also suffers from errors. In floating-point arithmetic, round-off noise sets the theoretical lower bound of simulation result error at the level of the machine epsilon, equal to approximately 2.2 × 10 −16 for double precision [26]. Truncation error plays a notable role if the order of the integration method is low. An additional source of error in synchronization between analog and numeric systems is contained in a coupling channel. The time-domain resolution of this channel is restricted by at most 1 GHz, but usually is much lower. Besides, the ADC and DAC have a limited level of resolution-e.g., 12 or 24 bit-which makes the use of a full double precision scale impossible.

Continuous System
Non-ideal components Analog noise Truncation error Round-off noise Under all of the abovementioned conditions, it is reasonable to ask whether the synchronization between analog and numerical systems is possible. Our previous works [22,23] provide some experimental evidence of its possibility but also report on valuable restrictions in error elimination. Theorem 1 establishes the theoretical bounds on synchronization accuracy in this case. Definition 1. The L 1 -volume of an attractor A is the L 1 -norm of the main diagonal of the hypercube circumscribed around the attractor. In other words, if the trajectory is denoted as x(t), then L 1volume is The concept of the L 1 -volume might be especially useful for estimating some attractor properties, as a conventional volume of a strange attractor may be equal to zero. Preliminary 1. Let F : P ⊕ X → X be an autonomous continuous chaotic system, defined by the equationẋ in the Banach space X, and g m ∈ P is a parameter vector. Denote this system as a master system. Preliminary 2. Let the slave system G : P ⊕ Y → Y be defined by the equationṡ Consider the particular case of diffused coupling and automated parameter matching: where g s ∈ P, C F : P ⊕ X ⊕ Y → Y is the synchronization law and C g : X ⊕ Y → P is the parameter adaptation law, and Y, P are also Banach spaces.
Consider the case when the systems F and G are completely synchronized; i.e., conditions (3) and (6) are satisfied. For simplicity, assume X = Y.
where p ∈ P. Let | · | be a metric in a function space. Then, |F −F| > 0.
A distorted master systemF can differ in both parameters and nonlinearities from F. Preliminary 4. Consider a discretization operator over Y L h : G(y(t)) → Γ(ŷ(n)) and P L h : C g (g s (t)) →Ĉ g (ĝ s (n)) with a local error e y (h, n): where h is a discretization period, and lim h→0 ε h (n) = 0.
The trajectory y(n) of the discrete system Γ : P ⊕ Y → Y mapped from G with the operator L h is defined only in discrete points t n and in those points is close to the trajectory of G, if the discretization error is low. A similar case is true for g s (n), which is a trajectory of the systemĈ g mapped from the system C g .
The discrete systems Γ andĈ g preserve all parameters and nonlinearities of tge original systems.

Preliminary 5.
There is a coupling channel with a key between the continuous master systemF and the discrete slave system Γ, which transforms the systemF continuous trajectoryx(t) into the discrete series {x(n)} with a discretization period h. For clarity, {x(n)} =x(t n ).
From these preliminaries, a conjecture follows. Conjecture 1. If |F −F| < δ and e y (h, n) < η and e gs (h, n) < ξ, where δ, η, ξ > 0 are certain small positive numbers, generalized synchronization betweenF and Γ occurs, and the trajectory of the error e(n) = (ŷ(n) −x(n);ĝ s (n) − g m ) lies within a bounded set A ⊂ Y ⊕ P. Its L 1 -volume V L 1 (A) satisfies the limit property: We propose that this would be true in case the rate of exponential decay of error dynamics e(n) between G and F is faster than the growth of errors caused by the dissimilarity between the original synchronized systems F and G and their perturbed counterpartsF and Γ. A rigorous proof is beyond the scope of this article.
The following corollaries are of practical importance.

Corollary 1.
Parameter identification via adaptive synchronization between two different systems is possible with the error e g (n).
In practice, while the maximal values of e g (n) can be rather high, the mean or median value might be rather small, allowing its use for identification and measurement.

Corollary 2.
The synchronization error between master and slave systems exŷ(n) depends on the integration method and is lower for more precise methods.
Both corollaries are illustrated in Section 3.

Adaptive Complete Synchronization between Numerical Models
As the sample model, we consider the Rössler chaotic system that is described by the following system of ordinary differential equations (ODE): where x, y, z are state variables and a, b, c are nonlinearity parameters.
In [27], Park proposed a technique to identify three unknown parameters of a system (10) using adaptive control. According to such a method, the master (m) and slave (s) systems are defined as and where u 1 , u 2 and u 3 are control laws andâ,b andĉ are estimated parameters. Based on the Lyapunov stability theory, synchronization of the master and slave systems is achieved under the control laws given by with the update rule for parameters described as follows: where e = (e 1 , e 2 , e 3 ) is the error vector between (11) and (12): To verify the considered approach for parameter estimation, we simulated the master and slave systems using the Dormand-Prince eighth-order integration method with a fixed integration step of h = 0.01 within 500 seconds. For systems (11) and (12), the initial conditions were chosen as x m = −1, y m = 0, z m = 1 and x s = 2, y s = −4, z s = −2 while the parameters were a = 0.2, b = 0.6 and c = 5.25. The initial values of three parameters for the slave system were equalâ =b =ĉ = 1. Figure 2 presents the results of adaptive synchronization between two computer models simulated with the eighth-order integration method. Parameters are fitted exactly, and the error reaches the machine epsilon in double precision. The numerical simulation also shows that the lowest reachable error in a complete adaptive synchronization depends on the accuracy of the integration method. Figure 3 presents the error e(t) dynamics between a master system implemented with an eighthorder Dormand-Prince method and a slave system implemented with various methods of orders 1 and 2 obtained in a numerical experiment. We considered eight integration methods listed below with the corresponding notations: the Euler method (Euler), Euler-Cromer method (Euler-Cromer), semi-implicit CD-method (CD), Heun's method (Heun), implicit midpoint method (IMP), explicit midpoint method (EMP) and trapezoidal rule (Trapezoidal).
The best results are achieved using second-order Heun and trapezoidal methods.

Adaptive Generalized Synchronization between Numerical Models
An interesting case is a generalized synchronization between two discrete systems, one of which (master) is simulated using the Dormand-Prince eighth-order method and the other using a low-order integration method with a much higher truncation error.
As expected, the minimal achieved error between slave systems is of the order of machine epsilon in double precision, while the error in master and slave systems is similar to the error obtained in the previous experiment.
As expected, there was no significant difference between integration methods in the resulting convergence rate. The final achieved accuracy had a steady state of machine epsilon order and below. Among all tested first and second-order methods, Heun's and trapezoidal methods tended to result in slightly more precise synchronization.

Circuit Implementation
The implementation of the chaotic circuit based on the Rössler system is shown in Figure 6. This circuit was designed with the classical principles of analog computer design. The circuit consists of three integrators made using the operational amplifiers U 1 , U 2 and U 4 , each corresponding to the line of the differential Equation (10). We start circuit synthesis with selection of R 0 and C 0 , defining the simulation time constant τ = R 0 C 0 . In our case, R 0 = 10 kΩ and C 0 = 10 nF.
Each state variable or constant term in the original equation is connected via a resistor to an integrator. If a state variable is multiplied by some constant ξ, then the corresponding resistor value decreases by ξ times; i.e., R = R 0 /ξ. If ξ = 1, the resistor R 0 is used. This principle allows the easy construction of circuit parts corresponding to the first and second lines of (10).
When considering the last line, one should keep in mind that the amplitudes of state variables require a decrease of M times compared to the original Rössler oscillator state variables. This is due to the fact that tge z variable in the Rössler system with standard parameters can exceed not only 10 but even 20, which in analog implementation means 20 V, while an allowable voltage range is commonly from −10 V to +10 V. Since the similar scale constant M is used for all state variables, only the last line of (10) needs to be changed. Let us write it as follows: where x = Mx and y = Mŷ are the state variables with no scaling, and V b and V c are voltages with values that are mathematically equal to tge values of appropriate parameters in numerical system: Recall also that multiplier AD633 (Analog Devices, Norwood, MA, US), the element U 5 in Figure 6, acts as follows: i.e., the output voltage W is the sum of the input voltage Z and the product of the input voltages (X 1 − X 2 ) and (Y 1 − Y 2 ) divided by 10. This scaling is necessary so that the multiplication operation does not lead to the appearance of impossible voltages in the circuit. Taking this into account, we rewrite the third line of (10):ẑ where the term ψ is the direct output from the U 5 (AD633) multiplier. From (18), it follows that V b requires division by 10M before application to the Z input of the multiplier. Furthermore, the factor 10 at the beginning of the line defines the resistor R 13 = R 0 /10 between the multiplier U 5 and the integrator U 4 . The exact equation of the whole circuit in Figure 6 is Let (19) may be convenient if different scales of state variables are required or capacitors of different values are to be used, in which cases it becomes difficult to apply the empirical principles described above. Furthermore, (19) may be helpful for more realistic circuit simulation.
All operational amplifiers in the circuit in Figure 6 work as inverting amplifiers. Therefore, the variable −ŷ at the output of U 2 is obtained, and an additional op-amp U 3 is required to change theŷ sign. Second, the circuit does not oscillate immediately after the power is applied, and zero initial conditions should be set for the oscillations to start. For this purpose, resistors R 7 and R 14 are used, which limit the op-amp gain in the absence of oscillations, meaning that a zero output signal is established. After achievinĝ x =ŷ =ẑ = 0 V, the resistor R 7 branch may be opened to force the circuit to go on the attractor. Resistor R 14 does not affect the chaotic oscillations. A resistor parallel to C 1 may also be introduced but was not required in our experiments.
Ultra-low distortion, low-noise operational amplifiers OPA2134 (Texas Instruments, Austin, TX, USA) were used to ensure the accuracy of analog computations. Parameters b and c were defined as voltages using National Instruments ELVIS III station with a high-precision DAC.

Circuit Operation
To validate the correct operation of the circuit, its phase portrait and bifurcation diagram were obtained. For the phase portrait of the analog circuit ( Figure 6) with parameters a = 0.2, b = 0.2, c = 5.7, refer to Figure 7a, which represents the recognizable chaotic scroll. In the bifurcation diagram (Figure 7b) which demonstrates chaotic attractor points as a function of parameter c, Feigenbaum sequences of bifurcations characterizing the Rössler system are also recognizable. Data acquisition and processing was performed using the NI ELVIS III station and the software developed in the LabVIEW 2019 environment.
However, the specific values of the nonlinearity parameters at which different oscillation modes are observed in the analog circuit differ from those in the numerical simulation. To ensure consistency between the diagrams for parameter c, the value of parameter b in the analog circuit has been adjusted. A feature of the Rössler system that makes it difficult to implement in the analog computer is the large ratio between the values of the variable z amplitude and parameter b. Therefore, in our implementation, the value b/40 = 0.005 V was applied to the summing input of U 5 multiplier. This value is presumably of the same order of magnitude as the weak output bias present in the multiplier and also is strongly distorted by analog noise. Thus, when specifying the value b, this offset must be taken into account. In our case, the bias had a value of about 0.0025 V. Any adjustments of the value of parameter b do not lead to a complete correspondence between bifurcation diagrams, since there are other nonlinearities and errors in specifying the values of components in the analog circuit that affect its operation; in particular, the amplitudes of the attractors differ.
Thus, system (19), which is a linear transformation of system (10), does not fully characterize the dynamics of the circuit. The exact identification of non-linear components is challenging. Therefore, complete synchronization between the circuit and the mathematical model is not possible. The same result was achieved in one of our previous works [22], where the Rössler system was also considered and error attractors in case of Pecora-Carroll synchronization [1] were obtained. Note that in [22], a slightly modified Rössler oscillator with an additional parameter d was studied; additionally, another circuit was designed with the parameter c defined as the ratio between two resistors, not as an external voltage.

Adaptive Generalized Synchronization between Circuit and Numerical Model with Respect to Numerical Method
The experimental study of adaptive generalized synchronization between the circuit and the computer model was conducted using the following parameters of the slave system: x s,1 = 2, y s,1 = −4, z s,1 = −2, initial parametersâ 1 = 0.2,b 1 = 0.6,ĉ 1 = 5.25. Heun's integration method was used, and the step-size was h = 0.0464. The auxiliary slave system was similar. Figure 8 presents the time series of the error between the circuit and computer model, taken as the L 2 -norm of the error in dynamics, and the time series of estimated parameters. One can clearly see that the estimated parameters are oscillating and do not converge to steady-state values. This behavior shows that the perturbations caused by the introduced errors are significant and cannot be tracked by the controller well. However, the overall stability of the solution is not violated.  Figure 9 presents the synchronization error between master and slave systems for different integration methods, smoothed with the median filter. There is a clear dependence of the averaged error on the integration method. It is interesting to note that the firstorder symplectic Euler-Cromer method presents an equally accurate synchronization as the second-order Heun's method. This result is significant as it allows the use of an algebraically simpler simulation algorithm without loss of accuracy. This is of great importance in embedded applications. A rough explanation of this phenomenon can be that the Rössler system phase space is approximately isomorphic to the harmonic oscillator phase space in regions in which |z| < ε, and ε > 0 is a small real number. However, this needs to be investigated further.  Figure 10 shows the accuracy of the master-slave and slav-auxiliary systems' synchronization achieved in our experiment. It can be compared with Figure 4, which shows that GS is preserved in the cases of analog and numeric chaotic systems.  Figure 11 illustrates one again the proposition that the rate of convergence between slave systems weakly depends on the integration method, and this can be compared with Figure 5. In the case of heterogeneous chaotic systems, the slave-auxiliary synchronization rate is approximately similar for all investigated methods. It is of interest that the averaged rate of all methods except for Euler and Trapezoidal methods is slightly superlinear. Summarizing the results, we can conclude that, even though the master and slave systems are heterogeneous, GS is experimentally verified, as well as CS. It is shown that the synchronization error in dynamics actually depends on the numerical method. Among all investigated methods, Heun's and Euler-Cromer methods look preferable.

Tracking an Additional Parameter in Adaptive Generalized Synchronization
Modify the Rössler system (10) as follows: where µ is the novel parameter, which is not tracked by the controller (13). This parameter is easily tuned in the circuit by changing the only resistor R 1 . In our experiment, the following values of R 1 were used: 5 kOhms, 8 kOhms, 9 kOhms, 10 kOhms, 10.5 kOhms, 11 kOhms, 12.8 kOhms and 20 kOhms. Different attractors of the Rössler system corresponding to these values are given in Figure 12. Parameter µ values differing from 10 kOhms introduce additional error in the implemented analog Rössler system. Attractors A e ⊂ R 3 of the error in dynamics e(n) between the circuit implementation of the master system (20) and computer model of the slave system (12) are given in blue in Figure 13.
Then, we equip the slave system with an additional adjustable parameterμ, which now reads as Modify the adaptive controller in the following way: with the following update rule for parameters: = −y m e 2 ,b = −e 3 ,ĉ = z m e 3 ,μ = y m e 1 .
(23) Figure 13. Error attractors between the circuit and the computer model: blue for the slave system (12) and red for the slave system (21).
The actual error in synchronization between the circuit and digital implementation of the system (20) using Heun's method with a step-size of h = 0.0518 s is given in red in Figure 12. This is substantially smaller in comparison with the error of the system (21). To obtain a quantitative measure of error attractors, we used the attractor L 1 -volume V L 1 (A e ) = ∑ j=x,y,z (max n (e j (n)) − min n (e j (n))). Its estimations are given in Figure 14.
One can see that additional system perturbation with the parameter µ increases the synchronization error when this parameter is not tracked, and the new value of µ sufficiently differs from 1. At the same time, tuning the corresponding parameterμ allows a notable reduction of the synchronization error.
To estimate the error in the parameter µ identification, the median value over 9000 samples ofμ(n) was computed. The estimated R 1 was calculated aŝ where M is the median and R 0 is 10 kOhms. Obtained results are given in Figure 15. Figure 14. Dependence of the error attractor L 1 -volume on the parameter µ value for slave systems without and with a controller onμ, corresponding to slave systems (12) and (21), respectively.
Although there was no complete synchronization between analog and computer systems, the relative error of parameter identification was on average less than 1.5%. This shows the practical applicability of this approach as a parameter identification technique. Figure 15. (a) Resistance R 1 identification results; (b) Evaluation of relative error in R 1 identification: the plot shows the actual R 1 value versus the relative error of its estimation | ∆R 1 R 1 |, where ∆R 1 =R 1 − R 1 .

Conclusions
For two chaotic systems of arbitrary natures that are close to each other in a reasonable sense, it is possible to achieve generalized synchronization. According to Conjecture 1, the synchronization error tends to vanish if the system difference becomes negligible. Once the driven system is equipped with the parameter estimation controller, one can fit uncertain parameters of the driving system even though there is no straight correspondence between the systems.
To experimentally verify Conjecture 1, we developed a novel design of a circuit modeling a Rössler chaotic system. Generalized synchronization between the circuit and computer models was established, and the synchronization error was estimated.
Several first and second-order integration methods were tested. We experimentally confirmed that the methods of a higher-accuracy order can improve the synchronization accuracy in computer simulations. The best result was obtained by Heun's method. It is of interest that the symplectic Euler-Cromer method gave almost similar results to Heun's method, with first-order accuracy and increased computationally efficiency.
We show that the use of adaptive generalized synchronization allows the parameters of an analog system to be tracked. We implemented a controller monitoring the parameter vector and found that the median value of one uncertain parameter fitted by this controller corresponds to the analog value with high accuracy. This opens the possibility of applying adaptive generalized synchronization for identification, secure communication and sensing using complete adaptive synchronization.
In future works, we will apply the developed approach to study other chaotic systems and also consider practical examples of applying adaptive generalized synchronization according to our findings. A more complete investigation of how symplectic methods can be used for the simulation of non-conservative chaotic systems will also be a topic of further investigation.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.