Solving a nonlinear ODE with derivative squared

Question

Travis 2019-2-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/442717-solving-a-nonlinear-ode-with-derivative-squared

评论： Bill Greene 2019-2-4

I'm trying to solve a nonlinear ODE which looks something like this:

. I know I can use the implicit solver ode15i but the problem is also stiff so I'd prefer to use ode15s. Is it possible to solve this type of nonlinear ode using ode15s? Any suggestions would be appreciated, thank you!

2 个评论
显示无隐藏无

Torsten 2019-2-4

As for all quadratic equations, there are two solutions for y'. Do you know which one you'll have to take ?

Bill Greene 2019-2-4

ode15i is based on backward differentiation formulas so I would expect it to be as effective as ode15s for stiff problems. That has also been my experience with the two solvers. Do you have an example stiff ODE where this is not the case?

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-2-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/442717-solving-a-nonlinear-ode-with-derivative-squared#answer_359105

编辑：Star Strider 2019-2-3

在 MATLAB Online 中打开

One approach:

syms a b c d y(t) T Y
Dy = diff(y);
DE = a*Dy^2 + b*Dy + c*y == d;
isoDE = isolate(DE,Dy)
[VF,Sbs] = odeToVectorField(isoDE)
odefcn = matlabFunction(VF, 'Vars',{T,Y,a b c d});
odefcn = @(T,Y,a,b,c,d)[((b+sqrt(a.*d.*4.0+b.^2-a.*c.*Y(1).*4.0)).*(-1.0./2.0))./a; ((b-sqrt(a.*d.*4.0+b.^2-a.*c.*Y(1).*4.0)).*(-1.0./2.0))./a]
a = 3;
b = 5;
c = 7;
d = 11;
[T,Y] = ode15s(@(T,Y)odefcn(T,Y,a,b,c,d), [0 5], [0;0]);
figure
plot(T, Y)
grid

It works!

2 个评论
显示无隐藏无

Travis 2019-2-3

Thanks for your help! Do you know if it's possible to solve this nonsymbollically? The actual equation I'm trying to solve is part of a system of PDEs which I'm discretizing in one domain to turn into a system of ODEs so I can use ode15s. So I'm unsure if I'd be able to do that symbolically.