返回信息流源程序如下:
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, ...
这究竟是什么问题,望各位纠正啊啊啊!!!
这是一条镜像帖。来源:北邮人论坛 / matlab / #11587同步于 2015/4/16
Matlab机器人发帖
关于一个matlab解微分方程组的问题,请教各位!!!
yanghuixuan
2015/4/16镜像同步0 回复
订阅后,新回复会通过你的通知中心匿名送达。
0 条回复
暂无回复 · 你可以订阅本帖等待新回复。