http://www.agnld.uni-potsdam.de/~marwan/matlab-tutorials/html/tisean.html

Using the TISEAN package within Matlab

Contents

For the study of nonlinear systems and for nonlinear time series analysis we will need some programmes or tools. Some of the algorithms are available as Matlab functions or can be simply programmed for Matlab. However, some other are rather complex or are not yet available as Matlab functions. But most of them are available in the TISEAN package, which is a collection of small commandline programms for nonlinear data analysis. We need to install TISEAN on our computer for some tasks. This software package as well as the installation instructions can be found at http://www.mpipks-dresden.mpg.de/~tisean/TISEAN_2.1/index.html

Matlab offers a way to use systems commands and, hence, the commandline programms of TISEAN. At first let us define a variable with the path of the TISEAN installation

tiseanPath = 'C:\Programme\MATLAB6p5\work\Tisean\';
tiseanPath = '/usr/users4/marwan/1.work/4.TeXte/matlab-tutorials/';
tiseanPath = '';

Now let us use the TISEAN command henon in order to compute the Henon map. In Matlab, a system call is performed using the command system

a = 1.4;
system([tiseanPath,'henon -B0.3 -A',num2str(a),' -l1000 -o']);
x = load('henon.dat');

plot(x(:,1),x(:,2),'.')

TISEAN 2.1 (C) R. Hegger, H. Kantz, T. Schreiber

henon: Henon map
opened henon.dat for output

The argument of the system command is a string with the same content what we would type on the systems commandline

[tiseanPath,'henon -B0.3 -A',num2str(a),' -l1000 -o']

ans =

henon -B0.3 -A1.4 -l1000 -o

We can commit arguments to the TISEAN command by converting numerical values to strings with num2str. For example, to plot the bifurcation diagram, we run through a range of values of a

x = [];
for a = 0:.02:2
    system([tiseanPath,'henon -B0 -A',num2str(a),' -l200 -o -V0']);
    henonData = load('henon.dat');
    x = [x, henonData(:,1)];
end
plot(0:.02:2, x', 'k.')

Let us consider again the Henon system (but now computed using Matlab)

a = 2; b = 0;
x(1,1) = .91; x(1,2) = 0;
for i = 2:10001
    x(i,1) = 1 - a * x(i-1,1)^2 + b * x(i-1,2);
    x(i,2) = x(i-1,1);
end
x(1,:) = [];
save henon.dat x -ascii -tabs

% In order to compute the Lyapunov exponent using TISEAN we
% call

system([tiseanPath,'lyap_r -s20 -o lyap.dat henon.dat']);

l = load('lyap.dat');
plot(l(:,1),l(:,2))

TISEAN 2.1 (C) R. Hegger, H. Kantz, T. Schreiber

lyap_r: Estimates the maximal Lyapunov exponent; Rosenstein et al.

Using henon.dat as datafile, reading column 1
Line in file too long. Increasing input size
Use 10000 lines.
Opened lyap.dat for writing
epsilon: 2.000000e-03 already found: 9979

Note that some TISEAN programmes add heading lines to the output files and that the comment mark is # instead of %. This requires that we change the # sign to the % sign in order to import the file into Matlab.

Exercise

Use the TISEAN command poincare and create a Poincare map of the Henon map (10000 data points). Plot it within Matlab.

Published with MATLAB® 7.3