分岔现象matlab,Matlab——混沌分岔图
经过近期的研究发现,目前对于系统单参数分岔图的计算共有以下的几种方法:
1)最大值法
即对系统微分方程(组)进行求解,对求解的结果用getmax函数进行取点,并绘图。
2)Poincare截面法
对系统参数的每一次取值,绘制其Poincare截面,进而得到其分岔图。
这种方法需要注意的是,自治系统的Poincare截面是选取一超平面,平面上点的分布即构成一Poincare截面,非自治系统的Poincare截面则是根据系统激励的频率进行取点并绘图。
本帖将以Lorenz系统为例,对这两种方法进行比较
首先对第二种方法进行阐述。
编程如下(matlab)
Lorenz系统:
function dy = Lorenz(t,y)
% Lorenz系统
% 系统微分方程:
%
dx/dt = -a(x-y)
%
dy/dt = x(r-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;
随r的分岔图求解程序:——按照x=y平面取截面
function Lorenz_bifur_r
Z=[];
for r=linspace(1,500,1000); %
舍弃前面迭带的结果,用后面的结果画图
[T,Y]=ode45('Lorenz',[0,1],[1;1;1;16;r;4]);
[T,Y]=ode45('Lorenz',[0,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)>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
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('r'),ylabel('|y|
where x=y')
getmax法——取最大值法
function [Xmax] = getmax(y)
a=length(y);
j=1;
for i=(a-1)/2:a
b=(y(i,1)-y(i-2,1))/2;
c=(y(i,1)+y(i-2,1))/2-y(i-1,1);
if
y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
Xmax(j)=y(i-1,1);
j=j+1;
elseif
y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)Xmax(j)=y(i-1,1)-b^2/(4*c);
j=j+1;
end
end
function Lorenz_bifur_r_getmax% 最大值法求解分岔图
clear all
t0=[0 100];%积分时间
%bifurcation
for r=linspace(1,500,1000); %r的变化精度
[t,y]=ode45('Lorenz',t0,[1;1;1;16;r;4]);
[Xmax]=getmax(y(:,1));
plot(r,Xmax,'b','markersize',1)
hold on
clear Xmax
end
分岔现象matlab,Matlab——混沌分岔图相关推荐
- hopf分岔的matlab仿真,对Hopf分岔现象及分析与的研究.doc
您所在位置:网站首页 > 海量文档  > 医药卫生 > 心理学理论/研究方法 对Hopf分岔现象及分析与的研究.doc ...
- matlab 电力系统分析 毕设,基于MATLAB的P-Q分解法电力系统潮流计算毕业设计
基于MATLAB的P-Q分解法电力系统潮流计算毕业设计 摘 要 电力系统潮流计算是研究电力系统稳定运行情况的一种重要的计算,在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用它来定量地分析比 ...
- 怎样在matlab q-q图上读出斜率,Q分解法潮流计算matlab小程序
<Q分解法潮流计算matlab小程序>由会员分享,可在线阅读,更多相关<Q分解法潮流计算matlab小程序(18页珍藏版)>请在人人文库网上搜索. 1.Q分解法潮流计算matl ...
- 毕肖普法matlab,瑞典条分法计算例题
确定安全系数 Fs T R TR Wi cos itgi Wi sin i cili 条分法是一种试算法,应选取 不同圆心位置和不同半径进行 计算,求最小的安全系数 三.例题分析 [...... 瑞典 ...
- 龙格现象 matlab,拉格朗日插值龙格现象的matlab实现
拉格朗日插值龙格现象的MATLAB实现 姓名:袁宽 学号:201430210065 专业:电气工程 题目:对于函数f(x) 1, 5 x 5进行拉格朗日插值.n 10,按等距节点求分21 x 段线性插 ...
- 基于MATLAB的混沌密码与数字图像加密应用、信息隐藏
基于MATLAB的混沌密码与数字图像加密应用.信息隐藏 实验目的:利用Matlab软件基于Logistic映射的加密算法对图像进行加密和解密,以及利用DEV C++软件.Excel表格绘制图像基于Lo ...
- matlab Lorenz混沌系统求解
1.内容简介 matlab Lorenz混沌系统求解 649-可以交流.咨询.答疑 2.内容说明 略 3.仿真分析 clc close all clear sigma = 10; alpha = 30 ...
- 基于matlab的智能天线波束方向图仿真,基于MATLAB的智能天线波束方向图仿真
第29卷第6期孝感学院学报V OL,基于M AT LA B的智能天线波束方向图仿真,汪 睿1,(1,3,摘 要:结合一种直线阵智能天线模型,关键词:智能天线,中图分类号:T N911,随着移动通信技术 ...
- Matlab plotyy画双纵坐标图实例
转载自:http://blog.sina.com.cn/s/blog_49d955150100lxoe.html Matlab plotyy画双纵坐标图实例 x = 0:0.01:20; y1 = 2 ...
- 单页双曲面 matlab,matlab绘制单叶双曲面
应用 MATLAB 绘制二次曲面图 1.用 surf 工 mesh 函数绘图 Surf 函数绘制的是三维表面图,mesh 函数绘制的是三维网格图,当二次曲面方程是标准方 程时,原方程式...... m ...
最新文章
- 云优后台提交显示parsererror_【Shopee新手开店】卖家后台管理及设置
- python简单代码画皮卡丘-实现童年宝可梦,教你用Python画一只属于自己的皮卡丘...
- python解释器哪一年_Python即Python解释器的发展史
- C# 空合并运算符 ??
- torch.tensor().permute(2,1,0)
- VTK:Utilities之ObserveError
- Spring 基于注解(annotation)的配置之@Autowired注解
- java so jnienv_JNI初步(五)jni ndk 一个.so文件依赖另一个.so文件的写法
- 6 CO配置-控制-一般控制-维护成本控制凭证的编号范围
- 【C语言】利用递归函数求n的阶乘
- pdf转换html乱码怎么办,pdf转word后乱码怎么办?
- X4扭曲字体或图形 coreldraw_cdrx4精简版下载|coreldraw x4 sp2 精简版下载增强版 15.2.3.1614 - 系统天堂...
- Python requests timeout 分析
- 台式计算机的打印机端口,打印机端口设置,高手教你如何搞定电脑打印机端口设置...
- 代码笔记源码php,读 PHP - Pimple 源码笔记(上)
- 小米路由器SSH资料攻略
- vue:单文件组件模板
- iuv_5g组网问题表
- 一台服务器,启动多个redis
- JDBC-----什么是JDBC
热门文章
- oracle-01441,ORA-01441:无法减小列长度,因为一些值过大
- WHmcs中的插件开发
- QTYX量化系统实战案例分享|箱底形态选股后潜伏介入之202209
- 魔兽作弊命令 只能单机用
- 报告 | “物联网+区块链”应用与发展白皮书(附下载)
- easypoi时间类型转化空格带字母T
- excel函数技巧:各种排名方式汇总
- Python百度云网盘搜索引擎源码及安装教程
- 微信小程序保存视频到相册wx.saveVideoToPhotosAlbum()
- python123第三周答案_python123平台第6周 函数练习