如何编写非自治非线性常微分方程的poincare映射
如何编写非自治非线性常微分方程的poincare映射%%%------------------------------------------------------------%%%
建议以后发贴最好是将方程直接贴上来,
这种相当模糊的问法,即使别人知道也不太愿意回答.
by xjzuo
%%%------------------------------------------------------------%%%
[ 本帖最后由 xjzuo 于 2007-1-24 17:00 编辑 ] 参考http://forum.vibunion.com/forum/viewthread.php?tid=12312的相关讨论
这个问题在论坛讨论了很多次了,建议搜索以下能够看到更多的讨论帖子
这样一个非自治非线性系统
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)];
但是现在需要非自治系统的poincare映射,拿位大侠能帮忙啊!急
编写非自治非线性常微分方程的poincare映射
如何编写回复
happy教授已经给了你提示,自己动动手先参考编一下,不要老想着让别人帮你编.另:请勿重复发没有意义的贴.
错误提示
??? Error using ===> odeargumentsInputs 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
=ode15s('shell',tspan,y0,opts); 原帖由 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 行
页:
[1]