5

我使用 ode45 来求解微分方程,但 tspan 始终必须使系统及时向前运行。我可以让 ode45 将系统向后运行到负 t 吗?

我遇到的具体问题使用了一个上升指数函数,即从 t=-infinity 到 t=0 的 a*exp(at),对于所有 t>0,该函数为零。如您所见,在我的问题中不可能翻译此函数并使用区间 t=0 到 t=infinity。所以,我希望 ode45 在负 t 下工作。

谢谢!

编辑:经过一番思考,我意识到我的问题措辞很糟糕,而且我的理解也很差,所以让我在这里重申一下。

我的 ODE 的形式为 x'=f(t)x,其中 f(t)=a*exp(at) 从 -infinity 到 0。系统从 t=-infinity 的一些初始条件开始。另一方面,ode45 要求 tspan 为正,并且还需要 t=0 时的初始条件。那么我怎样才能改变这两个标准呢?

谢谢 :)

4

1 回答 1

0

使用向后的 linspace 生成您的 t。

对于正向 (t),使用类似 t = 0:0.1:2 的值。对于反向 (rt),使用类似 rt = 2:0-.1:0 的值。

于 2013-02-14T04:37:30.250 回答