BBYR Achieve
返回信息流
这是一条镜像帖。来源:北邮人论坛 / matlab / #11587同步于 2015/4/16
Matlab机器人发帖

关于一个matlab解微分方程组的问题,请教各位!!!

yanghuixuan
2015/4/16镜像同步0 回复
源程序如下: function xdot=dox(t,x) global tao; m1=4;m2=2;g=9.8; c=[cos(x(1)) cos(x(2))]; s=[sin(x(1)) sin(x(2))]; c12=cos(x(1)+x(2)); s12=sin(x(1)+x(2)); xdot(1)=x(3); xdot(2)=x(4); D=[m1+2*m2+2*m2*c(2) m2+m2*c(2); m2*c(2)+m2 m2]; H=[-m2*s(2)*x(4).^2-2*m2*s(2)*x(3)*x(4)+m2*g*c12+(m1+m2)*g*c(1);m2*s(2)*x(3).^2+m2*g*c12]; xdot(3:4)=-inv(D)*H+inv(D)*tao; 主控制程序: global tao; tao=[7;12]; [t,x]=ode45('dox',[0 1],[3 -2 0 0]); plot(t,x(:,3),'w--') 结果报错为如下: 错误使用 odearguments (line 80) tspan 中的最后一个元素必须与第一个元素不同。 出错 ode45 (line 113) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... 这究竟是什么问题,望各位纠正啊啊啊!!!
订阅后,新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。