Improving the complexity of the Lorenz Dynamics

Hindawi Publishing Corporation Complexity Volume 2017, Article ID 3204073, 16 pages https://doi.org/10.1155/2017/3204073 Research Article Improving the Complexity of the Lorenz Dynamics María Pilar Mareca and Borja Bordel Department of Physic Electronics, Universidad Politécnica de Madrid, Avenida Complutense No. 30, 28040 Madrid, Spain Correspondence should be addressed to Marı́a Pilar Mareca; [email protected] Received 26 July 2016; Accepted 19 September 2016; Published 10 January 2017 Academic Editor: Michael Small Copyright © 2017 M. P. Mareca and B. Bordel. his 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. A new four-dimensional, hyperchaotic dynamic system, based on Lorenz dynamics, is presented. Besides, the most representative dynamics which may be found in this new system are located in the phase space and are analyzed here. he new system is especially designed to improve the complexity of Lorenz dynamics, which, despite being a paradigm to understand the chaotic dissipative lows, is a very simple example and shows great vulnerability when used in secure communications. Here, we demonstrate the vulnerability of the Lorenz system in a general way. he proposed 4D system increases the complexity of the Lorenz dynamics. he trajectories of the novel system include structures going from chaos to hyperchaos and chaotic-transient solutions. he symmetry and the stability of the proposed system are also studied. First return maps, Poincaré sections, and bifurcation diagrams allow characterizing the global system behavior and locating some coexisting structures. Numerical results about the irst return maps, Poincaré cross sections, Lyapunov spectrum, and Kaplan-Yorke dimension demonstrate the complexity of the proposed equations. 1. Introduction A chaotic system is a highly sensitive nonlinear system. he main characteristic of chaos is the sensibility to the initial conditions. his sensibility is sometimes known as “butterly efect” in honor of Lorenz [1] who described in 1963 the irst natural phenomenon showing a chaotic behavior (in this case, the atmospheric evolution) being represented by a diferential system (1) which generates a trajectory in the phase space similar to a butterly (see Figure 1): �̇ = � (� − �) , �̇ = �� − � − ��, �̇ = �� − ��. (1) he Lorenz system has been extensively employed: cyphers [2, 3], circuits [4], engineered systems [5], and so forth have been based on this dynamics. However, despite being the irst described chaotic system and, probably, the most used and the most popular one, it presents a complexity degree pretty low. For example, as we will see in Section 2, some encryptions based on Lorenz system may be broken by synchronizing the cypher with an external system (although its parameters were unknown). he underlying cause of this weakness is the great amount of redundant information present in the components of their dynamics that it is precisely associated with the low degree of complexity of the system. Even more, some works have shown that it is possible to recover the three systems’ components from only one of them [6]. In order to increase its degree of complexity, several works have tried to modify the Lorenz dynamics. Some of them try to increase the number of ixed points in the system [7]. Others are focused on including commutated terms which modify the number of wings in Lorenz’s attractor [8]. Some articles including new nonlinearities in the Lorenz system to improve its complexity can be also found [9]. As an extreme case, techniques to build � × � Lorenz-like attractors have been also reported [10]. Among all these previous works, a very important group of articles are those focused on the creation of Lorenz-like “hyperchaotic systems.” Traditionally, hyperchaos was deined [11] as a chaotic behavior characterized for more than one positive Lyapunov exponent; that is, the dynamics expand in more than one direction (therefore, hyperchaos can only be found in high order systems). However, in recent proposals, hyperchaos has been also applied to high-order 2 Complexity 2. Insecurity in Lorenz-Based Cyphers 300 It has been proven in many articles that Lorenz system can be self-synchronized [4]. In particular, some of the components of the dynamics contain so much information about the low (values of the parameters, topology being generated, etc.) in which the most weak synchronization schemes are useful when applied to the Lorenz dynamics. For example, even active-passive decomposition schemes [19], where only one of the terms in one of the equations is modiied (in contrast to other schemes [20], where one or various equations are totally removed), allow fast synchronization. 250 200 Z 150 100 50 −60 −40 −20 0 20 40 60 Figure 1: XZ composition of the Lorenz system solution for {� = 10, � = 178, and � = 8/3}. X systems which present only one positive Lyapunov exponent [12], besides the zero exponent, but whose value is much higher than the common values (e.g., for Lorenz system, the maximum value of the positive Lyapunov exponent is around � max = 5.2 when {� = 25, � = 420, and � = 8/3} and exhibits a Kaplan-Yorke dimension of �KY = 2.17). here have been proposed many hyperchaotic systems, with some designed to facilitate the electronic implementation [13] and others designed to create certain topologies in the phase space [14] and, even, new hyperchaotic dynamics have been described to be employed in diferent encryption schemes [15], with most of them being based on Lorenz equations [16– 18]. Nevertheless, these proposals do not cover the objective of increasing the complexity at all. While these new dynamics are more complex, that increased complexity is further supported by the fact of having at least one additional equation instead of the appearance of trajectories whose entropy is higher. Actually, the value of the maximum Lyapunov exponent in these cases is around � max = 2.4 and the rate of contraction is high, obtaining, thus, weak hyperchaos. herefore, the aim of this article is to describe a new hyperchaotic system, being able to generate highly complex and novel structures. For that, we seek to strengthen at least one direction of expansion and reduce the rate of contraction (but they should remain dominant in order to preserve the system as dissipative). he rest of the paper is organized as follows. In the next section, a study about the redundant data in the components of Lorenz’s dynamics is presented. In Section 3, the new hyperchaotic system is introduced and its basic topological properties are studied. In Section 4, the regular structures are located and analyzed. In Section 5, chaotic attractors are described. In Section 6, complex chaotic and hyperchaotic trajectories are analyzed, describing the special case of the chaotic-transient solutions. Section 7 concludes the paper. heorem 1. Consider two Lorenz systems (2), coupled using an active-passive decomposition scheme. If enough amount of time is elapsed, then both systems get completely synchronized (in amplitude and phase): ��̇ = �0 (�� − �� ) , ��̇ = �0 �� − �� − �� �� , ��̇ = �� �� − �0 �� , ��̇ = �0 (�� − �� ) , ��̇ = �0 �� − �� − �� �� , (2) ��̇ = �� �� − �0 �� . Proof. A new function, called error function, is deined as follows: �⃗ = (�1 , �2 , �3 ) = (�� , �� , �� ) − (�� , �� , �� ) . (3) hus, both coupled Lorenz systems get synchronized in amplitude and phase if �⃗ → 0⃗ when � → ∞. he error function satisies the following diferential system: �1̇ = −�0 �1 , �2̇ = �0 �1 − �2 − �� �� + �� �� , �3̇ = �� �� − �� �� − �0 �3 . (4) he irst diferential equation is decoupled, so it may be solved directly (5). Considering �0 > 0 (due to its physical meaning), �1 tends to zero (6) and the limit condition �� → �� is satisied: �1 (�) = ��−�0 � , lim � �→∞ 1 = lim ��−�0 � = 0. (5) (6) hen, the resulting diferential system (7) at � → ∞ may be used to prove the stability of the rest of the components: �→∞ �1 = 0, �2̇ = −�2 − �� �3 , �3̇ = �� �2 − �0 �3 . (7) Complexity 3 A Lyapunov function satisfying two conditions is proposed as follows: 1 � (�)⃗ = �⃗ ⋅ �.⃗ 2 (8) (1) �(�)⃗ presents a strict minimum at the origin, so that �(0)⃗ = 0 and �(�)⃗ > 0 ∀�⃗ ≠ 0.⃗ ⃗ (2) ��(�)/�� < 0 in any reduced neighborhood of the ⃗ origin. In fact, ��(�)/�� = �⋅⃗ �̇⃗ = �2 �2̇ +�3 �3̇ = �2 (−�2 − �� �3 ) + �3 (�� �2 − �0 �3 ) = −�2 �2 − �0 �3 �3 ∀�0 > 0. If �0 > 0, then Lyapunov’s theorem about the asymptotic stability guarantees that lim�→∞ �⃗ = 0⃗ and then (�� , �� , �� ) → (�� , �� , �� ) when � → ∞. heorem 1, in fact, enables the creation of cryptographic systems based on self-synchronized Lorenz system containing encrypted information. However, heorem 1 might be understood in another way. If so much information is carried in each one of Lorenz’s system components, are Lorenz-based cryptosystems safe? heorem 2. he information about the system in any of the Lorenz dynamic components allows any intruder system (9) not only to synchronize the temporal evolution but also to obtain the values of the employed parameters in the transmitter: ��̇ = �0 (�� − �� ) , ��̇ = � (�) �� − �� − �� �� , ��̇ = �� �� − � (�) �� . (9) Proof. Two irst diferential equations (10) determining the evolution laws of �(�) and �(�) are deined: �̇ (�) = −�2 �� , �̇ (�) = �3 �� . (10) Besides, both error functions for the temporal evolution of the systems and for the parameters values are also deined as follows: �⃗ = (�� , �� , �� ) − (�� , �� , �� ) , �⃗ = (� (�) , � (�)) − (�0 , �0 ) . (11) If it is satisied that �⃗ → 0⃗ and �⃗ → 0⃗ when � → ∞; then heorem 2 is veriied. For that, the expressions of the diferential systems for the error functions are deduced: −�0 �1 �̇⃗ = (� (�) �� − �0 �� − �2 − �� �� + �� �� ) , �� �� − � (�) �� − �� �� + �0 �� −�2 �� ̇ �⃗ = ( ). �3 �� (12) As it is said previously, the irst equation is decoupled and it can be proven that �1 tends to zero when � → ∞. hen, when � → ∞, the diferential system for the error functions may be simpliied: �̇⃗ = ( −�0 �1 �1 �� − �2 − �� �3 −� (�) �� + �� �2 + �0 �� −�2 �� ̇ �⃗ = ( ). �3 �� ), (13) A Lyapunov function is proposed (14), satisfying two conditions: ⃗ . ⃗ = 1 (�⃗ ⋅ �⃗ + �⃗ ⋅ �) � (�,⃗ �) 2 (14) (1) �(�,⃗ �)⃗ presents a strict minimum at the origin, so that �(0,⃗ 0)⃗ = 0 and �(�,⃗ �)⃗ > 0 ∀�,⃗ �⃗ ≠ 0.⃗ ⃗ (2) �(�,⃗ �)/�� < 0 in any neighborhood of the origin. In ̇ ⃗ fact, ��(�,⃗ �)/�� = �⃗ ⋅ �⃗̇ + �⃗ ⋅ �⃗ = �2 �2̇ + �3 �3̇ + �1 �1̇ + �2 �2̇ = �2 (�1 �� −�2 −�� �3 )+�3 (−�(�)�� +�� �2 +�0 �� )+ −�2 �� �1 + �3 �� �2 = −�22 + �3 �0 �� + �3 �� (�2 − �(�)) = (−�22 − �0 ⋅ �23 ) < 0 ∀�0 > 0. hen, Lyapunov’s theorem about the asymptotic stability guarantees that lim� →∞ �⃗ = 0⃗ and lim� →∞ �⃗ = 0;⃗ then (�� , �� , �� ) → (�� , �� , �� ) and (�(�), �(�)) → (�0 , �0 ) if � → ∞. hus, any intruder system may synchronize the temporal evolution of the transmitter, deduct the values of the parameters and, in conclusion, break the cryptosystem. For that, the same synchronization signal employed in the receiver is enough. In conclusion, Lorenz-based cryptosystems may be broken in a pretty easy way due to the low degree of complexity of that dynamics. Moreover, Orúe et al. [2] could prove how the Lorenz system is easily deciphered. herefore, new more complex dynamics are necessary in order to build secure cryptosystems. 3. A New Hyperchaotic Lorenz-Based System In order to get complex hyperchaos, we have considered the Lorenz system (1) and we have added a nonlinear controller � to the irst dynamical equation, including a diferential nonlinear evolution law for � as a new fourth component. We have sought to strengthen the rate of expansion and to reduce the rate of contraction of this dissipative system as explained below. From this procedure, we obtain the following new system: �̇ = � (� − �) + 2�, �̇ = 5� + �� − 4��, 4 Complexity 42 3.5 3.4 40 3.3 38 3.2 36 3.1 34 3 2.9 32 2.8 30 28 28 2.7 30 32 34 36 38 40 42 2.6 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 3.4 3.5 Figure 2: First return map for (a) Lorenz system {� = 10, � = 178, and � = 8/3} and (b) proposed system {� = 20, � = 2, � = 3, and � = 10}. (a) (b) �̇ = �� − 3�, �̇ = −�� − � (� − �) , where �, �, and � are real positive parameters. First of all, it must be noted that the proposed system still preserves the symmetry about �-axis of the Lorenz dynamics. Nevertheless, as can be seen in Figure 2, this symmetry is not preserved in the irst return map. One of the main discoveries of Lorenz in the system he investigated was the similarity of the irst return map for chaotic trajectories to the tent map [21]. he tent map is the most basic paradigm of discrete chaos, being mainly characterized for its symmetry. As can be seen in Figure 2(b), in the proposed system (for chaotic trajectories), the irst return map preserves some aspects of the Lorenz dynamics, but the resulting structure is not symmetric, as the cusp is splitted. his fact shows a irst sign of the degree of complexity reached. On the other hand, the proposed system presents three diferent ixed points (16) that maintain the symmetry about �-axis, inherited from the Lorenz system: (15) �0 (0, 0, 0, 0) , 5+� 3 3 , 0) , �1 (√ (5 + �), √ (5 + �), 4 4 4 5+� 3 3 �2 (−√ (5 + �), −√ (5 + �), , 0) . 4 4 4 (16) As can be seen, �0 (the coordinate origin) does not depend on any parameter, and �1 and �2 depend only on the parameter �. hen, as parameter � is always positive, the coordinates of �0, �1, and �2 are always real and diferent. he three points, which lie on the hyperplane � = 0, are getting closer as parameter � decreases (without ever colliding). herefore, the system exhibits no bifurcations involving collisions or splits of the ixed points such as the pitchfork bifurcation. Other important requirements for the novel system are as follows. (i) We maintain the two nonlinear terms of the Lorenz system but we introduce one new divergence in the second equation of system (15) with the aid of the parameter �. (ii) he system has dissipative structure and therefore the four dynamics parameters cannot be ixed at any particular value independently (17). his is a necessary condition for dissipative chaos that has to be fulilled. he volume occupied by the system in the phase space should decrease as time goes on. In dynamical systems theory, this is veriied by checking that the rate of the occupied volume per volume unit is negative (i.e., the rate of volume contraction): 1 �� (�) ( ) < 0. � �� (17) Expressing the dynamics as a matrix, we obtain the following: ( �̇ �̇ � (� − �) + 2� 5� + �� − 4�� ), ) = �⃗ = ( �� − 3� �̇ �̇ −�� − � (� − �) 1 �� (�) ⃗ = ( ��̇ + ��̇ + ��̇ + ��̇ ) ( ) = div (�) � �� �� �� �� �� (18) = � − (� + 3 + �) . So, the dynamics parameters must verify the following condition: � > � − (� + 3) . (19) As can be seen, the described dynamics depends on four parameters. However, in order to analyze the new system, it Complexity 5 Eigenvalues �0 Eigenvalues �1, �2 Interpretation Saddle point �1+ , �2+ , �3− , �4− order 1 Saddle point �1+ , �2+ , �3− , �4− order 1 Saddle spiral �1+ , �2+ , �3− , �4− order 1 Hopf bifurcation Saddle spiral �1− , �2− , �3− , �4− order 1 Table 1: Local stability study. 0<�<3 Parameters Analysis 3 < � < 15 System evolution when varying the parameter � for � = 10 15 < � < 25 � = 25 � > 25 ∀� System evolution when varying the parameter � for � = 2 0 < � < 11 � = 11 System evolution when varying the parameter � for � = 20 ��+ , real positive eigenvalue; 11 < � < 12 ��− , � > 12 real negative eigenvalue; �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− �1+ , �2− , �3− , �4− ��+ , Saddle spiral �1+ , �2+ , �3− , �4− order 1 Hopf bifurcation Saddle spiral �1− , �2− , �3− , �4− order 1 Saddle point �1− , �2− , �3− , �4− order 1 complex eigenvalue with positive real part; is important to choose a parameter coniguration which not only maintains the system’s complexity but also simpliies its analysis. herefore, we are taking ixed values for some system parameters. We ix the degree of divergence which controls the parameter � to � = 3 and inally the parameter �, which controls the convergence, is set to � = 2. If parameters � and � are ixed, the parameter � obeys > −2. For � < −2, the divergence is positive and the system tends to be not steady and diverges. When � is near zero, the system tends to converge slowly and generates hyperchaotic structures. In particular, when � = 2 and � = 38, we obtain a high value of the maximum Lyapunov exponent of approximately � max = 12 and a Kaplan-Yorke dimension close to �KY = 4. In contrast, when � > 0, if � becomes larger, system (15) shrinks strongly and will tend to be steady towards one of its equilibrium points. hus, the parameter � can be varied independently without change in the general behavior, which is determined essentially by the parameter �. hen, by the introduction of � and, especially, � as control parameters, we allow exploring the regular, chaotic, and hyperchaotic topologies of the system. According to these arguments and in order to study the maximum number of diferent representative structures and topologies, three diferent analyses are going to be performed: (1) System evolution when varying the parameter � for � = 10 (2) System evolution when varying the parameter � for �=2 (3) System evolution when varying the parameter � for � = 20 For all the three analyses described above, the eigenvalues associated with the matrix of the linearized system around the �1+ , �2+ , �3− , �4− Saddle spiral order 1 ��− , Interpretation Saddle spiral order 2 Saddle spiral order 2 Saddle spiral order 2 Spiral node Saddle spiral order 2 Saddle spiral order 2 Spiral node Spiral node complex eigenvalue with negative real part. equilibrium points have no analytical expression depending on the dynamics parameters. hus, to obtain their values and study the stability of the ixed points, possible bifurcations, and so forth, we made a numerical calculation. We take advantage in considering the symmetry of the system and we know, in general, that the linearized systems around the points �1 and �2 have the same eigenvalues, so their stability and behavior are identical. Table 1 studies the behavior of the ixed points for the three situations identiied above. From the analysis of the results shown in Table 1, two conclusions about the global behavior of the system may be deduced: (1) Taking � = 10, the evolution of the system is analyzed when varying the parameter �. hen a Hopf bifurcation appears at � = 25. At this point, the exterior ixed points change their behavior from instable to asymptotically stable. For values lower than but close to � = 25, the system admits (as we see in Section 4) two coexisting limit cycles which tend to be smaller when getting closer to � = 25 until they collapse over the ixed points. herefore, at {� = 25, � = 2, � = 3, and � = 10} the system presents a supercritical Hopf bifurcation. (2) On the other hand, with � = 20, the evolution of the system is analyzed when varying the parameter �. hen another supercritical Hopf bifurcation appears at � = 11. he system behavior is identical to what was described previously for the value � = 25. Various methods to determine the supercritical character of Hopf bifurcations have been described in the research literature. In particular, a Hopf bifurcation is considered subcritical if the maximum Lyapunov exponent in the bifurcation point is positive and supercritical in the opposite case. 6 Complexity Table 2: Lyapunov spectrum at the bifurcation points. Case {� = 25, � = 2, � = 3, and � = 10} {� = 20, � = 2, � = 3, and � = 11} 3 (−0.0419 −0.0503 −5.9441 −5.9613) Lyapunov spectrum (−0.0282 −0.0341 −6.4436 −6.4914) 2.5 Z 2 1.5 1 4 3 2 2 1 X 0 2 −2 W 0 −4 −2 0 4 Y −4 −1 (a) −2 3 −3 3 2.5 −1 2 Z 2.5 −2 1.5 −3 1 Y Z 2 −4 Figure 3: Projection over the subspace YZW of the limit cycle of new system, working in regular region (� = 25 and � = 10). Numerical algorithms allow determining the Lyapunov spectrum of the two cases of interest. Results are shown in Table 2. As can be seen, in both cases, the maximum Lyapunov exponent is negative, so, in fact, bifurcations are supercritical. In the proposed system, the complexity of the trajectories tends to increase for small values of parameter �, since the system is increasingly less convergent, as we will show later. hus, hyperchaotic and more complex structures are located when analyzing the system evolution along the parameter � for � = 2. In this range of low values of the parameter �, strong exponential growing (accompanied by strong contractions, since the system is dissipative) and instabilities appear, so we will do a careful study about the obtained results. 4. Regular Structures Regular structures may be found in any of the three studies proposed. In particular, two diferent regular types of trajectories are generated by the proposed system. On one hand, common limit cycles are generated in the neighborhoods of the Hopf bifurcations identiied in the previous section. On the other hand, more complex cycles are generated for diferent values of the parameter � when considering the value � = 10. In the next subsections, both cases are studied. 4.1. Limit Cycles and Hopf Bifurcation. Limit cycles (Figure 3) are mainly located in two ranges (both in the neighborhoods of the Hopf bifurcations): � ∈ (22.5, 25) for � = 10 and � ∈ (10.5, 11) for � = 20. In both cases, two diferent limit cycles coexist at the same time (as can be seen in Figure 4), with each one developing around of one of the exterior ixed points �1 and �2. 1.5 1 4 2 X 4 0 −2 −2 −4 0 2 Y −4 Figure 4: Coexisting limit cycles in the proposed system. (a) {� = 24, � = 2, � = 3, and � = 10} and (b) {� = 20, � = 2, � = 3, and � = 10.55}. Initial conditions: (1, 1, 1, 1) and (−2, −2, −2, −2). (b) his coexistence is maintained while the system remains in the regular region and disappears when chaotic trajectories appear (as these trajectories are developed around both ixed points). his fact is clearly shown in Figure 5, where a bifurcation diagram for two diferent initial conditions is represented. As can be seen, in the regular region, both diagrams (one represented in red color and the other in green color) evolve in parallel and independently. Nevertheless, when chaotic structures appear, both diagrams overlap. Finally, the existence of a Hopf bifurcation must be proven. In Figure 6, a composition of the evolution of the limit cycles when varying the control parameters is provided. As can be seen, limit cycles reduce their size when getting closer to the bifurcation point, until they collapse over the ixed points. his behavior is, in fact, corresponding to a supercritical Hopf bifurcation. 4.2. Other Regular Attractors and the Corresponding Bifurcation Diagram. If the system low evolution when varying the parameter � for � = 2 is considered, more complex structures are found (see Figure 7). A rich structure of regular trajectories appears within certain windows and miniature windows within the larger windows of parameter �. In this igure (Figure 7), two bifurcation diagrams are shown, representing Complexity 7 5 X −2 2 25 Figure 5: Bifurcation diagram using � as control parameter for two diferent initial conditions. In red color (1, 1, 1, 1) and in green color (−2, −2, −2, −2). Parameters: {� = 2, � = 3, and � = 10}. Control parameter, a 3 2.5 Z 2 1.5 1 4 3 30 2 X 28 1 26 a Figure 6: he evolution of the limit cycles with the parameter � in the neighborhood of Hopf bifurcation. Initial conditions: (1, 1, 1, 1). Parameters: {� = 2, � = 3, and � = 10}. 0 X 24 8 8 7 7 6 6 5 5 4 4 X 3 3 2 2 1 1 0 0 −1 −1 −2 30 40 50 60 70 Control parameter, a 80 90 100 −2 38.6 38.8 39 39.2 39.4 39.6 Control parameter, a Figure 7: Bifurcation diagram using � as control parameter. Parameters: {� = 2, � = 3, and � = 2}. (a) (b) 39.8 40 8 Complexity 30 40 20 20 W 10 W 0 0 −10 −20 −20 −40 40 −30 5 20 20 0 Y 10 0 −20 Y 0 X −10 −40 −20 0 −5 −5 −10 10 5 X Figure 8: Complex limit cycle for (a) {� = 10 and � = 3} and (b) {� = 50 and � = 2}. Projection over the subspace XYW. (b) (a) Z 3 3 2.5 2.5 2 Z 2 1.5 1.5 1 4 1 −1 3 2 Y 1 0 1 4 3 2 −2 −1 −3 X −4 −4 (a) (b) 100 300 200 100 50 Z −2 −3 Y X 0 0 Z −50 −100 200 100 0 Y −100 −200 −100 −50 50 0 0 −100 −200 −300 500 100 X 0 Y −500 −100 −50 50 0 100 X Figure 9: Coexisting limit cycles for {� = 40 and � = 7}. Initial conditions: (a) (150, 2, 1, 45), (b) (0.1, 0.1, 0.1, 0.1), (c) (100, 50, 50, 100), and (d) (200, 300, 200, 350). (c) diferent regular and chaotic regions in various ranges of the parameter � for � = 2. Several and interesting bifurcations may be identiied in these diagrams. For example, at � = 70, a period doubling bifurcation is exhibited. Additionally, around � = 50, another period doubling cascade gets stopped and reversed. Moreover, in the range � ∈ (39.1, 39.7), a period doubling cascade, a chaotic window, and a period doubling again end at � = 39.7 in chaos. As can be seen, at � = 39.7, the chaos appears abruptly by means of a tangent bifurcation. At this bifurcation parameter value, an underlying saddle point and an underlying stable node collide and annihilate one another. As a result, an orbit that has periods of chaos interspersed with periods of regular oscillation is generated (d) [22]. he same ideas shown before allow identifying another tangent bifurcation at � = 43. In some of the shown areas, limit cycles appear (e.g., in the range � ∈ (70, 100); see Figure 8(a)). However, other interesting regular structures are also generated, as can be seen in Figure 8(b). On the other hand, although the couple {� = 40 and � = 7} does not belong to any of the proposed analyses, the observed structures deserve to be commented on. In fact, in this range, the complexity of the system is enormous, while up to four diferent hidden limit cycles [23] can be generated by the dynamics for {� = 40 and � = 7} modifying the initial conditions (see Figure 9). Some of these cycles are regular limit cycles, but others exhibit a great complexity. Complexity 9 3.5 Table 3: Lyapunov exponents study. Topology � = 20; � = 7 � = 32; � = 5 � = 40; � = 2 3 2.5 2 Z 1.5 1 0.5 0 −5 0 5 Figure 10: Lorenz-like chaotic attractor in the proposed system {� = 20 and � = 10}. X 6 5 4 3 2 X 1 0 −1 −2 −3 3 4 5 6 7 8 Control parameter, d 9 10 11 Figure 11: Bifurcation diagram using � as control parameter. Parameters: {� = 20, � = 2, and � = 3}. 5. Chaotic Structures In the proposed system, the weakly chaotic structures maintain certain similarity with the classic Lorenz attractor (see Figure 10). In this igure, a pretty similar attractor to Lorenz’s butterly is shown and, as it was seen in Figure 2(b), the irst return map associated with this topology is also similar to that obtained in the Lorenz dynamics (apart from the changes in the symmetry). his chaotic region is well appreciated in the diagram of Figure 11, where parameter � is ixed to � = 20 and the parameter � is variable. We can appreciate that in the region � ∈ (3.8, 10.4) the dynamics is chaotic, and when � > 10.4 it starts the regular regime that inally converges to �0 ixed point by means of a Hopf bifurcation. Other interesting structures in the phase space, together with the associated irst return map, can be seen in Figure 12. he topology generated for {� = 20 and � = 7} maintains a great similarity to Lorenz’s butterly, but in the irst return Lyapunov exponents (0.62, 0, −4.7, −4.91) (0.43, 0, −3.2, −4.1) (1.11, 0, −2.3, −2.7) Kaplan-Yorke dimension 2.12 2.13 2.47 map there is a double cusp more complex than the Lorenz return map for �-component. he beginning of the symmetry breaking can also be seen. Similar topologies may be found using the bifurcation diagram provided in Figure 11. For {� = 32 and � = 5}, the variations are really evident. he irst return map still contains some information about the Lorenz dynamics but has only a slight inluence. At this point, complexity has started growing up. As we have said previously and as can be seen in the bifurcation diagram of Figure 7, the most complex structures are located in the system evolution when varying the parameter � for � = 2. In these ranges, highly complex chaotic structures and hyperchaotic trajectories may be located. In particular, for the couple of parameters {� = 40 and � = 2} a high complexity chaotic solution is generated. As can be seen, the trajectory in the phase space is much more complex than that obtained for the Lorenz dynamics. Moreover, the irst return map is an almost-random cloud of points which does not maintain any symmetry or regular structure as in the previous cases. Table 3 provides a study about the Lyapunov exponents and the Kaplan-Yorke dimension of the structures shown in Figure 12, which proves the increase of the complexity as the simpler Lorenz dynamics disappears. In particular, Figure 12(c) shows a high complex chaotic structure, which presents a Kaplan-Yorke dimension, �KY = 2.47. In the range � ∈ (38, 40), various additional complex trajectories may be found, including hyperchaotic solutions coexisting with high complexity chaotic trajectories. However, in the range � ∈ (40, 43), the proposed system behaves chaotically and various strange attractors with a low similarity to Lorenz butterly are found (see Figure 13). As can be seen and as was said in Section 3, the increase of the complexity is due to the explicit control of rates of convergence and divergence by means of the system parameters. As a inal result and with varying the parameters � and �, the hyperchaotic behavior appears. Next section is dedicated to these trajectories. 6. Hyperchaotic and Other Complex Structures In this section, diferent complex chaotic structures are located and analyzed. As we will see below, in the proposed system, the strengthening of the expansion dimensions has motivated the appearance of highly complex chaotic trajectories, sometimes coexisting with hyperchaotic ones. Other complex structures (such as the chaotic-transient solutions) and unbounded solutions, when the parameter � is small (in absolute value), also appear. hese structures are reviewed in the following subsections. 10 Complexity 3.5 3.6 3 3.4 3.2 2.5 3 2 2.8 Z 1.5 2.6 1 2.4 0.5 2.2 0 −6 −4 0 −2 2 4 6 2 2 2.2 2.4 2.6 3 2.8 3.2 3.4 3.6 X (a) 3.5 3.6 3.4 3 3.2 2.5 3 2 2.8 1.5 2.6 Z 2.4 1 2.2 0.5 0 −6 2 −4 −2 0 2 4 6 1.8 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 X (b) 4 4 3.5 3.5 3 3 2.5 2.5 Z 2 2 1.5 1.5 1 1 0.5 0 −8 −6 −4 −2 0 2 4 6 8 0.5 0.5 1 1.5 2 2.5 3 3.5 4 X Figure 12: Chaotic attractors of the XZ components generated by the proposed system and irst return maps obtained from �-component. Initial conditions: (0.1, 0.1, 0.1, 0.1). Parameters: (a) {� = 20 and � = 7}, (b) {� = 32 and � = 5}, and (c) {� = 40 and � = 2}. (c) Complexity 11 4 3 Z 2 1 0 10 5 Y 0 5 0 −5 10 4 3.5 3 2.5 2 Z 1.5 1 0.5 0 5 0 Y X −5 −10 −5 −10 −5 5 0 10 X Figure 13: Chaotic tridimensional attractors generated by the proposed system. Initial conditions: (0.1, 0.1, 0.1, 0.1). Parameters: (a) {� = 40 and � = 2} and (b) {� = 43 and � = 2}. (a) (b) 300 Table 4: Lyapunov exponents study for the hyperchaotic structure. Exponent �1 �2 �3 �4 200 100 Z 0 −100 −200 −300 500 0 100 0 Y −500 −100 −200 200 X Figure 14: Projection over the subspace XYZ of the chaotictransient solution. 6.1. Chaotic-Transient Solutions. As we have commented in the previous subsection, hyperchaotic and unbounded solutions coexist in most of the cases. However, another type of solutions which may be confused with hyperchaotic structures can exist. In some cases, the proposed system generates trajectories which seem to be hyperchaotic at irst. However, as time elapses, the components enter in a coherent state which does not correspond with hyperchaos. his situation, for example, appears for the couple of parameters {� = 10 and � = 2}. Figure 14 shows a projection of the cited trajectory over the subspace XYZ. As can be seen, a cylindrical topology appears, indicating an almost-regular behavior in some of the components. At the irst moments, this structure generates a Kaplan-Yorke dimension around �KY = 3.4, but, as time elapses, the system components reach a coherent state and complexity goes down rapidly. In fact, �-component and Y-component (as well as �component and �-component) evolve together with regular oscillations from � > 30 dimensionless units (see Figure 15). In conclusion, these topologies present a chaotic structure in a transitory way, but, as time elapses, they become more regular and lose complexity. Figure 16 shows irst two return maps, where the chaotic behavior can be seen on one hand and, on the other hand, Value 12.37 0 −0.036 −14.71 Kaplan-Yorke dimension 3.83 the regular behavior can be seen. As can be seen, the regular regime is clearly distinguished by the cloud of points over � = � straight line which is especially dominant in the return map of �-component (Figure 16(b)). On the other hand, the behavior of the chaos is characterized, in both igures, from a dense cloud of points outside of the straight line. In addition to these arguments, in Figure 15, it can be seen that the regular trajectory evolution at about � > 40 is dominant, indicating that, at the end, this regular trajectory is the cause of the loss of complexity. 6.2. Hyperchaotic Solutions. Although, probably, diferent hyperchaotic structures could be generated by the proposed system using diferent values of the control parameters, in this section, we will focus the study on the behavior in the range � ∈ (38, 38.5) taking the value � = 2. he irst hyperchaotic structure may be found at � = 38, although in the range � ∈ (38, 38.5) similar result could be obtained. Figure 17 shows the four diferent threedimensional projections of the four-dimensional speciied attractor. A basic proof of the presence of hyperchaos is obtained by means of the Lyapunov spectrum. Table 4 shows its values. As can be seen, in this system, like the generalized Chen system [11], only one exponent is positive and shows a high value. However, the third exponent is really close to zero (although it is negative). Using the ordered Lyapunov spectrum (�1 > � 2 > � 3 > � 4 ) we can observe a great divergence in the direction of the irst exponent together with a low contraction in the direction of the third exponent. his is the reason for the great complexity of the hyperchaotic solution being analyzed. he Kaplan-Yorke dimension �KY clearly shows this high 12 Complexity X Y Z 100 0 −100 0 10 20 30 Time 40 50 60 0 10 20 30 Time 40 50 60 0 10 20 30 40 50 60 40 50 60 200 0 −200 100 0 −100 Time W 200 0 −200 0 10 20 30 Time Figure 15: Temporal evolution in the chaotic-transient solutions. 800 10000 600 8000 400 6000 4000 200 2000 0 0 −200 −2000 −400 −4000 −600 −6000 −800 −800 −600 −400 −200 10000 8000 6000 4000 2000 0 800 −2000 600 −4000 400 −6000 200 −8000 −8000 0 Figure 16: First return map for the chaotic-transient solution. (a) �-component; (b) �-component. (a) complexity. It presents a value greater than three and close to four (�KY = 3.83). Also, we can observe that � 1 and � 4 , with opposite signs, compete with each other. Baier and homsen [24] performed an analysis of the complexity of the chaotic attractors in four dimensions and established the following: (a) For � 1 > 0 and (� 1 +� 3 ) < 0, there are ordinary chaos with 2 < �KY < 3. (b) For (� 1 + � 2 ) > 0 and (� 1 + � 3 ) > 0, there appears more complex chaos (hyperchaos) with 3 < �KY < 4. Our case belongs to the more complex chaotic behavior indicated in (b). One of the problems when working with dynamical systems with a high level of complexity is the low instability depending on the control parameters in the region of hyperchaos and transient chaos. his high level of complexity appears for low values of � parameter (� ≤ 2) and with � < 37, (b) and it is very sensitive on integration errors. In fact, precise integration algorithms and Lyapunov spectrum calculation routines are necessary in order to obtain accurate data. In the case of the result shown in Table 4, both routines, the integration one (with a Runge-Kutta four-order method and a ifth-order error control) and the Lyapunov spectrum one, have been implemented with a relative error between 10−7 and 10−8 and with a very long time propagation, so the convergence is reached, as can be seen in Figure 18. To obtain the Lyapunov spectrum, we use Wolf et al.’s algorithm [25] with a temporal step size between 10−3 and 10−4 . Finally, since the four-dimensional attractors cannot be represented in only one graphic, Poincaré sections are a traditional way of proving the complexity level and the presence of chaos. Although the Lyapunov exponents and the previous results prove the hyperchaotic behavior of the analyzed structure, Poincaré sections in Figure 19 demonstrate the existence of the hyperchaos in the system. hese Complexity 13 ×104 4 4000 2 2000 −2000 −2 ×10 5 0 X Z 0 −4 1 −4000 4 0.5 0 Y −0.5 −1 −4000 −2000 2000 0 4000 ×104 2 0 W X −2 Z −4 (b) (a) ×10 4 0 −2 −4 4 ×104 2 4 ×10 2.5 4 2 3 1.5 2 1 1 0.5 0 Z W 0 −0.5 −1 −1 −2 −1.5 −3 −2 −4 −4000 −3000 −2000 −1000 0 1000 2000 3000 4000 −2.5 −4 −3 −2 −1 0 1 Z X 2 Figure 17: Diferent projections of the hyperchaotic attractor. Initial conditions: (−3203, −26967, 32052, −5434). Kaplan Yorke dimension = 3.83; 15 Lyapunov exponents 10 5 0 −5 −10 −15 −20 0 50 100 150 200 Time 4 ×104 (d) (c) 20 3 250 300 350 400 Figure 18: Evolution of the Lyapunov exponents values versus the integration time. sections generate a hyperplane deined by the normal vector ⇀ �� = (0, 0, 0, 1) and the point �(0, 0, 0, 2) (see Figure 19). As can be seen in Figure 17 and as is shown more clearly in Figure 20, the hyperchaotic structures surround the origin, creating a hollow cavity. In this empty space, at � = 38, a second chaotic attractor coexists with the hyperchaotic one previously described which remains for that value of parameter �. In particular, this new topology exhibits chaotic behavior, with a great complexity, and develops an attractor (see Figure 21) very similar to that shown in Figure 13(a), although with a higher Kaplan-Yorke dimension, with a value of �KY = 2.52. his dimension is high but does not exceed the value of three. In Table 5, the obtained Lyapunov spectrum for this chaotic attractor is shown. he coexistence of two complex structures in certain areas of the phase space reinforces the arguments about the complexity of this new system. To our knowledge, no other chaotic system with this property has been reported. 7. Discussion and Conclusions In this article, a new hyperchaotic four-dimensional Lorenzbased system, especially designed to improve the complexity of traditional Lorenz dynamics, is numerically analyzed. 14 Complexity 5 8 4 6 3 4 2 2 1 0 0 −1 −2 −2 −4 −3 −6 −4 −5 −5 0 5 −8 −6 −4 0 −2 2 4 6 (b) (a) 4000 3000 2000 1000 0 −1000 −2000 −3000 −4000 −6000 −4000 −2000 0 2000 4000 6000 Figure 19: Poincaré sections of the system for diferent � and � values. Red color indicates that the low crosses the cut plane in the sense of the normal vector. Blue color indicates that the low crosses in the opposite sense. (a) Weak chaos with initial conditions: (0.1, 0.1, 0.1, 0.1){� = 20 and � = 10}. (b) Chaos with initial conditions: (0.1, 0.1, 0.1, 0.1){� = 38 and � = 2}. (c) Hyperchaos with initial conditions: (−3203, −26967, 32052, −5434){� = 38.4 and � = 2}. (c) 4000 Table 5: Lyapunov exponents study. Exponent �1 �2 �3 �4 3000 2000 1000 X 0 Value 1.2092 0 −2.3517 −2.8622 Kaplan-Yorke dimension 2.52 −1000 −2000 −3000 −4000 −2.5 −2 −1.5 −1 −0.5 0 W 0.5 1 1.5 2 2.5 ×104 Figure 20: Bidimensional XW projection of the hyperchaotic attractor for {� = 38 and � = 2}. Initial conditions: (−3203, −26967, 32052, −5434). First, a proof about the weakness of the Lorenz system, due to the amount of redundant information present in its components, is provided. Regular, chaotic, and hyperchaotic structures are located and analyzed by varying two parameters of the system. Analyses about the symmetry, stability, and bifurcation are also provided. Two Hopf bifurcations and other typical bifurcations such as period doubling bifurcations and a tangent bifurcation were observed when � and � control parameters were varied. Moreover, the Lyapunov exponents proved the increase in the trajectories complexity, Complexity 15 4 30 20 3 10 Z 2 0 W −10 1 −20 0 10 −30 4 5 0 Y −5 −5 −10 5 0 3 10 10 2 Z X 5 1 (a) W 0 0 −5 Y (b) 30 30 20 20 10 10 0 W 0 −10 −10 −20 −20 −30 4 −30 10 3 2 Z 5 0 1 −5 0 −10 10 5 Y X 5 0 0 −5 −10 −5 10 X (d) (c) Figure 21: Tridimensional projections of the highly complex chaotic attractor coexisting with the hyperchaotic one. Initial conditions: (0.1, 0.1, 0.1, 0.1). doubling in some cases the values obtained for the Lorenz dynamics. For small positive values of �, the third Lyapunov exponent is negative but small. his leads to a very weak contraction in this direction. On the other hand, the irst exponent, which is responsible for the divergence direction, is very high. Both results lead to a Kaplan-Yorke dimension of �KY = 3.83. A signiicant characteristic of this system is the coexistence of chaos and hyperchaos attractors as described in this paper. In addition to the chaos and hyperchaos obtained, complex chaotic-transient solutions occur in this novel system. Future work is addressed to study the coherence of this high complex system and also its application to chaos control, synchronization, and secure communications. Competing Interests he authors declare that there are no competing interests regarding the publication of this paper. Acknowledgments One of the authors, Borja Bordel, has received funding from the Ministry of Education through the FPU Program (Grant no. FPU15/03977), the Ministry of Economy and Competitiveness through SEMOLA project (TEC2015-68284-R) and from the Autonomous Region of Madrid through MOSIAGIL-CM project (Grant P2013/ICE-3019, cofunded by EU Structural Funds FSE and FEDER). he authors are grateful for discussions with Professor Vicente Alcober. References [1] E. N. Lorenz, “Deterministic nonperiodic low,” Journal of the Atmospheric Sciences, vol. 20, no. 2, pp. 130–141, 1963. [2] A. B. Orúe, V. Fernández Mármol, G. Álvarez Marañón et al., “Determination of the parameters for a Lorenz system and application to break the security of two-channel chaotic cryptosystems,” Physics Letters A, vol. 372, no. 34, pp. 5588–5592, 2008. [3] A. Ali-Pacha, N. Hadj-Said, A. M’Hamed, and A. Belgoraf, “Lorenz’s attractor applied to the stream cipher (Ali-Pacha generator),” Chaos, Solitons and Fractals, vol. 33, no. 5, pp. 1762–1766, 2007. [4] K. M. Cuomo, A. V. Oppenheim, and S. H. Strogatz, “Synchronization of Lorenz-based chaotic circuits with applications to communications,” IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, vol. 40, no. 10, pp. 626–633, 1993. [5] E. A. Lee, M. Niknami, T. S. Nouidui, and M. Wetter, “Modeling and simulating cyber-physical systems using CyPhySim,” in Proceedings of the 12th International Conference on Embedded Software, pp. 115–124, IEEE Press, October 2015. [6] L. M. Pecora, L. Moniz, J. Nichols, and T. L. Carroll, “A uniied approach to attractor reconstruction,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 17, no. 1, Article ID 013110, 2016. 16 [7] S. Yu and W. K. S. Tang, “Tetrapterous butterly attractors in modiied Lorenz systems,” Chaos, Solitons & Fractals, vol. 41, no. 4, pp. 1740–1749, 2009. [8] A. S. Elwakil, S. Ozoguz, and M. P. Kennedy, “Creation of a complex butterly attractor using a novel Lorenz-type system,” IEEE Transactions on Circuits and Systems I: Fundamental heory and Applications, vol. 49, no. 4, pp. 527–530, 2002. [9] C.-X. Liu and L. Liu, “A new three-dimensional autonomous chaotic oscillation system,” Journal of Physics: Conference Series, vol. 96, no. 1, Article ID 012173, 2008. [10] S. Yu, W. K. S. Tang, J. Lu, and G. Chen, “Generation of � × �wing lorenz-like attractors from a modiied Shimizu-Morioka model,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 55, no. 11, pp. 1168–1172, 2008. [11] O. E. Rössler, “Chaotic oscillations: an example of hyperchaos,” in Nonlinear Oscillations in Biology, vol. 17, pp. 141–156, American Mathematical Society, Providence, RI, USA, 1979. [12] X. Liu, X. S. Shen, and H. Zhang, “Multi-scroll chaotic and hyperchaotic attractors generated from Chen system,” International Journal of Bifurcation and Chaos, vol. 22, no. 2, Article ID 1250033, 15 pages, 2012. [13] C. L. Koliopanos, I. M. Kyprianidis, I. N. Stouboulos, A. N. Anagnostopoulos, and L. Magafas, “Chaotic behaviour of a fourthorder autonomous electric circuit,” Chaos, Solitons and Fractals, vol. 16, no. 2, pp. 173–182, 2003. [14] S. Cang, G. Qi, and Z. Chen, “A four-wing hyper-chaotic attractor and transient chaos generated from a new 4-D quadratic autonomous system,” Nonlinear Dynamics, vol. 59, no. 3, pp. 515–527, 2010. [15] C.-C. Yang, “Adaptive synchronization of Lü hyperchaotic system with uncertain parameters based on single-input controller,” Nonlinear Dynamics, vol. 63, no. 3, pp. 447–454, 2011. [16] C. Han, S. Yu, and G. Wang, “A sinusoidally driven lorenz system and circuit implementation,” Mathematical Problems in Engineering, vol. 2015, Article ID 706902, 11 pages, 2015. [17] K. Sun, X. Liu, C. Zhu, and J. C. Sprott, “Hyperchaos and hyperchaos control of the sinusoidally forced simpliied Lorenz system,” Nonlinear Dynamics, vol. 69, no. 3, pp. 1383–1391, 2012. [18] X. Wang and M. Wang, “A hyperchaos generated from Lorenz system,” Physica A. Statistical Mechanics and Its Applications, vol. 387, no. 14, pp. 3751–3758, 2008. [19] S. Boccaletti, J. Kurths, G. Osipov, D. L. Valladares, and C. S. Zhou, “he synchronization of chaotic systems,” Physics Report, vol. 366, no. 1-2, pp. 1–101, 2002. [20] L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990. [21] T. Yoshida, H. Mori, and H. Shigematsu, “Analytic study of chaos of the tent map: band structures, power spectra, and critical behaviors,” Journal of Statistical Physics, vol. 31, no. 2, pp. 279–308, 1983. [22] Y. Pomeau and P. Manneville, “Intermittent transition to turbulence in dissipative dynamical systems,” Communications in Mathematical Physics, vol. 74, no. 2, pp. 189–197, 1980. [23] C. Li and J. C. Sprott, “Coexisting hidden attractors in a 4-D simpliied lorenz system,” International Journal of Bifurcation and Chaos, vol. 24, no. 3, Article ID 1450034, 2014. [24] G. Baier and J. S. homsen, “Prototypes of attractors in four dimensions,” Physical Review E, vol. 48, no. 6, pp. R4172–R4174, 1993. [25] A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D: Nonlinear Phenomena, vol. 16, no. 3, pp. 285–317, 1985. Complexity Advances in Advances in Operations Research Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Applied Mathematics Algebra Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Probability and Statistics Volume 2014 The Scientiic World Journal Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 International Journal of Differential Equations Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014 Submit your manuscripts at https://www.hindawi.com International Journal of Advances in Combinatorics Hindawi Publishing Corporation http://www.hindawi.com Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of Journal of Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 International Journal of Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Discrete Mathematics Journal of Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Discrete Dynamics in Nature and Society Journal of Function Spaces Hindawi Publishing Corporation http://www.hindawi.com Abstract and Applied Analysis Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Journal of International Journal of Stochastic Analysis Optimization Hindawi Publishing Corporation http://www.hindawi.com Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 Volume 2014