Open access peer-reviewed chapter

Parrondian Games in Discrete Dynamic Systems

Written By

Steve A. Mendoza and Enrique Peacock-López

Submitted: 31 May 2018 Reviewed: 13 September 2018 Published: 19 November 2018

DOI: 10.5772/intechopen.81499

From the Edited Volume

Fractal Analysis

Edited by Sid-Ali Ouadfeul

Chapter metrics overview

1,050 Chapter Downloads

View Full Metrics


An interesting problem in nonlinear dynamics is the stabilization of chaotic trajectories, assuming that such chaotic behavior is undesirable. The method described in this chapter is based on the Parrondo’s paradox, where two losing games can be alternated, yielding a winning game. The idea of alternating parameter values has been used in chemical systems, but for these systems, the undesirable behavior is not chaotic. In contrast, ecological relevant map in one and two dimensions, most of the time, can sustain chaotic trajectories, which we consider as undesirable behaviors. Therefore, we analyze several of such ecological relevant maps by constructing bifurcation diagrams and finding intervals in parameter space that satisfy the conditions to yield a desirable behavior by alternating two undesirable behaviors. The relevance of the work relies on the apparent generality of method that establishes a dynamic pattern of behavior that allows us to state a simple conjecture for two-dimensional maps. Our results are applicable to models of seasonality for 2-D ecological maps, and it can also be used as a stabilization method to control chaotic dynamics.


  • chaos control
  • Parrondo’s paradox
  • switched dynamic systems
  • ecological maps
  • seasonality

1. Introduction

In population dynamics, discrete dynamic systems have been used to model the dynamics of ecological systems. One of the first maps used in ecology that suggested to study the new, X n + 1 , and the old, X n , non-overlapping populations is the logistic map. [1] Although a simple one-dimensional (1-D) map, the logistic map shows complex dynamics including chaos. Furthermore, the analyses of the logistic map gave us a better understanding of the properties of chaotic dynamics [2, 3, 4, 5, 6, 7].

In the case of 1-D discrete dynamics, for the last 18 years, alternate dynamics strategies have been the center of attention due to the so-called Parrondo paradox [8, 9, 10], where two losing games can be combined to yield a winning game. Furthermore, the idea that “lose + lose = win” has been extended to “chaos+chaos = periodic” in one-dimensional maps [11]. Just recently and for the first time, we were able to find the Parrondo dynamics in two 2-D maps [12]. In the contest of seasonality, we consider the alternation of undesirable dynamical behaviors yield a desirable behavior [13, 14]. So in the context of population dynamics we have considered cases where “undesirable + undesirable = desirable” dynamical behaviors occur as a result of a simple alternation of parameters [15, 16, 17, 18, 19, 20, 21].

In our present discussion, we extend our seasonality modeling strategy to several two-dimensional ecologically relevant maps and find that the “undesirable + undesirable = desirable”, the “chaos + chaos = periodic”, as well as, the “periodic + periodic = chaos” behaviors are not unique to 1-D maps. In Section 2, we consider a delayed logistic map, and in Section 3, we analyze a Lotka-Volterra map. In Section 4, we study a modified 2-D Ricker map, and in Section 5, we analyze the Beddington map. In Section 6, we discuss a modified Lotka-Volterra map, which includes a logistic prey growth. We conclude in Section 7 with a discussion and a summary of our results.


2. Delayed logistic equation

In our analysis of two-dimensional maps, we begin with the extended logistic map that incorporates a delay in population growth, defined by the following relation:

X n + 1 = Y n E1
Y n + 1 = C Y n 1 X n E2

where C is our bifurcation parameter. For the Lagged Logistic Equation, we consider C values from 0 to 2.27 in the original map, although with alternation, we can obtain a bifurcation diagram showing larger C values. Figure 1, shows the regular bifurcation map of the lagged logistic model. For all of the maps we study, both the X and Y graphs for a given show the same dynamics; for instance, parameters associated with chaotic dynamics in the X map are also associated with chaotic dynamics in the Y map; since we focus on a map’s dynamics, we only show the X function map.

Figure 1.

