matlab编写Lyapunov指数计算程序
matlab编写的Lyapunov指数计算程序 已有 2406 次阅读 2009-12-29 08:37 |个人分类:其它|系统分类:科普集锦|关键词:李雅普诺夫指数 一、计算连续方程Lyapunov指数的程序其中给出了两个例子:计算Rossler吸引子的Lyapunov指数 计算超混沌Rossler吸引子的Lyapunov指数/downloads39/sourcecode/math/detail132231.html二、recnstitution重构相空间,在非线性系统分析中有重要的作用function [Texp,Lexp]=lyapunov(n,tstart,stept,tend,ystart,ioutp);global DS;global P;global calculation_progress first_call;global driver_window;global TRJ_bufer Time_bufer bufer_i;%% Lyapunov exponent calcullation for ODE-system.%% The alogrithm employed in this m-file for determining Lyapunov% exponents was proposed in%% A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano,% "Determining Lyapunov Exponents from a Time Series," Physica D,% Vol. 16, pp. 285-317, 1985.%% For integrating ODE system can be used any MATLAB ODE-suite methods. % This function is a part of MATDS program - toolbox for dynamical system investigation% See: http://www.math.rsu.ru/mexmat/kvm/matds/%% Input parameters:% n - number of equation% rhs_ext_fcn - handle of function with right hand side of extended ODE-system.% This function must include RHS of ODE-system coupled with % variational equation (n items of linearized systems, see Example). % fcn_integrator - handle of ODE integrator function, for example: @ode45 % tstart - start values of independent value (time t)% stept - step on t-variable for Gram-Schmidt renormalization procedure.% tend - finish value of time% ystart - start point of trajectory of ODE system.% ioutp - step of print to MATLAB main window. ioutp==0 - no print, % if ioutp>0 then each ioutp-th point will be print.%% Output parameters:% Texp - time values% Lexp - Lyapunov exponents to each time value.%% Users have to write their own ODE functions for their specified% systems and use handle of this function as rhs_ext_fcn - parameter. %% Example. Lore