Iterative Maps - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : System : Information : Updates : Maple 2015 : updates/Maple2015/IterativeMaps

Iterative Maps

 

Bifurcation maps

Attractor maps

Escape maps

Bifurcation maps

The Bifurcation command generates two-dimensional images of an iterated map. Small changes in values of some parameters in the governing equations can give rise to qualitative changes in the evolving solution of the system, referred to as a bifurcation. For example, the evolving solution of a system of differential equations can switch between being periodic or chaotic in nature.

Images or bifurcation diagrams can illustrate the nature of the solution as the bifurcation parameters change value.

withIterativeMaps:withImageTools:

 

The Gauss map is an example of such a system. Its bifurcation diagram can be computed using the following iterated map of a single real variable.

 

xnew = r + ek x12

 

This map is sometimes referred to as the mouse map since for some parameter values the bifurcation diagram looks like a mouse.

 

In the following generated image the horizontal direction corresponds to r and the vertical direction corresponds to x.

GaussMapImage:=Bifurcationx,r+ⅇ4.8x12,0.5,                                                                   0.2,1.7,xmin=0.7,xmax=2.1,                                                                   iterations=10000 :

ArrayTools:-Dimensions GaussMapImage 

1..500,1..500,1..3

(1.1)

ColouringProcedures:-HueToRGBGaussMapImage:EmbedGaussMapImage

 

 

The Explore command can be used to obtain an application with which to interact with the bifurcation diagram.

This allows a convenient mechanism for exploring the behavior of the iterative process as the number of iterations or the parameter value change.

The following appliable module uses a local member img in order to take advantage of the in-place computational functionality of the Bifurcation command. The following procedure for GaussMap generates an image according to its arguments which specify the parameter k and the number of iterations of the map.

 

 

ExploreGaussMapk,logiter,                                         parameters=k=3.7..9.0,animate=false,logiter=0.0..4.0,label=`log[10](iterations)`,                                         initialvalues=k=4.96,logiter=4.0,                                         placement=left,orientation=vertical,                                        animate,title=Gauss map

Gauss map

 

k

Attractor maps

The Attractor command generates two-dimensional images of the iterated or evolving solution of a system of equations. The command gets its name from the fact that certain systems have the property that if the solution becomes close enough to attaining some special state then it remains close.

For example, the special state may be the situation of being close to a particular point, curve, or surface which is then referred to as an attractor.

 

The Rossler attractor appears in the solution of the following set of differential equations:

ⅆⅆ t xt = ytzt

ⅆⅆ t yt= xt + a yt

ⅆⅆ t zt= b+zt xtc

 

A diagram for this system can be computed using an iterative map which effectively approximates a solution using a fixed step-size dt.

A temporary variable tx is added to the list of variables, so that values for variables x, y, and z can be iterated in unison.

 

dt:=0.0015:

xnew:=x+dtyz:

ynew:=y+dttx+ay:

znew:=z+dtb+ztxc:

 

The behavior of this iterated map depends on the values of the parameters a, b, and c. Using the following values, the computed orbit exhibits chaotic behavior around a fixed attractor. For other values of the parameter c, the evolving solution becomes periodic in nature.

 

a,b,c:=0.1,0.1,14.0:

RosslerImage,xrng,yrngAttractorheight=400,width=400,                                                                                x,y,z,tx,                                                                                xnew,ynew,znew,x,                                                                                1,1,1,1,                                                                                xmin=24,xmax=24,ymin=24,ymax=22,fixview,                                                                                iterations=300000:

ArrayTools:-Dimensions RosslerImage 

1..400,1..400,1..3

(2.1)

ColouringProcedures:-HueToRGBRosslerImage:ImageTools:-EmbedRosslerImage:

 

The previous image has a view onto the x and y variables' space. An oblique view can be attained by introducing two additional variables and a projection. Here, a projection is made using a rotation about the first axis. A fixed viewing range is used.

 

The speed of computation is due to internal use of the Maple Compiler to apply the iterative map.

ang:=π8:RosslerImage,xrng,yrngCodeTools:-UsageAttractorheight=400,width=400,                                                                                                                         xview,yview,x,y,z,tx,                                                                                                                         x,28ysinang+zsinπ2ang,xnew,ynew,znew,x,                                                                                                                         1,1,1,1,1,1,                                                                                                                         xmin=24,xmax=24,ymin=20,ymax=28,fixview,                                                                                                                         iterations=500000:

memory used=3.69MiB, alloc change=3.66MiB, cpu time=30.00ms, real time=35.00ms, gc time=0ns

ColouringProcedures:-HueToRGBRosslerImage:ImageTools:-EmbedRosslerImage:

 

 

The Explore command can be used to obtain an application which can be used to interact with the Attractor diagram.

This allows a convenient mechanism for exploring the behavior of the solution as the number of iterations or the parameters change. The tilt angle can also be made interactive.

The following appliable module uses a local member img in order to take advantage of the in-place computational functionality of the Attractor command. The following procedure, Rossler, will generate an image according to its arguments which specify the parameter c, the number of iterations of the map, and the angle of tilt in the projected view.

 

 

The procedure Rossler can be passed to the Explore command to provide a convenient interface for varying the values of its arguments.

 

ExploreRosslerC,logiter,angle,                    parameters=C=4.0..14.0,animate,label='c',orientation=horizontal,logiter=4.0..6.7,label=`log[10](iterations)`,angle=π2..0,                    initialvalues=logiter=5.1,angle=π6,                    placement=left,orientation=vertical,                    animate=false,title=Rossler attractor

Rossler attractor

 

angle

Escape maps

The Escape command generates images from iterated maps involving several variables.

The dimensions of the resulting image correspond to the initial values of the first two specified variables of the map. The iterated computation can involve several variables, and operates by updating the value of each two-dimensional image using the previous value and (optionally) the values of any additional variables.

The iteration continues until either the supplied iteration limit is exceeded or a supplied escape condition is met.

 

One well-known application of such iterative maps is the generation of images for which the approximation of the boundaries of the regions that escape are fractal in nature.

 

expr:=zr+ziI3

expr:=zr+Izi3

(3.1)

fzi:=evalcℑexpr

fzi:=zi3+3zizr2

(3.2)

fzr:=evalcℜsubszi=zit,expr

fzr:=3zit2zr+zr3

(3.3)

mandelbroid:=Escapezi&comma;zr&comma;zit&comma;zrsqr&comma;zisqr&comma;                                                 fzi&plus;y&comma;fzr&plus;x&comma;zi&comma;zr2&comma;zi2&comma;                                                 y&comma;x&comma;y&comma;x2&comma;y2&comma;                                                 4<zrsqr&plus;zisqr&comma;                                                 1.2&comma;1.2&comma;1.2&comma;1.2&colon;

ArrayTools:-Dimensions mandelbroid 

1..500&comma;1..500&comma;1..3

(3.4)

ColouringProcedures:-HueToRGBmandelbroid&colon;Embedmandelbroid