Lyapunov spectrum for a user-defined dynamical system

Hello,
I have seen the nice documentation here about calculating the Lyapunov spectrum for a dynamical system defined in DynamicalSystems.jl . I have a dynamical system (which is a discretized PDE) which I time-integrate using ODEProblem and solve from DifferentialEquations.jl. Is there a way to calculate the Lyapunov spectrum using the solution generated from the latter call?

2 Likes

Hi,

I believe the intended way to do this would be to wrap your equation
into a ContinuousDynamicalSystem (see https://juliadynamics.github.io/DynamicalSystems.jl/latest/ds/general/ )
and then pass it to the lyapunovs function.

If you need more specific arguments for the solve function then have a look at the advanced docs
https://juliadynamics.github.io/DynamicalSystems.jl/latest/advanced/

1 Like

Hi ! @AhmedAlreweny
Would you mind sharing how exactly to do this? I have a custom-defined dynamical system and would like to calculate its lyapunov exponent without using the pre-defined lyapunovs function. Is there a way to do this within the DifferentialEquations.jl framework? Also I would like to do this for varying parmaters and initial conditions!