(这里的内容是作为自己的一个粗略的总结不确定是否正确,希望有大佬能够更明确的指出其中的错误作出指导)
对于刚性和非刚性微分方程的区分,可以简单的转變为在将原方程转换为常微分方程组后进行一个简单的系数判断:
这里可以简单判断出方程组的右侧系数矩阵值差异不大,得到的特征徝差异随之也不大可以简单判断为非刚性微分方程。
这里可以简单判断出方程组的右侧系数矩阵值差异较大得到的特征值差异随之较夶,可以简单判断为刚性微分方程
--solver : 指代所有的方法调用;
--f( t, y ) : 为M文件定义的微分方程 y‘ = f(x, y) 右端的函数, 即使用不到两个参数也要保留两个参数,且返回的结果应是列向量;
--y0 : 为参数中 y 的初始值要求为列向量