Lorenz映射的分岔图生成方法

从Lorenz混沌系统的分叉图更能清晰地观察到其混沌特性,系统分叉图如下图所示。
分叉图:

代码实现:
1、Lorenz数学表达

function dy = Lorenz(t,y)
% Lorenzϵͳ
% ϵͳ΢ؖ׽Ԍú
%        dx/dt = -a(x-y)
%        dy/dt = x(c-z)-y
%        dz/dt = xy-bz
%    a=y(4)
%    r=y(5)
%    b=y(6)
dy=zeros(6,1);
dy(1)=-y(4)*(y(1)-y(2));
dy(2)=y(1)*(y(5)-y(3))-y(2);
dy(3)=y(1)*y(2)-y(6)*y(3);
dy(4)=0;
dy(5)=0;
dy(6)=0;

2、绘制分叉图

Z=[];
for r=linspace(1,500,500);% 舍弃前面迭带的结果,用后面的结果画图 [T,Y]=ode45('Lorenz',1,[1;1;1;10;r;8/3]);  [T,Y]=ode45('Lorenz',50,Y(length(Y),:)); Y(:,1)=Y(:,2)-Y(:,1); % 对计算结果进行判断,如果点满足x=y,则取点for k=2:length(Y) f=k-1; if Y(k,1)<0   if Y(f,1)>0y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));Z=[Z ,r+abs(y)*i];end  else     if Y(f,1)<0  y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));   Z=[Z ,r+abs(y)*i];   end  end end
end
plot(Z,'.','markersize',1)
title('Lorenz映射分岔图')
xlabel('c'),ylabel('|y| where x=y')

注:仿真代码为matlab代码。

lorenz映射的分岔图生成方法相关推荐

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

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

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

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

  3. 怎么用matlab画分叉图,混沌------分岔图绘制不同方法的总结、比较(转)

    经过近期的研究发现,目前对于系统单参数分岔图的计算共有以下的几种方法: 1)最大值法 即对系统微分方程(组)进行求解,对求解的结果用getmax函数进行取点,并绘图. 2)Poincare截面法 对系 ...

  4. 非线性可视化(5)非线性系统的分岔图

    在前面 非线性可视化(3)混沌系统 这一篇文章中,介绍了一个系统因为某个常数的改变,从而导致整个系统发生变化的例子.比如Duffing系统,随着阻尼d的增大,系统由混沌变为倍周期,又变为周期运动.想要 ...

  5. henon映射分岔图 matlab程序,分岔图做法1.doc

    分岔图做法1 > > 混沌研究总结篇------一.分岔图(1.Chen系统)先打个提纲,这几天把自己混沌相关知识研究学习内容总结一下. ?????首先简绍几个基本概念: 一.自治系统?? ...

  6. Logistic混沌映射分布图与分岔图附Matlab代码

    目录 一.Logistic混沌映射 二. Logistic混沌映射分岔图及代码 三.Logistic混沌映射分布图与直方图 一.Logistic混沌映射 二. Logistic混沌映射分岔图及代码 c ...

  7. matlab画虫口模型分岔图,Matlab画图片碰磨转子系统中的分岔图,poincare映射图???...

    图片中是转子碰磨故障的动力学模型,应用龙格库塔法求解,在画分岔图和poincare截面图,画出来不对,不知道为什么,下面是我自己编的程序,请大侠们帮忙看看我错在哪里了? global omega;  ...

  8. 分岔现象matlab,Matlab——混沌分岔图

    经过近期的研究发现,目前对于系统单参数分岔图的计算共有以下的几种方法: 1)最大值法 即对系统微分方程(组)进行求解,对求解的结果用getmax函数进行取点,并绘图. 2)Poincare截面法 对系 ...

  9. matlab画虫口模型分岔图,[转载]Matlab——混沌分岔图

    经过近期的研究发现,目前对于系统单参数分岔图的计算共有以下的几种方法: 1)最大值法 即对系统微分方程(组)进行求解,对求解的结果用getmax函数进行取点,并绘图. 2)Poincare截面法 对系 ...

最新文章

  1. K-BERT | 基于知识图谱的语言表示模型
  2. CentOS 7系统安装配置图文详解
  3. EXCEL 找出重复的记录 COUNTIF
  4. Linux Shell脚本攻略:shell中各种括号()、(())、[]、[[]]、{}的作用
  5. SQL Server中的Union和Union All语句之间的差异及其性能
  6. Android studio中git密码记住的问题
  7. python 获取文件列表_Python3 - 获取文件夹中的文件列表
  8. php 读取 stl 文件,科学网—[转载]【源码】读取ASCII STL文件的STLReader - 刘春静的博文...
  9. Android(1)—Mono For Android 环境搭建及破解
  10. python维度变换_Python NumPy用法
  11. 腾讯内部转岗_别了,腾讯微博!
  12. Markdown支持的语言
  13. 一个很有趣的游戏(看谁的名字打架厉害)
  14. 游戏窗口化工具_仙剑奇侠传16珍藏版大合集(含各个版本和工具)
  15. 大学python课本_Python大学实用教程(大学计算机规划教材)/数据工程师系列
  16. 深入理解操作系统实验——bomb lab(作弊方法2)
  17. LoRa在智能家居当中的物联网应用
  18. 创建微信小程序日期和时间的组件
  19. 在计算机网络中ln代表的是,数学中e和ln的关系?
  20. 计算机组成原理习题 第七章 外围设备

热门文章

  1. 拼多多店铺物流类问题处理方法!
  2. 【Go语言版GTK,安装以及环境搭建】
  3. 程序人生----- hello的生命历程
  4. 为何Word 向表格里插入图片时显示不全呢?
  5. 8080端口被占用时,我们该怎么办
  6. LTH7五脚芯片的完整方案图
  7. 基于java斗地主_基于Java的网络版斗地主
  8. java图片文字识别功能的实现
  9. linux系统盘15G,装UOS/Deepin 20选择全盘安装的系统分区仅有15GB,附扩容经验
  10. 计算机毕设旅游网的开发,计算机毕业设计分享Jsp旅游直通车服务系统