Solving a nonlinear ODE with derivative squared
2 次查看(过去 30 天)
显示 更早的评论
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!
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/202374/image.png)
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?
采纳的回答
Star Strider
2019-2-1
编辑:Star Strider
2019-2-3
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 个评论
Star Strider
2019-2-3
As always, my pleasure!
I‘m not sure if it’s possible express systems of PDEs in the Symbolic Math Toolbox.
You most likelly need the Partial Differential Equation Toolbox (link). I haven’t used it recently, so I have no recent experience with it.
更多回答(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!