Gumowski-Mira-Attractor

Use this forum to post Vensim related questions.
josef
Junior Member
Posts: 2
Joined: Mon Aug 29, 2011 12:29 pm

Gumowski-Mira-Attractor

Postby josef » Mon Aug 29, 2011 12:33 pm

I like Vensim PLE and I am preparing a paper for a teachers’ inservice course using PLE and other tools.
I wonder if it is possible to model the Gumowski-Mira-attractor given by the following system:

There is a function g(x) defined (specially g(x)=-0.9x 2*(1+0.9)*x^2/(1+x^2)

Then

xn = yo + a(1-b*yo^2)*yo + g(xo)
yn = -xo + g(xn)

x0 = y0 = 0.1 (initial values)
xn = x new; xo = x old

yn = y new; yo = y old

Many thanks in advance

Josef Boehm

Administrator
Super Administrator
Posts: 2845
Joined: Wed Mar 05, 2003 3:10 am

Re: Gumowski-Mira-Attractor

Postby Administrator » Mon Aug 29, 2011 1:27 pm

I don't see why you cannot do it. To get the old values of x, y and G, use DELAY FIXED with a delay time of one time step.
Advice to posters seeking help (it really helps us to help you)
http://www.ventanasystems.co.uk/forum/viewtopic.php?f=2&t=4391

Units are important!
http://www.bbc.co.uk/news/magazine-27509559

tomfid
Administrator
Posts: 2041
Joined: Wed May 24, 2006 4:54 am
Location: Montana, USA
Contact:

Re: Gumowski-Mira-Attractor

Postby tomfid » Mon Aug 29, 2011 7:16 pm

Here's an example, which also shows in general how to handle discrete time math.

gumowski mira.mdl
(2.98 KiB) Downloaded 205 times


The behavior is really amazing.

Tom
/*
Advice to posters (it really helps us to help you)
viewtopic.php?f=2&t=4391

Blog: http://blog.metasd.com
Model library: http://models.metasd.com
Bookmarks: http://delicious.com/tomfid/SystemDynamics

*/

josef
Junior Member
Posts: 2
Joined: Mon Aug 29, 2011 12:29 pm

Re: Gumowski-Mira-Attractor

Postby josef » Wed Aug 31, 2011 6:13 am

Hi Tom,

many thanks. This helps a lot.
You are right the graphs are very amazing.
Thanks also for the link,

Regards
Josef

gwr
Senior Member
Posts: 205
Joined: Sun Oct 04, 2009 8:40 pm
Vensim version: DSS
Location: Germany
Contact:

Re: Gumowski-Mira-Attractor

Postby gwr » Wed Aug 31, 2011 8:27 pm

Hi Tom,

I do not see any reason for using a DELAY FIXED here. You can handle discrete systems - e.g. difference equations' systems - using the regular System Dynamics notation. What you need to do is:

1. Set the time step to 1.

2. Convert the Stock Equations into difference equations, e.g. Delta X = X(t+1) - X(t) = dx/dt with dt = 1.

3. In the case of the Gumowski-Mira-System a bit of algebra will convince you that no delays are needed.

See the model enclosed.

Cheers,

Guido
Attachments
gumowski mira.mdl
(5.84 KiB) Downloaded 146 times

tomfid
Administrator
Posts: 2041
Joined: Wed May 24, 2006 4:54 am
Location: Montana, USA
Contact:

Re: Gumowski-Mira-Attractor

Postby tomfid » Thu Sep 01, 2011 5:37 pm

Guido,

You're entirely correct. The INTEG notation that you used is probably the nicest way to do this for formal correspondence with discrete derivative notation. However, as long as TIME STEP=1 and Euler or Diff integration is used, the results will be identical whether INTEG, SMOOTH, or DELAY FIXED is used. (Diff integration is the same as Euler, but the rates and levels are stored differently, which makes it easier to see the initial values of the rates - sometimes useful for discrete systems like this.)

Tom
/*
Advice to posters (it really helps us to help you)
viewtopic.php?f=2&t=4391

Blog: http://blog.metasd.com
Model library: http://models.metasd.com
Bookmarks: http://delicious.com/tomfid/SystemDynamics

*/


Return to “Vensim”

Who is online

Users browsing this forum: No registered users and 3 guests