cam_1980 发表于 2008-8-25 21:51

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

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=;
    a=0.4;
    b=1.2;
    arfa=0.2;
    w=0.8;
    T=2*pi/w;
    tspan=;
    =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 发表于 2008-8-26 11:00

回复 楼主 cam_1980 的帖子

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

无水1324 发表于 2008-8-27 09:12

回复 沙发 cam_1980 的帖子

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

cam_1980 发表于 2008-8-27 12:43

好的，非常感谢无水兄的解答。
页: [1]
查看完整版本: 布鲁塞尔问题的poincare映射对不对？