Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2012-06-25
年龄: 33
帖子: 2
声望力: 0 ![]() |
![]()
题目: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 ********************************************************* 想知道第七八行是什么意思?如果第七八行是第一二行运行的结果,可不可以修改程序,不要第七八行而直接得出结果? 谢谢各位了! |
![]() |
![]() |