Lagged logistic map model, Eqs. (1,2), with C = 0 to C = 2.27 and the region from C = 1.90 to C = 2.27.

From Figure 1, we define our parameter value regions associated with complex or non-complex dynamics. The map on the left of the figure shows the whole bifurcation map; while the right magnifies the complex region. On the right hand figure, which is the magnified map, we can clearly see some periodic windows, but we pick parameter values associated with complex dynamics.

Next, we switch, or alternate, the parameter values between even and odd iterations through the following relation:

X n + 1 = f n X n Y n = Y n if n even f n X n Y n = Y n if n odd E3
Y n + 1 = g n X n Y n = C e Y n 1 X n if n even g n X n Y n = C o Y n 1 X n if n odd E4

The equation above describes our switching strategy in which we pick one parameter for every odd iteration, which we name Co, and use the even parameter, Ce, as our bifurcation parameter, for every even iteration. For the first type of behavior, we pick one parameter associated with complex dynamics as our Co value and switch it with our even parameter, Ce, in areas associated with chaotic dynamics. In our case, we see chaotic dynamics for C values greater than 2.0 when we construct the bifurcation diagram for Eqs. (1) and (2). In the resulting alternated, or switched, bifurcation diagram, represented by Eq. (3) for the current section, we look for regions that have periodic oscillations that are normally associated with chaos, thus resulting in the case “chaos + chaos = order.” Hence, for every switching map that we study, in the figures we also show the unswitched map for the same Ce parameter space. We make a note that the alternation of parameter values as defined by our switching strategy may result in an extension of C parameters yielding oscillations; that is we may see oscillations for C values greater than 2.27 as in the case for lagged logistic map. To compare our maps with Eqs. (1) and (2), we only study bifurcation maps using Eqs. (3) and (4) from C = 0 to 2.27.

For our analysis, we pick a Co as in the “chaos + chaos = order” case; however, we may choose a parameter within the periodic windows in the chaotic region, and when using the switched map, we focus on the Ce values that are less than the onset of chaos, which for the case of the lagged logistic pap is C = 2. We use the analysis discussed above for all of the cases in this paper.

For our first example of “chaos + chaos = periodic”, we consider the parameter value, C o = 2.10 and Eqs. (3) and (4). In our bifurcation diagram for Eqs. (3) and (4), in Figure 2, we consider C e greater than 2 and look for C e values that give us periodic oscillations. Figure 2 shows two maps at once, the left hand showing Eqs. (1) and (2), and the right hand graph shows Eqs. (1) and (2) with Co = 2.1. In this case, from Figure 2, we can see one region of periodicity from Ce = 2.26 to 2.27.

Figure 2.

Bifurcation diagram for Eqs. (1) and (2) and Eqs. (3) and (4), using Co 2.1.

Another combination of parameters yielding “chaos+chaos = periodic” uses Co = 2.15 and Eqs. (3) and (4), where Figure 3 shows a range of Ce values for which “chaos + chaos = periodic” holds, from Ce = 2.36 to 2.38. The same figure also shows other values for which the “chaos + chaos = periodic” relation holds, but these bands are not as prominent as the one we focus on. Through out the paper, we make a point that different Ce values give widely different behaviors and these differences in dynamic behaviors reveals the differences in the Ce values that give us desirable behaviors. For the rest of the paper, the approximate ranges of Ce will be given for one window that satisfies the “chaos + chaos = order” or “periodicity + periodicity = chaos”, since there are sometimes a variety of parameters meeting the relevant criteria for switching.

Figure 3.

Bifurcation diagram for Eqs. (1) and (2), and Eqs. (3) and (4), using C o = 2.15 .

We complete our analysis of delayed logistic map with one case in which “periodic + periodic = chaos”. As mentioned beforehand, we pick our value associated with periodic trajectories from the area associated with chaotic trajectories, and focus on Ce values less than the chaotic region for comparison. In particular, we choose C o = 2.19 as our periodic parameter for Eqs. (3) and (4); Figure 4 shows the corresponding bifurcation map, and, from the figure, we see one prominent example of “periodic+periodic = chaos” for Ce = 1.85 to 2.00.

Figure 4.

