声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3896|回复: 23

[分形与混沌] 无水大哥,请回答小弟的问题

[复制链接]
发表于 2008-1-7 20:33 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
附件是我画的Rossler系统的Poincare截面,问题:
1)从二)Poincare截面能否看出混沌?
2))画的Poincare截面对么?和二维相图几乎一样,为什么?
3)RosslerPoincare截面什么样?
小弟是初学者,诚心求教各位研兄.

下面是我画的Rossler系统的Poincare截面.doc

75 KB, 下载次数: 25

回复
分享到:

使用道具 举报

 楼主| 发表于 2008-1-7 20:35 | 显示全部楼层
:'( 怎么没人呀?
发表于 2008-1-7 20:57 | 显示全部楼层
你的图最好还是用插图形式直接贴出来,下载了再看确实蛮费劲的,呵呵!
发表于 2008-1-7 21:04 | 显示全部楼层
你画的图应该是有点问题,感觉可能是理解上有点偏差
把你的程序贴出来吧,大家帮你参考一下
 楼主| 发表于 2008-1-7 21:29 | 显示全部楼层
不好意思,我还不会操作
 楼主| 发表于 2008-1-7 21:34 | 显示全部楼层
我刚看了,octopussheng大哥的空间,里面有相似的讨论,发现确实,我的程序有问题.你们讨论的和我的程序相差很远.稍等,我试一下如何把自己的程序和图形发出.
 楼主| 发表于 2008-1-7 21:50 | 显示全部楼层
function xdot=rossler(t,x,flag,c)
a=0.2;   
b=0.2;
xdot=[-(x(2)+x(3));x(1)+a*x(2);  
x(3)*(x(1)-c)+b]
clear
options=odeset('RelTol',1e-10,'AbsTol',[1e-10 1e-10 1e-10 ]);
c=sqrt(3);  
[t,x]=ode45('rossler',[0:0.001:100],[-0.2 0.1 0.1],options,c);
figure(1);
subplot(2,1,1);
plot(t,x(:,1));
% axis([20 30 -0.1 0.1]);
subplot(2,1,2);
plot(t,x(:,2));
% axis([20 30 -1 1]);
figure(2);
plot(x(end-50000:end,1),x(end-50000:end,2)) figure(3);
final=fix(35*c/pi);
for i=1:final
g=(100001-7e4)+fix(2*pi*1000/c*i);
plot(x(g,1),x(g,2),'r*');  hold on
end
hold off
 楼主| 发表于 2008-1-7 21:51 | 显示全部楼层
两位研兄,请指教
 楼主| 发表于 2008-1-7 21:54 | 显示全部楼层
下面是我画的Rossler系统的Poincare截面:
一)        周期状态下的二维相图和Poincare截面:

                        二维相图

                                Pioncare截面
二)下面是混沌状态下的二维相图和Poincare截面:

                          二维吸引子

                               Pioncare截面
三)混沌状态下的二维吸引子和Pioncare截面:

                     二维相图

                   Pioncare截面
 楼主| 发表于 2008-1-7 21:59 | 显示全部楼层
octopussheng和咕噜噜研兄,请留下您的来年系方式吧,我这边快熄灯了.qq或邮箱均可.有时间小弟还要多多请教.
 楼主| 发表于 2008-1-7 22:02 | 显示全部楼层
我的邮箱: xiaoqiu810818@163.com.好郁闷这个图形我已经画了好长时间,也没理想的结果,万分感谢!小弟有礼了!
发表于 2008-1-7 22:27 | 显示全部楼层

回复 #11 xiaoqiu810818 的帖子

首先,感谢你的信任,我刚才才看到你的问题
第一个相图看出来应该是周期1解,因为只看到了一个明显的闭环,而你的Poincare好像是闭环上有间隔的点,这就有两种可能性:
1)是概周期的,Poincare组成闭环,所以你计算的时间长一点,看是不是形成闭环
2)是周期1的解,这时候可能是你的Poincare图做错了,(我现在没有时间仔细看,明天再看一下,我这了10:30准时关电的)

个人感觉你的Poincare截面取得有问题。
发表于 2008-1-8 08:35 | 显示全部楼层
个人认为你的poincare截取有问题
你的poincare截面理解上可以参考本论坛的一些帖子
 楼主| 发表于 2008-1-8 09:03 | 显示全部楼层
谢谢!等待ing.........
 楼主| 发表于 2008-1-8 10:40 | 显示全部楼层
这种自治系统截面如何选取?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-20 00:36 , Processed in 0.130501 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表