F=470.327*0.4536*9.8; m=527.436*0.4536; G=527.436*0.4536*9.8; b=0.08*0.4536*9.8/0.3048; ts=[0,2000]; x0=0; [t,x]=ode45('sudu',ts,x0,G,F,m,b); plot(t,x,'-'),grid; xlabel('t'); ylabel('v(t)'); title('速度——时间'); figure; v=x(1:21,:); h=[]; s=0; for i=1:21 if i<=20 s=s+(x(i)+x(i+1))/2; h(i)=s; s=s; else h(21)=h(20)+(x(20)+x(21))/2; end end dept=linspace(91.4400,91.4400,21); plot(v,h,v,dept,'r'),grid;

这是一段 MATLAB 代码，用于求解物体在重力和阻力作用下的速度随时间的变化，并绘制速度-时间和高度-时间图形。其中，F、m、b、x0、ts 分别表示物体受到的外力、质量、阻力系数、初始位置和时间范围。ode45 函数用于求解微分方程，并返回时间和速度的数组 t 和 x。接下来的 for 循环用于计算物体的高度随时间的变化，其中 dept 表示物体的初始高度，h 表示物体的高度随时间的变化。

你可以将代码复制粘贴到 MATLAB 环境中，运行该程序，得到速度-时间和高度-时间图形。