fhuabai2.m

来自「利用matlab进行动力学仿真之一」· M 代码 · 共 10 行

function ydot=fhuabai(t,y)
global l g ma mb
x1=[ma+mb mb*l*cos(y(2))
   cos(y(2)) l];
x2=[y(4)*y(4)*mb*l*sin(y(2)) 
   -g*sin(y(2))];
ydot=[y(3)
   y(4)
   inv(x1)*x2];