function dx=zhengzefun(t,x)
global w;
global f;
global j;
w=1;
f=0.5;
j=4.5;
dx=zeros(4,1);
%x(1)=p1 x(2)=q1 x(3)=p2 x(4)=q2
dx(1)=-w*x(2)-2*f*abs(((4*j-x(2)^2-x(1)^2)/4/j))^(1/2)*x(4)-f*x(2)^2*x(4)*abs(((4*j-x(2)^2-x(1)^2)/4/j))^(-1/2)/j;
dx(2)=w*x(1)-x(1)*x(2)*x(4)*abs(((4*j-x(2)^2-x(1)^2)/4/j))^(-1/2)/j;
dx(3)=-w*x(4)-2*f*abs(((4*j-x(2)^2-x(1)^2)/4/j))^(1/2)*x(2);
dx(4)=w*x(3);
上面这个是小弟最近在计算的非线性微分方程组,其中w f j。均为一些参数,当f=0.1时,我画出了该系统的poincare截面图,显示出是一系列的连续的闭合曲线,与理论以及别人的结果符合的很好。但是在选取f=0.5时,应该是在poincare界面上是一些不规则的点,但是发现x的取值出现了如下的情形(见附件),无法选取poincare截面了。种种时候该怎么做呐?请各位大虾指点一下,