如何在matlab中正确使用interp1来找到抛物线中的多个解决方案？

Question

我已经解决了一组常微分方程，以实现以下解决方案：

f = @(t,x) [(-B/m)*x(1);-g-(B/m)*x(2);x(1);x(2)];
[t,x_sol] = ode45(f,[0 5],[V_0*cos(alpha) V_0*sin(alpha) 0 0]);

并绘制了以下一种解决方案：

plot(t,x_sol(:,4));

我要做的就是找到Y（t）变为零的时间。但是，当我使用t_f = interp1(x_sol(:,4),t,0)时，我得到的第一个解为0。这可能有一个简单的解，但是我似乎无法获得另一个Y变为零的解决方案（距离图表约4秒）。这实际上是一个轨迹问题，我只需要做的就是找到轨迹时间。