- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 398293, 13 pages
Nonlinear Filtering Preserves Chaotic Synchronization via Master-Slave System
1Academia de Matemáticas, Universidad Politécnica de San Luis Potosí, Urbano Villalón 500, 78369 San Luis Potosí, SLP, Mexico
2División de Matemáticas Aplicadas, Instituto Potosino de Investigación Científica y Tecnológica, Camino a la Presa de San José 2055, 78216 San Luis Potosí, SLP, Mexico
3Departamento de Físico Matemáticas, Universidad Autónoma de San Luis Potosí, Alvaro Obregón 64, Col. Centro 78000 San Luis Potosí, SLP, Mexico
Received 17 December 2012; Revised 26 February 2013; Accepted 6 March 2013
Academic Editor: H. G. Enjieu Kadji
Copyright © 2013 J. S. González-Salas et al. 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.
We present a study on a class of interconnected nonlinear systems and give some criteria for them to behave like a filter. Some chaotic systems present this kind of interconnected nonlinear structure, which enables the synchronization of a master-slave system. Interconnected nonlinear filters have been defined in terms of interconnected nonlinear systems. Furthermore, their behaviors have been studied numerically and theoretically on different input signals.
In the study of nonlinear dynamics, much research was devoted to synchronization of chaotic systems, which have been unidirectional coupling , and one of the main applications is in secure communications systems . Under unidirectional coupling, there is a master system that forces a slave system. Let us begin by considering a master system whose temporal evolution is ruled by the following equation: where is the state vector, with defining a vector field . The slave system is given by where is the state vector, and the function describes the dynamics of the slave system and coupling. The master system behaves as an autonomous system and the slave system is a forced system that under certain conditions could behave as a nonautonomous system, where its dynamic is completely determined by the master system; that is, it means that the slave system acts in function of the master system. Hence, the slave system can be seen as a driven system by an external force, and then it can be studied as a filter when the input signal comes from a master system or any external signal unidirectional coupled.
There are several asymptotical behaviors reported based on master-slave configuration. Some cases are summarized as follows: Identical Synchronization (IS) implies coincidence of the corresponding states of the interacting systems . Lag Synchronization (LS) occurs when the trajectory of one oscillator is delayed by a specific time and is identical to the trajectory of the other oscillator . Phase Synchronization (PS) means the lock of chaotic oscillator phases, regardless of their amplitudes . Frequency Entrainment Synchronization (FES) occurs when two systems are oscillating with the same frequency . Generalized Synchronization (GS) is defined as the presence of some functional relationship between the states of the slave and master systems . Multimodal Generalized Synchronization (MGS) is presented when there are several basins of attraction for the slave system, and generalized synchronization is also presented [7, 8]. The aforementioned chaotic synchronization phenomena tell us that some types of synchronization are stronger than others. Another characteristic of these asymptotic behaviors is that they always can be described in terms of a functional relationship. For example, let be the GS function which relates the master system with the slave system; that is, . In general, the way how synchronization phenomena can be detected is by means of the auxiliary method ; thus the way one can know that unidirectional coupled systems present GS is when they present asymptotic behavior: where and are solutions of the slave and auxiliary systems, respectively, initialized with different initial condition . If this limit is satisfied, it indicates that the slave system has lost its sensitivity to initial conditions, thus losing their autonomy. Thence, the trajectory of the response system depends on the input signal from master system. If we consider the slave system as a forced system, then we will see that its behavior depends on an external forcing, like a filter does.
Forced Synchronization (FS) is defined as a phenomenon that occurs when oscillations of several chaotic systems , show correlated behavior because of an external signal applied to these oscillators . Meanwhile, the conditions in order that a chaotic system behaves as a filter for an external signal under a specific coupling are given in  where this phenomenon was called nonlinear filtering (NF). So, these three phenomena (GS, FS, and NF) can be recognized when the limit (3) is satisfied and can be studied from the viewpoint of forced chaotic systems. GS contains NF and MGS because these two phenomena present asymptotic behavior. Nevertheless, there exists a difference between NF and MGS. For MGS, the initial conditions determine the basin of attraction where the trajectory of the slave system converges, the fact despite that there exists a functional relationship between the slave and the master systems, the slave system does not behave like a filter. For NF, the slave system must act in function of the external force and never in function of its initial conditions.
The dynamic of at least one state of several chaotic systems has similar structure to low pass filter, for example, the first equation of the Lorenz system , where and can be seen as the output and input signals of a low pass linear filter, respectively. The second equation of the Lorenz system can be seen as a low pass nonlinear filter, , where is the output and and are the input signals of a low pass nonlinear filter. Without loss of generality, from the viewpoint of continuous dynamical systems, a filter can be seen as a forced dynamical system in which its response depends on its structure, given by their equations, and can be linear or nonlinear. Then a natural question would be the following: what is the effect on the response of forced nonlinear system given an input signal? In order to answer this question, we make a study of a chaotic system forced with different kinds of signals.
We are interested in explaining synchronization phenomenon of chaotic systems based on nonlinear interconnected structures, as these structures behave as a nonlinear filter that allows synchronization of master and slave systems. Thereby, the target is not to put forward a new nonlinear filter that can replace those designed with the specific purpose of preventing noise in the signal. Thus, the objective of this research is to show how the synchronization phenomenon underlies nonlinear interconnected structure. This kind of structure which is immersed in some chaotic systems could have some useful applications in filtering waves such as those originated from earthquakes and tsunamis, because they can be tuned to produce resonance at certain frequencies. So, the target of this work is to give features of the nonlinear filter as those given by the second and third equations of the Lorenz system in function of the parameters of the input signal , which gives us a complementary perspective with respect to the analysis made in , where it was only proved and showed that a forced-chaotic system filters constant, sinus, and random signals. To achieve this goal, we make a study on the effect of nonlinear filter based on interconnected nonlinear systems for different input signals , like sinusoidal, chaotic, and random signals. Several chaotic systems that have a similar structure to interconnected systems can be analyzed like a nonlinear filtering phenomenon.
The paper is organized as follows: Section 2 contains basic definitions of nonlinear filters based on -interconnected systems; Section 3 presents a way how tuning the parameters of an -interconnected system in order to the system behaves as a nonlinear filter; in Section 4 is the study of the response to amplitude and frequency of a nonlinear filter; Section 5 shows the response of the filter when the input signal is noise; Section 6 presents a study case of correlation coefficient analysis; in Section 7, we present the relation between nonlinear filters and chaotic systems; finally, conclusions are given in Section 8.
2. Nonlinear Filter Structure
In the theory of linear systems, it is very well known that a first-order low pass linear filter is given as follows: where is the output of the filter, and are parameters of the filter (4). These parameters control attenuation and amplitude of the input signal, respectively, and is the input signal to be filtered. Based on the configuration of low pass filter given by (4), a low pass nonlinear filter is defined as follows.
Definition 1. Let be an output signal, an input signal, and , , and parameters. Thus, a low pass nonlinear filter can be defined as follows: The parameter in the above definition is used to control the nonlinear term . Notice that when , the nonlinear term disappears, , and then the nonlinear filter (5) behaves as a linear filter. The structure of nonlinear filters (5) has been used to generate chaos, as it can be seen in the Lorenz system , but its states are interconnected. Therefore, a natural question emerges around nonlinear filtering: what are the characteristics of two nonlinear filters if they are interconnected and forced by the same signal ? An interconnected system via nonlinear filters is given as follows:
Definition 2. Let be output signals, an input signal, and , , , , , and parameters. So, let us define an interconnected system via low pass nonlinear filters as follows: Thus the system (6) is called two-interconnected systems.
Now, in Definition 2, the parameters and are coupling parameters for the outputs and , respectively. In general, an -interconnected system can be defined by coupled to the output with . This -interconnected system has a core based on low pass linear filters , where is a state vector, is a matrix, and is a constant vector. Thus, we can define -interconnected system with the following expression: where the nonlinear function is constituted by terms , , . An -interconnected system (7) is studied in  where the authors found that its model produces hyperbolic chaos when it is forced by a sinusoidal wave. On the other hand, system (7) is formed by a linear part , an input signal , and a nonlinear part , which can induce sensitivity to initial conditions . But a filter must act in function of the input signal and not in function of the initial condition. Therefore, it is important to find conditions in order to guarantee that an interconnected system behaves as a filter. Roughly speaking, the response of a filter only depends on the kind of input signal and if the interconnected systems depend on initial condition, this is classified as generalized forced synchronization phenomenon by forced systems . Therefore, according to context, a nonlinear filter is defined in the next way.
Definition 3. Let be a vector of output signals given by -interconnected system (7) with initial condition . System (7) is called an -interconnected nonlinear filter if it always presents asymptotic behavior:
3. Tuning of the Parameters
We select the entries of the matrices and in order for condition (8) to be satisfied. Let us start by considering , and then . Now, by using (7) we deduce that which can be rewritten as follows: In order to describe the asymptotic behavior of (9), we state the following theorem.
Theorem 4. If system (7), forced by signal , satisfyies the following conditions:(1)there exists a positive constant such that , for every eigenvalue of the linear part of system (7);(2) is a continuous Lipschitz; that is, there exists a positive function such that , then, system (7) is a filter for , provided that , where
Proof. Let , , and defined be before. We want to prove that implies . From (9) we have the following estimate: Using hypothesis , we obtain from the previous inequality the following estimate: If we multiply by , the last expression results in Using the hypothesis , it follows that Application of Gronwall's inequality yields which proves the theorem.
We can infer from the above theorem that all elements of the matrix must be negative; that is, for . On the other hand, due to the fact that (9) does not depend on matrix , then its terms are not conditioned and they can take any real value. These will be carried out in the remainder of this paper.
3.1. A Low Pass Nonlinear Filter
The analyzed case is the one-interconnected system given by (7); this case is confined to the particular case of the filter given by (5); where its elements are , , , and . According to the first condition of the above theorem, we have , and now considering the nonlinear part of the filter, if the Lipschitz condition is satisfied when , giving as a result , therefore the parameter can be any real value. Now, the interest is that the one-interconnected system behaves as a nonlinear filter; then Theorem 4 needs to be satisfied by . If the input signal is a constant, () and , then the response of the filter converges to zero; otherwise it diverges. Another case is when the input signal is , we have the following: which converges . Generally, condition (8) is always satisfied for oscillating functions with . However, for the case that is a polynomial of grade greater than , that is, for , condition (8) is not satisfied.
Due to the fact that our interest is to tune the value of the parameters , and , we make a numerical study on the effect of these parameters on the response of when the one-interconnected system (7) is being forced by the sinusoidal input signal . After a transient time, we calculate the maximal range of values of the one-interconnected system (7) by obtaining (), where and are the maximum and minimum values, respectively, of the response time series in function of the parameters . For example, if we want to calculate , then we fix the value of and simulate different time series of the one-interconnected system (7) for different values of . In general, the value of is calculated in function of the parameters which are varied and in each case condition (8) is verified if it is satisfied. Figures 1(a) and 1(b) show the graphs of for and for , respectively. These graphs show that the one-interconnected system (7) increments exponentially the amplitude of its response when the parameter increases its magnitude; meanwhile, the amplitude of is mildly incremented in a linear rate in function of the magnitude of the parameter .
3.2. Interconnected Nonlinear Filters
Now, the first case is a two-interconnected system (6); the matrices , , and the nonlinear function are assumed as follows:
System (6) is asked to behave as a nonlinear filter, so the conditions of Theorem 4 are used to select the value of the parameters . From the first condition of the theorem, we realize that , , and if , where and , then system (6) is a filter for the signal . In order to know more about the characteristics of the values of the parameters, we calculate the eigenvalues of the Jacobian of system (6), which results in the following:
Our purpose is that system (6) converges; then Re must be negative. Therefore, we need to guarantee that is true. Due to , then the parameters and need to be opposite sign, that is, either , or , . We choose arbitrarily the first relation of signs.
The effect of the parameters on the response of the interconnected nonlinear filters is computed by when system (6) is forced by . Several calculations were made and for each of them we verify that the nonlinear filtering condition (8) was satisfied. Figure 2(a) shows when the parameters of the filter (6) are fixed to , . Figure 2(b) shows for and .
In the remainder of this paper, we studied the effect of the parameters input signal on the response of the coupled filter (6); that is, we tune the parameters of the filter, vary the parameter's values of the input signal, and observe the effect of the response. Therefore, we need to fix (). Our interest is focus on the study of the response to the input when the linear filters , are coupled in a nonlinear way. Without loss of generality and seeking clarity in our study, we consider . The rest of the parameter values could vary in the intervals , , , , which is a rich variety of values where the filter works. We can select any value in these intervals and produce similar responses; nevertheless, we fix the values to , , , , , which will be used in the rest of the paper.
4. Response to the Amplitude and Frequency
In this section, we present a study on the effect of the filter's response (6) as function of parameter's values of the sinusoidal input signal . Figure 3(a) shows the output signals of the -interconnected nonlinear filter for the cases , and , . One can see that the orbits have a form of a Lissajous curve and the amplitude of the output signal as the frequency of the input signal is increased. In general, to see the effect of amplitude on the filter's response, we calculate the length of the Lissajous in this way where is the period of the orbit (, ). The effect of the input signal's frequency on the length of the orbit of the filter's response is shown in Figure 3(b). Thus, falls exponentially according to the increment of the frequency . System (6) presents a rejection to high frequencies because it is a low pass filter.
Now, fixing the frequency to and varying the amplitude , we observe that new frequency components appear in the Fourier Transform (FT) of the output signal . When the amplitude increases, then more peaks of frequency appear as multiples of the input signal's frequency. As is shown in Figures 4(a) and 4(b), and , respectively. This is a characteristic of a nonlinear filter that linear filters do not present.
The nonlinear filters display assorted behaviors, that is, contrary to the example shown in Figure 4, where the amplitude is fixed to and frequency is varied. Figures 5(a) and 5(b) show the FT of the response of the filter (6) for and , respectively. Comparing both plots of Figure 5, we see that frequency peaks appear at multiples of input signal's frequency .
4.1. A Theoretical Justification
For the purpose of giving a theoretical justification of the behavior of the filter (6) when it is forced by the sinusoidal signals, we can see in the interconnected nonlinear filter (6) that and are linear first-order differential equations which can be rewritten in the next form: Without loss of generality and for seeking clarity, we consider the parameter . Applying integration by parts to (20) and considering that the states and are functions that depend on time and frequency and that is a constant, we have the following: And for (21) we have
Note that in (22) and (23), the magnitude values of and decrease when the frequency increases. Furthermore, since the terms have factors , for , then we approximate (22) and (23) to where , are periodic functions with period given by the input signal. Substituting (25) in (26), we have Now, solving for from (27), it results that And for it results that Figure 6 shows the orbit of the solution of (28) and (29) and the orbit of the numerical solution of the filter (6) when and . The thin Lissajous curve is calculated with (6) and another with (28) and (29). We can see that both orbits oscillate in a Lissajous curve of one knot and the same symmetric form with respect to the -axis, but slightly different amplitude.
Now, we analyze the FT of the equations of system (6) which are given as follows: where is the FT of the signal . For the case when , we have Now, solving for and , we have Equations (33) and (34) show that variables , have the form of low pass filters (as we commented previously). Thus, the terms and in (33) are as follows: and using (33) and (35), we obtain where Equation (36) has the terms and , which are calculated by developing recursively (33), giving as a result
Equations (33) and (38) show that have terms which contain factors , , , and and if we develop recursively times (33), then new terms of the form for appear in the solution of , where is a so complex factor which depends on every parameter and variable of (32) and for the quotient . Therefore, when the amplitude is incremented, then the amplitude of the peaks does too. For this reason, the numerical evidence shows that if the parameter input signal is incremented, then apparently new components of frequency appear in the spectrum . These frequencies always form part of but with a very small amplitude. They are so small that they look like noise.
5. Response to Random Signals
Now, let , where is a random signal equidistributed on the interval and is its amplitude. When is the input signal of the filter (6), the phase space of the response is a limit cycle that has the form of a quarter moon which is shown in Figure 7(a). As a first glimpse the filter's response seems to be a periodic orbit like in the case of the sinusoidal input signal, but the time series and , which are respectively shown in Figures 7(b) and 7(c), are aperiodic signals.
The following remarks are pertinent. (i)The output signal of the filter has a normal (or Gaussian) distribution when the input signal is random with uniform distribution. For the case , Figures 8(a) and 8(b) show the histogram and the FT of the signal, respectively, where one can see that the FT of the response has an exponential decay.(ii)The filter presents saturation after a threshold value of the amplitude of input signal . Figure 9 shows the range of values of the response versus the amplitude of the parameter . To justify the exponential decay of the filter's response when it is forced with random signal, we consider , where is a constant which depends on , the amplitude of the random signal, and . Thus the following terms of (30) are and , where depends on too. So, the FT results are and for the magnitude of we have Figure 8(b) shows the form of the curve which corresponds to (40).
6. Correlation Coefficient
Because a filter acts in function of the input signal, another way to characterize the effect of the input signal on the filter's response is to compute the correlation between the input and output signals. If we have time series of data, the correlation coefficient of series and is defined as where , are the means and , are the standard deviations of the series and , respectively. We calculate correlation coefficients for different input signals and the output signals and which are given in Table 1. The first two cases in Table 1 correspond to a periodic and chaotic signals, and , and the last case considers a random signal as an input signal.
We make a further numerical study of the coefficient correlation in function of the parameters of the sinusoidal input signal. In Figure 10(a), we can see that the correlation between sinusoidal external signals is not constant and that it depends on both parameters and . On the contrary, Figure 10(b) shows that the absence of correlation between the response and their respective forces is independent of almost the whole range of values of the input signal parameters.
On the other hand, the last row in Table 1 shows that and . This means that the filter's response (6) does not have correlation when the input is a random signal. Despite the fact that the amplitude of the parameter is increased, the output signal of the filter does not have correlation with input signal (random), as shown in Figure 11. Therefore, the correlation between the input random signal and its filter's response cannot be induced by incrementing the value of its amplitude. A similar result occurs for the sinusoidal signal and its response .
7. Low Pass Filters in Chaotic Systems
The Lorenz system is a very well-known third-order chaotic system  which is defined as where is the state vector and , , and are parameters. Each component of system (42) is defined based on low pass filters. For example, the first equation is a linear low pass filter, where the state variables and are the input and output signals, respectively. The second and third equations comprise a -interconnected nonlinear filter whose outputs are and , and the input signal is . The -interconnected nonlinear filter of the Lorenz system (42) can be rewritten in terms of (7) as follows: where , , and . The projection of the Lorenz attractor onto the plane is shown in Figure 12(a) for the following parameter's values: , and , and . In order to study the response of the interconnected system (43) when the input signal is the state variable , we have calculated which is the maximal range of values of the time series . In Figure 12(b), we can see that for ; this means that the solution converges to a fixed point. For , there is a transition period in the behavior of the filter, and for the range of values of the filter's response is approximately constant.
When calculating the correlation coefficient between the components of the -interconnected system (43), it results that , , and . The high correlation between the state variables and occurs because of is a simple low pass lineal filter for the signal without any nonlinear term included in its equation (see in (42)). On the other hand, the time series has no correlation with which could be induced by the nonlinear term settles in the equation of system (42). Now, we put a parameter to control the amplitude of the input signal such that and calculate and as a function of the parameter . In Figure 13, it is shown that only is affected when the value of is incremented, while the absence of correlation between states and remains for all values of . For , the trajectories of system (43) diverge.
We have shown by means of numerical experiments in previous sections that the absence of correlation between the external force and the state variable of the filter response generally does not depend on the parameters of the input signal. Similar result occurs when the Lorenz autonomous system is considered like a -interconnected system which is forced by the signal . As shown in Figure 13, there exists a high correlation between the state variables and , and an absence of correlation between state variables and . But can the correlation between the state variables and be incremented by changing the value of the parameters of the Lorenz system? With the purpose of answering this question, we have calculated , as a function of the autonomous Lorenz system parameter keeping the value of the parameter constant and vice versa. Figure 14(a) shows that for approximately, correlations , have a transition period, and for , the correlations return to their previously showed behavior: and . Figure 14(b) shows that the parameter does not affect the correlation between and the responses , .
If Lorenz system parameters are tuned so that the Lorenz system behaves as a filter and can be used as a slave system, then generalized synchronization of master and slave systems always appears.
An -interconnected nonlinear system given by (6) always behaves as a filter if the condition given by Theorem 4 is satisfied. Despite of the nonlinearities in the structure of -interconnected nonlinear system, its dynamical behavior does not depend on the initial conditions but it does as a function of the input signal. The nonlinear filter's response to sinusoidal input signal presents several interesting phenomena such as frequency components at multiples of the input signal's frequency when the amplitude of the input signal increases, the number of peaks of the filter's response increases the amplitude of the filter's response falls exponentially as a function of the input signal's frequency.
The Lorenz system is described in terms of low pass filters which consists of a linear low pass filter and a -interconnected nonlinear low pass filter. This gives us the possibility to describe all the systems that conform the Lorenz family and others with similar structure in terms of low pass filters.
In several studies of chaos synchronization, specifically in forced systems , , it has been found that generalized synchronization of response systems occurs for specific external signals , but all these cases satisfy condition (8). Then, the general synchronization phenomenon and the phenomenon of nonlinear filter are seen to be the same. Nevertheless, condition (8), in generalized synchronization phenomenon, could not be satisfied, because of the trajectories generated with different initial conditions asymptotically could go to a different basin of attractions. This characteristic is different to the phenomenon of nonlinear filter which by definition condition (8) needs to be always satisfied if a system behaves like a nonlinear filter. Several chaotic systems have a similar structure to -interconnected system, so that we conjecture that this study can help to distinguish between generalized synchronization behavior and nonlinear filtering behavior of an -interconnected system.
The authors would like to thank Emanuel Rodríguez-Orozco for reading the paper and giving constructive comments. E. Jiménez-López is a doctoral fellow of CONACYT in the Graduate Program on Applied Science at UASLP, and IPICYT for the hospitality during his sojourn in DMAp-IPICYT. E. Campos-Cantón acknowledges CONACYT for the financial support through Project no. 181002.
- L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990.
- I. Campos-Cantón, E. Campos-Cantón, J. S. Murgúa-Ibarra, and M. E. Chavira-Rodríguez, “Secure communication system using chaotic signals,” Ingeniería, Investigación y Tecnología, vol. 10, no. 1, 2009.
- M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “From phase to lag synchronization in coupled chaotic oscillators,” Physical Review Letters, vol. 78, no. 22, pp. 4193–4196, 221997.
- M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “Phase synchronization of chaotic oscillators,” Physical Review Letters, vol. 76, no. 11, pp. 1804–1807, 1996.
- V. S. Anishchenko, T. E. Vadivasova, D. È. Postnov, and M. A. Safonova, “Synchronization of chaos,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 2, no. 3, pp. 633–644, 1992.
- H. D. I. Abarbanel, N. F. Rulkov, and M. M. Sushchik, “Generalized synchronization of chaos: the auxiliary system approach,” Physical Review E, vol. 53, no. 5, pp. 4528–4535, 1996.
- E. Campos, J. Urías, and N. F. Rulkov, “Multimodal synchronization of chaos,” Chaos, vol. 14, no. 1, pp. 48–54, 2004.
- L. Ontañón-García, E. Campos-Cantón, R. Femat, I. Campos-Cantón, and M. Bonilla-Marín, “Multivalued synchronization by poincar coupling,” Communications in Nonlinear Science and Numerical Simulation.
- J. S. González Salas, E. Campos Cantón, F. C. Ordaz Salazar, and I. Campos Cantón, “Forced synchronization of a self-sustained chaotic oscillator,” Chaos, vol. 18, no. 2, Article ID 023136, p. 9, 2008.
- E. C. Cantón, J. S. C. Salas, and J. Urías, “Filtering by nonlinear systems,” Chaos, vol. 18, no. 4, Article ID 043118, p. 4, 2008.
- E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of the Atmospheric Sciences, vol. 20, pp. 130–141, 1963.
- L. V. Turukina and A. Pikovsky, “Hyperbolic chaos in a system of resonantly coupled weakly nonlinear oscillators,” Physics Letters A, vol. 375, no. 11, pp. 1407–1411, 2011.