Synchronization of Chaotic Energy Resource System

Taking a four-dimensional energy resources demand-supply system between the East and West of China, this paper discusses its chaotic behavior, and via the unilateral coupling method we lead the system to synchronization successfully. But not all the values of coupling coefficient can lead to synchronization. The values of coupling coefficient have a range. By calculating the maximal relative Lyapunov exponents’ spectrum, we gained the value range of coupling coefficients. Within the value range, the two coupling systems can achieve synchronization, otherwise can’t. Further more, the values of coupling coefficient are in connection with the chaos synchronizing time. At last, we get the relationship of coupling coefficients and chaos synchronizing time.


Introduction
Generally, designing a controller to force a system to imitate the behavior of another chaotic system is called synchronization [1].Chaos synchronization and chaos control have been applied broadly in many fields such as in biological systems, chemical reactions, information processing, especially in secure communication area [2].There are plenty of methods have been proposed to achieve chaos synchronization, for example, sliding mode control [3], linear control [4], adaptive control [5], digital redesign control [6] and so on.
Nowadays, many countries attach importance to the exploitation and utilization of clear energies, and develop low carbon economy.Considering renewable energy resources, based on a three dimensional system, gained a four-dimensional energy resources demand-supply system between east and west in China by adding a new variable [7].This paper discusses the chaotic behavior of it, and uses the unilateral coupling method achieving chaos synchronization successfully.By means of calculating the maximal relative Lyapunov exponents, we get the value range of the coupling coefficient.Also, we get the relationship of coupling coefficient and chaos synchronizing time.

The Unilateral Coupling Method
In 1993, K. Pyragas proposed a way to control a nonlinear system, which is called the error variable negative feedback control also the unilateral coupling method [8].
Supposing there are two chaotic systems: , which are defined by following dynamic functions: where  is coupling coefficient, which are equal upon each variable and take positive values commonly, i.e.
When the parameters of the two systems are matched, as long as take an appropriate coupling coefficient K, the two systems can achieve synchronization.This method won't change the system's initial dynamic characteristic, because the coupling item K(X -Y) = 0 after synchronization.The characteristic of the method is we needn't analyze the system in advance.Further more, we can confirm the domain of the coupling coefficient by calculating.

Discussion of the Energy Resource Demand-Supply System between East and West in China
Adding a new variable to a three dimensional energy resource system, gained a four dimensional energy resources demand-supply system between east and west in China, which is defined by the functions below [7]: where x(t) is the energy resource shortage in A region, y(t) expresses the energy resources supply increment in B region, z(t) the energy resources import in A region, w(t) is the amount of renewable energy resources in A region;  , the system can generate complex chaotic attractor, the numerical simulation is shown in Figure 1, from which we can see the system has an abundant chaotic behavior.
The errors of corresponding variables between system (3) and system (4) are denoted as and which express the situation of chaos synchronization.When E i (t) and E equals to zero and then stable over time, we can say systems (3) and ( 4) have achieved synchronization, otherwise not.Randomly, chose c = 2, the errors E i (t) is shown as Figure 2, and the errors E is shown in Figure 3. Figure 2 (or Figure 3) indicates clearly that when c = 2, E i (t) (or E) tends to zero immediately and then stable, in other words, systems (3) and ( 4) achieved synchronization quickly.Thus, via the unilateral coupling method, we have made two systems achieve synchronization successfully.
But not all the values of coupling coefficient can lead to synchronization according to the following discussion.

Confirmation of Coupling Coefficients
Then let c = 0.01, the numerical simulation of the errors E i (t) is shown as From Figure 4 (or Figure 5), we can clearly see that when c = 0.01, E i (t) (or E) vibrates desultorily all the time, namely this value of coupling coefficient c won't      In a chaotic system, the maximal Lyapunov exponent is a quantity characterizes the rate of separation of infinitesimally close trajectories, also the quantity indicates the strength of butterfly effect.In coupled systems, when maximal relative Lyapunov exponent is less than zero, the two systems will achieve synchronization, otherwise won't [9].The definition of maximal relative Lyapunov exponent [9] is where D t is the distance of the two system's trajectories at the time of t, D 0 is the distance of the two system's trajectories at the initial time, E i (t) the error of corresponding variables between systems at t moment, E i (0) the error of corresponding variables between systems initially, N is dimension.The numerical relation between coupling coefficients and the maximal relative Lyapunov exponents is shown in Figure 7.
As shown in Figure 6, when c > 0.043, the maximal relative Lyapunov exponent is less than zero, thus systems ( 3) and ( 4) can achieve synchronization successfully, while c < 0.043, the maximal relative Lyapunov exponent is more than zero, the two systems can't achieve synchronization.So when c = 0.01, the maximal relative Lyapunov exponent is more than zero, the two systems haven't achieved synchronization, as shown in Figure 4 (or Figure 5).

