matlab庞加莱截面法画Lorenz系统分岔图附图Word格式文档下载.doc
- 文档编号:1449896
- 上传时间:2023-04-30
- 格式:DOC
- 页数:2
- 大小:34.50KB
matlab庞加莱截面法画Lorenz系统分岔图附图Word格式文档下载.doc
《matlab庞加莱截面法画Lorenz系统分岔图附图Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《matlab庞加莱截面法画Lorenz系统分岔图附图Word格式文档下载.doc(2页珍藏版)》请在冰点文库上搜索。
%舍弃前面迭带的结果,用后面的结果画图
[T,Y]=ode45(@Lorenz,[0,1],[1;
1;
16;
r;
4]);
[T,Y]=ode45(@Lorenz,[0,50],Y(length(Y),:
));
Y(:
1)=Y(:
2)-Y(:
1);
%对计算结果进行判断,如果点满足x=y,则取点
fork=2:
length(Y)
f=k-1;
ifY(k,1)<
0
ifY(f,1)>
y=Y(k,2)-Y(k,1)*(Y(f,2)-Y(k,2))/(Y(f,1)-Y(k,1));
Z=[Zr+abs(y)*i];
end
else
ifY(f,1)<
0
end
end
plot(Z,'
.'
'
markersize'
1)
title('
Lorenz映射分岔图'
)
xlabel('
r'
),ylabel('
|y|wherex=y'
functiondy=Lorenz(t,y)
%Lorenz系统
%系统微分方程:
%dx/dt=-a(x-y)
%dy/dt=x(r-z)-y
%dz/dt=xy-bz
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;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 庞加莱 截面 Lorenz 系统 分岔 附图