求助？？

yanzi

各位师兄师姐 帮忙看看这个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=[x(2);-2*betta*x(2)-ff+p*sin(v*t)];


function  poincareT1
global p
global v
global betta
global w1
global w2
global a1
global a2
global a3
global a4
betta=0.6;w1=231;w2=163.7;p=0.935;
a1=-0.01608;a2=0.03752;a3=-0.05628;a4=-0.00804;
v=26.1
T=2*pi/v;
tspan=[0:0.1*T:70000];                                                                                                                                                                                                                                                                                          
[t,x]=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)；


[ 本帖最后由 yanzi 于 2009-4-29 11:11 编辑 ]

yanzi

放大庞加莱图下面的黑点就是下面的图