Bifurcation diagram for Eqs. (1) and (2), and Eqs. (3) and (4), using Co 2.19.


3. Lotka-Volterra model

We begin our next section by discussing a discretized form of the Lotka-Volterra model. The Lotka-Volterra map describes predator prey interactions, assuming that the prey has a relatively high initial population, and that the predator’s growth rate is directly proportional to the prey’s growth rate.

The model follows a relation defined by the map below

X n + 1 = 1 + r X n rX n 2 C X n Y n E5
Y n + 1 = C X n Y n E6

In Figure 5, showing Eqs. (5) and (6), we look at the unswitched map, defined by showing the ranges of periodic and aperiodic behavior. As in the previous section, we use the unswitched bifurcation map as a comparison to the switched map when using certain parameters. For this section, we focus on the interval C = 0 to 2.8, and set r = 2 for this map and the rest of the maps that have an r parameter.

X n + 1 = f n X n Y n = r + 1 X n r X n 2 C e X n Y n if n even f n X n Y n = r + 1 X n r X n 2 C o X n Y n if n odd E7
Y n + 1 = g n X n Y n = C o X n Y n if n odd g n X n Y n = C e X n Y n if n even E8

Figure 5.

Bifurcation diagram for Eqs. (5) and (6), showing the interval studied, as well as a close up of the chaotic region.

As before, we pick a Co value associated with a chaotic trajectory and alternate with Ce, using Eq. (7), which we use as the bifurcation parameter, illustrated in Figure 6. For this figure, we use Co = 2.1, and we can easily find conditions in which “chaos + chaos = order.” In particular, we see this phenomena for parameter values of Ce = 2.33–2.40. Figure 7, shows another example of “chaos + chaos = order” using Eqs. (7) and (8) with a Co value of 2.22, and in the corresponding bifurcation diagram for roughly Ce = 2.58–2.65.

Figure 6.

Bifurcation diagram for Eqs. (5) and (6) and Eqs. (7) and (8) with the Co value 2.1.

Figure 7.

Bifurcation diagram for Eqs. (5) and (6) and Eqs. (7) and (8) with the Co value 2.22.

We conclude the present section with an example of “periodicity + periodicity = chaos”, using Eqs. (7) and (8) and Ce = 2.44. In this case, in Figure 8, we see some chaotic behavior for values of Ce in the interval 1.75,2.00 , a region that is periodic when using Eqs. (5) and (6).

Figure 8.

Bifurcation diagram for Eqs. (5) and (6) and Eqs. (7) and (8) with the Co value 2.44.


4. Modified 2-D Ricker map

Another interesting map includes an exponential term, describing the prey growth, with a simple predator–prey interaction term. The map is determined by the following equations:

