LyapSpec.zip:计算全谱 Lyapunov 指数或 k 个第一指数-matlab开发

该方法在第 11 届最优控制、动态博弈和阿姆斯特丹的非线性动力学，2010 年，Anton O. Беляков 函数 [LE,trJ,x] = LyapSpec(fun,T,x0,r,k) % 在二阶求解器的每个步骤中使用 Gram-Schmidt 归一化计算 Lyapunov 指数 (LE)，包括用于线性化系统的 Heun 方法和用于状态方程的 Leapfrog (Verlet) 方法。 % fun 是动态系统 (n) 及其雅可比矩阵 (nxn) 的函数% T 是时间间隔：T(1) 是开始时间，T(2) 是结束时间， % x0 是初始状态变量 (n) 的向量， % r 是采样率% k 是要计算的 LE 数，(1) % 回报： % LE 是李雅普诺夫指数 (k) 的向量， % trJ 是所有 LE (1) 的总和， % x 是最终状态变量 (n)。 %函数fun示例： 函数 [dz