数学建模，希望数学系的高手解答~~谢谢！！

设速度为v,时间为t，距离为h，加速度为 a
1. a=g v=at,dh/dt=v => h=1/2*a*t^2
故 H=1/2*g*T^2=1/2*9.8*5^2=122.5m

2. ma=mg-k1*v dv/dt=a =>
dv/dt=g-k1*v/m => v=(g/k1*m)*(1-exp(-k1/m*t))
dh/dt=v => h=g*m*(k1*t+m*exp(-k1/m*t)-m)/k1^2
故 H=g*m*(k1*t+m*exp(-k1/m*T)-m)/k1^2 代入数据

3. ma=mg-k1*v^2 dv/dt=a =>
dv/dt=g-k1*v^2/m =>
v=tanh((g*m*k1)^(1/2)*t/m)*(g*m*k1)^(1/2)/k1
dh/dt=v =>
h=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1
故 H=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*T/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1

考虑声速 设声速为vs
1.
T0=T-H/vs 下落时间
H=g*m*(k1*t+m*exp(-k1/m*T0)-m)/k1^2
求得H=
vs*(k1*T+(-lambertw(-g*m/k1/vs*exp(-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs))-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs)*m)/k1
代入数据

其他同理，想法很简单，但是计算比较复杂。
上述计算采用matlab完成。

设速度为v,时间为t，距离为h，加速度为 a
1. a=g v=at,dh/dt=v => h=1/2*a*t^2
故 H=1/2*g*T^2=1/2*9.8*5^2=122.5m

2. ma=mg-k1*v dv/dt=a =>
dv/dt=g-k1*v/m => v=(g/k1*m)*(1-exp(-k1/m*t))
dh/dt=v => h=g*m*(k1*t+m*exp(-k1/m*t)-m)/k1^2
故 H=g*m*(k1*t+m*exp(-k1/m*T)-m)/k1^2 代入数据

3. ma=mg-k1*v^2 dv/dt=a =>
dv/dt=g-k1*v^2/m =>
v=tanh((g*m*k1)^(1/2)*t/m)*(g*m*k1)^(1/2)/k1
dh/dt=v =>
h=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1
故 H=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*T/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1

考虑声速 设声速为vs
1.
T0=T-H/vs 下落时间
H=g*m*(k1*t+m*exp(-k1/m*T0)-m)/k1^2
求得H=
vs*(k1*T+(-lambertw(-g*m/k1/vs*exp(-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs))-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs)*m)/k1
代入数据

其他同理，想法很简单，但是计算比较复杂。
上述计算采用matlab完成。
设速度为v,时间为t，距离为h，加速度为 a
1. a=g v=at,dh/dt=v => h=1/2*a*t^2
故 H=1/2*g*T^2=1/2*9.8*5^2=122.5m

2. ma=mg-k1*v dv/dt=a =>
dv/dt=g-k1*v/m => v=(g/k1*m)*(1-exp(-k1/m*t))
dh/dt=v => h=g*m*(k1*t+m*exp(-k1/m*t)-m)/k1^2
故 H=g*m*(k1*t+m*exp(-k1/m*T)-m)/k1^2 代入数据

3. ma=mg-k1*v^2 dv/dt=a =>
dv/dt=g-k1*v^2/m =>
v=tanh((g*m*k1)^(1/2)*t/m)*(g*m*k1)^(1/2)/k1
dh/dt=v =>
h=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1
故 H=-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*T/m)-1)-1/2/k1*m*log(tanh((g*m*k1)^(1/2)*t/m)+1)+1/2*i*m*pi/k1

考虑声速 设声速为vs
1.
T0=T-H/vs 下落时间
H=g*m*(k1*t+m*exp(-k1/m*T0)-m)/k1^2
求得H=
vs*(k1*T+(-lambertw(-g*m/k1/vs*exp(-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs))-(-g*m*k1*t+g*m^2+k1^2*T*vs)/k1/m/vs)*m)/k1
代入数据

其他同理，想法很简单，但是计算比较复杂。
上述计算采用matlab完成。