我的Poincare截面程序

xiaoqiu810818

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

无水1324

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

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

咕噜噜

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

xiaoqiu810818

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

xiaoqiu810818

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

咕噜噜

你要画poincare截面还是映射？

无水1324

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

  真的不好意思，最近忙于一写琐碎的事情，很少来论坛，来了也是匆匆看一下。不知道你的问题解决没有。