X n + 1 = X n Exp [ r 1 X n Y n E9
Y n + 1 = C X n Y n E10

which is in essence modified and extended to 2-D Ricker-like map [22]. The corresponding switched map is defined below:

X n + 1 = f n X n Y n = X n Exp [ r 1 X n Y n if n even f n X n Y n = X n Exp [ r 1 X n Y n if n odd E11
Y n + 1 = g n Y n = C e X n Y n if n even g n Y n = C o X n Y n if n odd E12

Figure 9, showing Eqs. (9) and (10), considers the range of C values we focus on, from C = 0 to 2.8. We want to remark however, that this map also shows some interesting behavior beyond the interval of study, but we choose this interval to get a close up of the intervals of periodicity, since this interval is where we find our relevant behavior. For the X function we study, at higher values, the function stays at unity for values of C = 28 and higher, while the Y function stays at extinction, or Y = 0.

Figure 9.

Bifurcation map for Eqs. (9) and (10) showing the interval studied, as well as a close up of the chaotic region.

As in previous cases, we start with finding parameter values satisfying the “chaos + chaos = order” relation. To begin, we use Eqs. (11) and (12) with Co = 2.10, associated with aperiodic dynamics. Figure 10 zooms into the region for which “chaos + chaos = order” holds. From this diagram, we see a narrow region of periodicity from C = 2.775 to 2.790.

Figure 10.

Bifurcation diagram for Eqs. (9) and (10) and Eqs. (11) and (12) with the Co value 2.10.

We then use Eqs. (11) and (12) with Co = 2.26, for which Figure 11 hones in on the relevant Ce parameter values. The interval of Ce values is significantly wider in this case than the previous one, since we find “chaos + chaos = order” for 2.74–2.80.

Figure 11.

Bifurcation diagram for Eqs. (9) and (10) and Eqs. (11) and (12) with the Co value 2.26.

We finish this section by introducing one case in which “periodic + periodic = chaos”. We pick the periodic parameter Co = 2.333. In some maps, it is harder to find periodic windows, although they could usually be found sometimes but an extra significant figure is necessary such as in this case. We find chaos in this map from Ce = 1.71 to 2.00, as shown in Figure 12, periodic values when using Eqs. (11,12).

Figure 12.

Bifurcation diagram for Eqs. (9) and (10) and Eqs. (11) and (12) with the C o value 2.333.


5. Beddington model

Our next map is the Beddington 2-D map defined by the following equations:

X n + 1 = X n Exp ( r 1 X n Y n ) E13
Y n + 1 = C X n 1 Exp Y n E14

along with the corresponding alternation equation.

X n + 1 = f n X n Y n = X n Exp r 1 X n Y n if n even f n X n Y n = X n Exp r 1 X n Y n if n odd E15
Y n + 1 = g n X n Y n = C e X n 1 Exp Y n if n even g n X n Y n = C o X n 1 Exp Y n if n odd E16

Figure 13, showing Eqs. (13) and (14), shows the parameter range we use to analyze the map. We pick points between 0 and 14, and show the corresponding bifurcation diagrams within that range. We pick 14 as our maximum value because above that parameter, there are only steady state solutions.

Figure 13.

Bifurcation diagram for Eqs. (13) and (14) showing the interval studied, as well as a close up of the chaotic region.

We start with describing our first chaotic value, Co = 10, for Eqs. (15) and (16). Figure 14 shows the corresponding bifurcation diagram, and we see a relatively wide range of Ce values for which we have “chaos + chaos = order”. We find this behavior for most points of Ce between 4.54 and 4.7.

Figure 14.

Bifurcation diagram for Eqs. (13) and (14) and Eqs. (15) and (16) with the Co value 10.

We then use Eqs. (15) and (16), with Co = 4.0, and here we also see a relatively wide range of parameters in which we find that “chaos + chaos = periodicity”. Specifically, we see that alternating with Ce = 10.7–10.88 gives us the desired behavior, shown in Figure 15. Our last figure pertaining to this map, Figure 16, shows the “periodic + periodic = chaos” behavior, for Co = 6.0. Figure 16 shows the area of the map that is normally periodic, and shows characteristic chaotic behavior from Ce = 1.7–3.0, although this particular map shows some periodic windows than the other “periodic + periodic = chaos” maps.

Figure 15.

Bifurcation diagram for Eqs. (13) and (14) and Eqs. (15) and (16) with the Co value 4.0.

Figure 16.

Bifurcation diagram for Eqs. (13) and (14) and Eqs. (15) and (16) with the Co value 6.0.


6. Modified Lotka-Volterra map

Our last 2-D map considers a logistic growth, and an interaction term, and only a predation term for the predator. The dynamics of this map is considerably different than the previous two maps,

X n + 1 = 1 + r X n r X n 2 CX n Y n X n + h E17
Y n + 1 = CX n Y n X n + h E18

As before, the switched map is shown below.

X n + 1 = f n X n = 1 + r X n rX n 2 C e X n Y n X n + h if n even f n X n = X n r + 1 r X n 2 C o X n Y n X n + h if n odd E19
Y n + 1 = g n Y n = C o X n Y n X n + h if n odd g n Y n = C e X n Y n X n + h if n even E20

Aside from the r parameter, this map also has the h parameter, which we set equal to unity. Unlike the previous two maps we study that have relevant behaviors past C = 10, the max value of the unswitched map is C = 3.85, but chaos is only present above C = 3.0, as shown in Figure 17.

Figure 17.

Bifurcation diagram for Eqs. (17) and (18) showing the interval studied, as well as a close up of the chaotic region..

Our first chaotic point is Co = 3.3, and the corresponding bifurcation diagram is shown in Figure 18. There is a somewhat small region of periodicity from Ce = 3.704 to 3.724.

Figure 18.

Bifurcation diagram for Eqs. (17) and (18) and Eqs. (19) and (20) with the C o value 3.3.

The second to last figure, Figure 19 shows our final odd switching parameter, Co = 3.1 and the corresponding bifurcation diagram, which shows a similar range of periodic parameter values, specifically, Ce = 3.70–3.72.

Figure 19.

Bifurcation diagram for Eqs. (17) and (18) and Eqs. (19) and (20) with the C o value 3.1.

Our last figure, Figure 20, shows an example of “periodic + periodic = chaos”, where we switch with Co = 3.57 and see chaos for most values between Ce = 2.5 and 3.0.

Figure 20.

Bifurcation diagram for Eqs. (17) and (18) and Eqs. (19) and (20) with the C o value 3.57.


7. Discussion

In previous sections, we have analyzed five relevant ecological 2-D maps, setting a pattern of dynamic behavior similar to the well studied “chaos + chaos = periodic” in switched 1-D maps. Therefore, with the results discussed in this chapter, we can extend the 1-D maps conjecture to 2-D maps. The conjecture asserts that given a map with chaotic dynamics, we can find two parameters associated to chaotic trajectories that, when alternated yield a periodic trajectory. In general, we can consider these kinds of maps as nonautonomous maps because one of the parameters is a function of the iterations. In most case, we pick a parameter value for the even iterations and a different parameter for the odd iterations. But the connection with the Parrondo’s paradox is associated with the kind of alternating parameters, which in the conjecture are parameter associated with chaotic, or, in general, complex trajectories.

The case of “chaos + chaos = periodic” was presented for the first time by Almeida et al. [16] for simple 1-D maps, and just recently for 2-D maps by Mendoza et al. [12]. The implication of the so-called Parrondo’s dynamics has been used to model seasonality, but with the observation that, under the Parrondo dynamics, the case of “periodic + periodic = chaos” is also possible [15]. As generalization we have consider cases of “undesirable + undesirable = desirable” dynamics behaviors to analyze simple models of seasonality [23, 24, 25], which include migration or immigration [13, 14].

In the present analysis, we emphasize the use of bifurcation diagrams to find intervals of values in parameter space that could satisfy the “undesirable + undesirable = desirable” or “periodic + periodic = chaos” dynamics. Although we are interested in modeling ecological systems and in particular the effect of seasonality, one could use our results to look at the switched maps as a way to control chaotic dynamics. In particular an extension to continuous dynamic systems may be relevant or applicable to chemical and mechanical systems [26].

In summary, our approach of building bifurcation diagrams readily yield intervals of parameter values that can show the so-called Parrondian dynamics for 1-D and 2-D maps. We have concentrated on ecological relevant maps, but the approach applies to any kind of maps. In particular, we can easily find parameters that show desirable dynamics in switched maps, controlling complex or undesirable dynamics, with the by product that we can also avoid the alternation of desirable dynamics that could yield undesirable dynamics in switched maps. Finally, we believed that we have stablished a pattern of dynamic behavior that supports the conjecture described in previous paragraphs.


8. Conclusions

In previous sections, we have stablished a pattern of dynamic behavior for 2-D maps, which have been used to model ecological systems. The dynamic pattern allows to state that for any 2-D maps that shows chaotic dynamics for a set of parameters, we can always find two of such parameters that, when alternate, yield a periodic trajectory. This conjecture is an extension of the so-called Parrondo’s paradox, in the sense that two undesirable dynamics can be alternate to yield a desirable dynamics. In other words, we can always find a region in parameter space, where we can select a pair of such parameters. Therefore, we the developed methodology can be use, in general, as a chaos control approach, and, in particular, we can use it to model, in the case of ecological maps, seasonality. Although we interested in ecological relevant 2-D maps, we believed that our conjecture can be extended to other type of 1-D and 2-D maps. Finally, we consider that the major application of the methodology is in controlling chaotic dynamics.



The authors will like acknowledge the financial support of the National Science Foundation (CHE-0911380), and the Bronfman Science Center. One of us (EPL) would like to thank Professor Allen Rodgers, and the Chemistry Department of Cape Town University for their hospitality during my sabbatical leave.


  1. 1. May RM. Simple mathematical models with very complicated dynamics. Nature. 1976;261:459-467
  2. 2. May RM. Stability and Complexity in Model Ecosystems. Princeton: Princeton University Press; 1974
  3. 3. Kot M. Elements of Mathematical Ecology. Cambridge: Cambridge University press; 2001
  4. 4. Turchin P. Complex Population Dynamics. Princeton: Princeton University Press; 2003
  5. 5. Mangel M. The Theortical Biologist’s Toolbox. Cambridge: Cambridge University press; 2006
  6. 6. Allen LJS. An Introduction to Mathematical Biology. Upper Saddle River, NJ: Pearson Prentice Hall; 2007
  7. 7. May RM, McLean AR. Theoretical Ecology: Principles and Applications. Oxford: Oxford University Press; 2007
  8. 8. Harmer GP, Abbott D. Losing strategies can win by Parrondo’s paradox. Nature. 1999;402:864
  9. 9. Harmer GP, Abbott D. Parrondo’s paradox. Statistical Science. 1999;14:14
  10. 10. Harmer GP, Abbott D, Taylor PG. The paradox of Parrondo’s games. Proceedings of the Royal Society. 2000;456:247-259
  11. 11. Cánovas JS, Linero A, Peralta-Salas D. Dynamic Parrondo’s paradox. Physica D. 2006;218:177-184
  12. 12. Mendoza SA, Matt EW, Guimaraes-Blandon DR, Peacock-López E. Parrondo’s paradox or chaos control in discrete two-dimensional dynamic systems. Chaos, Solitons & Fractals. 2018;106:86-93
  13. 13. Peacock-López E. Seasonality as a Parondian game. Physics Letters A. 2011;375:3124-3129
  14. 14. Silva E, Peacock-López E. Seasonality and the logistic map. Chaos, Solitons & Fractals. 2017;95:152-156
  15. 15. Percus OE, Percus JK. Can two wrongs make a right? Coin-tossing games and Parrondo’s paradox. Mathematical Intelligencer. 2002;24:68-72
  16. 16. Almeida J, Peralta-Salas D, Romera M. Can two chaotic systems give rise to order? Physica D. 2005;200:124-132
  17. 17. Behrends E. Stochastic dynamics and Parrondo’s paradox. Physica D. 2008;237:198-206
  18. 18. Boyarsky A, Góra P, Aslam MS. Randomly chosen chaotic maps can give rise to nearly ordered behavior. Physica D. 2005;210:284-294
  19. 19. Amengual P, Meurs P, Cleuren B, Toral R. Reversal of chance in paradoxical games. Physica A. 2006;371:641-648
  20. 20. Levinohn EA, Mendoza SA, Peacock-López E. Switching induced complex dynamics in an extended logistic map. Chaos, Solitons & Fractals. 2012;45:426-432
  21. 21. Maier MPS, Peacock-López E. Switching induced oscillations in the logistic map. Physics Letters A. 2010;374:1028-1032
  22. 22. Ricker WE. Stock and recruitment. Journal of the Fisheries Research Board of Canada. 1954;11:559-623
  23. 23. Blasius B, Kurths J, Stone L, editors. Complex Population Dynamics. Singapore: World Scientific; 2007
  24. 24. Allen LJS. An Introduction to Mathematical Biology. Upper Saddle River: Prentice Hall; 2007
  25. 25. Kot M, Schaffer WM. The effects of seasonality on discrete models of population growth. Theoretical Population Biology. 1984;26:340-360
  26. 26. Kohar V, Ji P, Choudhary A, Sinha S, Kurths J. Synchronization in time-varying networks. Physical Review E. 2014;90:022812

Written By

Steve A. Mendoza and Enrique Peacock-López

Submitted: 31 May 2018 Reviewed: 13 September 2018 Published: 19 November 2018