d2x/dt2=0.435(U-dx/dt)[(U-dx/dt)^2+(dz/dt)^2]^(1/2)加速度方程
d2z/dt2=1.1856(-dz/dt)[(U-dx/dt)^2+(dz/dt)^2]^(1/2)-9.8加速度方程
U=20(z/10)^0.12
根据方程积分求运动轨迹,以及x、y方向的速度
t=0时x=0,z=500,v=0,时间间隔t=0.1s,求到z=0停止
谢谢,不胜感激
这个问题我之前回答过(编号2117613008218288267),当时因为题主给的数据有问题没能彻底解决。现在给的条件基本上一个没问题了,给你新版的程序:
function main和之前的差别,除了模型按照现在新的模型之外,最主要是增加了事件函数,以实现z到达0就能够结束。另外,微分方程的写法之前用匿名函数,现在改用具名函数,更清晰一些。仿真的最后一个点可能出现复数,所以绘图时删除最后一个点。