The Relation between Coupling Coefficients and Synchronizing Time
When c = 0.2, the errors E i (t) is shown as Figure 8, the errors E is shown as Figure 9.
Comparing Figure 7 (or Figure 8) and Figure 4 (or Figure 5), we see that when c = 2, the time achieving synchronization is shorter than the instance of c = 0.2.
In fact, sometime we need achieve synchronization immediately, while sometime we need transit a certain time and then achieve synchronization.So, the choice of a coupling coefficient's value is crucial, and it's neces-   sary to discuss the relation between coupling coefficients and the time achieving synchronization.
Figure 10 presents the relation between coupling coefficients and the synchronizing time.Though there comes some gurgitation, but generally, the bigger coupling coefficient is, the shorter two systems achieve synchronization.It is an important conclusion, in practice, we can choose an appropriate coupling coefficient according to demand.

Conclusion
Based on a four-dimensional energy resource demandsupply system between the East and West of China, via numerical simulation, this paper discusses its chaotic behavior, and uses the unilateral coupling method achiev- ing the chaos synchronization.By computing the maximal relative Lyapunov exponents' spectrum of systems ( 3) and ( 4), we gained the value range of coupling coefficients.When coupling coefficient c > 0.043, the maximal relative Lyapunov exponent is less than zero, thus the two systems can achieve synchronization sucessfully, while c < 0.043, the maximal relative Lyapunov exponent is more than zero, the two systems can't achieve synchronization.At last, we get the relationship between coupling coefficients and synchronizing time, i.e. the greater coupling coefficient is, the shorter systems achieve synchronization.
a i , b i , c i , d i and M, N are positive real constants [7].When M = 1.8, N = 1, a 1 = 0.1, a 2 = 0.15, b 1 = 0.06, b 2 = 0.082, b 3 = 0.07, c 1 = 0.2, c 2 = 0.5, d 1 = 0.1, d 2 = 0.06, d 3 = 0.07, and the initial condition y b y N y y c x y t y c c y y c c y c x y t y d y d y c x y t

Figure 4 ,
the errors E is shown as Figure 5.The Bifurcation diagram under various coupling coefficient from 0 -0.1 is shown as Figure 6.

Figure 2 .
Figure 2. When the coupling coefficient c = 2, the errors of corresponding variables between systems (3) and (4), E i (t), varies over time.

Figure 3 .
Figure 3.When the coupling coefficient c = 2, the errors of corresponding variables between systems (3) and (4), E, varies over time.

Figure 4 .
Figure 4.When the coupling coefficient c = 0.01, the errors of corresponding variables between systems (3) and (4), E i (t), varies over time.

Figure 5 .
Figure 5.When the coupling coefficient c = 0.01, the errors of corresponding variables between system (3) and system (4), E, varies over time.

Figure 6 .
Figure 6.The relation between coupling coefficients and x 2 of the system (3).
lead the two systems to synchronization.So, we conclude not all the values of coupling coefficient can lead the two systems to synchronization; Values of c have a domain, values in which can make the two systems achieve synchronization, values outside it are not appropriate.By calculating the maximal relative Lyapunov exponents' spectrum, we can get the value range of c.

Figure 7 .
Figure 7.The relation between coupling coefficients and the maximal relative Lyapunov exponents of systems (3) and (4).

Figure 8 .
Figure 8.When the coupling coefficient c = 0.2, the errors of corresponding variables between systems (3) and (4), E i (t), varies over time.

Figure 9 .
Figure 9.When the coupling coefficient c = 0.2, the errors of corresponding variables between systems (3) and (4), E, varies over time.

Figure 10 .
Figure 10.The relation between coupling coefficients and synchronizing time of systems (3) and (4).