lorenz画图matlab,Matlab画Lorenz系统的最大李雅普诺夫指数图

Lorenz 系统

文档分两个文件方程m文件和计算L指数m文件分开写，复制粘贴即可运行matlab2012a，改写方程文件和参数即可算自己的系统，其中最大L指数用的是经典的柏内庭(G.Benettin)计算方法，准确快速无误！附计算结果图！！

方程m文件：

function dX = Loren(t,X)

global a; % 变量不放入参数表中

global b;

global c;

x=X(1); y=X(2); z=X(3);

% Y的三个列向量为相互正交的单位向量

% 输出向量的初始化

dX = zeros(6,1);

% Lorenz吸引子

dX(1)=a*(y-x);

dX(2)=x*(b-z)-y;

dX(3)=x*y-c*z;

end

计算最大L指数文件

Z=[];

global a;

global b;

global c;

a=10;

c=8/3;

d0=1e-7;

for b=linspace(0,500,500)

lsum=0;

x=1;y=1;z=1;

x1=1;y1=1;z1=1+d0;

for i=1:100

[T1,Y1]=ode45('Loren',1,[x;y;z;16;b;4]);

[T2,Y2]=ode45('Loren',1,[x1;y1;z1;16;b;4]);

n1=length(Y1);n2=length(Y2);

x=Y1(n1,1);y=Y1(n1,2);z=Y1(n1,3);

x1=Y2(n2,1);y1=Y2(n2,2);z1=Y2(n2,3);

d1=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2);

x1=x+(d0/d1)*(x1-x);

y1=y+(d0/d1)*(y1-y);

z1=z+(d0/d1)*(z1-z);

if i>50

lorenz画图matlab,Matlab画Lorenz系统的最大李雅普诺夫指数图相关推荐

matlab 动态优化,基于Matlab的测控系统动态性能优化与仿真
随着测试技术的发展,人们采用传感器测控系统的动态性能指标来表征系统性能.描述传感器的主要动态性能指标是工作频带,系统的动态性能研究的重要一步是在辨识出合适的模型结构和模型参数的基础上,根据现有的工作频 ...
matlab中怎么画冲激函数,matlab怎么画冲激函数波形，这些知识你不一定知道
在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学习matlab的时候编写一些程序是非常头疼的.那么,matlab怎么画冲激函数波形呢?下面让我来给大家介绍下吧,希 ...
matlab怎么画周期性竖直线,matlab怎么画阶跃函数波形，这些经验不可多得
在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学习matlab的时候编写一些程序难度确实比较大,因为初学者对于很多函数的使用都不是很熟悉.为了帮助初学者学习mat ...
matlab表示函数波形,matlab如何画冲激函数波形 matlab画冲激函数波形的方法
今天给大家带来matlab如何画冲激函数波形,matlab画冲激函数波形的方法,让您轻松解决问题. 在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学 ...
matlab绘制庞加莱截面_matlab庞加莱截面法画Lorenz系统分岔图(附图).doc
利用庞加莱截面法画的Lorenz 系统的分岔图,复制改成其他系统即可运行,Matlab12a可以运行,附分岔图,见下页! function Lorenz_bifur_r Z=[]; for r=li ...
matlab绘制庞加莱截面_matlab庞加莱截面法画Lorenz系统分岔图(附图)
利用庞加莱截面法画的 Lorenz 系统的分岔图,复制改成其他系统即可运行, Matlab12a 可以运行,附分岔图,见下页! function Lorenz_bifur_r Z=[]; for ...
matlab怎么绘制零极点,matlab中画系统零极点的方法
写论文的时候由于需要画出系统的零极点图.但是之前不知道怎么用matlab画,今天研究了一下,拿出来和大家共享.所用到的matlab函数为zplane,matlab给出的解释如下: ZPLANE Z-p ...
matlab提取数据画图,matlab提取excel表格数据画图-如何将Excel中的数据导入MATLAB并在MATLAB中画出图......
在matlab上如何导入excel表格然后画图工具:matlab 2018b 1.打开matlab,点击主页下面导入数据,可以导入excel数据,在此将自己命名的huitushuju文件导入: 2. ...
【matlab】全面总结:MATLAB如何画出漂亮的图
在绘图命令中应该如何确定横坐标和纵坐标的标度尺寸(自己定范围) 如何在一个坐标中绘制多条曲线,并用不同的颜色区别. 希望得到答案,谢谢!!!!!!!!!!! help plot help axis ...

lorenz画图matlab,Matlab画Lorenz系统的最大李雅普诺夫指数图

lorenz画图matlab,Matlab画Lorenz系统的最大李雅普诺夫指数图相关推荐

最新文章

热门文章