Matlab 曲线达到特定x坐标需要多长时间？

如何计算曲线达到特定

x

坐标所需的时间（在Matlab中）。假设我们有：

---

dx/dt=x^2+y^2，dy/dt=5.x.y，曲线从点（a，b）开始。在ode45的帮助下，我得到了曲线图。我还需要计算曲线达到x=c（c>a）所需的时间。我被告知这可以通过插值来完成，但我不知道如何编写代码。

根据系统在

c

附近的行为，在输出上使用数据插值方法，如

interp1

，可能有效，也可能无效。解决这个问题的更严格的方法是使用事件（参见我的答案或），或者使用单结构输出参数形式的

ode45

，并结合常规的数据插值方法。这两种方法都使用设计用于底层ODE的方法。虽然更复杂，但事件可能是像您的案例那样准确确定交叉时间的最佳方式。

您看过吗？是的，但老实说，我没有放弃我应该编写的代码。