peter_de_jong_lyap.h

/*-----------------------------------------------------------------

ピーター・デ・ジョング写像のリアプノフスペクトラムの計算

calculating the Lyapunov spectrum of the Peter de Jong map

x[t+1] = sin( A * y[t] ) - cos( B * x[t] )

y[t+1] = sin( C * x[t] ) - cos( D * y[t] )

peter_de_jong_lyap.h

-----------------------------------------------------------------*/

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#define ITERATION 1000000

#define SKIP 10000

// parameters

// case 1

#define A -2.24

#define B  0.43

#define C -0.65

#define D -2.43

// case 2

/*

#define A  2.01

#define B -2.53

#define C  1.61

#define D -0.33

*/

// case 3

/*

#define A -2.0

#define B -2.0

#define C -1.2

#define D  2.0

*/

// 初期値

// initial values

#define X0 0.111

#define Y0 0.222

// 変数の数

// number of variables

#define DIM 2

/*------------------------------------------------------------------------*/

extern void Peter_de_Jong_eqs( double (*x), double (*y), double u[DIM][DIM] );

extern void gram_schmidt_orth( double u[DIM][DIM], double e[DIM][DIM] );

extern double calc_lyap_dim( double lambda[DIM] );

extern double log2( double x );