sdlmx 发表于 2007-7-27 15:21

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

方程是function si=sirs(t,x,b1,b3,p,c)
b1=0;
b3=5;
p=3;
c=2.5981;
si=;

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

我的程序是：
clc
clear
b1=0;
b3=5;
p=3;
c=2.5981;
options=odeset('RelTol',1e-7);
tt=0.01;
=ode45(@sirs,,[-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 发表于 2007-7-27 15:29

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

sdlmx 发表于 2007-7-27 15:36

提示有问题啊
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 编辑 ]

咕噜噜 发表于 2007-7-27 17:01

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

无水1324 发表于 2007-7-27 17:02

回复 #3 sdlmx 的帖子

检查你的参数吧
页: [1]
查看完整版本: 大侠帮忙看看这个方程的poincare map 是怎么画的啊？