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

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

有些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:
页: [1]
查看完整版本: 如何确定常微分方程是属于刚性还是非刚性的?