自己改写的matlab程序,画Lorenz系统的最大L指数图,判断混沌 运行无误,结果精确 附结果图!

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

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

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

    Lorenz 系统 文档分两个文件方程m文件和计算L指数m文件分开写,复制粘贴即可运行matlab2012a,改写方程文件和参数即可算自己的系统,其中最大L指数用的是经典的柏内庭(G.Benettin ...

  2. matlab绘制庞加莱截面_matlab庞加莱截面法画Lorenz系统分岔图(附图).doc

    利用庞加莱截面法 画的Lorenz 系统的分岔图,复制改成其他系统即可运行,Matlab12a可以运行,附分岔图,见下页! function Lorenz_bifur_r Z=[]; for r=li ...

  3. matlab绘制庞加莱截面_matlab庞加莱截面法画Lorenz系统分岔图(附图)

    利用庞加莱截面法 画的 Lorenz 系统的分岔图,复制改成其他系统即可运行, Matlab12a 可 以运行,附分岔图,见下页! function Lorenz_bifur_r Z=[]; for  ...

  4. 无法计算机您系统的体验指数,win7无法计算机您系统的windows体验指数如何解决...

    在win7系统中,自带有Windows体验指数测量功能,能够通过对硬件和软件进行检测,最后得出的分数可以反映计算机的运行性能,但是有不少win7旗舰版64位系统用户在进行分数评估的时候却弹出无法计算机 ...

  5. matlab 动态优化,基于Matlab的测控系统动态性能优化与仿真

    随着测试技术的发展,人们采用传感器测控系统的动态性能指标来表征系统性能.描述传感器的主要动态性能指标是工作频带,系统的动态性能研究的重要一步是在辨识出合适的模型结构和模型参数的基础上,根据现有的工作频 ...

  6. matlab中怎么画冲激函数,matlab怎么画冲激函数波形,这些知识你不一定知道

    在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学习matlab的时候编写一些程序是非常头疼的.那么,matlab怎么画冲激函数波形呢?下面让我来给大家介绍下吧,希 ...

  7. Lorenz系统Lyapunov指数图

    Lorenz系统Lyapunov指数图 Lyapunov指数是定量判断Lorenz系统混沌状态的重要指标,最大Lyapunov指数大于零则判定该系统是混沌系统. 结果图:

  8. matlab怎么画周期性竖直线,matlab怎么画阶跃函数波形,这些经验不可多得

    在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学习matlab的时候编写一些程序难度确实比较大,因为初学者对于很多函数的使用都不是很熟悉.为了帮助初学者学习mat ...

  9. matlab表示函数波形,matlab如何画冲激函数波形 matlab画冲激函数波形的方法

    今天给大家带来matlab如何画冲激函数波形,matlab画冲激函数波形的方法,让您轻松解决问题.        在学习理工科的过程中, matlab软件的使用是一个不可或缺的过程.对于初学者来说,学 ...

  10. matlab怎么绘制零极点,matlab中画系统零极点的方法

    写论文的时候由于需要画出系统的零极点图.但是之前不知道怎么用matlab画,今天研究了一下,拿出来和大家共享.所用到的matlab函数为zplane,matlab给出的解释如下: ZPLANE Z-p ...

最新文章

  1. 疫情之下,武汉女生在家中答辩,获得国外博士学位!
  2. NB-Iot烟感01:烟感探测器原理和规格
  3. Kotlin 简化Fragment使用的扩展方法
  4. 家里wifi网速越来越慢_wifi太慢怎么办?简单的操作几步就能让网速飞快
  5. blob转file对象_JavaScript Blob 对象解析
  6. 机器学习中的数学(二)--梯度下降法
  7. 敏捷世界里中层经理的角色
  8. java网页版斗地主_JAVA控制台版斗地主
  9. npn三种波形失真_RF测试笔记:三阶交调失真概述及测试
  10. 如何使用SPSS列联表分析,查看变量间的相关关系(上)
  11. 『交换机连交换机』交换机多种化连接方式
  12. 关于学习管理系统 LMS
  13. 丽思卡尔顿游轮正式加入万豪旅享家旅行计划
  14. 如何用计算机求极限,计算器的极限_500字
  15. python安装配置教程win10_win10怎么安装python_win10安装python的图文步骤
  16. Java之日志打印占位符
  17. VMware15安装
  18. CubeFS存储技术揭密(1) — 纠删码引擎系统设计
  19. 商城系统运营模式有哪些?有什么优缺点?
  20. 关于金融风险数据的ETL到分析

热门文章

  1. vue 下载pdf文件
  2. 主流的推荐系统算法总结
  3. 引起1月12日WIN10 Flash停用原因
  4. 在jsp代码添加背景图片
  5. 生意参谋活动复盘分析
  6. Clark变换和Park变换在三相系统和单相系统中的应用
  7. UCOSII系统移植详解
  8. 【2022PMP最全总结-敏捷专题课】
  9. android自定义图标下载,Android使用IconFont矢量图标库
  10. 安卓干货——安卓NFC通信