Ode23没有完成执行

Question

我是Matlab的新手，我正在尝试求解常数n的微分方程y'= - y / n。我在这样的脚本中定义函数：

pyproj

然后尝试解决这个问题：

Proj4

但是操作永远不会完成执行。它只是吃了我的RAM并在角落里说“忙”，直到我按ctrl + c来终止它。我在这里做错了什么？

Answer 1

你是下流的受害者。会发生的是没有无限数量的浮点数（见Is floating point math broken?），因此结果并不完全精确。当数字足够小时，计算机可能会将该数字视为0。 matlab用于ode23的过程基于有限差分，它涉及除法。考虑到下溢误差，可能会出现0问题的划分，或者可能溢出一个非常小的除数，趋向于无穷大并且不满足求解条件。