布鲁塞尔问题的poincare映射对不对？

cam_1980

function dy=blusel(t,x)
a=0.4;
b=1.2;
arfa=0.2;
w=0.8;
dy=[a-(b+1)*x(1)+x(1)^2*x(2)+arfa*cos(w*t);
b*x(1)-x(1)^2*x(2)];


function bluselpoincare
    x0=[0.43173;2.7369];
    a=0.4;
    b=1.2;
    arfa=0.2;
    w=0.8;
    T=2*pi/w;
    tspan=[0:T/100:5000*T];
    [t,x]=ode45(@blusel,tspan,x0);
    for i=1:1:5000
        m(i)=x(100*i,1);
        n(i)=x(100*i,2);
    end
    plot(m,n,'k.','markersize',5);

帮忙看看程序有没有错误，没有可对比的，不知道对错。

cam_1980

怎么没有人给看看呀，高手可以把代码运行一下看看有没有问题，弄出这么个图，真不知道对不对，不敢肯定。

无水1324

我看了非线性振动这本书上的Poincare图，差不多就是这样子的

cam_1980

好的，非常感谢无水兄的解答。