MATLAB Software


Dynamical Systems

2-dimensional Iterated Function Systems (IFS)
ifs.m
Plots the attractor of a generic 2-dimensional (linear) IFS by chaos game.
[X,h] = ifs(A, p, N, opt1)
Some examples of the use of ifs.m are here: ifs.zip

Self-organized criticality in a model of production and inventory dynamics
bcsw.m
A raw and simple program to numerically simulate the model developed by:
P. Bak, K. Chen, J. Scheinkman and M. Woodford (1993),
Aggregate Fluctuations from Independent Sectoral Shocks: Self-Organized Criticality
in a Model of Production and Inventory Dynamics, Ricerche Economiche, 47 (1), 3-30.
Click here to download a screenshot of the program.

Endogenous preferences and relational dynamics
biavati_sandri_zarri.m
The simulation software of the model presented in
Biavati M., Sandri M., Zarri L. (2002),
Preferenze endogene e dinamiche relazionali: un modello co-evolutivo,
in Sacco P.L. e Zamagni S. (a cura di),
Complessità relazionale e comportamento economico.
Verso un nuovo paradigma di razionalità, ll Mulino, Bologna



Plot of High Dimensional Data

Andrews Diagrams
andrews.m
Plots high-dimensional data using the method proposed by Andrews (1972).
Y = ANDREWS(X,N) plots each row of the matrix X as a trigonometric function defined in the interval (-pi,pi).
N is the number of points taken into account in the interval(-pi,pi). The default value for N is 100.
Y is the matrix containing the values f(t) of the trigonometric functions.
Reference: D.F. Andrews (1972), "Plots of High-Dimensional data", Biometrics, pp. 125-136
An example of the use of andrews.m is: testandr.m



Mathematica Packages


Nonlinear dynamical systems

Numerical Calculation of Lyapunov Exponents
The Lyapunov characteristic exponents play a crucial role in the description of the behavior
of dynamical systems. They measure the average rate of divergence or convergence of orbits
starting from nearby initial points. Therefore, they can be used to analyze the stability of limit
sets and to check sensitive dependence on initial conditions, that is, the presence of chaotic
attractors.
This package shows how to use Mathematica to compute the Lyapunov spectrum of a smooth
dynamical system.
The lce.m Mathematica package (for Mathematica < 5.2)
The lce.m Mathematica package (for Mathematica >7)
A notebook with all the (updated) code used in the paper (for Mathematica 12)



Time series analysis

MELISSA: a Mathematica package for Singular Spectrum Analysis
melissa.m

Singular spectrum analysis (SSA) is a relatively recent technique for time series analysis.
The idea behind SSA was originally purposed as a data adaptive method for choosing an optimal
embedding dimension for attractor reconstruction. Later the technique was developed as a "stand
alone" time series analysis technique.
During the last decade it has been very successful and has become a standard tool in many different
scientific fields, such as climatic, meteorological, geophysical, and astronomical time series analysis.
Download the tutorial for MELISSA (in italian)

---

Updated: 12 Oct 2022
Home