我需要找到导数,然后创建一个匿名函数v
,这是弹丸的垂直速度(即v
是高度函数的导数,即上面的函数)
y = @(t) (m*v0*sin(theta)/gamma + g*m^2/gamma^2) * (1 - exp(-gamma/m*(t-t0))) ...
- g*m/gamma * (t - t0) + y0;
答案 0 :(得分:1)
在Maxima中非常容易做,您可以在其中粘贴函数:
y(t) := (m*v0*sin(theta)/gamma + g*m^2/gamma^2) * (1 - exp(-gamma/m*(t-t0))) - g*m/gamma * (t - t0) + y0;
v(t) := diff(y(t),t);
display2d:false;
print(v(t));
输出
(((m*sin(theta)*v0)/gamma+(g*m^2)/gamma^2)*gamma*%e^-(((t-t0)*gamma)/m))/m -(g*m)/gamma
几乎可以直接粘贴到MATLAB中,只需将%e^
替换为exp()
:
v(t) = @(t)(((m*sin(theta)*v0)/gamma+(g*m^2)/gamma^2)*gamma*exp(-(((t-t0)*gamma)/m)))/m -(g*m)/gamma