关于#matlab#的问题,请各位专家解答!

考虑如下的运动方程

img

核心在于求解ode

b=1;
m=1;
g=9.8;
odefun = @(t,x) [x(2);  -b/m*sqrt(abs(x(2)))*sign(x(2)); x(4); -g-b/m*sqrt(abs(x(4)))*sign(x(4))];
x0 = [0;10;0;10];
[t,x] = ode45(odefun, [0,2], x0);
plot(x(:,1), x(:,3)); title('运动轨迹')