Using ode45 to solve system of ODEs with some initial conditions and some terminal conditions

3 次查看（过去 30 天）
Christian
Christian 2022-7-8
评论： Torsten 2022-7-12
I have a system of ODEs with , and with a mix of initial and terminal conditions on x. Specifically, for the first 8 components of x, we have the initial conditons, , while for the last 8 components, we have terminal conditions, . Mathematically, this should be enough to guarantee uniqueness of the solution, and I know that if you ever have only terminal or only initial conditions, ode45 handles it no problem.
The first thought would be to just separate the first 8 and last 8 components and just solve 2 system of ODEs, one with terminal conditions and the other with initial conditions, but unfortunately the equations governing both the first 8 components and the last 8 all involve all 16 components, so separating them out isn't possible.
Basically, I want to know if theres a way to corral ode45 into solving a system of ODEs such as this one, since I couldn't find a built-in way to do it when searching, or if there's some mathematical trick to where I can transform the system to only have initial conditions, or something like that.
Any help is appreciated!

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2022-7-8
This is a boundary-value problem.
Use BVP4C or BVP5C instead of ODE45.
  2 个评论
Christian
Christian 2022-7-12
编辑：Christian 2022-7-12
Awesome, reading the docs looks like that would work, thanks for the fast reply. The MATLAB docs and YouTube videos I tried to watch on bvp4c and bvp5c have only had very basic examples, so I've been having a hard time figuring out all the parameters that bvp4c and bvp5c require but which ode45 does not - specifically the boundary condition function, the bvpinit parameter, and the region parameter, which I guess is pretty much everything new. In my case, is piecewise smooth on 5 different regions, separated by , so the region parameter should be 5 as far as I can tell. But would you mind showing me an example of how I should go about coding f, the boundary conditon function, and the bvpinit to account for the regions and so that bvp4c or bvp5c would accept it? Thanks in advance again

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by