登录论坛

查看完整版本 : [MATLAB基础] matlab新手请教关于物体下落 的问题


晒月亮的狼
2012-06-25, 17:34
题目:75kg的物体下落,下落的阻力和速度的关系为F=-kv^2,(k=0,22kg/m),求下落过程中速度v(t)和路程x(t)的图像。
****************************************************
据题意有如下代码:
v=dsolve('m*D2v+2*k*v*Dv=0','v(0)=0,Dv(0)=g','t')
x=dsolve('m*g-k*(Dx)^2=m*D2x','x(0)=0,Dx(0)=0,D2x(0)=g','t')
m=75;
k=0.22;
g=9.8
t=0:0.5:30;
v =tanh(t.*(g*m*k)^(1/2)/m)*(g*m*k)^(1/2)/k
x=(m/k)*log(cosh(t.*(k*g/m)^1/2))
plot(t,v),gtext('v(t)')
hold on
plot(t,x),gtext('x(t)')
**************************************************
我不知道怎么解释第七八行。有人说是第一二行运行的结果,但是我单独运行一二行结果是这样的:
v =

tanh((g*m*k)^(1/2)*t/m)*(g*m*k)^(1/2)/k



x =

1/2*(-2*g^(1/2)*k^(1/2)*t*m^2-m^(5/2)*log(4)+m^(5/2)*log((exp(1/m^(1/2)*g^(1/2)*k^(1/2)*t)^4*g*m+2*exp(1/m^(1/2)*g^(1/2)*k^(1/2)*t)^2*g*m+g*m)/g/m))/m^(3/2)/k
*********************************************************
想知道第七八行是什么意思?如果第七八行是第一二行运行的结果,可不可以修改程序,不要第七八行而直接得出结果?
谢谢各位了!