如何编写非自治非线性常微分方程的poincare映射

bimhao

如何编写非自治非线性常微分方程的poincare映射

%%%------------------------------------------------------------%%%
建议以后发贴最好是将方程直接贴上来,
这种相当模糊的问法,即使别人知道也不太愿意回答.
by xjzuo
%%%------------------------------------------------------------%%%

[ 本帖最后由 xjzuo 于 2007-1-24 17:00 编辑 ]

happy

参考http://forum.vibunion.com/forum/viewthread.php?tid=12312的相关讨论

这个问题在论坛讨论了很多次了，建议搜索以下能够看到更多的讨论帖子

bimhao

xdot=[x(2);
     
-a0*x(1)-a1*x(2)-a2*x(1)*cos(Q2*t)-a3*x(1)^2-a4*x(3)^2-a5*x(1)*x(3)^2-a6*x(1)^3-a7*x(1)*x(3)-(a8*T11+a9*T12)*x(1)*cos(Q3*t)-(a10*T11+a11*T12)*x(3)*cos(Q3*t)-a12*T11*cos(Q3*t)+f1*cos(Q1*t);
     
x(4);

     -b0*x(3)-b1*x(4)-b2*x(3)*cos(Q2*t)-b3*x(1)*x(3)-b4*x(1)^2-b5*x(1)*x(3)^2-b6*x(1)^2*x(3)-b7*x(3)^3-(b8*T11+b9*T12)*x(1)*cos(Q3*t)-(b10*T11+b11*T12)*x(3)*cos(Q3*t)-b12*T12*cos(Q3*t)+f2*cos(Q1*t)];

bimhao

但是现在需要非自治系统的poincare映射,拿位大侠能帮忙啊!急

bimhao

如何编写

xjzuo

happy教授已经给了你提示，自己动动手先参考编一下，不要老想着让别人帮你编．
另：请勿重复发没有意义的贴．

bimhao

??? Error using ===> odearguments
Inputs to odearguments must be floats, namely single or double.

Error in ==> funfun\private\odearguments at 135
  dataType = superiorfloat(t0,y0,f0);

Error in ==> ode15s at 228
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in ==> fshell at 12
[t,y]=ode15s('shell',tspan,y0,opts);

eight

原帖由 bimhao 于 2007-2-6 19:46 发表
??? Error using ===> odearguments
Inputs to odearguments must be floats, namely single or double.

Error in ==> funfun\private\odearguments at 135
  dataType = superiorfloat(t0,y0,f0);
...

Inputs to odearguments must be floats, namely single or double

参数输入有问题，必须是浮点类型的，自己查看一下 ode15s 文件第 228 行