论坛 › 数学与非线性科学 › 大家帮忙看看这样的poincare图？

yanzi 发表于 2009-4-9 15:16

大家帮忙看看这样的poincare图？

这是我的程序:

function xdot=myff1(t,x,flag,v)
if abs(x(1))<=0.4*10^6
      ff=w1^2*(x(1))+a1;
elseif abs(x(1))<=0.3*10^6
      ff=0;
elseif abs(x(1))>0.4*10^6
      ff=w2^2*(x(1))+a2;
elseif abs(x(1))>=0.5*10^6
      ff=w1^2*(x(1))-a3;
else
      ff=w2^2*(x(1))-a4;
end   
xdot=;
      


functionpoincareT1
global p
global v
global betta
global w1
global w2
global a1
global a2
global a3
global a4
betta=0.36;w1=197.15;w2=200.52;p=0.935;
a1=-0.011661;a2=0.01527;a3=0.011329;a4=-0.012063;v=68
x0=;
T=2*pi/v;
tspan=;                                                                                                                                                                                                                                                                                          
=ode45('myff1',tspan,x0);
m=zeros(10000,1);
n=zeros(10000,1);
for i=1:10000
    m(i,1)=x(10*i,1);
    n(i,1)=x(10*i,2);
end
plot(m,n,'k.','markersize',1);
xlabel('x','fontsize',14);
ylabel('dx/dt','fontsize',14);

绘制出的poincare图按理来说应该就形如如下图所示，



但是我会出来的图是发散的 ，不知道哪里出了问题，还望各位师兄师姐赐教!!感谢之至

[ 本帖最后由 yanzi 于 2009-4-9 15:18 编辑 ]

查看完整版本: 大家帮忙看看这样的poincare图？