Matlab求解混沌系统最大李雅普诺夫指数

李雅普诺夫指数(Lyapunov)是一个较为典型的判断一个系统是否具有混沌特性以及混沌的程度分析方法。

李雅普诺夫指数：在相空间中初始时无限接近的两个轨道，随着时间的不断推移按指数收敛或发散的平均变化率，它可以定量描述混沌系统在局部范围里系统轨道间的分离程度。

假设在一维动力系统 Xn+1=F(Xn)，为了表示从整体上看相邻两状态分离的情况，必须对时间(或迭代次数)取平均，不妨设平均每次迭代所引起的指数分离中的指数为λ，原来相距为ε的两点经过n次迭代后[相距为]：

如果λ>0，运动轨道的局部不稳定，相邻点的轨道的终按指数方式分离，则在此作用下反复折叠，形成混沌吸引子。

在n维连续动力学系统中，将一个无穷小n维的球作为系统的初始条件，随着动力系统的演化向相空间的各个方向作伸展或收缩，球将变为椭球，将椭球的所有主轴按其长度顺序徘列，那么第i个李雅普诺夫指数根据第i个主轴的长度pi(t)的增加速率定义为:

将这n 个李雅普诺夫指数按照从大到小进行排序，得到李雅普诺夫指数谱：

λ1>=λ2>=λ3>=…>=λn。对于混沌系统，必须有一个正的李雅普诺夫数，所以通常采用计算最大李雅普诺夫指数的方法来判断系统是否是混沌的。

以sine函数为例，求其最大李指数

子函数：

function newx=sine_f1(x)newx=x(2)*sin(pi*x(1))/4;

主函数：

clc;cleard0=1e-7;Z=[]; r=0.5;for u=linspace(0,8,800)   le=0;   lsum=0;   x=[r;u];   x1=[r+d0;u];   for k=1:800       x=sine_f1(x);       x1=sine_f1(x1);       d1=sqrt((x(1)-x1(1))^2);%1范数       x1=x+(d0/d1)*(x1-x);       if k>100           lsum=lsum+log(d1/d0);       end       x=[x(1);u];       x1=[x1(1);u];   end   le=lsum/(k-100);   Z=[Z,u+le*1i]; endplot(Z,'-') xlabel('k','fontsize',12,'FontAngle','italic');ylabel('Lyapunov exponent','fontsize',12,'FontAngle','italic');grid on;

Matlab求解混沌系统最大李雅普诺夫指数相关推荐

matlab混沌指数的计算,matlab实现混沌系统最大李雅普诺夫指数
李雅普诺夫指数(Lyapunov)是一个较为典型的判断一个系统是否具有混沌特性以及混沌的程度分析方法. 李指数:在相空间中初始时无限接近的两个轨道,随着时间的不断推移按指数收敛或发散的平均变化率,它可 ...
matlab Lorenz混沌系统求解
1.内容简介 matlab Lorenz混沌系统求解 649-可以交流.咨询.答疑 2.内容说明略 3.仿真分析 clc close all clear sigma = 10; alpha = 30 ...
金融系统混沌分形维matlab,Matlab实现混沌系统的控制
基于MATLAB 的各类混沌系统的计算机模拟混沌是非线性系统所独有且广泛存在的一种非周期运动形式, 其覆盖面涉及到自然科学和社会科学的几乎每一个分支.1972年12月29日,美国麻省理工学院教授.混 ...
matlab系统稳态误差终值,Matlab求解给定系统稳态误差
对于如图线性系统,应用拉式变换终值定理,可以很容易求出系统给定稳态误差 ess=lim e(t) (t-无穷大)=limsE(s)(s-0)=limS(R(S)-C(S))= limS(R(S)-W( ...
Matlab求解给定系统稳态误差
对于如图线性系统,应用拉式变换终值定理,可以很容易求出系统给定稳态误差 ess=lim e(t) (t-无穷大)=limsE(s)(s-0)=limS(R(S)-C(S))= limS(R(S)-W( ...
MATLAB编程之混沌系统
一概述简单来说,混沌科学中,同一个动力学方程或者是方程组,给出的初始值相差甚微,长时间后,系统的运动轨迹却有很大差别:混沌系统广泛应用在加密等领域. 二 Matlab仿真混沌系统以che ...
基于混沌系统的文本加密算法研究（二）——经典混沌映射
基于混沌系统的文本加密算法研究(二)--经典混沌映射前言一.一维Logistic混沌映射二.二维Henon混沌映射三.三维Lorenz连续混沌映射总结代码 1.Logistic映射 2.H ...
信号与系统陈后金matlab,信号与系统(陈后金)_MATLAB.ppt
信号与系统(陈后金)_MATLAB 信号与系统 Signals and Systems 利用MATLAB进行信号与系统分析 MATLAB简介信号的MATLAB表示利用MATLAB进行系统的时域分析 ...
李雅普诺夫指数对应的特征方向
洛伦兹系统的李雅普诺夫指数的符号为 (+,0,−)(+,0,-)(+,0,−) Any continuous time-dependent dynamical system without a fix ...

Matlab求解混沌系统最大李雅普诺夫指数

Matlab求解混沌系统最大李雅普诺夫指数相关推荐

最新文章

热门文章