Govaerts, W.

Kouznetsov, Iouri Aleksandrovitsj

Meijer, Hil Gaétan Ellart

Neirynck, N.

English

MatcontM is a matlab toolbox for numerical analysis of bifurcations of fixed points and periodic orbits of maps. It computes codim 1 bifurcation curves and supports the computation of normal coefficients including branch switching from codim 2 points to secondary curves. Recently, the initialization and computation of connecting orbits was improved. Moreover, a graphical user interface was added enabling interactive control of all these computations. To further support these computations it allows to compute orbits of the map and its iterates and to represent them in 2D, 3D and numeric windows. We demonstrate the use of the toolbox in a study of Arnol'd tongues near a degenerate Neimark-Sacker (Chenciner) bifurcation. Here we illustrate the recent theory of [Baesens&amp;Mackay,2007] how resonance tongues interact with a quasi-periodic saddle-node bifurcation of invariant curves in maps. Using normal form coefficients we find evidence for one of their cases, but not the other. Actually, we find another unfolding, i.e. a third possibility. We also find a structure that resembles a quasi-periodic cusp bifurcation of invariant curves

University of Twente Research Information

ENOC 2011, 24-29 July 2011, Rome, ItalyA study of resonance tongues near a Chenciner bifurcation using MatcontMGovaerts W. 1, Kuznetsov Yu.A. 2, Meijer H.G.E.3 and Neirynck, N.11 Dept. of Applied Mathematics and Computer Science, Ghent University, Gent, Belgium2 Dept. of Mathematics, Utrecht University, Netherlands3 Department EEMCS, University of Twente, PO Box 217, 7500 AE Enschede, NetherlandsSummary. MatcontM is a matlab toolbox for numerical analysis of bifurcations of fixed points and periodic orbits of maps. Itcomputes codim 1 bifurcation curves and supports the computation of normal coefficients including branch switching from codim 2points to secondary curves. Recently, the initialization and computation of connecting orbits was improved. Moreover, a graphical userinterface was added enabling interactive control of all these computations. To further support these computations it allows to computeorbits of the map and its iterates and to represent them in 2D, 3D and numeric windows. We demonstrate the use of the toolbox in astudy of Arnol’d tongues near a degenerate Neimark-Sacker (Chenciner) bifurcation. Here we illustrate the recent theory of [4] howresonance tongues interact with a quasi-periodic saddle-node bifurcation of invariant curves in maps. Using normal form coefficientswe find evidence for one of their cases, but not the other. Actually, we find another unfolding, i.e. a third possibility. We also find astructure that resembles a quasi-periodic cusp bifurcation of invariant curves.MatContM: Numerical Bifurcation Analysis of MapsWe consider dynamical systems generated by (nonlinear) smooth mapsx 7→ f(x, α), x ∈ Rn, α ∈ Rm (1)with state variable x and free parameters α. The understanding of the dynamics of the map (1) requires most often numer-ical bifurcation analysis with the use of a dedicated software package. For this purpose we have written MATCONTM, aMATLAB continuation toolbox available at http://sourceforge.net/projects/matcont, [1]. It computes bifurcations of fixed(periodic) points in one parameter and continues the bifurcation curve in two parameters and detects codim 2 bifurcationpoints, similar to other packages. In addition, critical normal form coefficients are computed where for fixed points withhigh iterates we rely on symbolic derivatives or if these are not available, on automatic differentiation [3]. Also, we grow1D stable and unstable invariant manifolds of fixed points for n-dimensional maps, using generalizations of the standardalgorithms for n = 2. For planar systems the search of intersection points is automated allowing easy initialization ofcontinuation of homoclinic and heteroclinic orbits and also their tangencies. The most recent addition is a graphical userinterface that allows interactive use of the various routines, see Figure 1.Working with the command line versionThe command line version contains a continuer that computes bifurcation curves. These curves are defined by a system ofequations consisting of fixed point and bifurcation conditions. With one free parameter we have curves of fixed points andconnecting orbits. With two free system parameters we can compute bifurcation curves of limit point, period-doublingand Neimark-Sacker as well as tangencies of homoclinic and heteroclinic manifolds. For the actual continuation we dothe following steps.We first select an initial point. This point may be obtained from analysis, simulations or previous continuations. We theninitialize the continuer by setting certain variables. For instance, we set the map that we investigate, the system that definesthe bifurcation curve, the use of derivatives, the iteration number of a fixed point and so on. Next, we specify continueroptions such as tolerances and maximal stepsize. Then we call the continuer that computes the curve. Simultaneously,the continuer monitors the occurrence of zeros of test functions. Such a zero may correspond to a bifurcation of highercodimension. Finally, we can inspect the data and repeat this process.The organization of the GUIThe GUI has been developed using the object-oriented syntax added in Matlab since version 7.6 (2008a). There are twoimportant aspects of the GUI. First, there is a model that controls data and via the GUI the user can generate and managesuch data. Figure 2 illustrates the implementation of this idea. Second, we have built the GUI on top of the commandline (cl) version of MatcontM. That is, we keep the cl-version intact and call the continuer, the computational core, withcorrect settings from the GUI.The data consists of a current system and a current curve. The current system specifies the map currently analyzed andspecial user functions to be monitored in addition of standard test functions. There is also a current diagram to whichnewly computed curves are added. The current curve consists of an initial point type and a curvetype. This data can beinspected using the Data browser. We have the folllowing natural hierarchy in the data:• System This is the map that we currently analyze. A new system can be added or an existing one deleted.• Diagrams Computed curves are stored in an active diagram, i.e. directory. The user can select one active diagramamong multiple. Data can be transferred between diagrams as well.EN OC 2011 , 24-29 July 2011, Rome, ItalyENOC 2011, 24-29 July 2011, Rome, ItalyFigure 1: A screenshot of the GUI for MATCONTM. The main menu allows to select or edit the map, select initial point and curvetype, set various continuation settings, start the continuation. Output can be generated graphically with 2D and 3D plots or numericwindows.Cl MatcontMGUI session object Diagram managerBranch ManagerWindow ManagerPlot ManagerOptionsStarter DataContinuer DataCurrent SystemCurrent DiagramCurrent CurveSession dataSession managerFigure 2: Model view control scheme: The model controls the data via several managers.• Curves A diagram consists of computed curves. The automatic naming scheme adopted is as follows: dn PT CT(x).Here d stands for the direction along the curve (forward or backward), n is the iteration number, PT is the initialpoint type, CT is the curve type and x is the number of the curve among all curves of its particular type. A curvecan be loaded to start the continuation of a similar curve. We can also rename or delete curves.• (Special) Points A curve actually consists of a set of (ordered) points. We can browse these by loading them intothe MATLAB array editor or selecting special points. These special points correspond to bifurcations of highercodimension along the curve. The first and last point are also always given. We can check the settings of the starterand continuer by inspecting the Starterdata and Continuerdata, respectively.This data is managed through several functions, which we describe below.• Startdata This function controls the initial point, i.e. the state variables and parameters, and the iteration number.Here we select active parameters and adjust settings for derivatives. This is also called Starter.• Continuer This function controls the general parameters of the continuation, such as minimal and maximal step-sizes, the number of points to be computed, when to decide that the computation of a bifurcation point was success-ful, etc.• Branchmanager This function monitors the selection of the initializer currently chosen. There are certain codim 2bifurcations from which branches with higher iteration number (2, 3 or 4 times the original period) can be started.Therefore, the correct iteration number is managed by this function as well.EN OC 2011 , 24-29 July 2011, Rome, ItalyENOC 2011, 24-29 July 2011, Rome, ItalyFigure 3: A screenshot of the data browser for MATCONTM. One can specify new systems or edit existing ones. The main menu allowsto select or edit the map, select initial point and curve type, set various continuation settings, start the continuation. Output can begenerated graphically with 2D and 3D plots or numeric windows.• Curvemanager This function decides how a computed curve is named and where it is saved. It further decides(using Options) whether old curves are deleted. The user can move curves to different diagrams, rename or deletethem.• Windowmanager This function controls which windows are open, where they are placed on the screen and whichdata are represented.• Plotmanager This function interacts between the sessiondata and the windows, using the Options.• Options The user can specify preferences on the runtime options. One can pause the computation after each point,at special points only or never. Computed points are printed in a dialog window after a number of steps that can bespecified. The maximum number of computed curves of a particular type can be specified. This limits the size ofthe archive.New data can be generated by an actual continuation. This proceeds similar as when using the cl-version but automated.First, the continuer options are generated and some sanity checks are performed. We need to check whether we extendthe current curve or start a new curve. Command line data structures are set up. Next we call the initializer and then thecontinuer. During computation there is a simple output screen displaying the computed points and other messages. Thesemessages are start, end or failure of the continuation and special points together with their normal form coefficients. Theinformation of the current point can also be monitored using a numeric window. Finally, we can plot the computed dataduring the continuation as described above using the Plotmanager.We note that MatContM also allows simulation, i.e. the computation of orbits of the map and its iterates and theirrepresentation in 2D, 3D and numerical windows. This is useful both to find starting data for continuation and for thevalidation of the results obtained by continuation.Resonance bubbles near a Chenciner bifurcation: theoryWe use the software here to understand the precise bifurcation structure near a degenerate Neimark-Sacker (NS), i.e. aChenciner (CH), bifurcation. It is known that away from the CH point one can expect either stable or unstable invariantcurves to emanate from the NS bifurcation, see Figure 4(left). The boundary between the parameter region with twoinvariant curves and that with only a fixed point has a delicate bifurcation structure. These invariant curves may haveirrational enough, i.e. Diophantine, rotation number and then collide in a quasi-periodic saddle-node bifurcation of invari-ant curves. This happens on a set of positive measure but excludes a set of parameters that is the collection of diamondshaped regions. Inside these diamonds, the rotation number is rational, i.e. there is q-periodic behavior inside an Arnol’dtongue in parameter space. When the Arnol’d tongue passes through the diamond (in parameter space) a delicate bifur-cation structure arises as proposed by [4]. Their unfolding is shown schematically with local bifurcations only in Figure4(middle and right).Resonance bubbles near a Chenciner bifurcation: exampleWe want to see the structure described above in the following map [6]F : xyz 7→ ybx+ k + zyz − (bx+ k + zy − 1)ky/(c+ y2) (2)EN OC 2011 , 24-29 July 2011, Rome, ItalyENOC 2011, 24-29 July 2011, Rome, Italy213CHp:q tongue NS LPR1LPNSCHR1Figure 4: Unfolding near a Chenciner bifurcation with sub- and super-critical NS bifurcations on either side. Along a parabolic curvea stable and unstable invariant curve collide. This happens as a quasi-periodic saddle-node bifurcation or through a resonance bubble.Such a resonance occurs inside an Arnol’d tongue demarcated by two limit point curves. According to [4], this resonance bubbleconnects two R1 points via a NS curve that can either be simple (middle) or contain a figure eight (right). Labels: LP=limit point andR1=1:1 resonance.with c = 0.1 and k, b free parameters. A partial bifurcation diagram is shown in figure 5(left), other details may be foundin [6]. Arnol’d tongues of period q (4 < q < 101) are shown. Initial points were computed by simulation and checkedfor periodicity. Then the outer boundaries of the tongue were computed with continuation. Along these boundaries wedetected R1 points of fixed points of period q.Some of these R1 points come in pairs, one on each boundary of the tongue. This occurs for with−.53 < b < −.45. Thesepairs correspond to resonances near a quasi-periodic saddle-node bifurcation and their rough unfolding is sketched above.This occurs for with −.53 < b < −.45. MATCONTM automatically computes the normal form coefficient s yielding thelocal stability of the bifurcating invariant curve. For a simple resonance bubble, the sign of s on both sides of the tongue isequal, while in the other case it is not. First, we inspected the sign of s for each R1 point in a bubble. These coincided forall q < 100. Second, we checked the geometry of resonance bubbles. None of the NS-curves showed self-intersections.Summarizing, we find only the evidence for the simple unfolding. The situation is however more complicated. We detecttwo Chenciner bifurcations for some resonance bubbles of high periods, see Figure 5(right), yielding a third possibilityfor unfolding a bubble.0.5 0.7 0.9 1.1 1.3 1.6−0.54−0.52−0.5−0.48−0.46−0.44R4  CH  R3    kb1:51:61:71.29 1.33 1.37 1.41 1.45−0.53−0.5275−0.525−0.5225−0.52  kbCH1.196 1.1963 1.1965−0.5112−0.511−0.5108  kbR1R1NS LPCHFigure 5: Left: Partial bifurcation diagram of map (2). The black line indicates the NS bifurcation of the fixed point (x, y, z) =(1, 1, 1− b− k) with strong resonances (R3 and R4) and a Chenciner bifurcation (near b = −.5238, k = 1.30) indicated by red dots.Some weak resonance tongues are indicated. All other red dots correspond to R1 points of fixed point with period>4. The green dashedcurve indicates for which parameters approximately, two invariant curves collide. Middle: some tongues loop around and display tworesonance bubbles. Right: Resonance bubble within an Arnol’d tongue with period 15 with two Chenciner points.We take the resonance bubbles as an indication where the quasi-periodic saddle-node bifurcation occurs, which we haveplotted with a green dashed curve. Theoretically this curve should depart quadratically from the Chenciner point whichit appears to do indeed. Continuing we see that it turns around suggestive of a quasi-periodic cusp bifurcation. Indeed,also the Arnol’d tongues make a turn and show two resonance bubbles, see Figure 5(middle). We note however thatvery close to this point we have also found a homoclinic structure of a period 3 orbit starting at the R3 point. The fullbifurcation diagram is still under investigation. One feature that follows from the double folding of the Arnol’d tonguesis that three global invariant curves coexist for a small parameter region, with or without resonances. One emanates fromthe supercritical branch of the Neimark-Sacker bifurcation curve. It has small amplitude and is stable. There is a largestable invariant curve and one unstable in between, see Figure 6.References[1] W. Govaerts, R. Khoshsiar Ghaziani, Yu.A. Kuznetsov and H.G.E. Meijer, Numerical methods for two-parameter local bifurcation analysis of maps,SISC 29: 2644-2667, 2007.[2] W. Govaerts, R. Khoshsiar Ghaziani, Yu.A. Kuznetsov, and H.G.E. Meijer, Numerical continuation of connecting orbits of maps, Journal ofDifference Equations and Applications, 15 (8): 849-875, 2009.[3] J.D. Pryce, R. Khoshsiar Ghaziani, V. De Witte, W. Govaerts, Computation of normal form coefficients of cycle bifurcations of maps by algorithmicdifferentiation, Mathematics and Computers in Simulation 81: 109-119, 2010.EN OC 2011 , 24-29 July 2011, Rome, ItalyENOC 2011, 24-29 July 2011, Rome, Italy1.343 1.353 1.363 1.373−0.5255−0.5235−0.5215R1  R1  kb1.3496 1.35 1.3504−0.5252−0.5251−0.525LP10R1NS0.2 0.65 1.1 1.55 20.20.651.11.552xy0.2 0.65 1.1 1.55 20.20.651.11.552xy1 1.25 1.50.520.560.6Figure 6: Left: Bifurcation diagram showing an Arnol’d tongue of period 10 (LP 10-curves blue) emanating from the Neimark-Sackercurve (black). The tongue makes two turns accompanied by the occurrence of a resonance bubble with (closed) NS-curves of period10 (black). The inset shows an enlargement of the lower bubble. Middle: Phase portrait of (2) at (b, k) = (−.525, 1.35). Twostable invariant curves, there is an unstable one in between. Right: Phase portrait of (2) at (b, k) = (−.525115, 1.35015). The outerinvariant curve is close to 1:10 phase-locking. The two inner invariant curves have broken up in a resonance structure with saddles (x)and foci (o). In this bubble there is a supercritical Neimark-Sacker resulting in 10 small invariant curves that are periodically visitedunder iteration. The black lines indicate the unstable manifold of the saddle (for clarity only one is shown). The inset shows the smallinvariant curve and the unstable manifold more clearly.[4] Baesens C. and MacKay R.S., Resonances for weak coupling of the unfolding of a saddle-node periodic orbit with an oscillator, Nonlinearity 20:1283-1298, 2007[5] Chenciner A, Bifurcations de points fixes elliptiques: III. Orbites périodiques de petites périodes et élimination résonnante des couples de courbesinvariantes Publ. Math. IHES 66: 5-91, 1989[6] Frouzakis C.E., Adomaitis R.A. and Kevrekidis I.G., Resonance phenomena in an adaptively-controlled system, Internat. J. Bifur. Chaos Appl. Sci.Engrg. 1:, 83-106, 1991EN OC 2011 , 24-29 July 2011, Rome, Italy

A study of resonance tongues near a Chenciner bifurcation using MatcontM

https://ris.utwente.nl/ws/files/5367355/Govaerts_et_al_6page.pdf

A study of resonance tongues near a Chenciner bifurcation using MatcontM

Abstract

Similar works

Full text

Available Versions

University of Twente Research Information