Avoiding Pitfalls In Periodic Orbits of Maps With A Piecewise Approach

Numerical computation is widely recognized to be of great importance in many fields of science. Many conclusions in nonlinear science and complex systems have been drawn upon simulation on a digital computer. Although this importance is widely recognized, the way of computers work is, in general, not known in detail.

Let us start our discussion with the following example. Let be the following equation, known as a logistic map:

The idea behind this equation is very simple. We start a simulation at stage k = 0. In this stage, we have our initial condition. Suppose that r = 4, and It is easy to calculate that in the next stage, we will get ¼. After the third iteration, the correct answer is going to be stuck — fixed — with the value of ¾.

Using a commercial software, we proceed with this calculation, and we can get the results in Figure 1. The red line and the black line are produced by the equations F(xk) = 4xk(1 – xk), while G(xk) = 4xk – 4x2k. Although the equations are mathematically equivalent, the computer produces a stable and fixed result for F and a chaotic behavior for G. Nonlinear dynamic theory and a branch of math based only on calculus and differential equations cannot explain this properly.

Figure 1: Simulation of the logistic map from two different computer realizations. Image courtesy Erivelton Geraldo Nepomuceno

In our recent study entitled “Interval computing periodic orbits of maps using a piecewise approach,” we recognize this situation and developed an approach to, at least, decrease this situation. The basic idea relies on the fact that the majority of numerical computations based on floating-point are built to guarantee a good result, or, in technical terms, a result within a precision according to the number of bits used, only for basic arithmetical operations.

The example in Figure 1 presents two features that turn this example in a very challenging one. First, the computation is a set of basic arithmetic operations. Second and most important, there is a recurrence where a very tiny error in the beginning grows exponentially. In our paper, we examine some features of the function, which allows us to set bounds of the result with more consistency than other approaches. The function is seen as a piecewise function. In general, the proposed method has produced intervals that are significantly narrower than those obtained by Intlab (a Matlab toolbox) approach. However, it is also clear that our method requires a significantly larger number of iterates.

Related Study

In an accompanying article published in the same journal “Detecting unreliable computer simulations of recursive functions with interval extensions,” https://doi.org/10.1016/j.amc.2018.02.020, we used such ideas to show that some previous results that guarantee long simulations on a computer are not valid for current software and hardware. We also exploit the error propagation to compute the Lyapunov exponent, a usual feature to identify chaotic systems, in the article “Computation of the largest positive Lyapunov exponent using rounding mode and recursive least square algorithm”, https://doi.org/10.1016/j.chaos.2018.04.032.

Future Directions

We are confident that computers are incredible tools to develop science and technology. But we are also confident that we must understand the way that they work in order to be able to analyze and avoid pitfalls.

These findings are described in the article entitled Interval computing periodic orbits of maps using a piecewise approach, recently published in the journal Applied Mathematics and ComputationThis work was conducted by Erivelton G. Nepomuceno, Heitor M. Rodrigues Junior, and Samir A.M. Martins from the  Federal University of São João del-ReiMatjaž Perc from the University of Maribor and Beihang University, and Mitja Slavinec from the University of Maribor.

About The Author

Erivelton Geraldo Nepomuceno

Erivelton G. Nepomuceno received the BSc degree from UFSJ and the PhD in Electrical Engineering from the UFMG (Brazil). From 2013 to 2014 he was a Post-doctoral Research Fellow at Intelligent Systems and Networks Group, Imperial College London. He is currently Associate Professor at UFSJ and leader of Control and Modelling Group (GCOM). His research interests include chaos, complex systems, computer arithmetic, interval arithmetic, system identification. He is Associate Editor of the IEEE Latin America Transactions and Journal of Biomedical Research and Reviews. He is member of Brazilian Association of Automatica and IEEE Circuit and Systems Society (IEEE CAS). Since 2018, he has been appointed as member of the Technical Committee on Nonlinear Circuits and Systems (TC-NCAS) of the IEEE CAS

Speak Your Mind!

READ THIS NEXT

Investigating Mechanisms Of Cancer Cell-Intrinsic CYP Monooxygenases That Contribute To Tumor Progression

There has been great interest in finding structural targets in cancer for the widely prescribed type 2 diabetes drug metformin. Metformin is a synthetic product derived from galegine, a guanide natural product found in French lilac. Documented use of galegine for human disease dates back to 1620, when John Parkinson first described the use of […]

Lessons Learned From Serial Isolations Of Human Coronavirus 229E From Environmental Surfaces Of A Classroom During Influenza Season

The word “infection” has different connotations in the fields of virology and medicine. Generally, physicians refer to a condition wherein a patient displays symptoms of illness due to an infectious agent when they use the word “infection.” In virology, the word infection is not synonymous with outward signs of illness, because people can have a […]

3 Cases Of Energy Policy In Brazil

During the last few years, solutions to produce electricity in a decentralized manner have become increasingly attractive. In this context, can decentralized technologies enter centralized market arrangements, even if these market arrangements contain policies to facilitate investment in low carbon technologies? The answer may depend on the particular rules governing the sector. If, as expected, […]

What Type Of Star Is The Sun?

What type of star is the Sun, it is a G2V or a yellow dwarf in the main sequence of a star’s life cycle. Our star is not a particularly abnormal or unusual type, and experts agree it will be here for another 5 or so billion years. Our solar system is filled with a diverse […]

Every European Capital And Associated Country

One of the 7 continents in the world, Europe is located in the northern hemisphere, west of Asia and north of Africa. Europe is separated from Asia by the Caucus mountains and the Black Sea and separated from Africa by the Mediterranean Sea. Europe has a land area of ~4,000,000 sq. mi, occupying about 2% […]

First Ever Quantum-Encrypted Conference Call Completed

As research and quantum computing advances, more and more possible applications for the technology are coming to light. Quantum cryptography can be used to help people ensure that transmissions between them and other parties are confidential, thanks to the way elementary particles behave, and the quantum internet could serve as a specialized and secure portion […]

Irrigation Quantification From Space Exploiting Satellite Soil Moisture Products

Over 2 billion people are currently affected by water stress, a number that is expected to dramatically increase with population growth, economic development and through the effects of climate change (United Nations, 2018). Feeding the world’s growing and more affluent population will require a 25-70% increase in agricultural production by 2050 with a corresponding increase […]