资源描述:
《matlab编写Lyapunov指数计算程序汇总》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、matlab编写的Lyapunov指数计算程序汇总申明:以下各程序为个人在网络上收集的Lyapunov指数计算程序,未经过验证,不保证程序的正确性和计算结果的正确性,请大家见谅,也欢迎大家探讨!目录一、计算连续方程Lyapunov指数的程序二、wolf方法计算李雅普诺夫指数三、小数据量法计算Lyapunov指数的Matlab程序四、C-C方法计算时间延迟和嵌入维数计算Lyapunov指数计算关联维数混沌时间序列预测五、计算各种混沌系统李雅普洛夫指数的MATLAB源程序六、二维Henon映射的Lyapunov指数的计算
2、一、计算连续方程Lyapunov指数的程序,比较好用的其中给出了两个例子:计算Rossler吸引子的Lyapunov指数计算超混沌Rossler吸引子的Lyapunov指数http://www.pudn.com/downloads39/sourcecode/math/detail132231.html基于Volterra滤波器混沌时间序列多步预测----------------------------------------------参考文献:1、张家树.混沌时间序列的Volterra自适应预测.物理学报.2000
3、.032、ScottC.Douglas,TeresaH.-Y.Meng,NormalizedDataNonlinearitiesforLMSAdaptation.IEEETrans.Sign.Proc.Vol.421994----------------------------------------------文件说明:1、original_MultiStepPred_main.m程序主文件,直接运行此文件即可2、original_train.m训练函数3、original_test.m测试函数4、LorenzDa
4、ta.dll产生Lorenz离散序列5、normalize_1.m归一化6、PhaSpaRecon.m相空间重构7、PhaSpa2VoltCoef.m构造Volterra自适应FIR滤波器的输入信号矢量Un8、TrainTestSample_2.m将特征矩阵前train_num个为训练样本,其余为测试样本9、FIR_NLMS.dllNLMS自适应算法http://www.pudn.com/downloads45/sourcecode/math/detail150408.htmlrecnstitution重构相空间,在
5、非线性系统分析中有重要的作用MATLAB程序如下:1function[Texp,Lexp]=lyapunov(n,tstart,stept,tend,ystart,ioutp);2globalDS;3globalP;4globalcalculation_progressfirst_call;5globaldriver_window;6globalTRJ_buferTime_buferbufer_i;78%9%LyapunovexponentcalcullationforODE-system.10%11%Thealog
6、rithmemployedinthism-filefordeterminingLyapunov12%exponentswasproposedin13%14%A.Wolf,J.B.Swift,H.L.Swinney,andJ.A.Vastano,15%"DeterminingLyapunovExponentsfromaTimeSeries,"PhysicaD,16%Vol.16,pp.285-317,1985.17%18%ForintegratingODEsystemcanbeusedanyMATLABODE-suit
7、emethods.19%ThisfunctionisapartofMATDSprogram-toolboxfordynamicalsysteminvestigation20%See:http://www.math.rsu.ru/mexmat/kvm/matds/21%22%Inputparameters:23%n-numberofequation24%rhs_ext_fcn-handleoffunctionwithrighthandsideofextendedODE-system.25%Thisfunctionmus
8、tincludeRHSofODE-systemcoupledwith26%variationalequation(nitemsoflinearizedsystems,seeExample).27%fcn_integrator-handleofODEintegratorfunction,forexample:@ode4528%tstart-sta