xiaoqiu810818 发表于 2008-1-7 21:49

我的Poincare截面程序

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',);
c=sqrt(3);

=ode45('rossler',,[-0.2 0.1 0.1],options,c);
figure(1);
subplot(2,1,1);
plot(t,x(:,1));

% axis();
subplot(2,1,2);
plot(t,x(:,2));

% axis();
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
这是我修改别人的程序,望两位研兄指教,不胜感激.

无水1324 发表于 2008-1-7 22:28

回复 #1 xiaoqiu810818 的帖子

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*');

这个应该就是你计算Poincare截面点的程序把？

xiaoqiu810818 发表于 2008-1-9 08:45

是的无水大哥.怎么改一下呀?

咕噜噜 发表于 2008-1-9 09:05

现在你的程序要是修改需要大改，假如你做的是poincare截面
如果是poincare映射还好一点，但是也需要修改周期的选取

xiaoqiu810818 发表于 2008-1-9 17:34

院长大人,关键是怎么改呀?举个例子吧.小弟愚钝.

xiaoqiu810818 发表于 2008-1-9 17:45

我的邮箱:xiaoqiu810818@163.com.咕院长,欢迎来做客.

咕噜噜 发表于 2008-1-9 19:07

你要画poincare截面还是映射？

无水1324 发表于 2008-1-14 13:12

回复 6楼 的帖子

Xiaoqiu兄，你的问题我已经发到邮箱里面去了，你查收一下。

真的不好意思，最近忙于一写琐碎的事情，很少来论坛，来了也是匆匆看一下。不知道你的问题解决没有。
页: [1]
查看完整版本: 我的Poincare截面程序