求助？？
各位师兄师姐 帮忙看看这个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.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=;
=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 编辑 ] 放大庞加莱图下面的黑点就是下面的图
页:
[1]