yang95478 发表于 2006-11-2 22:52

请教求解刚性微分方程组的问题

现有一刚性比为2.5*10^4的线性微分方程组,请问有什么好的算法对其进行求解?即求解刚性微分方程组的方法有哪些?望高手赐教,谢谢!

dujunmin 发表于 2006-11-18 12:53

不知你是要自己编程序求解呢还是只想求出结果,如果是后者,直接利用Matlab中的ode类函数即可求解,方便的很。

heyong2002 发表于 2006-11-18 15:04

回复

求解刚性微分方程常用吉尔方法,他有标准的源程序,你只要编写自己的微分方程组就行了,他可以求出数值解

mxtfirst 发表于 2006-11-18 16:14

Malab的 ode15s ode23s0de23tb这三个可以试下,专门求解刚性ODE的

如果需要过程,可以选择DASSL求解器试下,它是用来求解刚性常微分方程和微分—代数方程的,
Fortran语言的源代码 http://www.netlib.org/ 有下载
页: [1]
查看完整版本: 请教求解刚性微分方程组的问题