Use pdepe to solve an ode

2 次查看（过去 30 天）
Guy Elisha
Guy Elisha 2022-6-30
回答： Torsten 2022-6-30
Hi everyone,
I am using the pdepe splver to solve a set of 2 pdes that arou coupled and it works great. However, I now want to couple another equation, which does not have a flux term. Therefore, by the pdepe notation,
m = 0, c = 1, f = 0, and s = u-1 (the s is just an example, not my actual s).
The problem is that pdepe discretization supports only parabolic and elliptic equations, with flux term involving spatial derivative.
How can I couple all three equations if two equations are pdes and one equation is ode.
Thank you!

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

回答（1 个）

Torsten
Torsten 2022-6-30
So your equation is not a PDE, but a simple ODE.
Either you try to artificially set du/dx = 0 at both ends and see whether the results make sense.
Or you use Bill Greene's code who claims that ODEs can be handled therein.
Or - as you wrote that there is another problem with your equations that pdepe cannot handle (sigmoid function) - discretize your equations in space and use ODE15S to solve the resulting system of ordinary differential equations in the grid points. Look up "method-of-lines" for more details.

类别

在 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