论坛 › 数学与非线性科学 › 大家看一下此Poincare图是否反应混沌运动？

octopussheng 发表于 2007-12-17 14:21

做映射图时，选点的步长是保证图正确与否的标准啊！不能随便选择的！

quiet1016 发表于 2007-12-17 22:50

去掉瞬态响应，得出，看样子是混沌啊

quiet1016 发表于 2007-12-17 22:52

晕，图片传不上来

quiet1016 发表于 2007-12-18 17:53

步长随意选取（最好选整数）

function ff
options = odeset('RelTol',1e-5,'AbsTol',);
sol= ode45(@ischao,,,options);% =
step=7.5;    %步长
x = 0:0.1:1000;
y1 = deval(sol,x,1);
y2 = deval(sol,x,2);
%plot(y1(100:step:end),y2(100:step:end),'*')   %y1 vs y2
subplot(1,2,1)
plot(y1(1:step:end-step),y1(1+step:step:end),'*');%y1(t-1) vs y1(t)
subplot(1,2,2)
plot(y2(1:step:end-step),y2(1+step:step:end),'*');%y2(t-1) vs y2(t)
function dy = ischao(t,y)
dy = zeros(2,1);    % a column vector
dy=;

quiet1016 发表于 2007-12-18 17:55

step=7

quiet1016 发表于 2007-12-18 18:00

楼主是不是把公式贴错了呀？

好像怎么都得不到楼主贴的那个POINCARE图

octopussheng 发表于 2007-12-18 20:32

呵呵，确实啊！我也没算到！

16443 发表于 2007-12-19 09:44

回复 #19 quiet1016 的帖子

subplot(1,2,1)
plot(y1(1:step:end-step),y1(1+step:step:end),'*');%y1(t-1) vs y1(t)
subplot(1,2,2)
plot(y2(1:step:end-step),y2(1+step:step:end),'*');%y2(t-1) vs y2(t)


这样做相图是什么意思啊 ？

无水1324 发表于 2007-12-19 10:14

回复 #19 quiet1016 的帖子

你这个步长选择的依据是什么？
Poincare图示根据什么原理做出来的？

16443 发表于 2007-12-19 10:22

就是么，给介绍一下1

无水1324 发表于 2007-12-19 10:24

回复 #25 16443 的帖子

根据现有的两种做Poincare方法，他这个都不能很好的解释。难道是一种新的理论算法，请quiet1016不吝赐教！

quiet1016 发表于 2007-12-19 12:24

随便取的步长

不是做Poincare，是y(t-x)与y(t)的关系，其中x＝7.5,我就是想看看楼主之前给的图究竟是怎么得到的

lzhyxcx 发表于 2007-12-20 16:31

回复 #21 quiet1016 的帖子

对不起，我的方程与第一个Poincare图不符，其中一个参数不一样。谢谢各位的讨论！

无水1324 发表于 2007-12-20 18:12

回复 #28 lzhyxcx 的帖子

不会吧。:@):@L

zhuoer 发表于 2007-12-21 15:11

建议将系统的Lyapunov指数求出来，更能准确的判断是否进入了混沌

查看完整版本: 大家看一下此Poincare图是否反应混沌运动？