考虑如下的运动方程
核心在于求解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('运动轨迹')