scholarly journals Dynamics of Numerics of Nonautonomous Equations with Periodic Solutions: Introducing the Numerical Floquet Theory

2013 ◽  
Vol 2013 ◽  
pp. 1-11
Author(s):  
Melusi Khumalo

Nonautonomous systems with periodic solutions are encountered frequently in applications. In this paper, we will consider simple systems whose solutions are periodic with a known period. Their transformation under linearized collocation methods is investigated, using a technique called stroboscopic sampling, a discrete version of the well-known Poincaré map. It is shown that there is an inextricable relationship between AN stability (or BN stability) of the numerical methods and the correct qualitative behaviour of solutions.

2007 ◽  
Vol 17 (03) ◽  
pp. 837-850 ◽  
Author(s):  
SHIGEKI TSUJI ◽  
TETSUSHI UETA ◽  
HIROSHI KAWAKAMI

The Bonhöffer–van der Pol (BVP) oscillator is a simple circuit implementation describing neuronal dynamics. Lately the diffusive coupling structure of neurons attracts much attention since the existence of the gap-junctional coupling has been confirmed in the brain. Such coupling is easily realized by linear resistors for the circuit implementation, however, there are not enough investigations about diffusively coupled BVP oscillators, even a couple of BVP oscillators. We have considered several types of coupling structure between two BVP oscillators, and discussed their dynamical behavior in preceding works. In this paper, we treat a simple structure called current coupling and study their dynamical properties by the bifurcation theory. We investigate various bifurcation phenomena by computing some bifurcation diagrams in two cases, symmetrically and asymmetrically coupled systems. In symmetrically coupled systems, although all internal elements of two oscillators are the same, we obtain in-phase, anti-phase solution and some chaotic attractors. Moreover, we show that two quasi-periodic solutions disappear simultaneously by the homoclinic bifurcation on the Poincaré map, and that a large quasi-periodic solution is generated by the coalescence of these quasi-periodic solutions, but it disappears by the heteroclinic bifurcation on the Poincaré map. In the other case, we confirm the existence a conspicuous chaotic attractor in the laboratory experiments.


1980 ◽  
Vol 47 (3) ◽  
pp. 645-651 ◽  
Author(s):  
L. A. Month ◽  
R. H. Rand

The stability of periodic motions (nonlinear normal modes) in a nonlinear two-degree-of-freedom Hamiltonian system is studied by deriving an approximation for the Poincare´ map via the Birkhoff-Gustavson canonical transofrmation. This method is presented as an alternative to the usual linearized stability analysis based on Floquet theory. An example is given for which the Floquet theory approach fails to predict stability but for which the Poincare´ map approach succeeds.


2020 ◽  
Vol 54 (1) ◽  
pp. 335-358 ◽  
Author(s):  
Dongling Wang ◽  
Aiguo Xiao ◽  
Jun Zou

In this work, we study the long time behavior, including asymptotic contractivity and dissipativity, of the solutions to several numerical methods for fractional ordinary differential equations (F-ODEs). The existing algebraic contractivity and dissipativity rates of the solutions to the scalar F-ODEs are first improved. In order to study the long time behavior of numerical solutions to fractional backward differential formulas (F-BDFs), two crucial analytical techniques are developed, with the first one for the discrete version of the fractional generalization of the traditional Leibniz rule, and the other for the algebraic decay rate of the solution to a linear Volterra difference equation. By means of these auxiliary tools and some natural conditions, the solutions to F-BDFs are shown to be contractive and dissipative, and also preserve the exact contractivity rate of the continuous solutions. Two typical F-BDFs, based on the Grünwald–Letnikov formula and L1 method respectively, are studied. For high order F-BDFs, including convolution quadrature schemes based on classical second order BDF and product integration schemes based on quadratic interpolation approximation, their numerical contractivity and dissipativity are also developed under some slightly stronger conditions. Numerical experiments are presented to validate the long time qualitative characteristics of the solutions to F-BDFs, revealing very different decay rates of the numerical solutions in terms of the the initial values between F-ODEs and integer ODEs and demonstrating the superiority of the structure-preserving numerical methods.


10.26524/cm41 ◽  
2018 ◽  
Vol 2 (2) ◽  
pp. 80-88
Author(s):  
Gokulakrishnan S ◽  
Chandrasekar V

In this paper, we derive the discrete version of the Bernoulli’s formula according to the generalized α- difference operator for negative `l ,and to find the sum of several type of arithmetic series in the field of Numerical Methods. Suitable example are provided to illustrate the main results.


1968 ◽  
Vol 11 (5) ◽  
pp. 743-745 ◽  
Author(s):  
T.A. Burton

In [1] Demidovic considered a system of linear differential equationswith A(t) continuous, T-periodic, odd, and skew symmetric. He proved that all solutions of (1) are either T-periodic or 2T-periodic0 In [2] Epstein used Floquet theory to prove that all solutions of (1) are T-periodic without the skew symmetric hypothesis. Epstein's results were then generalized by Muldowney in [7] using Floquet theory. Much of the above work can also be interpreted as being part of the general framework of autosynartetic systems discussed by Lewis in [5] and [6]. According to private correspondence with Lewis it seems that he was aware of these results well before they were published. However, it appears that these theorems were neither stated nor suggested in the papers by Lewis.


2010 ◽  
Vol 51 (3) ◽  
pp. 350-368 ◽  
Author(s):  
LIE-HUI ZHANG ◽  
YONG WANG

AbstractCriteria for guaranteeing the existence, uniqueness and asymptotic stability (in the sense of Liapunov) of periodic solutions of a forced Liénard-type equation under certain assumptions are presented. These criteria are obtained by application of the Manásevich–Mawhin continuation theorem, Floquet theory, Liapunov stability theory and some analysis techniques. An example is provided to demonstrate the applicability of our results.


2012 ◽  
Vol 2012 ◽  
pp. 1-17 ◽  
Author(s):  
Hunki Baek

A Lotka-Volterra-type predator-prey system with state-dependent feedback control is investigated in both theoretical and numerical ways. Using the Poincaré map and the analogue of the Poincaré criterion, the sufficient conditions for the existence and stability of semitrivial periodic solutions and positive periodic solutions are obtained. In addition, we show that there is no positive periodic solution with period greater than and equal to three under some conditions. The qualitative analysis shows that the positive period-one solution bifurcates from the semitrivial solution through a fold bifurcation. Numerical simulations to substantiate our theoretical results are provided. Also, the bifurcation diagrams of solutions are illustrated by using the Poincaré map, and it is shown that the chaotic solutions take place via a cascade of period-doubling bifurcations.


1973 ◽  
Vol 40 (4) ◽  
pp. 1091-1096 ◽  
Author(s):  
G. A. Thurston

Application of Newton’s method to nonlinear vibration problems can lead to a sequence of nonhomogeneous ordinary differential equations with periodic coefficients. The form of the complementary solutions are known from Floquet theory. This paper suggests a method for avoiding “secular terms” that grow with time in the particular solution. The method consists of finding a single periodic solution of the complementary solutions and its adjoint. If the periodic solution exists, a frequency correction can be computed that eliminates secular terms. After the frequency correction, the rest of the particular solution is periodic and can be computed by the infinite determinant method or other numerical methods. In oversimplified terms, the procedure is to find the improved approximation to the period by variation of parameters and the next approximation to the amplitudes by undetermined coefficients which is a simpler computation than variation of parameters.


Sign in / Sign up

Export Citation Format

Share Document