matlab力学题目求解,求代码一质量为1kg的物体从静止开始下落,假设受到的空气阻力大小与速度成正比,方向相反：f=-v(t),假设初始时刻,t=0时,x(0)=0,v(0)=0以向下为正方向,建立此系统的常微分方程,

Question

matlab力学题目求解,求代码
一质量为1kg的物体从静止开始下落,假设受到的空气阻力大小与速度成正比,方向相反：f=-v(t),假设初始时刻,t=0时,x(0)=0,v(0)=0以向下为正方向,建立此系统的常微分方程,并用ode23求解出0到5秒之间,物体的运动轨迹以及速度大小变化

Answer 1

f=@(t,x) [-1;x(1);x(2)];%x(1)=a,x(2)=v,x(3)=s
[t x]=ode45(f,[0,5],[9.8,0,0]);
plot(t,x(:,2),t,x(:,3))