如何确定常微分方程是属于刚性还是非刚性的?

Anaiad · 发表于 2007-7-4 16:32

您需要登录才可以下载或查看，没有账号？我要加入

x

有些Matlab资料中提到求解器ode45适用与求解非刚性的常微分方程，而ode23S适用与求解刚性的微分常微分方程，这里的刚性是怎么定义的呢?望高手赐教！

octopussheng · 发表于 2007-7-5 09:33

我的办法是试算，呵呵，如果用ode45计算不收敛的话就换个算法试试，再不行的话那就是自己的参数有问题了

appleseed05 · 发表于 2007-7-5 15:42

刚性方程的定义：
特征值的实部为负数，实部最大的的特征值的实部与实部最小的的特征值的实部的比值远大于1
你可以求一下特征值试试

特别感谢gghhjj，前两天他告诉我的:lol :lol

Anaiad · 发表于 2007-7-5 17:16

不太明白，实部既然是负数怎么会远远大于1呢？

[ 本帖最后由 Anaiad 于 2007-7-5 17:22 编辑 ]

appleseed05 · 发表于 2007-7-5 18:28

不是说是比值嘛，比如-500/-10=50

Anaiad · 发表于 2007-7-5 18:54

不好意思，看错了:loveliness:
明白了，谢谢你的帮助:loveliness:

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 如何确定常微分方程是属于刚性还是非刚性的?