我在matlab脚本中有一个微分方程组的解决方案。我将我的函数定义如下:
syms V_x(t) V_y(t) X(t) Y(t);
ode_V_x = diff(V_x,t) == -(B/m)*V_x;
ode_V_y = diff(V_y,t) == -g - (B/m)*V_y;
ode_X = diff(X,t) == V_x;
ode_Y = diff(Y,t) == V_y;
然后我用适当的初始条件解决它们如下:
[V_xSol(t), V_ySol(t), XSol(t), YSol(t)] = dsolve(odes,conds);
例如,Y(t)的解决方案是:
YSol(t) = exp(-t/5)*((981*exp(t/5))/4 - 12891678040772023/35184372088832) - (981*t)/20 + 4262710785985975/35184372088832
现在我需要找到Y(t) = 0
的时间值。我想通过执行:interp1
在matlab中使用t_f = interp1([0,5],YSol,0);
函数,但它不起作用。错误如下:Values V must be of type double or single
。在matlab中使用interp1在symfun
类型函数上的正确方法是什么?附:自变量't'尚未定义为先前persay的任何向量(我理解matlab喜欢向量)。
我不确定interpl
功能是否适合这项工作。您可以使用fsolve
来达到您想要的效果。
fsolve
的用法:
fsolve(function handle, initial guess)
这是针对您在matlab中的问题实现的:
t = fsolve(@(t)exp(-t/5)*((981*exp(t/5))/4 - 12891678040772023/35184372088832) - (981*t)/20 + 4262710785985975/35184372088832, 5)
在这种情况下,答案是4.3243
。请注意,答案取决于您给求解器的初始猜测。在你的情况下,我注意到任何高于2的初始值都应该有效。否则任何东西都会收敛到不同/错误的答案。
你可以在这里阅读更多关于fsolve
的信息:https://www.mathworks.com/help/optim/ug/fsolve.html