A Study of the Dynamics of a new Piecewise Smooth Map
Dhrubajyoti Biswas, Soumyajit Seth, Mita Bor

TL;DR
This paper investigates a novel 1D piecewise map combining tent and logistic maps, analyzing its complex dynamics through Lyapunov spectra, bifurcation diagrams, and numerical simulations, including a Simulink implementation.
Contribution
It introduces a new piecewise smooth map, explores its dynamics analytically and numerically, and demonstrates its implementation in Simulink.
Findings
The map exhibits both smooth and non-smooth behaviors.
Complex bifurcation structures are observed.
Unique dynamical phenomena are numerically demonstrated.
Abstract
In this article, we have studied a 1D map, which is formed by combining the two well-known maps i.e. the tent and the logistic maps in the unit interval i.e. [0, 1]. The proposed map can behave as the piecewise smooth or non-smooth maps (depending on the behaviour of the map just before and after the border) and then the dynamics of the map has been studied using analytical tools and numerical simulations. Characterization has been done by primarily studying the Lyapunov spectra and the corresponding bifurcation diagrams. Some peculiar dynamics of this map have been shown numerically. Finally, a Simulink implementation of the proposed map has been demonstrated.
| Remarks | ||
| Line of fixed points for | ||
| Line of fixed points for | ||
| and the characteristic fixed point of the logistic map. | ||
| Only one fixed point at | ||
| Fixed point at and the | ||
| characteristic fixed point of logistic map | ||
| Two fixed points of characteristics tent | ||
| and logistic map along with (3 fixed points overall) | ||
| Characteristic fixed point of logistic map | ||
| along with | ||
| Characteristic fixed point of tent map | ||
| along with |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
[1]Dhrubajyoti Biswas
D. Biswas, S. Seth, M. Bor
A Study of the Dynamics of a new Piecewise Smooth Map
Soumyajit Seth
Nonlinear Dynamics Laboratory, Department of Physical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal-741246. Email: [email protected]
Mita Bor
Department of Physics, St. Xavier's College (Autonomous), Kolkata-700016. Email: [email protected]
Abstract
In this article, we have studied a D map, which is formed by combining the two well-known maps i.e. the tent and the logistic maps in the unit interval i.e. . The proposed map can behave as the piecewise smooth or non-smooth maps (depending on the behaviour of the map just before and after the border) and then the dynamics of the map has been studied using analytical tools and numerical simulations. Characterization has been done by primarily studying the Lyapunov spectra and the corresponding bifurcation diagrams. Some peculiar dynamics of this map have been shown numerically. Finally, a Simulink implementation of the proposed map has been demonstrated.
1 The Mathematical Setup
The well known logistic and the tent maps (traditionally defined on the unit interval ) in the study of the chaotic dynamics [13] are given by
[TABLE]
and
[TABLE]
respectively. The logisitc map is a quadratic continous map on the interval , whereas the tent map is a piecewise linear map on the same interval. The representations $L_{r}$'' and '' have been used to represent them throughout the article. The control parameters ( for the logistic map and for the tent map) determine the dynamics exhibited by both of these. The Tent Map gives a fixed point in the parameter interval . As, the parameter is increased more, the fixed point looses it's stability and a bifurcation happens where a two piece chaotic orbit is born which gradually merges to give rise to robust chaos' [[5](#bib.bib5), [12](#bib.bib12)] upto the parameter value $2$. As in the parameter range $[1,2]$, only the chaotic orbit exists without any periodic attractor or co-existing attractors, the chaotic orbit is robust' in the mentioned parameter range which is shown in Figure. (LABEL:fig17a). For the case of the logistic map, it gives the period doubling bifurcation which gradually goes to chaos [16]. The most interesting feature is that the periodic windows exist in-between chaotic orbits in the bifurcation parameter range, which is shown in Figure. (LABEL:fig17b). These maps have been widely studied along with their applications in population modelling, stock market profiling, encryption-decryption systems etc. in various literature over the year [15, 10, 7, 9].
This article deals with a map born out of an amalgamation of these two well known maps, aptly named as the Mixed Map'' (MM), denoted by M'', and is defined as follows:
[TABLE]
where, is the nth iteration. It is clear from the above expression that the MM is the tent map with the parameter value in the region and it is the logistic map with the parameter value in the region . The map has a discontinuity at the point . Figure LABEL:fig1 shows the structure of the map for various values of for a fixed set of .
As it is clearly evident from Figure. (LABEL:fig1), there is a discontinuity at the point of transition of the map i.e. at which is known as `Border'. As the system is D, the border will be a point at . Depending upon the values of and parameters and , the map can behave as Piecewise Continuous or Piecewise Discontinuous. The map is continuous on each of the regions before and after border, but is discontinuous at border. In case of a Piecewise Discontinuous Map, there is a borderline in the Poincare section such that two arbitrarily close points on the two sides of the border land far apart at the next observation instant [2, 8]. Otherwise, the map is known as the Piecewise Continuous Map [11, 3]. The value of the discontinuity at border would be given by the difference in the functional values of and . We can define a quantity as follows:
[TABLE]
which can quantify the amount of discontinuity in the map at border, as a function of the value of . Owing to the fact that is defined differently for and , the expression of would change for two regimes and , which is as follows:
[TABLE]
For, , the map will have one discontinuity at . Therefore the system will have one border for . But when , the system will have two borders i.e. one at and at where discontinuities will occur. Therefore, the map, which we are taking here, can have single or multiple borders depending on the value of .
Quite clearly, as seen in Figures. (LABEL:fig1b) & (LABEL:fig2),
[TABLE]
At this time the map will be piecewise continuous map with discontinuity at border .
It is quite easy to see, both mathematically and intuitively, that
[TABLE]
This is because or means that the existence of the Mixed Map concept is lost and it is either fully tent-like or fully logistic-like respectively.
The MM, as defined in Equation. LABEL:eqn3 has three parameters which determines the dynamics of the system. While this presents a very wide range of possibilities, it becomes a hard problem to keep track of various parameters and determine their effect on the dynamics of the system when they are varied together. Thus, to simplify the problem, we define a subset of the MM, called the `Reduced Mixed Map' (RMM), in a reduced parameter space of two variables and , where the new parameter is related to the old ones parameters and as follows:
[TABLE]
In the definition of the RMM, the value of remains the same as that of the original Mixed Map and it denotes the Borders' of the map. The value of $r$ is bounded in the interval $[0,2]$ so as to keep the iterates of the map bounded to $[0,1]$. A convenient notation for the RMM would be to reuse the M' notation from Equation. LABEL:eqn3 as follows:
[TABLE]
The structure of the map is shown, for various values of , in Figure. (LABEL:fig3). We can similarly define the amount of discontinuity at for the RMM by plugging in Equation. LABEL:eqn8 into Equation. LABEL:eqn5, which gives us:
[TABLE]
Equations. LABEL:eqn6 and LABEL:eqn7 follows from the previous discussions as well. A 3D plot of Equation. LABEL:eqn10 is shown in Figure (LABEL:fig4). Here, we have shown that the behaviour of if both and are varied in diagram.
2 Location and Stability of the fixed points of the RMM
The RMM can have atmost fixed points whose existence, stability and expression depends on the values of and , with one of them always being at .
The existence of the fixed point is always guaranteed because whatever be the value of and , the line always intersects the map at the origin. The stability of the fixed point is determined by the value of both the parameters. If , that means the map is totally logistic-like, the slope at is given by . Then, by condition of stability,
[TABLE]
Here, we are taking .
But if the value of , then near to the map is tent-like, and thus its slope at that point is given by . Thus, by condition of stability,
[TABLE]
The other two fixed points do not always exist. If , then the map is fully logistic-like, and therefore, the only other fixed point is
[TABLE]
The slope at this point is given by , and therefore, the condition of stability is
[TABLE]
Quite clearly, the limiting case of results in a line of fixed points from till if . This is demonstrated in Figures. (LABEL:fig5a), (LABEL:fig5b) and (LABEL:fig5c). But if the value of is less than , and , then there is no other fixed point other than , which is stable, having an attracting basin as that of the whole of the unit interval , as demonstrated in Figures. (LABEL:fig5d), (LABEL:fig5e) and (LABEL:fig5f).
The case of and is slightly nontrivial and there exists a fixed point other than , given by for . But if , the fixed point again vanishes leaving behind only . This is clearly demonstrated in Figures. (LABEL:fig6a), (LABEL:fig6b) and (LABEL:fig6c). This fixed point has the same stability range as discussed in Equation. LABEL:eqn14. Again, similar to the previous discussion, if and , there again exists a fixed point given by , as shown in Figure. (LABEL:fig7a), (LABEL:fig7b) and (LABEL:fig7c), which is super-stable [1]. This is evident from the fact that \big{|}\frac{d}{dx}M_{r}^{\alpha}(0.5)\big{|}=0.
The fixed points formed on the line till for the case of and has no defined stability in the usual sense because perturbations from the fixed point neither diverge nor converge. But for initial conditions , we see the iterates converge to a fixed point rapidly (in a single step), and the fixed point is given by:
[TABLE]
We see that there exists a region/basin of attraction in both cases where there exists fixed points. For the case and , we see that the fixed point formed at (see Figure.(LABEL:fig7)) has a basin of attraction given by . For the case and , the stable fixed point , which only exists for , has a basin of attraction as same as the previous case (i.e. for case).
The last set of non-trivial fixed points are generated when . For the case , the map again has only two fixed points, (which is unstable) and , which is globally stable upto (similar to Equation. LABEL:eqn14). This is described in Figures. (LABEL:fig8a) to (LABEL:fig8f). The other possibility is that , which results in a peculiar case where three fixed points can exist for a certain range of . For , the fixed point is the characteristic tent map fixed point and for , the fixed point is the characteristic logistic map fixed point. This has been demonstrated in Figures. (LABEL:fig9a) to (LABEL:fig9f). The values of and is given by
[TABLE]
[TABLE]
These can be obtained easily by finding the first intersection point of the tent map and the logistic map respectively with the line. The fixed points themselves are the characteristic fixed points of the tent and the logistic maps and their stability has been discussed before for various values of . The existance and stabilities of all the fixed points have been shown in Table LABEL:table1.
3 Numerical Results
In this section, we will present the various numerical results obtained for the RMM for various values of and . Bifurcation diagrams form the main interest in this section because it clearly shows how the state variable of the map changes with the change of parameters.
The bifurcation diagrams of the parent maps of the RMM is shown in Figures. (LABEL:fig10a) and (LABEL:fig10b). We can see that the bifurcation diagram of the logistic map is distinctively different from that of the tent map, in terms existence of different dynamical phenomena like periodic doubling bifurcations, periodic windows in between chaos, existance of chaos without any other periodic attractors or co-existing attractors. In the case of the RMM, there are two parameters, namely, and , and both of them are suitable candidates for plotting the bifurcation diagrams in the various domains as discussed in the previous section. In our article, we will keep as a constant, and vary and observe the bifurcation diagrams. This would be done for various different values of to understand the overall behaviour of the map for various values of and .
If we keep the value of , i.e. the map is piecewise smooth as well as continuous at border, the bifurcation diagram and corresponding it's lyapunov exponents are shown in Figure. (LABEL:fig11). We see the bifurcation diagram, is in a sense, a mixture of that of the logistic map and the tent map's bifurcation diagrams. As seen in Figure. (LABEL:fig10b) in the case of the tent map, zero is a stable fixed point till . After that, the fixed point jumps to the nontrivial value and the bifurcation diagram has a structure similar to that of a logistic map where after a normal period doubling bifurcation periodic windows exist in between chaos. The corresponding lyapunov exponent shows the existence of the fixed points, periodic orbits and the values of which periodic windows exist in-between chaos. This behaviour is characterized by negative lyapunov exponent for the fixed points and the periodic orbits, zero for the quasi-periodic orbits and positive values for the chaotic orbits.
Even more interesting bifurcation diagrams emerge when we break the continuity of the map at border in case of ) i.e. the piecewise discontinuous map with single or multiple borders. If we use , i.e. there will be two borders, one is at and another one is at , the bifurcation diagram shows a very interesting features. There exists many interesting dynamical structure - there are period-, period- upto period- orbits in this case. The time series waveforms have been shown in Figure. (LABEL:fig12.1). The period- orbit occurs near the value of whereas the distinct period- orbit exists near (see Figure. (LABEL:fig12e)). The period- orbit cannot be seen in case of both logistic and tent map bifurcation diagram. In Figure. (LABEL:fig12a), there exists a -piece chaotic orbits just after the bifurcation near . After that a period- orbit emerges, which gradually gives a periodic orbit as we change the parameter further. The zoomed bifurcation diagram and corresponding it's lyapunov exponents in between to parameter values have been shown in Figure.(LABEL:fig12g) and Figure.(LABEL:fig12h) respectively. As the parameter is varied more, another bifurcation occurs, where a normal period doubling bifurcation gives a period- orbit, which gradudally goes to - piece chaotic orbit. Also, in between to , in the parameter range, a interior crisis [4] happens where a chaotic orbit suddenly expands it's shape. This happens when a chaotic attractor just overlaps with the co-existing unstable chaotic orbit and the main chaotic orbit suddenly expands. Apart from all these, this map also has the normal periodic points and period- orbits as well. Figure. (LABEL:fig12f) demonstrates `sensitive dependence on initial conditions' which is a strong indicator for chaos [6].
Further changing the value of gives rise to other interesting phenomena. Also, leads to approach to fully logistic-like bifurcation diagram whereas leads to approach to fully tent-like bifurcation diagram. These have been demonstrated in Figure. (LABEL:fig13). The zoomed versions of the bifurcation diagrams in Figures (LABEL:fig13a, LABEL:fig13b, and LABEL:fig13c) have been shown in Figure. LABEL:fig13.1. For the sake of completion, bifurcation diagrams of the case , and are also shown in Figures. (LABEL:fig14a),(LABEL:fig14b) and (LABEL:fig14c). The lyapunov exponents are shown in Figures (LABEL:fig13.2 and LABEL:fig14.1) respectively. In the Figures (LABEL:fig13j and LABEL:fig13k), one can see the existences of the two lyapunov exponents in the parameter values, but if one can look closely in the two diagrams, it can be said that the values of the lyapunov exponents toggle between the two lyapunov exponent values in that parameter ranges as the parameter changes gradually and as they have plotted closely, it looks like the existance of the two lyapunov exponents in the same parameter values. The bifurcation diagrams for the cases and have been shown in Figure (LABEL:fig14.2) for the sake of clarity in parameter range to .
From Figure.(LABEL:fig2), it can be said that the amount of discontinuities are same in case of and . But for the first case, as there is only one border which is at , and for the next case, there are two borders, one is at and another is at although the amount of discontinuities are same for the two cases, the bifurcations are different due to different borders. The different bifurcation diagrams have been shown in Figure.(LABEL:fig18).
Upto this, we have taken and . In this consideration, the map becomes continuous across the border in case of . Now, we are taking , which makes the map piecewise discontinuous with discontinuity at .
The bifurcation in the Figure.(LABEL:fig18.1) shows a reverse period adding bifurcation upto . This phenomenon is obversed generally in case of a D piecewise discontinuous map [8]. As the parameter is evolved more. we get a normal period doubling bifurcations which goes to chaotic region. The periodic windows also exist in between chaos here as well.
The zoomed figures of all these bifurcations and corresponding it's lyapunov exponents have been shown in Figure.(LABEL:fig19).
4 A Simulink based implementation of the Mixed Map
In this section, we have shown a Simulink based implementation of the MM (not to be confused with the RMM, which was discussed untill now). The implementation is shown below in Figure. (LABEL:fig15).
The implementation as shown in Figure. (LABEL:fig15) assume the presence of various blocks (which are readily available in Simulink) like the product (multiplier), amplifier, minimum finder, controllable switch and a delay. A basic circuit diagram for implementing discrete maps is shown in [14]. Following a similar method, first, the individual maps are synthesized from the available blocks. The logistic map, which is defined continuously in the interval can be easily synthesized by using a subtracter, product (multiplier) and an amplifier. The tent map, which is piecewise continuous in is implemented by using a subtracter, a minimum finder and an amplifier. The minimum finder takes two inputs and outputs the signal which is lesser of the two. The switching condition is implemented by using a controllable switch which is driven by the output of the circuit. Finally, a unit delay block is used to implement the discrete nature of our system. An initial condition block is used to set the initial condition i.e. , of our system. The values of the switching parameter, the initial condition and the gains of the amplifiers (which set the control parameters and as discussed in Equation. LABEL:eqn3) can be changed manually for each simulation run.
Some sample outputs are shown by running the simulation: Figures. (LABEL:fig16a) and (LABEL:fig16b) shows the period-11 orbit and a aperiodic orbit which were previously observed using numerical techniques. Figure. (LABEL:fig16c) exhibits Sensitive dependence on Initial Conditions (SIC) for two very nearby initial conditions for a pair of identical maps.
5 Conclusion
In this article, we have proposed a novel chaotic map, called the Mixed Map (MM), which is formed from the amalgamation of two well known maps i.e. logistic map and the tent map. It has two parameters and , which are the control parameter and the transition parameter respectively. This map was found to be both piecewise discontinous and continous as well, depending upon the values of the parameters. The amount of discontinuity was investigated and plotted. To reduce the complexity of the map, a Reduced Mixed Map (RMM) was defined for further studies. All the fixed points for the RMM for the various cases of the parameter values was found and their stability was classified. The presence of a three fixed point map was observed for certain values of the parameter . The results have been summarized in Table. LABEL:table1. Numerical studies were done on the proposed RMM and bifurcation diagrams, lyapunov exponents and orbits were plotted for all the interesting cases. A stable period-11 orbit was observed for case. It was also shown that for the bifurcation structure becomes similar to that of the logistic bifurcation structure (as because the map becomes almost totally logistic-like). For , the structure becomes similar to that of the tent-map bifurcation structure. For the case of and has been studied here as a special case. There, we have got a period adding cascade in reverse direction in bifurcation diagram. It can also be shown the dynamics of the map for different with . We can address these issues as our future work. Finally, a simple simulink implementation was presented for the proposed MM, which was used to verify some of the time series diagrams and properties we discussed in the previous sections.
@noitemerr
DB acknowledges the HTRA scholarship provided by Indian Institute of Technology Madras, Chennai. SS acknowledges the support of DST-INSPIRE, Government of India (Ref. No: IF150667). MB acknowledges the support provided by St. Xavier's College (Autonomous), Kolkata.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Kathleen T Alligood, Tim D Sauer, and James A Yorke. Chaos . Springer, 1996.
- 2[2] Viktor Avrutin, Michael Schanz, and Soumitro Banerjee. Multi-parametric bifurcations in a piecewise–linear discontinuous map. Nonlinearity , 19(8):1875, 2006.
- 3[3] Soumitro Banerjee, MS Karthik, Guohui Yuan, and James A Yorke. Bifurcations in one-dimensional piecewise smooth maps-theory and applications in switching circuits. IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications , 47(3):389–394, 2000.
- 4[4] Soumitro Banerjee and George C Verghese. Nonlinear phenomena in power electronics . IEEE, 1999.
- 5[5] Soumitro Banerjee, James A Yorke, and Celso Grebogi. Robust chaos. Physical Review Letters , 80(14):3049, 1998.
- 6[6] Eli Glasner and Benjamin Weiss. Sensitive dependence on initial conditions. Nonlinearity , 6(6):1067–1075, 1993.
- 7[7] Toshiki Habutsu, Yoshifumi Nishio, Iwao Sasase, and Shinsaku Mori. A secret key cryptosystem by iterating a chaotic map. In Workshop on the Theory and Application of of Cryptographic Techniques , pages 127–140. Springer, 1991.
- 8[8] Parag Jain and Soumitro Banerjee. Border-collision bifurcations in one-dimensional discontinuous maps. International Journal of Bifurcation and Chaos , 13(11):3341–3351, 2003.
