Variable precision arithmetic in ode45 and others?

18 次查看（过去 30 天）

显示更早的评论

Michael Croucher 2020-10-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630869-variable-precision-arithmetic-in-ode45-and-others

评论： Walter Roberson 2022-6-23

采纳的回答： Walter Roberson

Is it possible to use vpa in ordinary differential equation solvers like ode45?

采纳的回答

Walter Roberson 2020-10-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630869-variable-precision-arithmetic-in-ode45-and-others#answer_528469

编辑：Walter Roberson 2020-10-30

Sort of. The internal symbolic engine has two numeric ode routines. At the moment I do not know if there is a matlab level interface to those routines.

You might need to use feval(symengine) to access the routines.

http://www.cfm.brown.edu/people/dobrush/am33/MuPad/MuPad22.html

9 个评论
显示 7更早的评论隐藏 7更早的评论

FastCar 2022-6-23

Any luck?

Walter Roberson 2022-6-23

Mathworks said that those routines are gone and will not be back.

请先登录，再进行评论。

更多回答（2 个）

John D'Errico 2020-10-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630869-variable-precision-arithmetic-in-ode45-and-others#answer_528454

Short answer - no.

Long answer - nnnnoooooo.

ODE45 uses only singles or doubles. Nothing symbolic. vpa is a symbolic tool.

2 个评论
显示无隐藏无

Michael Croucher 2020-10-29

Thanks John. I thought as much but wanted to check.

John D'Errico 2020-10-29

I suppose it may happen one day. For example, vpaintegral appeared a few years ago. They could provide a tool to do the same for an ODE. But I don't see it happening soon.

请先登录，再进行评论。

Steven Lord 2020-10-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/630869-variable-precision-arithmetic-in-ode45-and-others#answer_529584

What types of differential equations are you trying to solve that you need to use arbitrary precision arithmetic?

If you're trying to solve the ODE symbolically rather than numerically, take a look at the dsolve function in Symbolic Math Toolbox.

If you have computed a symbolic equation for the right-hand side of your system of ODEs, you can convert that into a form that the numeric ODE solvers in MATLAB (like ode45) can accept that operates on double arrays using odeFunction also in Symbolic Math Toolbox.

This documentation page may be of interest.