Re: the circle map
- To: mathgroup at smc.vnet.net
- Subject: [mg52267] Re: the circle map
- From: p-valko at tamu.edu (Peter Valko)
- Date: Sat, 20 Nov 2004 03:41:36 -0500 (EST)
- References: <cneuml$qbc$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Roger,
Can you tell me why is it that in the following code of yours:
Clear[x, y, n];
a0 = 0.41209;
x[n_] := x[n] = Mod[-a0*x[n - 1] - y[n - 1], 1];
y[n_] := y[n] = Mod[x[n - 1], 1] ;
x[0] = 0.7;
y[0] = .65;
a = Table[{x[n], y[n]}, {n, 0, 10000}];
ListPlot[a, PlotRange -> All] ;
we get a fractal-like pic, but changing to a0 = 0.41208 we do not?
(And why is that replacing Mod[-,1] by FractionalPart[-] in the above
code will not give the same phenomenon?
Peter
Roger Bagula <tftn at earthlink.net> wrote in message news:<cneuml$qbc$1 at smc.vnet.net>...
> I've done a lot of searches on chaos
> and Mathematica and have never seem this.
> It is sensative chaos , in both the angle based a0 and the
> initial starting point.
> The circle was used by Chua as a starting point in his lectures on Chaos.
>
> Clear[x,y,a,b,s,g,a0]
> (* circle map: from Chaos in Digital Filters ,Chua,Lin,
> IEEE transactions on Circuits and Systems,vol 35 no 6 June 1988*)
> (* very sensitive to intial conditions*)
> a0=Cos[Pi/6]/2;
> x[n_]:=x[n]=Mod[-a0*x[n-1]-y[n-1],1]
> y[n_]:=y[n]=Mod[x[n-1],1]
> x[0]=0.7;y[0]=.65;
> a=Table[{x[n],y[n]},{n,0, 10000}];
> ListPlot[a, PlotRange->All]
>
> Respectfully, Roger L. Bagula
>
> tftn at earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
> alternative email: rlbtftn at netscape.net
> URL : http://home.earthlink.net/~tftn
- Follow-Ups:
- Re: Re: the circle map
- From: DrBob <drbob@bigfoot.com>
- Re: Re: the circle map