大侠帮忙看看这个方程的poincare map 是怎么画的啊？

sdlmx

方程是function si=sirs(t,x,b1,b3,p,c)
b1=0;
b3=5;
p=3;
c=2.5981;
si=[x(2);-c^2+x(1)^2+0.00002*(b1*x(2)+((2*p-2)/p)*x(1)*x(2)+((p^3*b3)/(2-2*p))*sin(2*t))];

我为什么不能画出下面给出的图呢？

我的程序是：
clc
clear
b1=0;
b3=5;
p=3;
c=2.5981;
options=odeset('RelTol',1e-7);
tt=0.01;
[t,x]=ode45(@sirs,[0:tt/100:5*tt],[-5.19,0],[],b1,b3,p,c);
plot(x(300:end,1),x(300:end,2))
%pojialai
hold on
for i=100:10:500
plot(x(i,1),x(i,2),'r.')
end

图形是：


[ 本帖最后由 sdlmx 于 2007-7-27 15:26 编辑 ]

无水1324

clc
clear
b1=0;
b3=5;
p=3;
c=2.5981;
options=odeset('RelTol',1e-7);
T=pi;
tt=100*T;
[t,x]=ode45(@sirs,[0:T/100:tt],[-5.19,0],[],b1,b3,p,c);
plot(x(8000:100:end,1),x(8000:100:end,2),'k.')

sdlmx

提示有问题啊
Warning: Failure at t=1.027494e+001.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.842171e-014) at time t.
> In ode45 at 355

我是按照文章上给的参数值来画的啊

[ 本帖最后由 sdlmx 于 2007-7-27 15:39 编辑 ]

咕噜噜

你确定你给的系统和图没有错？我也得不到你的图啊

无水1324

检查你的参数吧