matlab编写的Lyapunov指数计算程序精选.docVIP

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 ioutp0 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_ex